Search

In this study, wind speed sensitivity and calibration bias were first determined for Surface Water and Ocean Topography (SWOT) satellite Ka-band Radar Interferometer (KaRIn) Normalized Radar Backscatter Cross Section (NRCS) data at VV and HH polarizations. Here, the calibration bias was estimated by comparing the KaRIn NRCS with collocated simulations from a model developed using Global Precipitation Measurement (GPM) satellite Dual-frequency Precipitation Radar (DPR) data. To recalibrate the bias, the correlation coefficient between the KaRIn data and the simulations was estimated, and the data with the corresponding top 10% correlation coefficients were used to estimate the recalibration coefficients. After recalibration, a Ka-band NRCS model was developed from the KaRIn data to retrieve ocean surface wind speeds. Finally, wind speed retrievals were evaluated using the collocated European Center for Medium-Range Weather Forecasts (ECMWF) reanalysis winds, Haiyang-2C scatterometer (HY2C-SCAT) winds and National Data Buoy Center (NDBC) and Tropical Atmosphere Ocean (TAO) buoy winds. Evaluation results show that the Root Mean Square Error (RMSE) at both polarizations is less than 1.52 m/s, 1.34 m/s and 1.57 m/s, respectively, when compared to ECMWF, HY2C-SCAT and buoy collocated winds. Moreover, both the bias and RMSE were constant with the incidence angles and polarizations. This indicates that the winds from the SWOT KaRIn data are capable of correcting the sea state bias for sea surface height products. [ABSTRACT FROM AUTHOR]

Temperature and water vapor play crucial roles in the Earth's climate system, and it is important to understand and monitor the variation in the thermodynamic profile within the lower troposphere. Among various observation platforms for understanding the vertical structure of temperature and humidity, ground-based Fourier-transform infrared (FTIR) can provide detailed information about the lower troposphere by complementing the limitations of radiosonde or satellite methods. However, these ground-based systems have limitations in terms of cost, operation, and mobility. Herein, we introduce a cost-effective and easily deployable FTIR observation system designed to enhance monitoring capabilities for atmospheric conditions. The atmospheric downwelling radiance spectrum of sky is measured by applying a real-time radiative calibration using a blackbody. From the observed radiance spectrum, the thermodynamic profile (temperature and the water vapor mixing ratio) of the lower troposphere was retrieved using an algorithm based on the optimal estimation method (OEM). The retrieved vertical structure results in the lower troposphere were similar to the fifth-generation reanalysis database (ERA-5) of the European Center for Medium-range Weather Forecasts (ECMWF) and the National Centers for Environmental Prediction final analysis (NCEP FNL). This provides a potential possibility for monitoring atmospheric conditions by a compact FTIR system. [ABSTRACT FROM AUTHOR]

Based on the half-hourly Integrated Multi-Satellite Retrievals for Global Precipitation Measurement (IMERG) data product (0.1°), the fifth-generation European Center for Medium-Range Weather Forecasting atmospheric reanalysis dataset (ERA5), sounding data, and the Weather Research and Forecasting Model (WRF-ARW), this study explored the developmental process of a typical extreme precipitation event in Liaoning Province on 2 June 2021. This study focused on the impact of Changbai Mountain on this precipitation process and its corresponding physical mechanisms. The research findings revealed that Changbai Mountain significantly affected the precipitation event in three main aspects: blocking drag, forcing uplift, and leeside convergence. The blocking drag caused by the mountain topography led to an extension in the duration of heavy rainfall. The dynamic lifting and leeside convergence associated with the mountainous terrain also substantially increased the amount of precipitation. Furthermore, the topography hindered the movement of the Bohai Bay cold pool and enhanced the intensity of the cold pool, contributing to the sustained extreme precipitation in Liaoning Province. Lastly, the terrain sensitivity experiment demonstrated that when the height of Changbai Mountain was reduced, the convergence uplift, moisture condensation, and cold pool intensity were weakened, leading to significant changes in precipitation intensity and spatial distribution. These findings further confirm the crucial role of Changbai Mountain in the occurrence and development of local precipitation in Liaoning Province. [ABSTRACT FROM AUTHOR]

The observed radiation data from the second-generation Hyperspectral Infrared Atmospheric Sounder (HIRAS-II) on the Fengyun-3E (FY-3E) satellite contain useful vertical atmosphere information which can distinguish and retrieve vertical profiles of atmospheric gas components including ozone (O3), carbon monoxide (CO), and methane (CH4). This paper utilizes FY-3E/HIRAS-II observational data to optimize each gas channel using the improved Optimal Sensitivity Profile method (OSP) channel algorithm and establishes a typical convolutional neural network model (CNN) and a representative U-shaped network model (UNET) with deep features and shallow feature links to perform atmospheric profile retrieval calculations of O3, CO, and CH4. We chose the clear sky data of the Indian and its southern seas in December 2021 and January 2022, with reanalysis data from European Center for Medium-Range Weather Forecasts Reanalysis v5 (ERA5) and European Center for Medium-Range Weather Forecasts Atmospheric Composition Reanalysis v4 (EAC4) serving as the reference values. The retrieval outcomes were then compared against advanced numerical forecast models including the Whole Atmosphere Community Climate Model (WACCM), Global Forecast System (GFS), and satellite products from an Atmospheric Infrared Sounder (AIRS) and Infrared Atmospheric Sounding Interferometer (IASI). Experimental results show that the generalization ability and retrieval accuracy of CNN are slightly higher compared with UNET. For O3 profile retrieval, the mean percentage error (MPE) of the whole layers for CNN and UNET data in relation to ERA5 data was less than 8%, while the root-mean-square error (RMSE) was below 1.5 × 10−7 kg/kg; for CH4 profile retrieval, the MPE of the whole layers for CNN and UNET data in relation to EAC4 data was less than 0.7%, while the RMSE was below 1.5 × 10−8 kg/kg. The retrieval of O3 and CH4 are resulted in a significant improvement compared to the forecast data and satellite products in most pressure levels; for CO profile retrieval, the MPE of the whole layers for CNN and UNET data in relation to EAC4 data was less than 11%, while the RMSE was below 4 × 10−8 kg/kg. The error of the CO retrieval results was higher than that of the forecast data at the pressure level of 200~500 hPa and lower than that of similar satellite products with most pressure levels. The experiments indicated that the neural network method effectively determines the atmospheric gas profiles using infrared hyperspectral data, exhibiting a positive performance in accuracy and retrieval speed. [ABSTRACT FROM AUTHOR]

Sea surface temperature (SST) has important impacts on the global ecology, and having a good understanding of the predictability, i.e., the possibility of achieving accurate prediction, of SST can help us monitor the marine environment and climate change, and guide the selection and design of SST prediction methods. However, existing studies for analyzing SST mostly measure the rising or falling trends of SST. To address this issue, we introduce a temporal-correlated entropy to quantify the predictability of SST series from both global coarse-grained and local fine-grained aspects, and make SST prediction with multiple deep learning models to prove the effectiveness of such predictability evaluation method. In addition, we explore the dynamics of SST predictability by dividing the time range of interest into consecutive time periods, evaluating the corresponding predictability of SST for each time period, and analyzing the stability of the predictability of SST over time. According to the experiments, the SST predictability values near the poles and equator are really high. The average SST predictability values of the East China Sea, Bohai Sea, and Antarctic Ocean are 0.719, 0.706, and 0.886, respectively, and the size relationship of the SST predictability in the three local sea areas is consistent with our prediction results using multiple representative SST prediction methods, which corroborates the reliability of the predictability evaluation method. In addition, we found that the SST predictability in the Antarctic Ocean changes more dramatically over time than in the East China Sea and the Bohai Sea. The results of SST predictability and its dynamic analysis indicate that global warming, ocean currents, and human activities all have significant impacts on the predictability of SST. [ABSTRACT FROM AUTHOR]

Accurate precipitation forecasting is challenging, especially on the sub-seasonal to seasonal scale (14–90 days) which mandates the bias correction. Quantile mapping (QM) has been employed as a universal method of precipitation bias correction as it is effective in correcting the distribution attributes of mean and variance, but neglects the correlation between the model and observation data and has computing inefficiency in large-scale applications. In this study, a quantile mapping of matching precipitation threshold by time series (MPTT-QM) method was proposed to tackle these problems. The MPTT-QM method was applied to correct the FGOALS precipitation forecasts on the 14-day to 90-day lead times for the Pearl River Basin (PRB), taking the IMERG-final product as the observation. MPTT-QM was justified by comparing it with the original QM method in terms of precipitation accumulation and hydrological simulations. The results show that MPTT-QM not only improves the spatial distribution of precipitation but also effectively preserves the temporal change, with a better precipitation detection ability. Moreover, the MPTT-QM-corrected hydrological modeling has better performance in runoff simulations than the QM-corrected modeling, with significantly increased KGE metrics ranging from 0.050 to 0.693. MPTT-QM shows promising values in improving the hydrological utilities of various lead time precipitation forecasts. [ABSTRACT FROM AUTHOR]

Data-driven machine learning technology can learn and extract features, a factor which is well recognized to be powerful in the warning and prediction of severe weather. With the large-scale deployment of the radar wind profile (RWP) observational network in China, dynamical variables with higher temporal and spatial resolution in the vertical become strong supports for machine-learning-based severe convection prediction. Based on the RWP mesonet that has been deployed in Beijing, this study uses the measurements from four triangles composed of six RWP stations to determine the profiles of divergence, vorticity, and vertical velocity before rainfall onsets. These dynamic feature variables, combined with cloud properties from Himawari-8 and ERA-5 reanalysis, serve as key input parameters for two rainfall forecast models based on the random forest (RF) classification algorithm. One is for the rainfall/non-rainfall forecast and another for the rainfall grade forecast. The roles of dynamic features such as divergence, vorticity, and vertical velocity are examined from ERA-5 reanalysis data and RWP measurements. The contribution of each feature variable to the performance of the RF model in independent tests is also discussed here. The results show that the usage of RWP observational data as the RF model input tends to result in better performance in rainfall/non-rainfall forecast 30 min in advance of rainfall onset than using the ERA-5 data as inputs. For the rainfall grade forecast, the divergence and vorticity that were estimated from the RWP measurements at 800 hPa show importance in improving the model performance in heavy and moderate rain forecasts. This indicates that the atmospheric dynamic variable measurements from RWP have great potential to improve the prediction skill of convection with the aid of a machine learning model. [ABSTRACT FROM AUTHOR]

