scholarly journals Challenges of Increased Resolution for the Local Ensemble Tangent Linear Model

2020 ◽  
Vol 148 (6) ◽  
pp. 2549-2566
Author(s):  
Douglas R. Allen ◽  
Sergey Frolov ◽  
Rolf Langland ◽  
Craig H. Bishop ◽  
Karl W. Hoppel ◽  
...  
Keyword(s):  
Linear Model ◽  
Fast Moving ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Forecast Model ◽  
Nonlocal Operators ◽  
Upper Troposphere ◽  
Global Environmental ◽  
Nonlinear Normal ◽  
Normal Mode Initialization

Abstract An ensemble-based linearized forecast model has been developed for data assimilation applications for numerical weather prediction. Previous studies applied this local ensemble tangent linear model (LETLM) to various models, from simple one-dimensional models to a low-resolution (~2.5°) version of the Navy Global Environmental Model (NAVGEM) atmospheric forecast model. This paper applies the LETLM to NAVGEM at higher resolution (~1°), which required overcoming challenges including 1) balancing the computational stencil size with the ensemble size, and 2) propagating fast-moving gravity modes in the upper atmosphere. The first challenge is addressed by introducing a modified local influence volume, introducing computations on a thin grid, and using smaller time steps. The second challenge is addressed by applying nonlinear normal mode initialization, which damps spurious fast-moving modes and improves the LETLM errors above ~100 hPa. Compared to a semi-Lagrangian tangent linear model (TLM), the LETLM has superior skill in the lower troposphere (below 700 hPa), which is attributed to better representation of moist physics in the LETLM. The LETLM skill slightly lags in the upper troposphere and stratosphere (700–2 hPa), which is attributed to nonlocal aspects of the TLM including spectral operators converting from winds to vorticity and divergence. Several ways forward are suggested, including integrating the LETLM in a hybrid 4D variational solver for a realistic atmosphere, combining a physics LETLM with a conventional TLM for the dynamics, and separating the LETLM into a sequence of local and nonlocal operators.

scholarly journals GPS radio occultation with TerraSAR-X and TanDEM-X: sensitivity of lower troposphere sounding to the Open-Loop Doppler model

2014 ◽  
Vol 7 (12) ◽  
pp. 12719-12733 ◽  
Author(s):  
F. Zus ◽  
G. Beyerle ◽  
S. Heise ◽  
T. Schmidt ◽  
J. Wickert
Keyword(s):  
Radio Occultation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Systematic Difference ◽  
Open Loop ◽  
Negative Bias ◽  
Upper Troposphere ◽  
Gps Radio Occultation ◽  
Global Positioning ◽  
Data Source

Abstract. The Global Positioning System (GPS) radio occultation (RO) technique provides valuable input for numerical weather prediction and is considered as a data source for climate related research. Numerous studies outline the high precision and accuracy of RO atmospheric soundings in the upper troposphere and lower stratosphere. In this altitude region (8–25 km) RO atmospheric soundings are considered to be free of any systematic error. In the tropical (30° S–30° N) Lower (<8 km) Troposphere (LT), this is not the case; systematic differences with respect to independent data sources exist and are still not completely understood. To date only little attention has been paid to the Open Loop (OL) Doppler model. Here we report on a RO experiment carried out on-board of the twin satellite configuration TerraSAR-X and TanDEM-X which possibly explains to some extent biases in the tropical LT. In two sessions we altered the OL Doppler model aboard TanDEM-X by not more than ±5 Hz with respect to TerraSAR-X and compare collocated atmospheric refractivity profiles. We find a systematic difference in the retrieved refractivity. The bias mainly stems from the tropical LT; there the bias reaches up to ±1%. Hence, we conclude that the negative bias (several Hz) of the OL Doppler model aboard TerraSAR-X introduces a negative bias (in addition to the negative bias which is primarily caused by critical refraction) in our retrieved refractivity in the tropical LT.

