Integrated Forecasting System Research Articles

Evaluation of isoprene emissions from the coupled model SURFEX–MEGANv2.1

Abstract. Isoprene, a key biogenic volatile organic compound, plays a pivotal role in atmospheric chemistry. Due to its high reactivity, this compound contributes significantly to the production of tropospheric ozone in polluted areas and to the formation of secondary organic aerosols. The assessment of biogenic emissions is of great importance for regional and global air quality evaluation. In this study, we implemented the biogenic emission model MEGANv2.1 (Model of Emissions of Gases and Aerosols from Nature, version 2.1) in the surface model SURFEXv8.1 (SURface EXternalisée in French, version 8.1). This coupling aims to improve the estimation of biogenic emissions using the detailed vegetation-type-dependent treatment included in the SURFEX vegetation ISBA (Interaction between Soil Biosphere and Atmosphere) scheme. This scheme provides vegetation-dependent parameters such as leaf area index and soil moisture to MEGAN. This approach enables a more accurate estimation of biogenic fluxes compared to the stand-alone MEGAN model, which relies on average input values for all vegetation types. The present study focuses on the assessment of the SURFEX–MEGAN model isoprene emissions. An evaluation of the coupled SURFEX–MEGAN model results was carried out by conducting a global isoprene emission simulation in 2019 and by comparing the simulation results with other MEGAN-based isoprene inventories. The coupled model estimates a total global isoprene emission of 443 Tg in 2019. The estimated isoprene is within the range of results obtained with other MEGAN-based isoprene inventories, ranging from 311 to 637 Tg. The spatial distribution of SURFEX–MEGAN isoprene is consistent with other studies, with some differences located in low-isoprene-emission regions. Several sensitivity tests were conducted to quantify the impact of different model inputs and configurations on isoprene emissions. Using different meteorological forcings resulted in a ±5 % change in isoprene emissions using MERRA (Modern-Era Retrospective analysis for Research and Applications) and IFS (Integrated Forecasting System) compared with ERA5. The impact of using different emission factor data was also investigated. The use of PFT (plant functional type) spatial coverage and PFT-dependent emission potential data resulted in a 12 % reduction compared to using the isoprene emission potential gridded map. A significant reduction of around 38 % in global isoprene emissions was observed in the third sensitivity analysis, which applied a parameterization of soil moisture deficit, particularly in certain regions of Australia, Africa, and South America. The significance of coupling the SURFEX and MEGAN models lies particularly in the ability of the coupled model to be forced with meteorological data from any period. This means, for instance, that this system can be used to predict biogenic emissions in the future. This aspect of our work is significant given the changes that biogenic organic compounds are expected to undergo as a result of changes in their climatic factors.

Evaluation of downward and upward solar irradiances simulated by the Integrated Forecasting System of ECMWF using airborne observations above Arctic low-level clouds

Abstract. The simulations of upward and downward irradiances by the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts are compared with broadband solar irradiance measurements from the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign. For this purpose, offline radiative transfer simulations were performed with the ecRad radiation scheme using the operational IFS output. The simulations of the downward solar irradiance agree within the measurement uncertainty. However, the IFS underestimates the reflected solar irradiances above sea ice significantly by −35 W m−2. Above open ocean, the agreement is closer, with an overestimation of 28 W m−2. A sensitivity study using measured surface and cloud properties is performed with ecRad to quantify the contributions of the surface albedo, cloud fraction, ice and liquid water path and cloud droplet number concentration to the observed bias. It shows that the IFS sea ice albedo climatology underestimates the observed sea ice albedo, causing more than 50 % of the bias. Considering the higher variability of in situ observations in the parameterization of the cloud droplet number concentration leads to a smaller bias of −27 W m−2 above sea ice and a larger bias of 48 W m−2 above open ocean by increasing the range from 36–69 to 36–200 cm−3. Above sea ice, realistic surface albedos, cloud droplet number concentrations and liquid water paths contribute most to the bias improvement. Above open ocean, realistic cloud fractions and liquid water paths are most important for reducing the model–observation differences.