The incident and emitted radiative energy data for the top of the atmosphere (TOA) are essential in climate research. Since East Asia (11–61°N, 80–175°E) is complexly composed of land and ocean, real-time satellite data are used importantly for analyzing the detailed energy budget or climate characteristics of this region. Therefore, in this study, the radiative energy budget for East Asia, during the year 2021, was analyzed using GEO-KOMPSAT-2A/Advanced Metrological Imager (GK-2A/AMI) and the European Centre for Medium-range Weather Forecasts reanalysis (ERA5) data. The results showed that the net fluxes for the TOA and surface were −4.09 W·m−2 and −8.24 W·m−2, respectively. Thus, the net flux difference of 4.15 W·m−2 between TOA and surface implied atmospheric warming. These results, produced by GK-2A/AMI, were well-matched with the ERA5 data. However, they varied with surface characteristics; the atmosphere over ocean areas warmed because of the large amounts of longwave radiation emitted from surfaces, while the atmosphere over the plain area was relatively balanced and the atmosphere over the mountain area was cooled because large amount of longwave radiation was emitted to space. Although the GK2A/AMI radiative products used for this study have not yet been sufficiently compared with surface observation data, and the period of data used was only one year, they were highly correlated with the CERES (Clouds and the Earth's Radiant Energy System of USA), HIMAWARI/AHI (Geostationary Satellite of Japan), and ERA5 data. Therefore, if more GK-2A/AMI data are accumulated and analyzed, it could be used for the analysis of radiant energy budget and climate research for East Asia, and it will be an opportunity to greatly increase the utilization of total meteorological products of 52 types, including radiative products. [ABSTRACT FROM AUTHOR]

Ocean primary productivity generated by phytoplankton is critical for ocean ecosystems and the global carbon cycle. Accurate ocean primary productivity forecasting months in advance is beneficial for marine management. Previous persistence-based prediction studies ignore the temporal memories of multiple relevant factors and the seasonal forecasting skill drops quickly with increasing lead time. On the other hand, the emerging ensemble climate forecasts are not well considered as new predictability sources of ocean conditions. Here we proposed a joint forecasting model by combining the seasonal climate predictions from ten heterogeneous models and the temporal memories of relevant factors to examine the monthly predictability of ocean productivity from 0.5- to 11.5-month lead times. The results indicate that a total of ~90% and ~20% productive oceans are expected to be skillfully predicted by the combination of seasonal SST predictions and local memory at 0.5- and 4.5-month leads, respectively. The joint forecasting model improves by 10% of the skillfully predicted areas at 6.5-month lead relative to the prediction by productivity persistence. The hybrid data-driven and model-driven forecasting approach improves the predictability of ocean productivity relative to individual predictions, of which the seasonal climate predictions contribute largely to the skill improvement over the equatorial Pacific and Indian Ocean. These findings highlight the advantages of the integration of climate predictions and temporal memory for ocean productivity forecasting and may provide useful seasonal forecasting information for ocean ecosystem management. [ABSTRACT FROM AUTHOR]

We combined observations of ocean surface winds from Indian SCATterometer SATellite-1 (SCATSAT-1) with a background wind field from a numerical weather prediction (NWP) model available at National Centre for Medium-Range Weather Forecast (NCMRWF) to generate a 6-hourly gridded hybrid wind product. A distinctive feature of the study is to produce a global gridded wind field from SCATSAT-1 scatterometer passes with spatio-temporal data gaps at regular synoptic hours relevant for forcing models and other NWP studies. We are following the concept from the modern particle filter technique, which does not represent the model probability density function (PDF) as Gaussian. We generated the 6-hourly hybrid winds for 2018 and validated them using the wind speed from daily gridded level-4 SCATSAT-1 winds (L4AW), Cross Calibrated Multi-Platform (CCMP) dataset and global buoy data from National Data Buoy Centre (NDBC). The results suggest the potential of the technique to produce scatterometer winds at the desired temporal frequency with significantly less noise and bias along the swath. The study shows that the generated hybrid winds are of prime quality compared with the already existing daily products available from Indian Space Research Organization (ISRO). [ABSTRACT FROM AUTHOR]

High-resolution meteorological datasets are urgently needed for understanding the hydrological cycle of the Tibetan Plateau (TP), where ground-based meteorological stations are sparse. Rapid advances in remote sensing create possibilities to represent spatiotemporal properties of precipitation at a high resolution. In this study, the hourly precipitation characteristics over the TP from two gridded precipitation products, one from global reanalysis (the fifth generation of the European Center for Medium-Range Weather Forecasts atmospheric reanalysis of the global climate; ERA5) and the other is simulated by Global-to-Regional Integrated forecast SysTem (GRIST) global nonhydrostatic model, are compared against satellite-gauge merged precipitation analysis (China Merged Precipitation Analysis; CMPA) from 27 July to 31 August 2014, and a satellite-retrieved precipitation estimate from the Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement (IMERG) is also evolved. Two aspects are mainly focused on: the spatial distribution and the elevation dependence of hourly precipitation characteristics (including precipitation amount, frequency, intensity, diurnal variations, and frequency–intensity structure). Results indicate that: (1) The precipitation amount, frequency, and intensity of CMPA and IMERG decrease with altitude in the Yarlung Tsangpo river valley (YTRV), but increase at first and then decrease with altitude (except for intensity) in the eastern periphery of TP (EPTP). ERA5 performed well on the variation of precipitation amount with altitude (especially in EPTP), but poorly on the frequency and intensity. GRIST is the antithesis of ERA5, but they all overestimate (underestimate) the frequency (intensity) at all heights; (2) With increasing altitude, the diurnal phase of precipitation of CMPA and IMERG shifted from night to evening in the two sub-regions. IMERG's diurnal phase is 1 to 3 h earlier than CMPA's, and the discrepancy decreases (increases) as the altitude increases in YTRV (EPTP). The diurnal phase of precipitation amount and frequency in ERA5 and GRIST is significantly earlier than CMPA, and the frequency peaks around midday except in the basin. GRIST's simulation of the diurnal variation in intensity at various altitudes is consistent with CMPA; (3) ERA5 and GRIST overestimate (underestimate) the frequency of weak (intense) precipitation, with ERA5's deviance being the most severe. The deviations increased with altitude. These findings provide intensive metrics to evaluate precipitation in complex terrain and are helpful for deepening the understanding of simulated biases for further improving performance in high-resolution simulation. [ABSTRACT FROM AUTHOR]

Satellite-based and reanalysis precipitation products have experienced increasing popularity in agricultural, hydrological and meteorological applications, but their accuracy is still uncertain in different areas. In this study, six frequently used high-resolution daily precipitation products, including Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), Global Satellite Mapping of Precipitation (GSMaP), Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG), Multi-Source Weighted-Ensemble Precipitation (MSWEP), Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Cloud Classification System-Climate Data Record (PERSIANN-CCS-CDR) and European Center for Medium-Range Weather Forecasts Reanalysis V5-Land (ERA5-Land), were comprehensively evaluated and compared in nine regions of mainland China between 2015 and 2019. The results reveal that, in general, GSMaP is the best precipitation product in different agricultural regions, especially based on the Pearson correlation coefficient (CC) and critical success index (CSI). ERA5-Land and MSWEP tend to have the highest probability of detection (POD) values, and MSWEP tends to have the smallest relative root mean squared error (RRMSE) values. GSMaP performs better at almost all precipitation levels and in most agricultural regions in each season, while MSWEP has the best performance for capturing the time series of mean daily precipitation. In addition, all precipitation products perform better in summer and worse in winter, and they are more accurate in the eastern region. The findings of this study will contribute to understanding the uncertainties of precipitation products, improving product quality and guiding product selection. [ABSTRACT FROM AUTHOR]

In this study, the polarization characteristics of the newly launched Fengyun-3E (FY-3E) microwave sounding instruments are discussed, and its data quality is also assessed using one month of observation by the double-difference method. By comparison with the equivalent channels onboard Fengyun-3D (FY-3D) and advanced technology microwave sounder (ATMS), the data quality of FY-3E Microwave Humidity Sounder-II (MWHS-II) is improved and comparable to ATMS, while the data of FY-3E Microwave Temperature Sounder-III (MWTS-III) are slightly worse than data of FY-3D. The data of FY-3E MWTS-III are more susceptible to the early-morning orbit than the data of MWHS-II. In addition, striping noise is still present in channels 5–10 of MWTS-III. After the assessments, FY-3E microwave data are preprocessed and assimilated in the global forecast system for the Chinese Meteorology Administration (CMA-GFS). A total of six individual experiments over the period from 16 July to 15 August 2021 were conducted and the impact was evaluated with the composite score used in operation. It is shown that not only the forecasts for the southern hemisphere and tropics are improved significantly, but also the predictions for the northern hemisphere show some improvements in an overall neutral change from adding FY-3E microwave sounding instruments. The impact of FY-3E microwave radiance is equivalent to ATMS as they are assimilated individually. Furthermore, we note that the forecast impact is affected by the cloud detection scheme to a large extent. [ABSTRACT FROM AUTHOR]