A Hybrid Differential-Ensemble Linear Forecast Model for 4D-Var

2021 ◽  
Vol 149 (1) ◽  
pp. 3-19
Author(s):  
T. J. Payne
Keyword(s):  
Data Assimilation ◽  
Linear Model ◽  
Numerical Weather Prediction ◽  
Weather Prediction ◽  
Forecast Model ◽  
Regression Technique ◽  
Time Step ◽  
Dynamical Core ◽  
Numerical Weather ◽  
Better Than

AbstractA key component of the 4D-Var data assimilation method used widely for numerical weather prediction is the linear forecast model, which is approximately tangent linear to the forecast model. Traditionally this has been based on differentiating the forecast model, though recently some authors have experimented with an ensemble regression technique, the localized ensemble tangent linear model (LETLM). We propose a hybrid of the two, in which a simplified conventional tangent-linear model (e.g., just the dynamical core) is used together with an LETLM-like adjustment every time step to account for the remaining processes (in this example, the parameterized physics). This is much cheaper than the LETLM, and in tests using the Met Office’s linear model performs considerably better than either a pure LETLM (with a very large ensemble) or the existing linear model.

scholarly journals The Pan-Canadian High Resolution (2.5 km) Deterministic Prediction System

2016 ◽  
Vol 31 (6) ◽  
pp. 1791-1816 ◽  
Author(s):  
Jason A. Milbrandt ◽  
Stéphane Bélair ◽  
Manon Faucher ◽  
Marcel Vallée ◽  
Marco L. Carrera ◽  
...  
Keyword(s):  
High Resolution ◽  
Weather Prediction ◽  
Forecast Model ◽  
Added Value ◽  
Precipitation Forecast ◽  
Initial Surface ◽  
Prediction System ◽  
Global Environmental ◽  
Fractional Cloud ◽  
Land Data Assimilation System

Abstract Since November 2014, the Meteorological Services of Canada (MSC) has been running a real-time numerical weather prediction system that provides deterministic forecasts on a regional domain with a 2.5-km horizontal grid spacing covering a large portion of Canada using the Global Environmental Multiscale (GEM) forecast model. This system, referred to as the High Resolution Deterministic Prediction System (HRDPS), is currently downscaled from MSC’s operational 10-km GEM-based regional system but uses initial surface fields from a high-resolution (2.5 km) land data assimilation system coupled to the HRDPS and initial hydrometeor fields from the forecast of a 2.5-km cycle, which reduces the spinup time for clouds and precipitation. Forecast runs of 48 h are provided four times daily. The HRDPS was tested and compared to the operational 10-km system. Model runs from the two systems were evaluated against surface observations for common weather elements (temperature, humidity, winds, and precipitation), fractional cloud cover, and also against upper-air soundings, all using standard metrics. Although the predictions of some fields were degraded in some specific regions, the HRDPS generally outperformed the operational system for a majority of the scores. The evaluation illustrates the added value of the 2.5-km model and the potential for improved numerical guidance for the prediction of high-impact weather.

Impact of Infrared, Microwave, and Radio Occultation Satellite Observations on Operational Numerical Weather Prediction

2014 ◽  
Vol 142 (11) ◽  
pp. 4164-4186 ◽  
Author(s):  
L. Cucurull ◽  
R. A. Anthes
Keyword(s):  
Radio Occultation ◽  
Positive Impact ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Correlation Coefficients ◽  
Forecast Model ◽  
Satellite Observation ◽  
Satellite Observations ◽  
Satellite Systems ◽  
The Impact

