TOMOREF operator for assimilation of GNSS tomography wet refractivity fields in WRF DA system

Author(s):  
Natalia Hanna ◽  
Estera Trzcina ◽  
Maciej Kryza ◽  
Witold Rohm

<p>The amount of water vapor in the atmosphere is highly variable and not easy to measure. One of the methods to provide reliable information about the amount and distribution of the humidity in the troposphere is GNSS (Global Navigation Satellite Systems) tomography. The GNSS tomography uses the observations of signal delays between satellites and ground-based receivers over the field covered by a GNSS network. This method enables deriving the 3D distribution of wet refractivity at a low cost in all weather conditions, with high temporal and spatial resolution.</p><p>The first applications of the GNSS tomography data in the Weather Research and Forecasting Data Assimilation (WRF DA) system were performed by the adaptation of the GPSREF observation operator. In this study, we present a new tool, namely the TOMOREF observation operator, which consists of three parts: forward, tangent linear, and adjoint operators. As the input data in the assimilation process, the wet refractivity fields from two tomographic models (TUW, WUELS) are used. The analysis is carried out for a 2-week long period (May 29 – June 14, 2013) in Central Europe when severe weather conditions occurred, including heavy precipitation events. The data assimilation results are verified against radiosonde observations, synoptic data, and ERA5 reanalysis. Moreover, the performance of the TOMOREF and GPSREF operators is examined. For the forecasts of relative humidity (RH) at a pressure level of 300 hPa, the implementation of the TOMOREF operator vanishes the negative impact caused by the GPSREF operator. Additionally, the improvement of the root mean square error of the forecasts of RH up to 0.5% is observed. Comparing to the assimilation of Zenith Total Delay observations, the application of the tomographic data has overall a greater influence on the WRF model. Consequently, the GNSS tomography data can be valuable in operational weather forecasting.</p>

2019 ◽  
Author(s):  
Natalia Hanna ◽  
Estera Trzcina ◽  
Gregor Möller ◽  
Witold Rohm ◽  
Robert Weber

Abstract. From Global Navigation Satellite Systems (GNSS) signals, accurate and high-frequency atmospheric parameters can be determined in all-weather conditions. GNSS tomography is a novel technique that takes advantage of these parameters, especially of slant troposphere observations between GNSS receivers and satellites, traces these signals through a 3D grid of voxels and estimates by an inversion process the refractivity of the water vapour content within each voxel. In the last years, the GNSS tomography development focused on numerical methods to stabilize the solution, which has been achieved to a great extent. Currently, we are facing new challenges and possibilities in the application of GNSS tomography in numerical weather forecasting – the main research objective of this paper. In the first instance, refractivity fields were estimated using two different GNSS tomography models (TUW, WUELS), which cover the area of Central Europe during the period of 29 May–14 June 2013, when heavy precipitation events were observed. For both models, Slant Wet Delays (SWD) were calculated based on estimates of Zenith Total Delay (ZTD) and horizontal gradients, provided for 72 GNSS sites by Geodetic Observatory Pecny (GOP). In total, three sets of SWD observations were tested (set0 without compensation for hydrostatic anisotropic effects, set1 with compensation of this effect, set2 cleaned by wet delays outside the inner voxel model). The GNSS tomography outputs have been assimilated into the nested (12- and 36-km horizontal resolution) Weather Research and Forecasting (WRF) model, using its three-dimensional variational data assimilation (WRFDA 3DVar) system, in particular its radio occultation observations operator (GPSREF). As only total refractivity is assimilated in GPSREF, it was calculated as the sum of the hydrostatic part derived from the ALADIN-CZ model and the wet part from the GNSS tomography. We compared the results of the GNSS tomography data assimilation to the radiosonde (RS) observations. The validation shows the improvement in the weather forecasting of relative humidity (bias, standard deviation) and temperature (standard deviation) during heavy precipitation events. Future improvements to the assimilation method are also discussed.