The Arctic stratospheric polar vortex is a key component of the climate system, which has significant impacts on surface temperatures in the mid-latitudes and polar regions. Therefore, understanding polar vortex variability is helpful for extended-range weather forecasting. The present study evaluates long-term changes in the position and strength of the polar vortex in the Arctic lower stratosphere during the winters from 1980/81 to 2013/14. Simulations of the Coupled Model Intercomparison Project Phase 6 (CMIP6) models are compared with Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA2) reanalysis dataset. Overall, the CMIP6 models well capture the spatial characteristics of the polar vortex with spatial correlation coefficients between the potential vorticity (PV) in the lower stratosphere from simulations and MERRA2 products generally greater than 0.85 for all CMIP6 models during winter. There is a good agreement in the position and shape of the polar vortex between the CMIP6 multi-model mean and MERRA2, although there exist differences between simulations of individual CMIP6 models. However, most CMIP6 models underestimate the strength of polar vortex in the lower stratosphere, with the largest negative bias up to about −20%. The present study further reveals that there is an anticorrelation between the polar vortex strength bias and area bias simulated by CMIP6 models. In addition, there is a positive correlation between the trend of EP-flux divergence for wavenumber one accumulated in early winter and the trend in zonal mean zonal wind averaged in late winter. As for the long-term change in polar vortex position, CanESM5, IPSL-CM5A2-INCA, UKESM1-0-LL, and IPSL-CM6A-LR well capture the persistent shift of polar vortex towards the Eurasian continent and away from North America in February, which has been reported in observations. These models reproduce the positive trend of wavenumber-1 planetary waves since the 1980s seen in the MERRA2 dataset. This suggests that realistic wave activity processes in CMIP6 models play a key role not only in the simulation of the strength of the stratospheric polar vortex but also in the simulation of the polar vortex position shift. [ABSTRACT FROM AUTHOR]

Data fusion is an effective method to obtain high-precision and high-spatiotemporal-resolution precipitable water vapor (PWV) products, which play an important role in understanding climate change and meteorological monitoring. However, existing fusion methods have some shortcomings, such as ignoring the applicability of the model space or the high complexity of model operation. In this study, the high-precision and high-temporal-resolution Global Navigation Satellite System (GNSS) PWV was used to calibrate and optimize the ERA5 PWV product of the European Center for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) with high spatial resolution to improve its accuracy, and its applicability was verified at the spatiotemporal scale. First, this study obtained accurate GNSS PWV from meteorological data from stations and used it as the true value to analyze the distribution of the ERA5 PWV in mainland China. The results showed that the ERA5 PWV showed significant spatial and temporal differences. Then, a backpropagation neural network (BPNN) fusion correction model with additional constraints was established. The correction results showed that the bias of the ERA5 PWV mainly fluctuated near 0, the correlation between the ERA5 PWV and GNSS PWV was increased to 0.99, and the positive improvement rate of the root-mean-square error (RMSE) was 95%. In the temporal scale validation, the RMSE of the ERA5 PWV decreased from 2.05 mm to 1.67 mm, an improvement of 18.54%. In the spatial scale validation, the RMSE of the four seasons decreased by 0.26–80% (spring), 0.28–70.71% (summer), 0.28–45.23% (autumn), and 0.30–40.75% (winter). Especially in the summer and plateau mountainous areas where the ERA5 PWV performance was poor, the model showed suitable stability. Finally, the fusion model was used to generate a new PWV product, which improved the accuracy of ERA5 PWV on the basis of ensuring the spatial resolution. [ABSTRACT FROM AUTHOR]

In this study, we investigate sea state estimation from spaceborne GNSS-R. Due to the complex scattering of electromagnetic waves on the rough sea surface, the neural network approach is adopted to develop an algorithm to derive significant wave height (SWH) from CYGNSS data. Eighty-nine million pieces of CYGNSS data from September to November 2020 and the co-located ECMWF data are employed to train a three-hidden-layer neural network. Ten variables are considered as the input parameters of the neural network. Without the auxiliary of the wind speed, the SWH retrieved using the trained neural network exhibits a bias and an RMSE of −0.13 and 0.59 m with respect to ECMWF data. When considering wind speed as the input, the bias and RMSE were reduced to −0.09 and 0.49 m, respectively. When the incidence angle ranges from 35 ° to 65 ° and the SNR is above 7 dB, the retrieval performance is better than that obtained using other values. The measurements derived from the "Block III" satellite offer worse results than those derived from other satellites. When the distance is considered as an input parameter, the retrieval performances for the areas near the coast are significantly improved. A soft data filter is used to synchronously improve the precision and ensure the desired sample number. The RMSEs of the retrieved SWH are reduced to 0.45 m and 0.41 m from 0.59 m and 0.49 m, and only 16.0% and 14.9% of the samples are removed. The retrieved SWH also shows a clear agreement with the co-located buoy and Jason-3 altimeter data. [ABSTRACT FROM AUTHOR]

This paper reconstructs, for the first time, the motion dynamics of an eruptive cloud formed during the catastrophic eruption of the Sheveluch volcano in November 1964 (Volcanic Explosivity Index 4+). This became possible due to the public availability of atmospheric reanalysis data from the ERA-40 archive of the European Center for Medium-Range Weather Forecasts (ECMWF) and the development of numerical modeling of volcanic ash cloud propagation. The simulation of the eruptive cloud motion process, which was carried out using the FALL3D and PUFF models, made it possible to clarify the sequence of events of this eruption (destruction of extrusive domes in the crater and the formation of an eruptive column and pyroclastic flows), which lasted only 1 h 12 min. During the eruption, the ash cloud consisted of two parts: the main eruptive cloud that rose up to 15,000 m above sea level (a.s.l.), and the co-ignimbrite cloud that formed above the moving pyroclastic flows. The ashfall in Ust-Kamchatsk (Kamchatka) first occurred out of the eruptive cloud moving at a higher speed, then out of the co-ignimbrite cloud. In Nikolskoye (Bering Island, Commander Islands), ash fell only out of the co-ignimbrite cloud. Under the turbulent diffusion, the forefront of the main eruptive cloud rose slowly in the atmosphere and reached 16,500 m a.s.l. by 04:07 UTC on November 12. Three days after the eruption began, the eruptive cloud stretched for 3000 km over the territories of the countries of Russia, Canada, the USA, Mexico, and over both the Bering Sea and the Pacific Ocean. It is assumed that the well-known long-term decrease in the solar radiation intensity in the northern latitudes from 1963–1966, which was established according to the world remote sensing data, was associated with the spread of aerosol clouds formed not only by the Agung volcano, but those formed during the 1964 Sheveluch volcano catastrophic eruption. [ABSTRACT FROM AUTHOR]

Coupled general circulation models (CGCMs), as tools of predicting climate variability, are constantly being improved due to their immense value in a host of theoretical and practical, real-world problems. Consequently, four new parameterization schemes are introduced in the First Institute of Oceanography Earth System Model (FIO-ESM), and a new climate prediction System (CPS) is built up based on modified and original FIO-ESM. Here, turbulence from the sea surface to deep ocean were fully described, and seasonal forecasts of El Niño-Southern Oscillation (ENSO) and year-to-year prediction of Pacific Decadal Oscillation (PDO) were made with both the modified and original FIO-ESM-CPS. The results illustrate that the anomaly correlation coefficient (ACC) of the Niño 3.4 index significantly increased, and the root mean square error (RMSE) significantly decreased, respectively, in the modified FIO-ESM-CPS as compared to the original. The RMSE is improved by over 20% at 4- and 5-month lead times. Over longer leads, and in the modified FIO-ESM-CPS, forecast ENSO amplitudes are far closer to observations than the original CGCM, which significantly overestimates amplitudes. PDO prediction skill is also improved in the modified FIO-ESM-CPS with ACC improving by 36% at the 4-year lead time and RMSE decreasing by 21% at the 3-year lead time. [ABSTRACT FROM AUTHOR]

In recent years, Global Navigation Satellite System (GNSS) radio occultation (RO) has become a critical observation system for global operational numerical weather prediction. Constellation Observing System for Meteorology, Ionosphere, Climate (COSMIC) 2 (COSMIC-2) has been a backbone RO mission for NOAA. NOAA also began to purchase RO data from commercial sources in 2020. To ensure the consistent quality of RO data from different sources, NOAA Center for Satellite Applications and Research (STAR) has developed capabilities to process all available RO data from different missions. This paper describes the STAR RO processing systems which convert the pseudo-range and carrier phase observations to excess phases and bending angles (BAs). We compared our COSMIC-2 data products with those processed by the University Corporation for Atmospheric Research (UCAR) COSMIC Data Analysis and Archive Center (CDAAC). We processed more than twelve thousand COSMIC-2 occultation profiles. Our results show that the excess phase difference between UCAR and STAR is within a few centimeters at high altitudes, although the difference increases towards the lower atmosphere. The BA profiles derived from the excess phase are consistent with UCAR. The mean relative BA differences at impact height from 10 to 30 km are less than 0.1% for GLObal NAvigation Satellite System (GLONASS) L2C signals and Global Positioning System (GPS) L2C and L2P signals. The standard deviations are 1.15%, 1.15%, and 1.32% for GLONASS L2C signal and for GPS L2C and L2P signals, respectively. The BA profiles agree with those derived from European Center for Medium-range Weather Forecast (ECMWF) reanalysis version 5 (ERA5). The Signal-to-Noise-Ratio (SNR) plays an essential role in the processing. The STAR BA profiles with higher L1 SNRs (L1 at 80 km) tend to yield more consistent results than those from UCAR, with a negligible difference and a smaller deviation than lower SNR profiles. Profiles with lower SNR values tend to show a more significant standard deviation towards the surface during the open-loop stage in the lower troposphere than those of higher SNR. We also found that the different COSMIC-2 clock solutions could contribute to the significant relative BA difference at high altitudes; however, it has little effect on the lower troposphere comparisons given larger BA values. [ABSTRACT FROM AUTHOR]