Abstract A comparison of the impact of infrared (IR), microwave (MW), and radio occultation (RO) observations on NCEP’s operational global forecast model over the month of March 2013 is presented. Analyses and forecasts with only IR, MW, and RO observations are compared with analyses and forecasts with no satellite data and with each other. Overall, the patterns of the impact of the different satellite systems are similar, with the MW observations producing the largest impact on the analyses and RO producing the smallest. Without RO observations, satellite radiances are over– or under–bias corrected and RO acts as an anchor observation, reducing the forecast biases globally. Positive correlation coefficients of temperature impacts are generally found between the different satellite observation analyses, indicating that the three satellite systems are affecting the global temperatures in a similar way. However, the correlation in the lower troposphere among all three systems is surprisingly small. Correlations for the moisture field tend to be small in the lower troposphere between the different satellite analyses. The impact of the satellite observations on the 500-hPa geopotential height forecasts is much different in the Northern and Southern Hemispheres. In the Northern Hemisphere, all the satellite observations together make a small positive impact compared to the base (no satellite) forecasts. The IR and MW, but not the RO, make a small positive impact when assimilated alone. The situation is considerably different in the Southern Hemisphere, where all the satellite observations together make a much larger positive impact, and all three observation types (IR, MW, and RO) make similar and significant impacts.

scholarly journals Towards All-sky Assimilation of Microwave Temperature Sounding Channels in Environment Canada’s Global Deterministic Weather Prediction System

2021 ◽  
Author(s):  
Maziar Bani Shahabadi ◽  
Mark Buehner
Keyword(s):  
Standard Deviation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Forecast Model ◽  
Specific Humidity ◽  
Horizontal Wind ◽  
Observation Error ◽  
Horizontal Wind Speed ◽  
Observation Operator ◽  
Background State

AbstractThe all-sky assimilation of radiances from microwave instruments is developed in the 4D-EnVar analysis system at Environment and Climate Change Canada (ECCC). Assimilation of cloud-affected radiances from Advanced Microwave Sounding Unit A (AMSUA) temperature sounding channels 4 and 5 for non-precipitating scenes over the ocean surface is the focus of this study. Cloud-affected radiances are discarded in the ECCC operational data assimilation system due to the limitations of forecast model physics, radiative transfer models, and the strong non-linearity of the observation operator. In addition to using symmetric estimate of innovation standard deviation for quality control, a state-dependent observation error inflation is employed at the analysis stage. The background state clouds are scaled by a factor of 0.5 to compensate for a systematic overestimation by the forecast model, before being used in the observation operator. The changes in the fit of the background state to observations show mixed results. The number of AMSUA channels 4 and 5 assimilated observations in the all-sky experiment is 5-12% higher than in the operational system. The all-sky approach improves temperature analysis when verified against ECMWF operational analysis in the areas where the extra cloud-affected observations were assimilated. Statistically significant reductions in error standard deviation by 1-4% for the analysis and forecasts of temperature, specific humidity, and horizontal wind speed up to maximum 4 days were achieved in the all-sky experiment in the lower troposphere. These improvements result mainly from the use of cloud information for computing the observation-minus-background departures. The operational implementation of all-sky assimilation is planned for Fall 2021.

scholarly journals Global 3D Features of Error Variances of GPS Radio Occultation and Radiosonde Observations

2020 ◽  
Vol 13 (1) ◽  
pp. 1
Author(s):  
Xu Xu ◽  
Xiaolei Zou
Keyword(s):  
Radio Occultation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Error Variance ◽  
Specific Humidity ◽  
Temperature Error ◽  
Observation Error ◽  
Bending Angle ◽  
Gps Radio Occultation ◽  
Low Latitudes