A new implementation of FLEXPART with Enviro-HIRLAM meteorological input, and a case study during a heavy air pollution event

ABSTRACT We integrated Enviro-HIRLAM (Environment-High Resolution Limited Area Model) meteorological output into FLEXPART (FLEXible PARTicle dispersion model). A FLEXPART simulation requires meteorological input from a numerical weather prediction (NWP) model. The publicly available version of FLEXPART can utilize either ECMWF (European Centre for Medium-range Weather Forecasts) Integrated Forecast System (IFS) forecast or reanalysis NWP data, or NCEP (U.S. National Center for Environmental Prediction) Global Forecast System (GFS) forecast or reanalysis NWP data. The primary benefits of using Enviro-HIRLAM are that it runs at a higher resolution and accounts for aerosol effects in meteorological fields. We compared backward trajectories generated with FLEXPART using Enviro-HIRLAM (both with and without aerosol effects) to trajectories generated using NCEP GFS and ECMWF IFS meteorological inputs, for a case study of a heavy haze event which occurred in Beijing, China in November 2018. We found that results from FLEXPART were considerably different when using different meteorological inputs. When aerosol effects were included in the NWP, there was a small but noticeable difference in calculated trajectories. Moreover, when looking at potential emission sensitivity instead of simply expressing trajectories as lines, additional information, which may have been missed when looking only at trajectories as lines, can be inferred.

Impact of Saharan dust outbreaks on short‐range weather forecast errors in Europe

AbstractMineral dust, the most abundant atmospheric aerosol by mass, interacts with radiation directly and alters cloud properties indirectly. Many operational numerical weather prediction models account for aerosol direct effects by using climatological mean concentrations and neglect indirect effects. This simplification may lead to shortcomings in model forecasts during outbreaks of Saharan dust towards Europe, when climatological mean dust concentrations deviate strongly from actual concentrations. This study investigates errors in model analyses and short‐range forecasts during such events. We investigate a pronounced dust event in March 2021 using the pre‐operational ICOsahedral Nonhydrostatic weather and climate model with Aerosols and Reactive Trace gases (ICON‐ART) with prognostic calculation of dust and the operational European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) model, which deploys a dust climatology. We compare model analysis and forecast with measurements from satellite and in situ instruments. We find that inclusion of prognostic aerosol and direct radiative effects from dust improves forecasts of surface radiation during clear‐sky conditions. However, dust‐induced cirrus clouds are strongly underestimated, highlighting the importance of representing indirect effects adequately. These findings are corroborated by systematic quantification of forecast errors against satellite measurements. For this we construct an event catalogue with 49 dust days over Central Europe between January 2018 and March 2022. We classify model cells by simulated and observed cloudiness and simulated dustiness in the total atmospheric column. We find significant overestimations of brightness temperature for cases with dust compared with cases without dust. For surface shortwave radiation, we find median overestimations of 6.2% during cloudy conditions with dust optical depth greater than 0.1, however these are not significant compared with cloudy conditions without dust. Our findings show that the pre‐operational ICON‐ART and the operational IFS model still do not reproduce cloudiness adequately during events with Saharan dust over Central Europe. Missing implementations of prognostic dust, particularly of indirect effects on cloud formation, lead to significant underestimations of cloudiness and potentially overestimations of surface radiation.

Deep learning‐based postprocessing for hourly temperature forecasting