The interferometric imaging radar altimeter (InIRA) adopts a short baseline along with small incidence angles to acquire interferometric signals from the sea surface with high accuracy, thus the wide-swath sea surface height (SSH) and backscattering coefficient ( σ 0 ) can be obtained simultaneously. This work presents an approach to jointly retrieve the wind speed and significant wave height (SWH) for the Chinese Tiangong-2 interferometric imaging radar altimeter (TG2-InIRA). This approach utilizes a multilayer perceptron (MLP) joint retrieval model based on σ 0 and SSH data. By comparing with the European Center for Medium-Range Weather Forecasts (ECMWF) reanalysis data, the root mean square errors (RMSEs) of the retrieved wind speed and the SWH are 1.27 m/s and 0.36 m, respectively. Based on the retrieved SWH, two enhanced wind speed retrieval models are developed for high sea states and low sea states, respectively. The results show that the RMSE of the retrieved wind speed is 1.12 m/s when the SWHs < 4 m; the RMSE is 0.73 m/s when the SWHs ≥ 4 m. Similarly, two enhanced SWH retrieval models for relatively larger and relatively smaller wind speed regions are developed based on the retrieved wind speed with corresponding RMSEs of 0.19 m and 0.16 m, respectively. The comparison between the retrieved results and the buoy data shows that they are highly consistent. The results show that the additional information of SWH can be used to improve the accuracy of wind speed retrieval at small incidence angles, and also the additional information of wind speed can be used to improve the SWH retrieval. The stronger the correlation between wind speed and SWH, the greater the improvement of the retrieved results. The proposed method can achieve joint retrieval of wind speed and SWH accurately, which complements the existing wind speed and SWH retrieval methods for InIRA. [ABSTRACT FROM AUTHOR]

For more than a decade, the European Centre for Medium-Range Weather Forecasts (ECMWF) has used in-situ observations of 2 m temperature and 2 m relative humidity to operationally constrain the temporal evolution of model soil moisture. These observations are not available everywhere and they are indirectly linked to the state of the surface, so under various circumstances, such as weak radiative forcing or strong advection, they cannot be used as a proxy for soil moisture reinitialization in numerical weather prediction. Recently, the ECMWF soil moisture analysis has been updated to be able to account for the information provided by microwave brightness temperatures from the Soil Moisture and Ocean Salinity (SMOS) mission of the European Space Agency (ESA). This is the first time that ECMWF uses direct information of the soil emission from passive microwave data to globally adjust the estimation of soil moisture by a land-surface model. This paper presents a novel version of the ECMWF Extended Kalman Filter soil moisture analysis to account for remotely sensed passive microwave data. It also discusses the advantages of assimilating direct satellite radiances compared to current soil moisture products, with a view to an operational implementation. A simple assimilation case study at global scale highlights the potential benefits and obstacles of using this new type of information in a global coupled land-atmospheric model. [ABSTRACT FROM AUTHOR]

Unlike most altimetric missions, CryoSat-2 is not equipped with an onboard microwave radiometer (MWR) to provide wet tropospheric correction (WTC) to radar altimeter measurements, thus, relying on a model-based one provided by the European Center for Medium-range Weather Forecasts (ECMWF). In the ambit of ESA funded project CP4O, an improved WTC for CryoSat-2 data over ocean is under development, based on a data combination algorithm (DComb) through objective analysis of WTC values derived from all existing global-scale data types. The scope of this study is the analysis and inter-calibration of the large dataset of total column water vapor (TCWV) products from scanning MWR aboard Remote Sensing (RS) missions for use in the WTC computation for CryoSat-2. The main issues regarding the computation of the WTC from all TCWV products are discussed. The analysis of the orbital parameters of CryoSat-2 and all other considered RS missions, their sensor characteristics and inter-calibration is presented, providing an insight into the expected impact of these datasets on the WTC estimation. The most suitable approach for calculating the WTC from TCWV is investigated. For this type of application, after calibration with respect to an appropriate reference, two approaches were found to give very similar results, with root mean square differences of 2 mm. [ABSTRACT FROM AUTHOR]

Aerosols in the atmosphere play an essential role in the radiative transfer process due to their scattering, absorption, and emission. Moreover, they interrupt the retrieval of atmospheric properties from ground-based and satellite remote sensing. Thus, accurate aerosol information needs to be obtained. Herein, we developed an optimal-estimation-based aerosol optical depth (AOD) retrieval algorithm using the hyperspectral infrared downwelling emitted radiance of the Atmospheric Emitted Radiance Interferometer (AERI). The proposed algorithm is based on the phenomena that the thermal infrared radiance measured by a ground-based remote sensor is sensitive to the thermodynamic profile and degree of the turbid aerosol in the atmosphere. To assess the performance of algorithm, AERI observations, measured throughout the day on 21 October 2010 at Anmyeon, South Korea, were used. The derived thermodynamic profiles and AODs were compared with those of the European center for a reanalysis of medium-range weather forecasts version 5 and global atmosphere watch precision-filter radiometer (GAW-PFR), respectively. The radiances simulated with aerosol information were more suitable for the AERI-observed radiance than those without aerosol (i.e., clear sky). The temporal variation trend of the retrieved AOD matched that of GAW-PFR well, although small discrepancies were present at high aerosol concentrations. This provides a potential possibility for the retrieval of nighttime AOD. [ABSTRACT FROM AUTHOR]

In this study, we present a nationwide machine learning model for hourly PM 2.5 estimation for the continental United States (US) using high temporal resolution Geostationary Operational Environmental Satellites (GOES-16) Aerosol Optical Depth (AOD) data, meteorological variables from the European Center for Medium Range Weather Forecasting (ECMWF) and ancillary data collected between May 2017 and December 2020. A model sensitivity analysis was conducted on predictor variables to determine the optimal model. It turns out that GOES16 AOD, variables from ECMWF, and ancillary data are effective variables in PM 2.5 estimation and historical reconstruction, which achieves an average mean absolute error (MAE) of 3.0 μ g/m 3 , and a root mean square error (RMSE) of 5.8 μ g/m 3 . This study also found that the model performance as well as the site measured PM 2.5 concentrations demonstrate strong spatial and temporal patterns. Specifically, in the temporal scale, the model performed best between 8:00 p.m. and 11:00 p.m. (UTC TIME) and had the highest coefficient of determination (R 2 ) in Autumn and the lowest MAE and RMSE in Spring. In the spatial scale, the analysis results based on ancillary data show that the R 2 scores correlate positively with the mean measured PM 2.5 concentration at monitoring sites. Mean measured PM 2.5 concentrations are positively correlated with population density and negatively correlated with elevation. Water, forests, and wetlands are associated with low PM 2.5 concentrations, whereas developed, cultivated crops, shrubs, and grass are associated with high PM 2.5 concentrations. In addition, the reconstructed PM 2.5 surfaces serve as an important data source for pollution event tracking and PM 2.5 analysis. For this purpose, from May 2017 to December 2020, hourly PM 2.5 estimates were made for 10 km by 10 km and the PM 2.5 estimates from August through November 2020 during the period of California Santa Clara Unite (SCU) Lightning Complex fires are presented. Based on the quantitative and visualization results, this study reveals that a number of large wildfires in California had a profound impact on the value and spatial-temporal distributions of PM 2.5 concentrations. [ABSTRACT FROM AUTHOR]

From August to October 2020, a serious wildfire occurred in California, USA, which produced a large number of particulate matter and harmful gases, resulting in huge economic losses and environmental pollution. Particulate matter delays the GNSS signal, which affects the like precipitable water vapor (LPWV) derived by the GNSS non-hydrostatic delay. Most of the information of GNSS-derived LPWV is caused by water vapor, and a small part of the information is caused by particulate matter. A new method based on the difference (ΔPWV) between the PWV of virtual radiosonde stations network and GNSS-derived LPWV is proposed to detect the changes of particulate matter in the atmosphere during the 2020 California wildfires. There are few radiosonde stations in the experimental area and they are far away from the GNSS station. In order to solve this problem, we propose to use the multilayer perceptron (MLP) neural network method to establish the virtual radiosonde network in the experimental area. The PWV derived by the fifth-generation European center for medium-range weather forecasts reanalysis model (PWVERA5) is used as the input data of machine learning. The PWV derived by radiosonde data (PWVRAD) is used as the training target data of machine learning. The ΔPWV is obtained based on PWV derived by the virtual radiosonde station network and GNSS in the experimental area. In order to further reduce the influence of noise and other factors on ΔPWV, this paper attempts to decompose ΔPWV time series by using the singular spectrum analysis method, and obtain its principal components, subsequently, analyzing the relationship between the principal components of ΔPWV with particulate matter. The results indicate that the accuracy of PWV predicted by the virtual radiosonde network is significantly better than the fifth-generation European center for the medium-range weather forecast reanalysis model, and the change trend of ΔPWV is basically consistent with the change law of particulate matter in which the value of ΔPWV in the case of fire is significantly higher than that before and after the fire. The mean of correlation coefficients between ΔPWV and PM10 at each GNSS station before, during and after wildfires are 0.068, 0.397 and 0.065, respectively, which show the evident enhancement of the correlation between ΔPWV and particulate matter during wildfires. It is concluded that because of the high sensitiveness of ΔPWV to the change of particulate matter, the GNSS technique can be used as an effective new approach to detect the change of particulate matter and, then, to detect wildfires effectively. [ABSTRACT FROM AUTHOR]

Water vapor is one of the most important greenhouse gases in the world. There are many techniques that can measure water vapor directly or remotely. In this work, we first study the Global Positioning System (GPS)- and the Global Navigation Satellite System (GLONASS)-derived Zenith Wet Delay (ZWD) time series based on 11 years of the second reprocessing campaign of International Global Navigation Satellite Systems (GNSS) Service (IGS) using 320 globally distributed stations. The amount of measurement, the local environment, and the antenna radome are shown to be the main factors that affect the GNSS ZWDs and the corresponding a posteriori formal errors. Furthermore, antenna radome is able to effectively reduce the systematic bias of ZWDs and a posteriori formal errors between the GPS- and GLONASS-based solutions. With the development of the GLONASS, the ZWD differences between the GPS- and the GLONASS-based solutions have gradually decreased to sub-mm-level after GLONASS was fully operated. As the GPS-based Precipitable Water Vapor (PWV) is usually used as the reference to evaluate the other PWV products, the PWV consistency among several common techniques is evaluated, including GNSSs, spaceborne sensors, and numerical products from the European Center for Medium-Range Weather Forecasts (ECMWF). As an example of the results from a detailed comparison analysis, the long-term global analysis shows that the PWV obtained from the GNSS and the ECMWF have great intra-agreements. Based on the global distribution of the magnitude of the PWV and the PWV drift, most of the techniques showed superior agreement and proved their ability to do climate research. With a detailed study performed for the ZWDs and PWV on a long-term global scale, this contribution provides a useful supplement for future research on the GNSS ZWD and PWV. [ABSTRACT FROM AUTHOR]