Global Positioning System (GPS) radio occultation (RO) and radiosonde (RS) observations are two major types of observations assimilated in numerical weather prediction (NWP) systems. Observation error variances are required input that determines the weightings given to observations in data assimilation. This study estimates the error variances of global GPS RO refractivity and bending angle and RS temperature and humidity observations at 521 selected RS stations using the three-cornered hat method with additional ERA-Interim reanalysis and Global Forecast System forecast data available from 1 January 2016 to 31 August 2019. The global distributions, of both RO and RS observation error variances, are analyzed in terms of vertical and latitudinal variations. Error variances of RO refractivity and bending angle and RS specific humidity in the lower troposphere, such as at 850 hPa (3.5 km impact height for the bending angle), all increase with decreasing latitude. The error variances of RO refractivity and bending angle and RS specific humidity can reach about 30 N-unit2, 3 × 10−6 rad2, and 2 (g kg−1)2, respectively. There is also a good symmetry of the error variances of both RO refractivity and bending angle with respect to the equator between the Northern and Southern Hemispheres at all vertical levels. In this study, we provide the mean error variances of refractivity and bending angle in every 5°-latitude band between the equator and 60°N, as well as every interval of 10 hPa pressure or 0.2 km impact height. The RS temperature error variance distribution differs from those of refractivity, bending angle, and humidity, which, at low latitudes, are smaller (less than 1 K2) than those in the midlatitudes (more than 3 K2). In the midlatitudes, the RS temperature error variances in North America are larger than those in East Asia and Europe, which may arise from different radiosonde types among the above three regions.

scholarly journals The Sentinel-3 SLSTR Atmospheric Motion Vectors Product at EUMETSAT

2021 ◽  
Vol 13 (9) ◽  
pp. 1702
Author(s):  
Kévin Barbieux ◽  
Olivier Hautecoeur ◽  
Maurizio De Bartolomei ◽  
Manuel Carranza ◽  
Régis Borde
Keyword(s):  
Land Surface ◽  
Weather Prediction ◽  
Forecast Model ◽  
Low Earth Orbit ◽  
Reference Image ◽  
Motion Vectors ◽  
Infrared Imager ◽  
Atmospheric Motion ◽  
Atmospheric Motion Vectors ◽  
Leo Satellite

Atmospheric Motion Vectors (AMVs) are an important input to many Numerical Weather Prediction (NWP) models. EUMETSAT derives AMVs from several of its orbiting satellites, including the geostationary satellites (Meteosat), and its Low-Earth Orbit (LEO) satellites. The algorithm extracting the AMVs uses pairs or triplets of images, and tracks the motion of clouds or water vapour features from one image to another. Currently, EUMETSAT LEO satellite AMVs are retrieved from georeferenced images from the Advanced Very-High-Resolution Radiometer (AVHRR) on board the Metop satellites. EUMETSAT is currently preparing the operational release of an AMV product from the Sea and Land Surface Temperature Radiometer (SLSTR) on board the Sentinel-3 satellites. The main innovation in the processing, compared with AVHRR AMVs, lies in the co-registration of pairs of images: the images are first projected on an equal-area grid, before applying the AMV extraction algorithm. This approach has multiple advantages. First, individual pixels represent areas of equal sizes, which is crucial to ensure that the tracking is consistent throughout the processed image, and from one image to another. Second, this allows features that would otherwise leave the frame of the reference image to be tracked, thereby allowing more AMVs to be derived. Third, the same framework could be used for every LEO satellite, allowing an overall consistency of EUMETSAT AMV products. In this work, we present the results of this method for SLSTR by comparing the AMVs to the forecast model. We validate our results against AMVs currently derived from AVHRR and the Spinning Enhanced Visible and InfraRed Imager (SEVIRI). The release of the operational SLSTR AMV product is expected in 2022.

scholarly journals Application of Adjoint-Derived Forecast Sensitivities to the 24–25 January 2000 U.S. East Coast Snowstorm

2005 ◽  
Vol 133 (11) ◽  
pp. 3148-3175 ◽  
Author(s):  
Daryl T. Kleist ◽  
Michael C. Morgan
Keyword(s):  
Response Function ◽  
Initial Conditions ◽  
Forecast Error ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Vertical Motion ◽  
Initial Condition ◽  
Response Functions ◽  
Eastern United States ◽  
Surface Cyclone