AbstractIn this article, a prediction model based on spatiotemporal stacked ResNet (Res‐STS) for hourly temperature prediction is designed. On the timescale, the Res‐STS removes the gate structure of the long short‐term memory (LSTM) model, and the data of multiple consecutive time nodes are stacked together to preserve all temporal characteristics of the data. A point‐to‐point data mapping relationship is developed on the spatial scale to maximize the impact of large‐scale environmental background field characteristics on a single grid point. Based on the historical gridded data from the China Meteorological Administration land data assimilation system (CLDAS) and the optimal factor dataset of the European Centre for Medium‐Range Weather Forecasts Integrated Forecasting System (ECMWF‐IFS) from 2017 to 2020, hourly temperature prediction models based on convolutional long short‐term memory (ConvLSTM) and Res‐STS model are developed, respectively. Furthermore, the prediction results of the two models in 2021 are compared with the ECMWF‐IFS. The results show that the root mean square error (RMSE) of the prediction results by ConvLSTM and Res‐STS models are both smaller than that of ECMWF‐IFS. Specially, the Res‐STS model performs best: it reduces the RMSE by 20.8% (24.5%) compared with the ConvLSTM (ECMWF‐IFS). Specifically, the RMSE peaks in the afternoon when the daily maximum temperature occurs, while it is relatively smaller at night. Res‐STS demonstrates a significant improvement in forecast performance compared with ECMWF‐IFS, while ConvLSTM's correction during the period of maximum temperature occurrence has been enhanced. Moreover, the forecast performance of the Res‐STS model is least affected by terrain compared with those of the ConvLSTM and ECMWF‐IFS. For the regions with terrain height greater than 1 km, the model Res‐STS evidently improves the RMSE.

The Role of Non‐Local Effects on Surface Sensible Heat Flux Under Different Types of Thermal Structures Over the Arctic Sea‐Ice Surface

AbstractThe effects of atmospheric thermal structure on the surface energy flux are poorly understood over the Arctic sea‐ice surface. Here, we explore the mechanism of sensible heat exchange under different types of thermal structures over the Arctic sea‐ice surface by using data collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. The quadrant analysis indicates that strong surface temperature inversions below 100 m enhance non‐local effects on the positive (upward) sensible heat flux () through entrainment of large eddies from the convective boundary layer aloft. However, strong surface inversions restrict the contributions of large eddies to the negative (downward) due to intensified surface stability. By inspecting the existing parameterization schemes, we found that the European Center for Medium‐Range Weather Forecasts Integrated Forecasting System scheme fails to predict the impacts of non‐local processes on the positive , and an adjustment term to correct the bias of parameterized is proposed.

Convective gravity wave events during summer near 54° N, present in both AIRS and Rayleigh–Mie–Raman (RMR) lidar observations

Abstract. We connect tropospheric deep convective events over western Europe, as measured by the 8.1 µm radiance observations from the Atmospheric Infrared Sounder (AIRS) on NASA's Aqua satellite, to horizontal brightness temperature variance in the 4.3 µm AIRS channel (maximum sensitivity at around 40 km) and temperature perturbations in vertical lidar profiles (between 33-43 km) over Kühlungsborn, Germany (54.12∘ N, 11.77∘ E). Although the lidar and AIRS are sensitive to different parts of the gravity wave spectrum, they both capture the same peaks in gravity wave activity tied to convection. This suggests that a broad range of vertical wavelengths is present in the convective gravity waves. To account for wave propagation conditions from the troposphere to the stratosphere, we also consider the horizontal winds in the troposphere and stratosphere using the ECMWF Integrated Forecasting System (IFS) operational analysis. In this work, we highlight sporadic peaks in gravity wave activity in summer greatly exceeding those typical of summer, which is generally a season with lower wave activity compared to winter. Although these events are present in roughly half of the years (between 2003 and 2019), we focus our study on two case study years (2014 and 2015). These case study years were chosen because of the high cadence of lidar soundings close in time to the convective events. These events, while sporadic, could contribute significantly to the zonal mean momentum budget and are not accounted for in weather and climate models.

Twelve Times Faster yet Accurate: A New State‐Of‐The‐Art in Radiation Schemes via Performance and Spectral Optimization