Solar radiation is the most important source of energy on the Earth. The Gobi area in the eastern Xinjiang region, due to its geographic location and climate characteristics, has abundant solar energy resources. In order to provide detailed scientific data supporting solar energy development in this area, we used ground-based data to evaluate the applicability of the five reanalysis data sources: the Clouds and the Earth's Radiant Energy System (CERES), the European Center for Medium-Range Weather Forecasts Reanalysis version 5 (ERA5), the Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA2), and the Japanese 55-year Reanalysis (JRA-55). Our results indicated that the CERES data show underestimated short-wave radiation and overestimated long-wave radiation. The correlation coefficients (r) between the ERA5 dataset and the net long-wave and short-wave radiation in observation were 0.92 and 0.91, respectively, and the r between the MERRA2 dataset and the net long-wave and short-wave radiation in observation were both 0.88. The JRA-55 dataset overestimated the long-wave radiation flux and underestimated the short-wave radiation flux. The clearness index (kt) of all datasets was poor during autumn and winter, the ERA5 estimates were cloudy when the actual condition was sunny, while the JRA-55 estimates were sunny when the actual condition was cloudy. Overall, the radiation flux in the ERA5 dataset had the highest applicability in the Gobi region. [ABSTRACT FROM AUTHOR]

This study applies the Advanced Radiative Transfer Modeling System (ARMS), which was developed to accelerate the uses of Fengyun satellite data in weather, climate, and environmental applications in China, to characterize the biases of seven infrared (IR) bands of the Advanced Geosynchronous Radiation Imager (AGRI) onboard the Chinese geostationary meteorological satellite, Fengyun–4A. The AGRI data are quality controlled to eliminate the observations affected by clouds and contaminated by stray lights during the mid–night from 1600 to 1800 UTC during spring and autumn. The mean biases, computed from AGRI IR observations and ARMS simulations from the National Center for Environmental Prediction (NCEP) Final analysis data (FNL) as input, are within −0.7–1.1 K (0.12–0.75 K) for all seven IR bands over the oceans (land) under clear–sky conditions. The biases show seasonal variation in spatial distributions at bands 11–13, as well as a strong dependence on scene temperatures at bands 8–14 and on satellite zenith angles at absorption bands 9, 10, and 14. The discrepancies between biases estimated using FNL and the European Center for Medium–Range Weather Forecasts Reanalysis–5 (ERA5) are also discussed. The biases from water vapor absorption bands 9 and 10, estimated using ERA5 over ocean, are smaller than those from FNL. Such discrepancies arise from the fact that the FNL data are colder (wetter) than the ERA5 in the middle troposphere (upper–troposphere). [ABSTRACT FROM AUTHOR]

Tropospheric delay is a major error source in the Global Navigation Satellite System (GNSS), and the weighted mean temperature (Tm) is a key parameter in precipitable water vapor (PWV) retrieval. Although reanalysis products like the National Centers for Environmental Prediction (NCEP) and the European Center for Medium-Range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERA-Interim) data have been used to calculate and model the tropospheric delay, Tm, and PWV, the limitations of the temporal and spatial resolutions of the reanalysis data have affected their performance. The release of the fifth-generation accurate global atmospheric reanalysis (ERA5) and the second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) provide the opportunity to overcome these limitations. The performances of the zenith tropospheric delay (ZTD), zenith wet delay (ZWD), Tm, and zenith hydrostatic delay (ZHD) of ERA5 and MERRA-2 data from 2016 to 2017 were evaluated in this work using GNSS ZTD and radiosonde data over the globe. Taking GNSS ZTD as a reference, the ZTD calculated from MERRA-2 and ERA5 pressure-level data were evaluated in temporal and spatial scales, with an annual mean bias and root mean square (RMS) of 2.3 and 10.9 mm for ERA5 and 4.5 and 13.1 mm for MERRA-2, respectively. Compared to radiosonde data, the ZHD, ZWD, and Tm derived from ERA5 and MERRA-2 data were also evaluated on temporal and spatial scales, with annual mean bias values of 1.1, 1.7 mm, and 0.14 K for ERA5 and 0.5, 4.8 mm, and –0.08 K for MERRA-2, respectively. Meanwhile, the annual mean RMS was 4.5, 10.5 mm, and 1.03 K for ERA5 and 4.4, 13.6 mm, and 1.17 K for MERRA-2, respectively. Tropospheric parameters derived from MERRA-2 and ERA5, with improved temporal and spatial resolutions, can provide a reference for GNSS positioning and PWV retrieval. [ABSTRACT FROM AUTHOR]

The number of global precipitation datasets (PPs) is on the rise and they are commonly used for hydrological applications. A comprehensive evaluation on their performance in hydrological modeling is required to improve their performance. This study comprehensively evaluates the performance of eight widely used PPs in hydrological modeling by comparing with gauge-observed precipitation for a large number of catchments. These PPs include the Global Precipitation Climatology Centre (GPCC), Climate Hazards Group Infrared Precipitation with Station dataset (CHIRPS) V2.0, Climate Prediction Center Morphing Gauge Blended dataset (CMORPH BLD), Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks Climate Data Record (PERSIANN CDR), Tropical Rainfall Measuring Mission multi-satellite Precipitation Analysis 3B42RT (TMPA 3B42RT), Multi-Source Weighted-Ensemble Precipitation (MSWEP V2.0), European Center for Medium-range Weather Forecast Reanalysis 5 (ERA5) and WATCH Forcing Data methodology applied to ERA-Interim Data (WFDEI). Specifically, the evaluation is conducted over 1382 catchments in China, Europe and North America for the 1998-2015 period at a daily temporal scale. The reliabilities of PPs in hydrological modeling are evaluated with a calibrated hydrological model using rain gauge observations. The effectiveness of PPs-specific calibration and bias correction in hydrological modeling performances are also investigated for all PPs. The results show that: (1) compared with the rain gauge observations, GPCC provides the best performance overall, followed by MSWEP V2.0; (2) among the eight PPs, the ones incorporating daily gauge data (MSWEP V2.0 and CMORPH BLD) provide superior hydrological performance, followed by those incorporating 5-day (CHIRPS V2.0) and monthly (TMPA 3B42RT, WFDEI, and PERSIANN CDR) gauge data. MSWEP V2.0 and CMORPH BLD perform better than GPCC, underscoring the effectiveness of merging multiple satellite and reanalysis datasets; (3) regionally, all PPs exhibit better performances in temperate regions than in arid or topographically complex mountainous regions; and (4) PPs-specific calibration and bias correction both can improve the streamflow simulations for all eight PPs in terms of the Nash and Sutcliffe efficiency and the absolute bias. This study provides insights on the reliabilities of PPs in hydrological modeling and the approaches to improve their performance, which is expected to provide a reference for the applications of global precipitation datasets. [ABSTRACT FROM AUTHOR]

Inland open water bodies often pose a systematic error source in the passive remote sensing retrievals of soil moisture. Water temperature is a necessary variable used to compute water emissions that is required to be subtracted from satellite observation to yield actual emissions from the land portion, which in turn generates accurate soil moisture retrievals. Therefore, overestimation of soil moisture can often be corrected using concurrent water temperature data in the overall mitigation procedure. In recent years, several data sets of lake water temperature have become available, but their specifications and accuracy have rarely been investigated in the context of passive soil moisture remote sensing on a global scale. For this reason, three lake temperature products were evaluated against in-situ measurements from 2007 to 2011. The data sets include the lake surface water temperature (LSWT) from Global Observatory of Lake Responses to Environmental Change (GloboLakes), the Copernicus Global Land Operations Cryosphere and Water (C-GLOPS), as well as the lake mix-layer temperature (LMLT) from the European Centers for Medium-Range Weather Forecast (ECMWF) ERA5 Land Reanalysis. GloboLakes, C-GLOPS, and ERA5 Land have overall comparable performance with Pearson correlations (R) of 0.87, 0.92 and 0.88 in comparison with in-situ measurements. LSWT products exhibit negative median biases of −0.27 K (GloboLakes) and −0.31 K (C-GLOPS), whereas the median bias of LMLT is 1.56 K. When mapped from their respective native resolutions to a common 9 km Equal-Area Scalable Earth (EASE) Grid 2.0 projection, similar relative performance was observed. LMLT and LSWT data are closer in performance over the 9 km grid cells that exhibit a small range of lake cover fractions (0.05–0.5). Despite comparable relative performance, ERA5 Land shows great advantages in spatial coverage and temporal resolution. In summary, an integrated evaluation on data accuracy, long-term availability, global coverage, temporal resolution, and regular forward processing with modest data latency led us to conclude that LMLT from the ERA5 Land Reanalysis product represents the most optimal path for use in the development of a long-term soil moisture product. [ABSTRACT FROM AUTHOR]