2021 ◽  
Author(s):  
Natalia Hanna ◽  
Estera Trzcina ◽  
Maciej Kryza ◽  
Witold Rohm

<p>The numerical weather model starts from the initial state of the Earth's atmosphere in a given place and time. The initial state is created by blending the previous forecast runs (first-guess), together with observations from different platforms. The better the initial state, the better the forecast; hence, it is worthy to combine new observation types. The GNSS tomography technique, developed in recent years, provides a 3-D field of humidity in the troposphere. This technique shows positive results in the monitoring of severe weather events. However, to assimilate the tomographic outputs to the numerical weather model, the proper observation operator needs to be built.</p><p>This study demonstrates the TOMOREF operator dedicated to the assimilation of the GNSS tomography‐derived 3‐D fields of wet refractivity in a Weather Research and Forecasting (WRF) Data Assimilation (DA) system. The new tool has been tested based on wet refractivity fields derived during a very intense precipitation event. The results were validated using radiosonde observations, synoptic data, ERA5 reanalysis, and radar data. In the presented experiment, a positive impact of the GNSS tomography data assimilation on the forecast of relative humidity (RH) was noticed (an improvement of root‐mean‐square error up to 0.5%). Moreover, within 1 hour after assimilation, the GNSS data reduced the bias of precipitation up to 0.1 mm. Additionally, the assimilation of GNSS tomography data had more influence on the WRF model than the Zenith Total Delay (ZTD) observations, which confirms the potential of the GNSS tomography data for weather forecasting.</p>


2019 ◽  
Vol 12 (9) ◽  
pp. 4829-4848 ◽  
Author(s):  
Natalia Hanna ◽  
Estera Trzcina ◽  
Gregor Möller ◽  
Witold Rohm ◽  
Robert Weber

Abstract. From Global Navigation Satellite Systems (GNSS) signals, accurate and high-frequency atmospheric parameters can be determined in all-weather conditions. GNSS tomography is a technique that takes advantage of these parameters, especially of slant troposphere observations between GNSS receivers and satellites, traces these signals through a 3-D grid of voxels, and estimates by an inversion process the refractivity of the water vapour content within each voxel. In the last years, the GNSS tomography development focused on numerical methods to stabilize the solution, which has been achieved to a great extent. Currently, we are facing new challenges and possibilities in the application of GNSS tomography in numerical weather forecasting, the main research objective of this paper. In the first instance, refractivity fields were estimated using two different GNSS tomography models (TUW, WUELS), which cover the area of central Europe during the period of 29 May–14 June 2013, when heavy-precipitation events were observed. For both models, slant wet delays (SWDs) were calculated based on estimates of zenith total delay (ZTD) and horizontal gradients, provided for 88 GNSS sites by Geodetic Observatory Pecny (GOP). In total, three sets of SWD observations were tested (set0 without compensation for hydrostatic anisotropic effects, set1 with compensation of this effect, set2 cleaned by wet delays outside the inner voxel model), in order to assess the impact of different factors on the tomographic solution. The GNSS tomography outputs have been assimilated into the nested (12 and 36 km horizontal resolution) Weather Research and Forecasting (WRF) model, using its three-dimensional variational data assimilation (WRFDA 3D-Var) system, in particular, its radio occultation observation operator (GPSREF). As only total refractivity is assimilated in GPSREF, it was calculated as the sum of the hydrostatic part derived from the ALADIN-CZ model and the wet part from the GNSS tomography. We compared the results of the GNSS tomography data assimilation to the radiosonde (RS) observations. The validation shows the improvement in the weather forecasting of relative humidity (bias, standard deviation) and temperature (standard deviation) during heavy-precipitation events. Future improvements to the assimilation method are also discussed.


2018 ◽  
Vol 33 (2) ◽  
pp. 561-582 ◽  
Author(s):  
Kuan-Jen Lin ◽  
Shu-Chih Yang ◽  
Shuyi S. Chen