AbstractRadiation schemes are critical components of Earth system models that need to be both efficient and accurate. Despite the use of approximations such as 1D radiative transfer, radiation can account for a large share of the runtime of expensive climate simulations. Here we seek a new state‐of‐the‐art in speed and accuracy by combining code optimization with improved algorithms. To fully benefit from new spectrally reduced gas optics schemes, we restructure code to avoid short vectorized loops where possible by collapsing the spectral and vertical dimensions. Our main focus is the ecRad radiation scheme, where this requires batching of adjacent cloudy layers, trading some simplicity for improved vectorization and instruction‐level parallelism. When combined with common optimization techniques for serial code and porting widely used two‐stream kernels fully to single precision, we find that ecRad with the TripleClouds solver becomes 12 times faster than the operational radiation scheme in ECMWF's Integrated Forecast System (IFS) cycle 47r3, which uses a less accurate gas optics model (RRMTG) and a more noisy solver (McICA). After applying the spectral reduction and extensive optimizations to the more sophisticated SPARTACUS solver, we find that it’s 2.5 times faster than IFS cy47r3 radiation, making cloud 3D radiative effects affordable to compute in large‐scale models. The code optimization itself gave a threefold speedup for both solvers. While SPARTACUS is still under development, preliminary experiments show slightly improved medium‐range forecasts of 2‐m temperature in the tropics, and in year‐long coupled atmosphere‐ocean simulations the 3D effects warm the surface substantially.

Development of an Optimal Postprocessing Model Using the Microgenetic Algorithm to Improve Precipitation Forecasting in South Korea

Abstract We developed an advanced postprocessing model for precipitation forecasting using a microgenetic algorithm (MGA). The algorithm determines the optimal combination of three general circulation models: the Korean Integrated Model, the Unified Model, and the Integrated Forecast System model. To measure model accuracy, including the critical success index (CSI), probability of detection (POD), and frequency bias index, the MGA calculates optimal weights for individual models based on a fitness function that considers various indices. Our optimized multimodel yielded up to 13% and 10% improvement in CSI and POD performance compared to each individual model, respectively. Notably, when applied to an operational definition that considers precipitation thresholds from three models and averages the precipitation amount from the satisfactory models, our optimized multimodel outperformed the current operational model used by the Korea Meteorological Administration by up to 1.0% and 6.8% in terms of CSI and false alarm ratio performance, respectively. This study highlights the effectiveness of a weighted combination of global models to enhance the forecasting accuracy for regional precipitation. By utilizing the MGA for the fine-tuning of model weights, we achieved superior precipitation prediction compared to that of individual models and existing standard postprocessing operations. This approach can significantly improve the accuracy of precipitation forecasts. Significance Statement We developed an optimized multimodel for predicting precipitation occurrence using advanced techniques. By integrating various weather models with their optimized weights, our approach outperforms the method of using an arithmetic average of all models. This study underscores the potential to enhance regional precipitation forecasts, thereby facilitating more precise weather predictions for the public.

Deep learning for quality control of surface physiographic fields using satellite Earth observations

Abstract. A purposely built deep learning algorithm for the Verification of Earth System ParametERization (VESPER) is used to assess recent upgrades to the global physiographic datasets underpinning the quality of the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF), which is used in both numerical weather prediction and climate reanalyses. A neural network regression model is trained to learn the mapping between the surface physiographic dataset, plus the main meteorologic fields from ERA5, and the MODIS satellite skin temperature observations. Once trained, this tool is applied to rapidly assess the quality of upgrades to the physiographic fields used by land surface schemes. Upgrades which improve the prediction accuracy of the machine learning tool indicate a reduction in the errors in the surface fields used as input to the surface parameterization schemes. Conversely, incorrect specifications of the surface fields decrease the accuracy with which VESPER can make predictions. We apply VESPER to assess the accuracy of recent upgrades to the permanent lake and glacier covers, as well as of planned upgrades to represent seasonally varying water bodies (i.e. ephemeral lakes). We show that, for grid cells where the lake fields have been updated, the prediction accuracy of VESPER in the land surface temperature (as quantified by the mean absolute error) improves by 0.37 K on average, whilst for the subset of points where the lakes have been completely removed and replaced with bare ground, the improvement is 0.83 K. We also show that updates to the glacier cover improve the prediction accuracy by 0.22 K. We highlight how neural networks such as VESPER can assist the research and development of surface parameterizations and their input physiography to better represent Earth’s surface coupled processes in weather and climate models.