The radio occultation technique provides stable atmospheric measurements that can work as a benchmark for calibrating and validating satellite-sounding data. Launched on 25 June 2019, the Constellation Observing System for Meteorology, Ionosphere, and Climate 2 and Formosa Satellite Mission 7 (COSMIC-2/FORMOSAT-7) are expected to produce about 5000 high-quality RO observations daily over the tropics and subtropics. COSMIC-2 constellation consists of 6 Low Earth Orbit (LEO) satellites in 24° inclination orbits at 720 km altitude and distributed mainly between 45°N to 45°S. The COSMIC-2 observations have uniform temporal coverage between 30°N to 30°S. This paper presents an independent inversion algorithm to invert COSMIC-2 geometry and phase data to bending angle and refractivity. We also investigate the quality of Global Navigation Satellite System (GNSS) and LEO position vectors derived from the UCAR COSMIC Data Analysis and Archive Center (CDAAC). The GNSS and LEO position vectors are stable with LEO position variations < 1.4 mm/s. The signal-to-noise ratio (SNR) on the L1 band ranges from 300–2600 v/v with a mean of 1600 v/v. The inversion algorithm developed at NOAA Center for Satellite Applications and Research (STAR) uses the Full Spectrum Inversion (FSI) method to invert COSMIC-2 geometry and phase data to bending angle and refractivity profiles. The STAR COSMIC-2 bending angle and refractivity profiles are compared with in situ radiosonde, the current COSMIC-2 products derived from CDAAC, and the collocated European Center for Medium-Range Weather Forecasts (ECMWF) climate reanalysis data ERA5. The mean bias at 8–40 km altitude among the UCAR, ERA5, and STAR is <0.1% for both bending and refractivity, with a standard deviation in the range of 1.4–2.3 and 0.9–1.1% for bending angles refractivity, respectively. In the lowest 2 km, the RO bias relative to ERA-5 shows a strong latitudinal and SNR dependence. [ABSTRACT FROM AUTHOR]

The meteorological reanalysis data has been widely applied to derive zenith tropospheric delay (ZTD) with a high spatial and temporal resolution. With the rapid development of artificial intelligence, machine learning also begins as a high-efficiency tool to be employed in modeling and predicting ZTD. In this paper, we develop three new regional ZTD models based on the least squares support vector machine (LSSVM), using both the International GNSS Service (IGS)-ZTD products and European Centre for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) data over Europe throughout 2018. Among them, the ERA5 data is extended to ERA5S-ZTD and ERA5P-ZTD as the background data by the model method and integral method, respectively. Depending on different background data, three schemes are designed to construct ZTD models based on the LSSVM algorithm, including the without background data, with the ERA5S-ZTD, and with the ERA5P-ZTD. To investigate the advantage and feasibility of the proposed ZTD models, we evaluate the accuracy of two background data and three schemes by segmental comparison with the IGS-ZTD of 85 IGS stations in Europe. The results show that the overall average Root Mean Square Errors (RMSE) value of all sites is 30.1 mm for the ERA5S-ZTD, and 10.7 mm for the ERA5P-ZTD. The overall average RMSE is 25.8 mm, 22.9 mm, and 9 mm for the three schemes, respectively. Moreover, the overall improvement rate is 19.1% and 1.6% for the ZTD model with ERA5S-ZTD and ERA5P-ZTD, respectively. In order to explore the reason of the lower improvement for the ZTD model with ERA5P-ZTD, the loop verification is performed by estimating the ZTD values of each available IGS station. In actuality, the monthly improvement rate of estimated ZTD is positive for most stations, and the biggest improvement rate can even reach about 40%. The negative rate mainly comes from specific stations, these stations are located on the edge of the region, near the coast, as well as the lower similarity between the individual verified station and training stations. [ABSTRACT FROM AUTHOR]

The present study aims to quantify the potential of hyperspectral thermal infrared sounders such as the Infrared Atmospheric Sounding Interferometer (IASI) and the future IASI next generation (IASI-NG) for retrieving the ice cloud layer altitude and thickness together with the ice water path. We employed the radiative transfer model Radiative Transfer for TOVS (RTTOV) to simulate cloudy radiances using parameterized ice cloud optical properties. The radiances have been computed from an ice cloud profile database coming from global operational short-range forecasts at the European Center for Medium-range Weather Forecasts (ECMWF) which encloses the normal conditions, typical variability, and extremes of the atmospheric properties over one year (Eresmaa and McNally (2014)). We performed an information content analysis based on Shannon's formalism to determine the amount and spectral distribution of the information about ice cloud properties. Based on this analysis, a retrieval algorithm has been developed and tested on the profile database. We considered the signal-to-noise ratio of each specific instrument and the non-retrieved atmospheric and surface parameter errors. This study brings evidence that the observing system provides information on the ice water path ( I W P ) as well as on the layer altitude and thickness with a convergence rate up to 95% and expected errors that decrease with cloud opacity until the signal saturation is reached (satisfying retrievals are achieved for clouds whose I W P is between about 1 and 300 g / m 2 ). [ABSTRACT FROM AUTHOR]

Satellite-based precipitation estimates with high quality and spatial-temporal resolutions play a vital role in forcing global or regional meteorological, hydrological, and agricultural models, which are especially useful over large poorly gauged regions. In this study, we apply various statistical indicators to comprehensively analyze the quality and compare the performance of five newly released satellite and reanalysis precipitation products against China Merged Precipitation Analysis (CMPA) rain gauge data, respectively, with 0.1° × 0.1° spatial resolution and two temporal scales (daily and hourly) over southern China from June to August in 2019. These include Precipitation Estimates from Remotely Sensed Information using Artificial Neural Networks Cloud Classification System (PERSIANN-CCS), European Center for Medium-Range Weather Forecasts Reanalysis v5 (ERA5-Land), Fengyun-4 (FY-4A), Global Satellite Mapping of Precipitation (GSMaP), and Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG). Results indicate that: (1) all five products overestimate the accumulated rainfall in the summer, with FY-4A being the most severe; additionally, FY-4A cannot capture the spatial and temporal distribution characteristics of precipitation over southern China. (2) IMERG and GSMaP perform better than the other three datasets at both daily and hourly scales; IMERG correlates slightly better than GSMaP against CMPA data, while it performs worse than GSMaP in terms of probability of detection (POD). (3) ERA5-Land performs better than PERSIANN-CCS and FY-4A at daily scale but shows the worst correlation coefficient (CC), false alarm ratio (FAR), and equitable threat score (ETS) of all precipitation products at hourly scale. (4) The rankings of overall performance on precipitation estimations for this region are IMERG, GSMaP, ERA5-Land, PERSIANN-CCS, and FY-4A at daily scale; and IMERG, GSMaP, PERSIANN-CCS, FY-4A, and ERA5-Land at hourly scale. These findings will provide valuable feedback for improving the current satellite-based precipitation retrieval algorithms and also provide preliminary references for flood forecasting and natural disaster early warning. [ABSTRACT FROM AUTHOR]

Recent global warming has been accompanied by high water temperatures (HWTs) in coastal areas of Korea, resulting in huge economic losses in the marine fishery industry due to disease outbreaks in aquaculture. To mitigate these losses, it is necessary to predict such outbreaks to prevent or respond to them as early as possible. In the present study, we propose an HWT prediction method that applies sea surface temperatures (SSTs) and deep-learning technology in a long short-term memory (LSTM) model based on a recurrent neural network (RNN). The LSTM model is used to predict time series data for the target areas, including the coastal area from Goheung to Yeosu, Jeollanam-do, Korea, which has experienced frequent HWT occurrences in recent years. To evaluate the performance of the SST prediction model, we compared and analyzed the results of an existing SST prediction model for the SST data, and additional external meteorological data. The proposed model outperformed the existing model in predicting SSTs and HWTs. Although the performance of the proposed model decreased as the prediction interval increased, it consistently showed better performance than the European Center for Medium-Range Weather Forecast (ECMWF) prediction model. Therefore, the method proposed in this study may be applied to prevent future damage to the aquaculture industry. [ABSTRACT FROM AUTHOR]

The dynamical tropopause is the interface between the stratosphere and the troposphere, whose variation gives indication of weather and climate changes. In the past, the dynamical tropopause height determination mainly depends on analysis and diagnose methods. While, due to the high computational cost, it is difficult to obtain tropopause structures with high spatiotemporal resolution in real time by these methods. To solve this problem, the statistical method is used to establish the dynamical tropopause pressure retrieval model based on Fengyun-4A geostationary meteorological satellite observations. Four regression schemes including random forest (RF) regression are evaluated. By comparison with GEOS-5 (the Goddard Earth Observing System Model of version 5) and ERA-Interim (European Center for Medium-Range Weather Forecasts Reanalysis-Interim) reanalysis, it is found that among the four schemes, the RF-based retrieval model is most accurate and reliable (RMSEs (root mean square errors) are 25.99 hPa and 43.05 hPa, respectively, as compared to GEOS-5 and ERA-Interim reanalysis). A series of sensitivity experiments are performed to investigate the contributions of the predictors in the RF-based model. Results suggest that 6.25 μm channel information representing the distributions of the potential vorticity and water vapor in upper troposphere has the greatest contribution, while 10.8 and 12 μm channels information have relatively weak influences. Therefore, a simplified model without involving a brightness temperature of 10.8 and 12 μm can be adopted to improve the calculation efficiency. [ABSTRACT FROM AUTHOR]

Recently, surface diffuse solar radiation (Rdif) has been attracting a growing interest in view of its function in improving plant productivity, thus promoting global carbon uptake, and its impacts on solar energy utilization. To date, very few radiation products provide estimates of Rdif, and systematic validation and evaluation are even more scare. In this study, Rdif estimates from Reanalysis Fifth Generation (ERA5) of European Center for Medium-Range Weather Forecasts and satellite-based retrieval (called JiEA) are evaluated over East Asia using ground measurements at 39 stations from World Radiation Data Center (WRDC) and China Meteorological Administration (CMA). The results show that JiEA agrees well with measurements, while ERA5 underestimates Rdif significantly. Both datasets perform better at monthly mean scale than at daily mean and hourly scale. The mean bias error and root-mean-square error of daily mean estimates are −1.21 W/m2 and 20.06 W/m2 for JiEA and −17.18 W/m2 and 32.42 W/m2 for ERA5, respectively. Regardless of over- or underestimation, correlations of estimated time series of ERA5 and JiEA show high similarity. JiEA reveals a slight decreasing trend at regional scale, but ERA5 shows no significant trend, and neither of them reproduces temporal variability of ground measurements. Data accuracy of ERA5 is more robust than JiEA in time but less in space. Latitudinal dependency is noted for ERA5 while not for JiEA. In addition, spatial distributions of Rdif from ERA5 and JiEA show pronounced discrepancy. Neglect of adjacency effects caused by horizontal photon transport is the main cause for Rdif underestimation of ERA5. Spatial analysis calls for improvements to the representation of clouds, aerosols and water vapor for reproducing fine spatial distribution and seasonal variations of Rdif. [ABSTRACT FROM AUTHOR]