Abstract Ensemble-based data assimilation (EDA) has been used for tropical cyclone (TC) analysis and prediction with some success. However, the TC position spread determines the structure of the TC-related background error covariance and affects the performance of EDA. With an idealized experiment and a real TC case study, it is demonstrated that observations in the core region cannot be optimally assimilated when the TC position spread is large. To minimize the negative impact from large position uncertainty, a TC-centered EDA approach is implemented in the Weather Research and Forecasting (WRF) Model–local ensemble transform Kalman filter (WRF-LETKF) assimilation system. The impact of TC-centered EDA on TC analysis and prediction of Typhoon Fanapi (2010) is evaluated. Using WRF Model nested grids with 4-km grid spacing in the innermost domain, the focus is on EDA using dropsonde data from the Impact of Typhoons on the Ocean in the Pacific field campaign. The results show that the TC structure in the background mean state is improved and that unrealistically large ensemble spread can be alleviated. The characteristic horizontal scale in the background error covariance is smaller and narrower compared to those derived from the conventional EDA approach. Storm-scale corrections are improved using dropsonde data, which is more favorable for TC development. The analysis using the TC-centered EDA is in better agreement with independent observations. The improved analysis ameliorates model shock and improves the track forecast during the first 12 h and landfall at 72 h. The impact on intensity prediction is mixed with a better minimum sea level pressure and overestimated peak winds.


2009 ◽  
Vol 137 (1) ◽  
pp. 299-314 ◽  
Author(s):  
Xiang-Yu Huang ◽  
Qingnong Xiao ◽  
Dale M. Barker ◽  
Xin Zhang ◽  
John Michalakes ◽  
...  

Abstract The Weather Research and Forecasting (WRF) model–based variational data assimilation system (WRF-Var) has been extended from three- to four-dimensional variational data assimilation (WRF 4D-Var) to meet the increasing demand for improving initial model states in multiscale numerical simulations and forecasts. The initial goals of this development include operational applications and support to the research community. The formulation of WRF 4D-Var is described in this paper. WRF 4D-Var uses the WRF model as a constraint to impose a dynamic balance on the assimilation. It is shown to implicitly evolve the background error covariance and to produce the flow-dependent nature of the analysis increments. Preliminary results from real-data 4D-Var experiments in a quasi-operational setting are presented and the potential of WRF 4D-Var in research and operational applications are demonstrated. A wider distribution of the system to the research community will further develop its capabilities and to encourage testing under different weather conditions and model configurations.


2016 ◽  
Author(s):  
Tzvetan Simeonov ◽  
Dmitry Sidorov ◽  
Felix Norman Teferle ◽  
Georgi Milev ◽  
Guergana Guerova

Abstract. Global Navigation Satellite Systems (GNSS) meteorology is an established operational service providing hourly updated GNSS tropospheric products to the National Meteorologic Services (NMS) in Europe. In the last decade through the ground-based GNSS network densification and new processing strategies like Precise Point Positioning (PPP) it has become possible to obtain sub-hourly tropospheric products for monitoring severe weather events. In this work one year (January–December 2013) of sub-hourly GNSS tropospheric products (Zenith Total Delay) are computed using the PPP strategy for seven stations in Bulgaria. In order to take advantage of the sub-hourly GNSS data to derive Integrated Water Vapour (IWV) surface pressure and temperature with similar temporal resolution is required. As the surface observations are on 3 hourly basis the first step is to compare the surface pressure and temperature from numerical weather prediction model Weather Forecasting and Research (WRF) with observations at three synoptic stations in Bulgaria. The mean difference between the two data-sets for 1) surface pressure is less than 0.5 hPa and the correlation is over 0.989, 2) temperature the largest mean difference is 1.1 °C and the correlation coefficient is over 0.957 and 3) IWV mean difference is in range 0.1–1.1 mm. The evaluation of WRF on annual bases shows IWV underestimation between 0.5 and 1.5 mm at five stations and overestimation at Varna and Rozhen. Varna and Rozhen have also much smaller correlation 0.9 and 0.76. The study of the monthly IWV variation shows that at those locations the GNSS IWV has unexpected drop in April and March respectively. The reason for this drop is likely problems with station raw data. At the remaining 5 stations a very good agreement between GNSS and WRF is observed with high correlation during the cold part of 2013 i.e. March, October and December (0.95) and low correlation during the warm part of 2013 i.e. April to August (below 0.9). The diurnal cycle of the WRF model shows a dry bias in the range of 0.5-1.5 mm. Between 00 and 01 UTC the GNSS IWV tends to be underestimate IWV which is likely due to the processing window used. The precipitation efficiency from GNSS and WRF show very good agreement on monthly bases with a maximum in May-June and minimum in August–September. The annual precipitation efficiency in 2013 at Lovech and Burgas is about 6 %.