A deep-learning real-time bias correction method for significant wave height forecasts in the Western North Pacific

Significant wave height (SWH) is one of the most important parameters characterizing ocean waves, and accurate numerical ocean wave forecasting is crucial for coastal protection and shipping. However, due to the randomness and nonlinearity of the wind fields that generate ocean waves and the complex interaction between wave and wind fields, current forecasts of numerical ocean waves have biases. In this study, a spatiotemporal deep-learning method was employed to correct gridded SWH forecasts from the European Centre for Medium-Range Weather Forecasting System Integrated Forecast System Global Model (ECMWF-IFS). This method was built on the trajectory gated recurrent unit deep neural network, and it conducts real-time rolling correction for the 0–240-h SWH forecasts from ECMWF-IFS. The correction model is co-driven by wave and wind fields, providing better results than those based on wave fields alone. A novel pixel-switch loss function was developed. The pixel-switch loss function can dynamically fine-tune the pre-trained correction model, focusing on pixels with large biases in SWH forecasts. According to the seasonal characteristics of SWH, four correction models were constructed separately, for spring, summer, autumn, and winter. The experimental results show that, compared with the original ECMWF SWH predictions, the correction was most effective in spring, when the mean absolute error decreased by 12.972–46.237%. Although winter had the worst performance, the mean absolute error decreased by 13.794–38.953%. The corrected results improved the original ECMWF SWH forecasts under both normal and extreme weather conditions, indicating that our SWH correction model is robust and generalizable.

A neural-network-based method for generating synthetic 1.6 µm near-infrared satellite images

Abstract. In combination with observations from visible satellite channels, near-infrared channels can provide valuable additional cloud information, e.g. on cloud phase and particle sizes, which is also complementary to the information content of thermal infrared channels. Exploiting near-infrared channels for operational data assimilation and model evaluation requires a sufficiently fast and accurate forward operator. This study presents an extension to the method for fast satellite image synthesis (MFASIS) that allows for simulating reflectances of the 1.6 µm near-infrared channel based on a computationally efficient neural network with the same accuracy that has already been achieved for visible channels. For this purpose, it is important to better represent vertical variations in effective cloud particle radii, as well as mixed-phase clouds and molecular absorption in the idealized profiles used to train the neural network. A new approach employing a two-layer model of water, ice and mixed-phase clouds is described, and the relative importance of the different input parameters characterizing the idealized profiles is analysed. A comprehensive data set sampled from Integrated Forecasting System (IFS) forecasts together with different parameterizations of the effective water and ice particle radii is used for the development and evaluation of the method. Further evaluation uses a month of ICOsahedral Non-hydrostatic development based on version 2.6.1 (ICON-D2) hindcasts with effective radii directly determined by the two-moment microphysics scheme of the model. In all cases, the mean absolute reflectance error achieved is about 0.01 or smaller, which is an order of magnitude smaller than typical differences between reflectance observations and corresponding model values. The errors related to the imperfect training of the neural networks present only a small contribution to the total error, and evaluating the networks takes less than a microsecond per column on standard CPUs. The method is also applicable for many other visible and near-infrared channels with weak water vapour sensitivity.

Mediterranean tropical-like cyclone forecasts and analysis using the ECMWF ensemble forecasting system with physical parameterization perturbations