The latest reanalysis of the European Center for Medium-Range Weather Forecasts (ECMWF), ERA5, can provide atmospheric data for calculating Zenith Tropospheric Delay (ZTD) with hourly temporal resolution, which is a key factor in Global Navigation Satellite System (GNSS) high-precision application. This paper is aimed at evaluating the performance of ZTD derived from ERA5 reanalysis data over China using 219 GNSS stations of the Crustal Movement Observation Network of China (CMONOC) covering the period from 2015 to 2016. The site-specific hourly ZTD at these stations is obtained by integration method and Saastamoinen model method on ERA5 pressure-level and surface-level reanalysis data with the temporal resolution of 1 h and the spatial resolution of 0.25° × 0.25°. Firstly, the atmospheric temperature and pressure that derived from ERA5 are compared with temperature and pressure obtained from meteorological sensors available at 193 GNSS stations. The biases are 2.31 °C and 1.26 mbar implying the accuracy and feasibility of ERA5 pressure and temperature for calculating ZTD over China. Secondly, the performance of ERA5 ZTD is systematically evaluated using ZTD from 219 GNSS sites. The average bias and Root Mean Square (RMS) of ERA5 pressure-level ZTD at all test stations in integration method are approximately 2.97 mm and 11.49 mm respectively, and those of ERA5 surface-level ZTD in model method are 7.97 mm, 39.25 mm, which indicates that ERA5 pressure-level ZTD has a higher accuracy over China. Further analysis indicates that the accuracies of ZTD derived from ERA5 pressure-level and surface-level data are approximately 13.8% and 10.9% higher than those from of ERA-Interim pressure-level and surface-level data. Moreover, ERA5 is able to accurately capture the short-term (hourly) variation of ZTD, which further indicates the better performance of ERA5. Thirdly, the temporal and spatial variation characteristics of ERA5 ZTD accuracy are further analyzed over China. The results show that the ZTD in the southern region has the lower accuracy compared with that in the northern region over China due to the influence of latitude and altitude. Furthermore, it is found that the ERA5 ZTD over China has obvious seasonality, with higher accuracy in winter and lower accuracy in summer. [ABSTRACT FROM AUTHOR]

Southeast China, a non-core region influenced by the El Niño–Southern Oscillation (ENSO), has been seldom investigated before. However, the occurrence of ENSO will affect the redistribution of precipitation and the temperature (T) spatial pattern on a global scale. This condition will further lead to flood or drought disasters in Southeast China. Therefore, the method of monitoring the occurrence of ENSO is important and is the focus of this paper. The spatiotemporal characteristics of precipitable water vapor (PWV) and T are first analyzed during ENSO using the empirical orthogonal function (EOF). The results showed that a high correlation spatiotemporal consistency exist between PWV and T. The response thresholds of PWV and T to ENSO are determined by moving the window correlation analysis (MWCA). If the sea surface temperature anomaly (SSTA) at the Niño 3.4 region exceeded the ranges of (−1.17°C, 1.04°C) and (−1.15°C, 1.09°C), it could cause the anomalous change of PWV and T in Southeast China. Multichannel singular spectral analysis (MSSA) is introduced to analyze the multi-type signals (tendency, period, and anomaly) of PWV and T over the period of 1979–2017. The results showed that the annual abnormal signal and envelope line fluctuation of PWV and T agreed well in most cases with the change in SSTA. Therefore, a standard PWV and T index (SPTI) is proposed on the basis of the results to monitor ENSO events. The PWV and T data derived from the grid-based European Center for Medium-Range Weather Forecasting (ECMWF) reanalysis products and GNSS/RS stations in Southeast China were used to validate the performance of the proposed SPTI. Experimental results revealed that the time series of average SPTI calculated in Southeast China corresponded well to that of SSTA with a correlation coefficient of 0.66 over the period of 1979–2017. The PWV values derived from the Global Navigation Satellite System (GNSS) and radiosonde data at two specific stations (WUHN and 45004) were also used to calculate the SPTI. The results showed that the correlation coefficients between SPTI and SSTA were 0.73 and 0.71, respectively. Such results indicate the capacity of the proposed SPTI to monitor the ENSO in Southeast China. [ABSTRACT FROM AUTHOR]

Launched on 5 June 2019, the BuFeng-1 A/B twin satellites were part of the first Chinese global navigation satellite system reflectometry (GNSS-R) satellite mission. In this letter, a brief introduction of the BF-1 mission and its preliminary results of sea surface wind retrieval are presented. Empirical fully developed sea (FDS) geophysical model functions (GMFs) relating the normalized bistatic radar cross-section to the sea surface wind speed are proposed for the BF-1 GNSS-R instruments. The FDS GMFs are derived from the collocated BF-1 observations, the European Center for Medium-Range Weather Forecasts (ECMWF) reanalysis data, and the advanced scatterometer (ASCAT) satellite observations. The preliminary tests reveal that the root-mean-square error (RMSE) between the derived wind speed and the reanalysis is 2.63 m/s for wind speeds in the range of 0.5–40.5 m/s. Further comparisons with the ASCAT observations and mooring buoys show that the RMSEs are 2.04 m/s and 1.77 m/s, respectively, at low-to-moderate wind speeds. This study demonstrates the effectiveness of BF-1 and provides a basis for the future GMF development of the BF-1 A/B mission. [ABSTRACT FROM AUTHOR]

Some seasonal natural floods can be attributed to typhoons that bring a large amount of atmospheric water vapor, and variations in atmospheric water vapor can be reflected in the precipitable water vapor (PWV). Therefore, monitoring typhoons based on the anomalous variations of the PWV is the focus of this paper. The anomalous variations of ERA5(fifth-generation reanalysis dataset of the European Centre for Medium-range Weather Forecasting)-derived PWV with other atmospheric parameters related to typhoons, such as precipitation, pressure, and wind, were first analyzed during typhoon periods. After that, a typhoon-monitoring method with and without considering the typhoon's acceleration was proposed according to the time of the maximum value of the PWV during the typhoon period in this paper. Corresponding experiments based on the measured and simulated data were performed to evaluate the proposed method. The experimental measurement of Typhoon Hato revealed that the velocity of the typhoon's movement estimated by the proposed method was close to the observed value, and the maximum difference between the estimated and observed values was less than 3 km/h. A simulated experiment was also carried out in which the acceleration of the typhoon's movement was also considered. The simulated results verified the reliability and feasibility of the proposed method. The estimated velocity and acceleration of the typhoon's movement were almost equal to the true values under the cases of using different numbers of stations and selecting various typhoon locations. Such results obtained above indicate that the method proposed in this paper has a significant potential application value for typhoon monitoring. [ABSTRACT FROM AUTHOR]

Precipitation monitoring and early warning systems are required to reduce negative flood impacts. In this study, the performance of ensemble precipitation forecasts of three numerical weather prediction (NWP) models within the THORPEX interactive grand global ensemble (TIGGE) as well as the integrated multi-satellite retrievals for global precipitation measurement (GPM), namely IMERG, for precipitation estimates were evaluated in recent severe floods in Iran over the March–April 2019 period. The evaluations were conducted in three aspects: spatial distribution of precipitation, mean areal precipitation in three major basins hard hit by the floods, and the dichotomous evaluation in four precipitation thresholds (25, 50, 75, and 100 mm per day). The results showed that the United Kingdom Met Office (UKMO) model, in terms of spatial coverage and satellite estimates as well as the precipitation amount, were closer to the observations. Moreover, with regard to mean precipitation at the basin scale, UKMO and European Center for Medium-Range Weather Forecasts (ECMWF) models in the Gorganrud Basin, ECMWF in the Karkheh Basin and UKMO in the Karun Basin performed better than others in flood forecasting. The National Centers for Environmental Forecast (NCEP) model performed well at low precipitation thresholds, while at high thresholds, its performance decreased significantly. On the contrary, the accuracy of IMERG improved when the precipitation threshold increased. The UKMO had better forecasts than the other models at the 100 mm/day precipitation threshold, whereas the ECMWF had acceptable forecasts in all thresholds and was able to forecast precipitation events with a lower false alarm ratio and better detection when compared to other models. [ABSTRACT FROM AUTHOR]

Land surface temperature (LST) is a key variable in surface-atmosphere energy and water exchanges. The main goals of this study are to (i) evaluate the LST of the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim and ERA5 reanalyses over Iberian Peninsula using the Satellite Application Facility on Land Surface Analysis (LSA-SAF) product and to (ii) understand the main drivers of the LST errors in the reanalysis. Simulations with the ECMWF land-surface model in offline mode (uncoupled) were carried out over the Iberian Peninsula and compared with the reanalysis data. Several sensitivity simulations were performed in a confined domain centered in Southern Portugal to investigate potential sources of the LST errors. The Copernicus Global Land Service (CGLS) fraction of green vegetation cover (FCover) and the European Space Agency's Climate Change Initiative (ESA-CCI) Land Cover dataset were explored. We found a general underestimation of daytime LST and slightly overestimation at night-time. The results indicate that there is still room for improvement in the simulation of LST in ECMWF products. Still, ERA5 presents an overall higher quality product in relation to ERA-Interim. Our analysis suggested a relation between the large daytime cold bias and vegetation cover differences between (ERA5 and CGLS FCocver) with a correlation of −0.45. The replacement of the low and high vegetation cover by those of ESA-CCI provided an overall reduction of the large Tmax biases during summer. The increased vertical resolution of the soil at the surface, has a positive impact, but much smaller when compared with the vegetation changes. The sensitivity of the vegetation density parameter, that currently depends on the vegetation type, provided further proof for a needed revision of the vegetation in the model, as there is a reasonable correlation between this parameter and the Tmax mean errors when using the ESA-CCI vegetation cover (while the same correlation cannot be reproduced with the original model vegetation). Our results support the hypothesis that vegetation cover is one of the main drivers of the LST summertime cold bias in ERA5 over Iberian Peninsula. [ABSTRACT FROM AUTHOR]