Water ◽  
2020 ◽  
Vol 12 (7) ◽  
pp. 1918
Author(s):  
Jiaying Zhang ◽  
Liao-Fan Lin ◽  
Rafael L. Bras

Precipitation estimates from numerical weather prediction (NWP) models are uncertain. The uncertainties can be reduced by integrating precipitation observations into NWP models. This study assimilates Version 04 Integrated Multi-satellite Retrievals for the Global Precipitation Measurement (GPM) (IMERG) Final Run into the Weather Research and Forecasting (WRF) model data assimilation (WRFDA) system using a four-dimensional variational (4D-Var) method. Three synoptic-scale convective precipitation events over the central United States during 2015–2017 are used as case studies. To investigate the effect of logarithmically transformed IMERG precipitation in the WRFDA system, this study reports on several experiments with six-hour and hourly assimilation windows, regular (nontransformed) and logarithmically transformed observations, and a constant observation error in regular and logarithmic spaces. Results show that hourly assimilation windows improve precipitation simulations significantly compared to six-hour windows. Logarithmically transformed precipitation does not improve precipitation estimations relative to nontransformed precipitation. However, better predictions of heavy precipitation can be achieved with a constant error in the logarithmic space (corresponding to a linearly increasing error in the regular space), which modifies the threshold of rejecting observations, and thus utilizes more observations. This study provides a cost function with logarithmically transformed observations for the 4D-Var method in the WRFDA system for future investigations.


2006 ◽  
Vol 134 (1) ◽  
pp. 113-133 ◽  
Author(s):  
Teddy R. Holt ◽  
Dev Niyogi ◽  
Fei Chen ◽  
Kevin Manning ◽  
Margaret A. LeMone ◽  
...  

Abstract Numerical simulations are conducted using the Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) to investigate the impact of land–vegetation processes on the prediction of mesoscale convection observed on 24–25 May 2002 during the International H2O Project (IHOP_2002). The control COAMPS configuration uses the Weather Research and Forecasting (WRF) model version of the Noah land surface model (LSM) initialized using a high-resolution land surface data assimilation system (HRLDAS). Physically consistent surface fields are ensured by an 18-month spinup time for HRLDAS, and physically consistent mesoscale fields are ensured by a 2-day data assimilation spinup for COAMPS. Sensitivity simulations are performed to assess the impact of land–vegetative processes by 1) replacing the Noah LSM with a simple slab soil model (SLAB), 2) adding a photosynthesis, canopy resistance/transpiration scheme [the gas exchange/photosynthesis-based evapotranspiration model (GEM)] to the Noah LSM, and 3) replacing the HRLDAS soil moisture with the National Centers for Environmental Prediction (NCEP) 40-km Eta Data Assimilation (EDAS) operational soil fields. CONTROL, EDAS, and GEM develop convection along the dryline and frontal boundaries 2–3 h after observed, with synoptic-scale forcing determining the location and timing. SLAB convection along the boundaries is further delayed, indicating that detailed surface parameterization is necessary for a realistic model forecast. EDAS soils are generally drier and warmer than HRLDAS, resulting in more extensive development of convection along the dryline than for CONTROL. The inclusion of photosynthesis-based evapotranspiration (GEM) improves predictive skill for both air temperature and moisture. Biases in soil moisture and temperature (as well as air temperature and moisture during the prefrontal period) are larger for EDAS than HRLDAS, indicating land–vegetative processes in EDAS are forced by anomalously warmer and drier conditions than observed. Of the four simulations, the errors in SLAB predictions of these quantities are generally the largest. By adding a sophisticated transpiration model, the atmospheric model is able to better respond to the more detailed representation of soil moisture and temperature. The sensitivity of the synoptically forced convection to soil and vegetative processes including transpiration indicates that detailed representation of land surface processes should be included in weather forecasting models, particularly for severe storm forecasting where local-scale information is important.