Abstract. Mediterranean tropical-like cyclones, called medicanes, present a multi-scale nature, and their track and intensity have been recognized as highly sensitive to large-scale atmospheric forcing and diabatic heating as represented by the physical parameterizations in numerical weather prediction. Here, we analyse the structure and investigate the predictability of medicanes with the aid of the European Centre for Medium-Range Weather Forecast (ECMWF) Integrated Forecast System (IFS) ensemble forecasting system with 25 perturbed members at 9 km horizontal resolution (compared with the 16 km operational resolution). The IFS ensemble system includes the representation of initial uncertainties from the ensemble data assimilation (EDA) and a recently developed uncertainty representation of the model physics with perturbed parameters (stochastically perturbed parameterizations, SPP). The focus is on three medicanes, Ianos, Zorbas, and Trixie, among the strongest in recent years. In particular, we have carried out separate ensemble simulations with initial perturbations, full physics SPP, with a reduced set of SPP, where only convection is perturbed to highlight the convective nature of medicanes and an operational ensemble combining the SPP and the initial perturbations. It is found that compared with the operational analysis and satellite rainfall data, the forecasts reproduce the tropical-like features of these cyclones. Furthermore, the SPP simulations compare to the initial-condition perturbation ensemble in terms of tracking, intensity, precipitation, and, more generally, in terms of ensemble skill and spread. Moreover, the study confirms that similar processes are at play in the development of the investigated three medicanes, in that the predictability of these cyclones is linked not only to the prediction of the precursor events (namely the deep cutoff low) but also to the interaction of the upper-level advected potential vorticity (PV) streamer with the tropospheric PV anomaly that is diabatically produced by latent heat.

Assessment of Typhoon Precipitation Forecasts Based on Topographic Factors

For this paper, a new global atmospheric model (Global-to-Regional Integrated forecast SysTem; GRIST) with improved sub-grid scale orographic parameterization was verified and assessed, with an emphasis on the precipitation caused by typhoons. Four typical typhoon cases were selected for the verification of the model. The results indicate that, compared to the control experiments, the sensitivity experiments consistently simulated the trends in the three-hour cumulative precipitation changes and the high-value regions of total precipitation better. However, the improved experiments only had an ameliorating effect on the cumulative precipitation modelling biases for Typhoon LEKIMA and Typhoon HAGUPIT, not all of them. Precipitation bias is smaller on flat land than that on mountainous land, but the precipitation bias on windward/leeward slopes depends on the typhoon case. Precipitation modelling accuracy varies considerably between flat and mountainous terrain but very little between windward and leeward slopes. The precipitation simulation is poor for all terrains, with large precipitation thresholds in three typhoon cases, but for Typhoon HOTA, after improving the terrain, the model has the ability to forecast the heavy rainfall scenarios of the mountainous terrain.

The role of in situ ocean data assimilation in ECMWF subseasonal forecasts of sea‐surface temperature and mixed‐layer depth over the tropical Pacific ocean

AbstractThe tropical Pacific plays an important role in modulating the global climate through its prevailing sea‐surface temperature spatial structure and dominant climate modes like El Niño–Southern Oscillation (ENSO), the Madden–Julian Oscillation (MJO), and their teleconnections. These modes of variability, including their oceanic anomalies, are considered to provide sources of prediction skill on subseasonal timescales in the Tropics. Therefore, this study aims to examine how assimilating in situ ocean observations influences the initial ocean sea‐surface temperature (SST) and mixed‐layer depth (MLD) and their subseasonal forecasts. We analyze two subseasonal forecast systems generated with the European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), where the ocean states were initialized using two Observing‐System Experiment (OSE) reanalyses. We find that the SST differences between forecasts with and without ocean data assimilation grow with time, resulting in a reduced cold‐tongue bias when assimilating ocean observations. Two mechanisms related to air–sea coupling are considered to contribute to this growth of SST differences. One is a positive feedback between the zonal SST gradient, pressure gradient, and surface wind. The other is the difference in Ekman suction and mixing at the Equator due to surface wind‐speed differences. While the initial mixed‐layer depth (MLD) can be improved through ocean data assimilation, this improvement is not maintained in the forecasts. Instead, the MLD in both experiments shoals rapidly at the beginning of the forecast. These results emphasize how initialization and model biases influence air–sea interaction and the accuracy of subseasonal forecasts in the tropical Pacific.