Plastic pollution in the world's natural waters is of growing concern and currently receiving significant attention. However, remote sensing of marine plastic litter is still in the developmental stage. Most progress has been made in spectral remote sensing using visible to short-wave infrared wavelengths where optical physics applies. Thermal infrared (TIR) sensing could potentially monitor plastic water pollution but has not been studied in detail. We applied radiative transfer theory to predict TIR sensitivity to changes in the surface fraction of water covered by plastic litter and found that the temperature difference between the water surface and the surroundings controls the TIR signal. Hence, we mapped this difference for various months and times of the day using global SST (sea surface temperature) and t2m (temperature at 2 m height) hourly estimates from the European Centre for Medium-Range Weather Forecasts (ECMWF), ERA5. The maps show how SST-t2m difference varied, altering the anticipated effectivity of TIR floating plastic litter remote sensing. We selected several locations of interest to predict the effectivity of TIR sensing of the plastic surface fraction. TIR remote sensing has promising potential and is expected to be more effective in areas with a high air–sea temperature difference. [ABSTRACT FROM AUTHOR]

Satellite-based normalized difference vegetation index (NDVI) data are widely used for estimating vegetation greenness. Seasonal climate predictions of spring (April–May–June) NDVI over Eurasia are explored by applying the year-to-year increment approach. The prediction models were developed based on the coupled modes of singular value decomposition (SVD) analyses between Eurasian NDVI and climate factors. One synchronous predictor, the spring surface air temperature from the NCEP's Climate Forecast System (SAT-CFS), and three previous-season predictors (winter (December–January–February) sea-ice cover over the Barents Sea (SICBS), winter sea surface temperature over the equatorial Pacific (SSTP), and winter North Atlantic Oscillation (NAO) were chosen to develop four single-predictor schemes: the SAT-CFS scheme, SICBS scheme, SSTP scheme, and NAO scheme. Meanwhile, a statistical scheme that involves the three previous-season predictors (i.e., SICBS, SSTP, and NAO) and a hybrid scheme that includes all four predictors are also proposed. To evaluate the prediction skills of the schemes, one-year-out cross-validation and independent hindcast results are analyzed, revealing the hybrid scheme as having the best prediction skill. The results indicate that the temporal correlation coefficients at 92% of grid points over Eurasia are significant at the 5% significance level in the hybrid scheme, which is the best among all the schemes. Furthermore, spatial correlation coefficients (SCCs) of the six schemes are significant at the 1% significance level in most years during 1983–2015, with the averaged SCC of the hybrid scheme being the highest (0.60). The grid-averaged root-mean-square-error of the hybrid scheme is 0.04. By comparing the satellite-based NDVI value with the independent hindcast results during 2010–2015, it can be concluded that the hybrid scheme shows high prediction skill in terms of both the spatial pattern and the temporal variability of spring Eurasian NDVI. [ABSTRACT FROM AUTHOR]

Precise modeling of tropospheric delay and weighted mean temperature (Tm) is critical for Global Navigation Satellite System (GNSS) positioning and meteorology. However, the model data in previous models cover a limited time span, which limits the accuracy of these models. Besides, the vertical variations of tropospheric delay and Tm are not perfectly modeled in previous studies, which affects the performance of height corrections. In this study, we used the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis from 1979 to 2017 to build a new empirical model. We first carefully modeled the lapse rates of tropospheric delay and Tm. Then we considered the temporal variations by linear trends, annual, and semi-annual variations and the spatial variations by grids. This new model can provide zenith hydrostatic delay (ZHD), zenith wet delay (ZWD), and Tm worldwide with a spatial resolution of 1° × 1°. We used the ECMWF ERA-Interim data and the radiosonde data in 2018 to validate this new model in comparison with the canonical GPT2w model. The results show that the new model has higher accuracies than the GPT2w model in all parameters. Particularly, this new model largely improves the accuracy in estimating ZHD and Tm at high-altitude (relative to the grid point height) regions. [ABSTRACT FROM AUTHOR]

Snowfall data are vital in calculating the surface mass balance of the Antarctic Ice Sheet (AIS), where in-situ and satellite measurements are sparse at synoptic timescales. CloudSat data are used to construct Antarctic snowfall data at synoptic timescales to compensate for the sparseness of synoptic snowfall data on the AIS and to better understand its surface mass balance. Synoptic CloudSat snowfall data are evaluated by comparison with daily snow accumulation measurements from ten automatic weather stations (AWSs) and the fifth generation of the European Centre for Medium-Range Weather Forecasts climate reanalysis (ERA5) snowfall. Synoptic snowfall data were constructed based on the CloudSat measurements within a radius of 1.41°. The results show that reconstructed CloudSat snowfall at daily and two-day resolutions cover about 28% and 29% of the area of the AIS, respectively. Daily CloudSat snowfall and AWS snow accumulation have similar trends at all stations. While influenced by stronger winds, >73.3% of extreme snow accumulation events correspond to snowfall at eight stations. Even if the CloudSat snowfall data have not been assimilated into the ERA5 dataset, the synoptic CloudSat snowfall data are almost identical to the daily ERA5 snowfall with only small biases (average root mean square error and mean absolute error < 3.9 mm/day). Agreement among the three datasets suggests that the CloudSat data can provide reliable synoptic snowfall data in most areas of the AIS. The ERA5 dataset captures a large number of extreme snowfall events at all AWSs, with capture rates varying from 56% to 88%. There are still high uncertainties in ERA5. Nevertheless, the result suggests that ERA5 can be used to represent actual snowfall events on the AIS at synoptic timescale. [ABSTRACT FROM AUTHOR]

Haiyang-2A (HY-2A) has been working in-flight for over seven years, and the accuracy of HY-2A calibration microwave radiometer (CMR) data is extremely important for the wet troposphere delay correction (WTC) in sea surface height (SSH) determination. We present a comprehensive evaluation of the HY-2A CMR observation using the numerical weather model (NWM) for all the data available period from October 2011 to February 2018, including the WTC and the precipitable water vapor (PWV). The ERA(ECMWF Re-Analysis)-Interim products from European Centre for Medium-Range Weather Forecasts (ECMWF) are used for the validation of HY-2A WTC and PWV products. In general, a global agreement of root-mean-square (RMS) of 2.3 cm in WTC and 3.6 mm in PWV are demonstrated between HY-2A observation and ERA-Interim products. Systematic biases are revealed where before 2014 there was a positive WTC/PWV bias and after that, a negative one. Spatially, HY-2A CMR products show a larger bias in polar regions compared with mid-latitude regions and tropical regions and agree better in the Antarctic than in the Arctic with NWM. Moreover, HY-2A CMR products have larger biases in the coastal area, which are all caused by the brightness temperature (TB) contamination from land or sea ice. Temporally, the WTC/PWV biases increase from October 2011 to March 2014 with a systematic bias over 1 cm in WTC and 2 mm in PWV, and the maximum RMS values of 4.62 cm in WTC and 7.61 mm in PWV occur in August 2013, which is because of the unsuitable retrieval coefficients and systematic TB measurements biases from 37 GHz band. After April 2014, the TB bias is corrected, HY-2A CMR products agree very well with NWM from April 2014 to May 2017 with the average RMS of 1.68 cm in WTC and 2.65 mm in PWV. However, since June 2017, TB measurements from the 18.7 GHz band become unstable, which led to the huge differences between HY-2A CMR products and the NWM with an average RMS of 2.62 cm in WTC and 4.33 mm in PWV. HY-2A CMR shows high accuracy when three bands work normally and further calibration for HY-2A CMR is in urgent need. Furtherly, 137 global coastal radiosonde stations were used to validate HY-2A CMR. The validation based on radiosonde data shows the same variation trend in time of HY-2A CMR compared to the results from ECMWF, which verifies the results from ECMWF. [ABSTRACT FROM AUTHOR]

In-situ observation, climate reanalyses, and satellite remote sensing are used to study the annual cycle of turbulent latent heat flux (LHF) in the Agulhas Current system. We assess if the datasets do represent the intense exchange of moisture that occurs above the Agulhas Current and the Retroflection region, especially the new reanalyses as the former, the National Centers for Environmental Prediction Reanalysis 2 (NCEP2) and the European Centre for Medium-Range Weather Forecast (ECMWF) reanalysis second-generation reanalysis (ERA-40) have lower sea and less distinct surface temperature (SST) in the Agulhas Current system due to their low spatial resolution thus do not adequately represent the Agulhas Current LHF. We use monthly fields of LHF, SST, surface wind speed, saturated specific humidity at the sea surface (Qss), and specific humidity at 10 m (Qa). The Climate Forecast System Reanalysis (CFSR), the European Centre for Medium-Range Weather Forecast fifth generation (ERA-5), and the Modern-Era Retrospective analysis for Research and Applications version-2 (MERRA-2) are similar to the air–sea turbulent fluxes (SEAFLUX) and do represent the signature of the Agulhas Current. ERA-Interim underestimates the LHF due to lower surface wind speeds than other datasets. The observation-based National Oceanography Center Southampton (NOCS) dataset is different from all other datasets. The highest LHF of 250 W/m2 is found in the Retroflection in winter. The lowest LHF (~100 W/m2) is off Port Elizabeth in summer. East of the Agulhas Current, Qss-Qa is the main driver of the amplitude of the annual cycle of LHF, while it is the wind speed in the Retroflection and both Qss-Qa and wind speed in between. The difference in LHF between product are due to differences in Qss-Qa wind speed and resolution of datasets. [ABSTRACT FROM AUTHOR]