Abstract The 24–25 January 2000 eastern United States snowstorm was noteworthy as operational numerical weather prediction (NWP) guidance was poor for lead times as short as 36 h. Despite improvements in the forecast of the surface cyclone position and intensity at 1200 UTC 25 January 2000 with decreasing lead time, NWP guidance placed the westward extent of the midtropospheric, frontogenetically forced precipitation shield too far to the east. To assess the influence of initial condition uncertainties on the forecast of this event, an adjoint model is used to evaluate forecast sensitivities for 36- and 48-h forecasts valid at 1200 UTC 25 January 2000 using as response functions the energy-weighted forecast error, lower-tropospheric circulation about a box surrounding the surface cyclone, 750-hPa frontogenesis, and vertical motion. The sensitivities with respect to the initial conditions for these response functions are in general very similar: geographically isolated, maximized in the middle and lower troposphere, and possessing an upshear vertical tilt. The sensitivities are maximized in a region of enhanced low-level baroclinicity in the vicinity of the surface cyclone’s precursor upper trough. However, differences in the phase and structure of the gradients for the four response functions are evident, which suggests that perturbations could be constructed to alter one response function but not necessarily the others. Gradients of the forecast error response function with respect to the initial conditions are used in an iterative procedure to construct initial condition perturbations that reduce the forecast error. These initial condition perturbations were small in terms of both spatial scale and magnitude. Those initial condition perturbations that were confined primarily to the midtroposphere grew rapidly into much larger amplitude upper-and-lower tropospheric perturbations. The perturbed forecasts were not only characterized by reduced final time forecast error, but also had a synoptic evolution that more closely followed analyses and observations.

scholarly journals Estimation and evaluation of COSMIC radio occultation excess phase using undifferenced measurements

2017 ◽  
Vol 10 (5) ◽  
pp. 1813-1821
Author(s):  
Pengfei Xia ◽  
Shirong Ye ◽  
Kecai Jiang ◽  
Dezhong Chen
Keyword(s):  
Radio Occultation ◽  
Weather Prediction ◽  
Lower Troposphere ◽  
Double Difference ◽  
Cosmic Radio ◽  
Processing Strategy ◽  
Single Difference ◽  
Excess Phase ◽  
Mean Square Errors ◽  
Better Than

Abstract. In the GPS radio occultation technique, the atmospheric excess phase (AEP) can be used to derive the refractivity, which is an important quantity in numerical weather prediction. The AEP is conventionally estimated based on GPS double-difference or single-difference techniques. These two techniques, however, rely on the reference data in the data processing, increasing the complexity of computation. In this study, an undifferenced (ND) processing strategy is proposed to estimate the AEP. To begin with, we use PANDA (Positioning and Navigation Data Analyst) software to perform the precise orbit determination (POD) for the purpose of acquiring the position and velocity of the mass centre of the COSMIC (The Constellation Observing System for Meteorology, Ionosphere and Climate) satellites and the corresponding receiver clock offset. The bending angles, refractivity and dry temperature profiles are derived from the estimated AEP using Radio Occultation Processing Package (ROPP) software. The ND method is validated by the COSMIC products in typical rising and setting occultation events. Results indicate that rms (root mean square) errors of relative refractivity differences between undifferenced and atmospheric profiles (atmPrf) provided by UCAR/CDAAC (University Corporation for Atmospheric Research/COSMIC Data Analysis and Archive Centre) are better than 4 and 3 % in rising and setting occultation events respectively. In addition, we also compare the relative refractivity bias between ND-derived methods and atmPrf profiles of globally distributed 200 COSMIC occultation events on 12 December 2013. The statistical results indicate that the average rms relative refractivity deviation between ND-derived and COSMIC profiles is better than 2 % in the rising occultation event and better than 1.7 % in the setting occultation event. Moreover, the observed COSMIC refractivity profiles from ND processing strategy are further validated using European Centre for Medium-Range Weather Forecasts (ECMWF) analysis data, and the results indicate that the undifferenced method reduces the noise level on the excess phase paths in the lower troposphere compared to the single-difference processing strategy.

Sign in / Sign up

Export Citation Format

Share Document