Impact of assimilating NOAA VIIRS aerosol optical depth (AOD) observations on global AOD analysis from the Copernicus Atmosphere Monitoring Service (CAMS)

Abstract. Global monitoring of aerosols is required to analyse the impacts of aerosols on air quality and to understand their role in modulating the climate variability. The Copernicus Atmosphere Monitoring Service (CAMS) provides near-real-time forecasts and reanalyses of aerosols using the ECMWF Integrated Forecasting System (IFS), constrained by the assimilation of MODIS and Polar Multi-Sensor Aerosol Optical Properties (PMAp) aerosol optical depth (AOD). Given the potential end of lifetime of MODIS AOD, implementing new AOD observations in the CAMS operational suite is a priority to ensure the continuity of the CAMS forecast performances. The objective of this work is to test the assimilation of the NOAA VIIRS AOD product from S-NPP and NOAA20 satellites in the IFS model. Simulation experiments assimilating VIIRS on top or in place of MODIS were carried out from June 2021 to November 2021 to evaluate the impacts on the AOD analysis. For maritime aerosol background, the assimilation of VIIRS and the use of VIIRS from NOAA20 as an anchor reduce the analysis AOD values compared to MODIS-based experiments, in which the analysis values were too high due to the positive bias of MODIS/Terra over ocean. Over land, the assimilation of VIIRS induces a large increase in the analysis over biomass burning regions where VIIRS shows larger AOD than MODIS due to differences in the aerosol models and cloud filtering between MODIS and VIIRS retrieval algorithms. For dust source regions, the analysis is reduced when VIIRS is assimilated on top of or in place of MODIS, particularly over the Sahara, the Arabian Peninsula and a few places in Asia in the July–August period. The assimilation of VIIRS leads to an overall reduction of the bias in AOD analysis evaluated against AERONET measurements, with the largest bias reduction over Europe and desert and maritime sites.

Forecast evaluation of the North Pacific jet stream using AR Recon dropwindsondes

AbstractThe term jet stream generally refers to a narrow region of intense winds near the top of the midlatitude or subtropical troposphere. It is in the midlatitude jet stream where instabilities and waves may develop into synoptic‐scale systems, which in turn makes accurately resolving the structure of the jet stream and associated features critical for atmospheric development, predictability, and impacts, such as extreme precipitation and winds. Using dropwindsonde observations collected during the Atmospheric River Reconnaissance (AR Recon) campaign from 2020 to 2022, this study assesses the North Pacific jet stream structure in the European Centre for Medium‐Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS). Results show that the IFS has a slow‐wind bias on the lead times assessed, with the strongest winds (≥50 m·s−1) having a bias of up to −1.88 m·s−1 on forecast day 4. Also, the IFS cannot resolve the sharp potential vorticity (PV) gradient across the jet stream and tropopause, and this PV gradient weakens with forecast lead time. Cases with larger wind biases are characterized by higher PV biases and PV biases tend to be larger for cases with a higher horizontal PV gradient. These results suggest that further model‐based experiments are needed to identify and address these biases, which could ultimately yield increased forecast accuracy.

Evaluating the Precipitation Biases over the Western Periphery of the Sichuan Basin with the ECMWF Operational Forecast Model