2021 ◽  
Author(s):  
Florian Zus ◽  
Galina Dick ◽  
Jens Wickert

<p>Global Navigation Satellite Systems (GNSS) have revolutionized positioning, navigation, and timing, becoming a common part of our everyday life.  A geophysical key application is atmospheric water vapor monitoring using GNSS ground station data. GNSS water vapor data, derived from regional ground networks hereby close gaps in the established meteorological observing systems. No other observing system provides data with such high temporal and spatial resolution. The data from European GNSS networks are therefore already widely operationally used to improve regional weather forecasts in several countries. However, the impact of the currently provided data products to the forecast systems is still limited due to the limited atmospheric information content, which is provided by the currently used Zenith Total Delay (ZTD) data.</p><p>In this talk we introduce the new project EGMAP (Exploitation of GNSS tropospheric gradients for severe weather Monitoring And Prediction). This project will pioneer the development and usage of next generation data products; tropospheric gradients. The new data products, developed and provided within the project, are expected to improve the impact of the currently provided GNSS data to weather forecast systems. The main innovations, which will be addressed by the project are: (1) Developments to provide high quality ZTDs and tropospheric gradients in near-real-time for the German SAPOS network; (2) Developments to make use of ZTDs and tropospheric gradients in numerical weather prediction, i.e., implement operators in the variational/ensemble data assimilation system of the Weather Research and Forecasting (WRF) model; (3) Impact studies with the state of the art numerical weather model. In this talk we provide an overview and the current status of the project.</p>


2019 ◽  
Vol 11 (20) ◽  
pp. 2457
Author(s):  
Guangxin He ◽  
Juanzhen Sun ◽  
Zhuming Ying ◽  
Lejian Zhang

Automated and accurate radar dealiasing algorithms are very important for their assimilation into operational numerical weather forecasting models. A radar radial velocity dealiasing algorithm aimed at radar data assimilation is introduced and assessed using from several S-band and C-band radar observations under the severe weather conditions of hurricanes, typhoons, and deep continental convection in this paper. This dealiasing algorithm, named automated dealiasing for data assimilation (ADDA), is a further development of the dealiasing algorithm named the China radar network (CINRAD) improved dealiasing algorithm (CIDA), originally developed for China’s CINRAD (China Next Generation Weather Radar) radar network. The improved scheme contains five modules employed to remove noisy data, select the suitable first radial, preserve the convective regions, execute multipass dealiasing in both azimuthal and radial directions and conduct the final local dealiasing with an error check. This new dealiasing algorithm was applied to two hurricane cases, two typhoon cases, and three intense-convection cases that were observed from the CINRAD of China, Taiwan‘s radar network, and NEXRAD (Next Generation Weather Radar) of the U.S. with a continuous period of more than 12 h for each case. The dealiasing results demonstrated that ADDA performed better than CIDA for all selected cases. This algorithm not only produced a high success rate for the S-band radar, but also a reasonable performance for the C-band radar.


Sign in / Sign up

Export Citation Format

Share Document