Abstract Based on the hourly merged precipitation product, the performance of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecasting System (IFS) in simulating the diurnal variations of precipitation during warm season over the western periphery of the Sichuan basin (SCB) has been evaluated, and the underlying physical causes associated with the wet biases have also been investigated. The results show that the IFS well reproduces the spatial distributions of precipitation amount, frequency, and intensity over the SCB, as well as their diurnal variations, but the simulated precipitation peaks earlier than the observation with notable wet biases over the western periphery of the SCB. In addition, the strong wet biases exhibit notable regional differences over the western periphery of the SCB. The simulated wet biases over the southwestern periphery of the SCB expand westward to higher altitudes along the windward slope, with the maximum wet biases occurring at night. The westward expansion of the simulated stronger upward motions results in a westward shift of precipitation. However, the simulated precipitation over the northwestern periphery of the SCB has small difference in terms of the location; hence, the overestimated precipitation is associated with the stronger atmospheric instability, resulting from the higher potential temperature and the larger specific humidity near the surface. The findings revealed in this study indicate that the ECMWF forecast shows distinct uncertainties over the different complex terrain, and thus offers a promising way forward for improvements of model physical processes.

Parametric Postprocessing of Dual-Resolution Precipitation Forecasts

Abstract All major weather centers issue ensemble forecasts, which differ both in ensemble size and spatial resolution, even while covering the same domain. These parameters directly determine both the forecast skill of the prediction and the computation cost. In the past few years, the plans for upgrading the configuration of the Integrated Forecast System of the European Centre for Medium-Range Weather Forecasts (ECMWF) from a single forecast with 9-km resolution and a 51-member ensemble with 18-km resolution induced an extensive study of the forecast skill of both raw and postprocessed dual-resolution predictions comprising ensemble members of different horizontal resolutions. We investigate the predictive performance of the censored shifted gamma (CSG) ensemble model output statistic (EMOS) approach for statistical postprocessing with the help of dual-resolution, 24-h, precipitation accumulation ensemble forecasts over Europe with various forecast horizons. We consider the operational 50-member ECMWF ensemble as of high resolution and extend it with a low-resolution (29-km grid), 200-member experimental forecast. The investigated dual-resolution combinations consist of subsets of these two forecast ensembles with equal computational cost, which is equivalent to the cost of the operational ensemble. Our case study verifies that, compared with the raw ensemble combinations, EMOS postprocessing results in a significant improvement in forecast skill and that skill is statistically indistinguishable between any of the analyzed mixtures of dual-resolution combinations. Furthermore, the semilocally trained CSG EMOS provides an efficient alternative to the state-of-the-art quantile mapping without requiring additional historical data.

A mountain ridge model for quantifying oblique mountain wave propagation and distribution

Abstract. Following the current understanding of gravity waves (GWs) and especially mountain waves (MWs), they have a high potential for horizontal propagation from their source. This horizontal propagation and therefore the transport of energy is usually not well represented in MW parameterizations of numerical weather prediction and general circulation models. In this study, we present a mountain wave model (MWM) for the quantification of horizontal propagation of orographic gravity waves. This model determines MW source locations from topography data and estimates MW parameters from a fit of idealized Gaussian-shaped mountains to the elevation. Propagation and refraction of these MWs in the atmosphere are modeled using the Gravity-wave Regional Or Global Ray Tracer (GROGRAT). Ray tracing of each MW individually allows for an estimation of momentum transport due to both vertical and horizontal propagation. The MWM is a capable tool for the analysis of MW propagation and global MW activity and can support the understanding of observations and improvement of MW parameterizations in GCMs. This study presents the model itself and gives validations of MW-induced temperature perturbations to ECMWF Integrated Forecast System (IFS) numerical weather prediction data and estimations of GW momentum flux (GWMF) compared to HIgh Resolution Dynamics Limb Sounder (HIRDLS) satellite observations. The MWM is capable of reproducing the general features and amplitudes of both of these data sets and, in addition, is used to explain some observational features by investigating MW parameters along their trajectories.