European Centre For Medium-Range Weather Forecasts Research Articles

The Influence of Typhoon-Induced Wave on the Mesoscale Eddy

The strong wind-induced current and sea level have influences on the wave distribution in a tropical cyclone (TC). In particular, the wave–current interaction is significant in the period in which the TC passed the mesoscale eddy. In this study, the wave fields of Typhoon Chan-hom (2015) are hindcastly simulated using a coupled oceanic model that utilizes a nested triangle grid, i.e., the finite-volume community ocean model-simulating waves nearshore (FVCOM-SWAVE) model. The forcing wind field is composited from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis data and the simulation using a parametric Holland model, denoted as H-E. The open boundary fields include tide data from TPOX.5 and the hybrid coordinate ocean model (HYCOM) global datasets, including sea surface temperature (SST), sea surface salinity, sea surface current, and sea level data. The simulated oceanic parameters (e.g., the significant wave height, SWH) are validated against the measurements from the Jason-2 altimeter, yielding a root mean square error (RMSE) of 0.58 m for the SWH, a correlation (COR) coefficient of 0.94, and a scatter index (SI) of 0.23. Similarly, the simulated SSTs are compared with the remote sensing products of the remote sensing system (REMSS) and the measurements from Argos, yielding an RMSE of <0.8 °C, a COR of >0.95, and an SI of <0.04. The significant zonal asymmetry of the wave distribution along the typhoon track is observed. The Stokes drift is calculated from the FVCOM-SWAVE simulation results, and then the contribution of the Stokes transport is estimated using the Ekman–Stokes numbers. It is found that the ratio of the Stokes transport to the total net transport can reach >80% near the typhoon center, and the ratio is reduced to approximately <20% away from the typhoon center, indicating that Stokes transport is an essential aspect in the water mixing during a TC. The mesoscale eddies are detected by the sea level anomalies (SLA) fusion data from AVISO. It is found that the significant wave heights, Stokes drift, and Stokes transport inside the eddy area were higher than those outside the eddy area. These parameters inside the cold mesoscale eddies were higher than t inside the warm mesoscale eddies. Otherwise, the SST mainly increased within the cold mesoscale eddies area, while decreased within the warm mesoscale eddies area. The influence of mesoscale eddies on the SST was in proportion to the eddy radius and eddy EKE.

Comparison of Cloud/Rain Band Structures of Typhoon Muifa (2022) Revealed in FY‐3E MWHS‐2 Observations With All‐Sky Simulations

AbstractAll‐sky simulations of brightness temperature (TB) at Microwave Humidity Sounder‐2 (MWHS‐2) channel 15 from Fengyun‐3E (FY‐3E) are generated based on the European Center for Medium‐Range Weather Forecasts (ECMWF) reanalysis version 5 (ERA5) and the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) analysis for four typhoon cases. For Typhoon Muifa (2022), the MWHS‐2 TB observations have a strongest symmetric component within about 210‐km radial distances. Although much weaker, the symmetric component in the GFS all‐sky simulations compares more favorably with that of observations than the ERA5 all‐sky simulations. Due to ice scattering effects at 183.31 ± 7 GHz, spatial distribution patterns of TB observations and all‐sky simulations are quite similar to those of ice water path. The more cloud ice, the lower the TB. Deep convective systems within Muifa are captured by TB observations but not by simulations. Using an objective azimuthal‐spectral‐analysis based center‐fixing algorithm, the center positions of Typhoon Muifa can be well determined from both MWHS‐2 observations and GFS all‐sky simulations, but not from the ERA5 all‐sky TB simulations. Compared with the best track, the average center position errors for Muifa are 16.2, 24.9 and 74.4 km based on FY‐3E MWHS‐2 TB observations, the GFS all‐sky TB simulations, and the ERA5 all‐sky simulations, respectively. The conclusions of TC cloud/rain band structures and center positions of all‐sky simulations using the NCEP‐GFS analysis being better matching with those of FY‐3E MWHS‐2 observations than the ERA5 all‐sky simulations for Typhoon Muifa also hold true for three other typhoons.

Diversity, quality, and quantity of real ship data on the black-box and gray-box prediction models of ship fuel consumption

Establishing an efficient and reliable fuel consumption prediction model by using big data from ships facilitates optimized decision-making and ensures the green and intelligent development of ships. As data are crucial in model construction, this study presents a general multi-source data processing method for obtaining high-quality training data. A long short-term memory (LSTM) neural network, suitable for time-series data, was used to develop the fuel consumption black-box model. This was combined with the ship theory to establish a fuel consumption theoretical model, thereby generating the LSTM based gray-box model. We explored the impact of data diversity, quality, and quantity on black-box and gray-box models. Analysis of the navigation data of a passenger ship and meteorological data collected from the European Centre for Medium-Range Weather Forecasts (ECMWF) and MeteoBlue indicated that the combination of variables obtained via the feature selection of the least absolute shrinkage and selection operator (LASSO) statistical method yielded the best overall prediction performance. Moreover, the gray-box model was relatively stable in terms of the changes in effective variables. An analysis of data quality revealed that the systematic processing of outliers, which improves the accuracy of both models by 6.19% compared with direct deletion. Furthermore, the gray-box models use less amounts of data than the black-box models to achieve higher accuracy.

The Relationship of CLWC and Rainfall to the Synoptic Cases of Two Case Studies Over Iraq

Climate change and the improvement of the environment is one of the most important topics that preoccupied the world, and the main factor for it is precipitation. To understand how to improve the environment, precipitation must be increased by understanding the mechanisms of cloud formation and benefiting from the water content of the liquid cloud, which is the basic building block of precipitation. There is no precipitation without liquid water content. This study examined the effect of depressions on rain and liquid cloud water content, in addition to studying the relationship of liquid water content and precipitation, and classifying clouds according to liquid water content over Iraq. We found that the liquid water content and the amount of precipitation are affected by depressions. Hourly data has been taken cloud liquid water content were used for four primary times, (00, 06, 12, and 18 UTC), at all pressure levels from (1000 to 175 hpa) for the study Selected days in Winter 2019 Taken from the generation of European Center for Medium-Range Weather Forecasting (ECMWF) from ERA5 it is the fifth generation ECMWF atmospheric reanalysis of the global climate covering the period from January 1950 to present. ERA5 is produced by the Copernicus Climate Change Service at ECMWF, sea level pressure maps and geopotential height maps were taken from ECMWF ERA5 Reanalysis. And low, medium, and high clouds and the amount of precipitation were taken from the (Iraqi Meteorological Organization and Seismic Monitoring). Two programs (Matlab) and (Origin) were used to analyze and plot the data.

Impacts of Dropsonde Observations on Forecasts of Atmospheric Rivers and Associated Precipitation in the NCEP GFS and ECMWF IFS Models

Abstract Atmospheric River Reconnaissance has held field campaigns during cool seasons since 2016. These campaigns have provided thousands of dropsonde data profiles, which are assimilated into multiple global operational numerical weather prediction models. Data denial experiments, conducted by running a parallel set of forecasts that exclude the dropsonde information, allow testing of the impact of the dropsonde data on model analyses and the subsequent forecasts. Here, we investigate the differences in skill between the control forecasts (with dropsonde data assimilated) and denial forecasts (without dropsonde data assimilated) in terms of both precipitation and integrated vapor transport (IVT) at multiple thresholds. The differences are considered in the times and locations where there is a reasonable expectation of influence of an intensive observation period (IOP). Results for 2019 and 2020 from both the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Centers for Environmental Prediction (NCEP) global model show improvements with the added information from the dropsondes. In particular, significant improvements in the control forecast IVT generally occur in both models, especially at higher values. Significant improvements in the control forecast precipitation also generally occur in both models, but the improvements vary depending on the lead time and metrics used. Significance Statement Atmospheric River Reconnaissance is a program that uses targeted aircraft flights over the northeast Pacific to take measurements of meteorological fields. These data are then ingested into global weather models with the intent of improving the initial conditions and resulting forecasts along the U.S. West Coast. The impacts of these observations on two global numerical weather models were investigated to determine their influence on the forecasts. The integrated vapor transport, a measure of both wind and humidity, saw significant improvements in both models with the additional observations. Precipitation forecasts were also improved, but with differing results between the two models.

Forecasting large hail and lightning using additive logistic regression models and the ECMWF reforecasts

Abstract. Additive logistic regression models for lightning (ARlig) and large hail (ARhail) were developed using convective parameters from the ERA5 reanalysis, hail reports from the European Severe Weather Database (ESWD), and lightning observations from the Met Office Arrival Time Difference network (ATDnet). The models yield the probability of lightning and large hail in a given timeframe over a particular grid point. To explore the value of this approach to medium-range forecasting, the models were applied to the European Centre for Medium Range Weather Forecasts (ECMWF) reforecasts to reconstruct probabilistic lightning and large hail forecasts for 11 ensemble members, from 2008 to 2019 and for lead times up to 228 h. The lightning and large hail models were based on different predictor parameters: most unstable convective available potential energy (CAPE), 925–500 hPa bulk shear, mixed layer mixing ratio, wet bulb zero height (for large hail), most unstable lifted index, mean relative humidity between 850 and 500 hPa, 1 hourly accumulated convective precipitation and specific humidity at 925 hPa (for lightning). First, we compared the lightning and hail ensemble forecasts for different lead times with observed lightning and hail focusing on a recent hail outbreak. Second, we evaluated the predictive skill of the model as a function of forecast lead time using the area under the ROC curve (AUC) as a validation score. This analysis showed that ARhail has a very high predictive skill (AUC > 0.95) for a lead time up to 60 h. ARhail retains a high predictive skill even for extended forecasts (AUC = 0.86 at 180 h lead time). Although ARlig exhibits a lower predictive skill than ARhail, lightning forecasts are also skilful both in the short term (AUC = 0.92 at 60 h) and in the medium range (AUC = 0.82 at 180 h). Finally, we compared the performance of the 4-dimensional hail model with that of composite parameters such as the significant hail parameter (SHP) or the product of CAPE and the 925–500 hPa bulk shear (CAPESHEAR). Results show that ARhail outperforms CAPESHEAR at all lead times and SHP at short-to-medium lead times. These findings suggests that the combination of additive logistic regression models and ECMWF ensemble forecasts can create highly skilful medium-range hail and lightning forecasts for Europe.

Evaluation of Tropopause Height from Sentinel-6 GNSS Radio Occultation Using Different Methods

The tropopause is described as the distinction between the troposphere and the stratosphere, and the tropopause height (TPH) is an indicator of climate change. GNSS Radio Occultation (RO) can monitor the atmosphere globally under all weather conditions with a high vertical resolution. In this study, four different techniques for identifying the TPH were investigated. The lapse rate tropopause (LRT) and cold point tropopause (CPT) methods are the traditional methods for determining the TPH based on temperature profiles according to the World Meteorological Organization (WMO) definition. Two advanced methods based on the covariance transform (CT) method are used to estimate the TPH from the refractivity (TPHN) and the TPH from the bending angle (TPHα). Data from the Sentinel-6 satellite were used to evaluate the different algorithms for the identification of the TPH. The analysis shows that the CPT height is greater than the LRT height and that the CPT is only valid in tropical regions. The CPT height, TPHN, and TPHα were compared with the LRT height. In general, the TPHα had the largest value, followed by the TPHN, and the LRT had the lowest value of TPH at and near the equator. In the equatorial region, the maximum TPH results from the TPHα (approximately 17.5 km), and at the poles, the minimum TPH results from the LRT (approximately 9 km). The results were also compared with the European Center for Medium-Range Weather Forecasts (ECMWF), and there was a strong correlation of 0.999 between the different approaches for identifying the TPH from RO and the ECMWF model. The identification of the TPH is critical for the transfer of mass, water, and trace gases between the troposphere and stratosphere.

Comprehensive evaluation of typical planetary boundary layer (PBL) parameterization schemes in China – Part 2: Influence of uncertainty factors

Abstract. This study focuses on the uncertainties that influence numerical simulation results of meteorological fields (horizontal resolution: 75, 15, and 3 km; vertical resolution: 48 and 62 levels; near-surface (N-S) scheme: MM5 and Eta schemes; initial and boundary conditions: Final (FNL) and European Center for Medium-Range Weather Forecasting (ECMWF) reanalysis data; underlying surface update: model default and latest updates; model version: version 3.6.1 and 3.9.1). By further evaluating and analyzing the uncertainty factors, it is expected to provide relevance for those scholars devoted to factor analysis in order to make the results closer to the observed values. In this study, a total of 12 experiments are set up to analyze the effects of the uncertainties mentioned above, and the following conclusions are drawn: (1) horizontal resolution has a greater effect than vertical resolution; (2) the simulated effects of temperature and wind speed in the N-S scheme are smaller than those in the planetary boundary layer (PBL) scheme; (3) the initial and boundary conditions of different products have the most remarkable effect on relative humidity, while the simulation results of ECMWF data are the best; (4) the updates with urban and water bodies as the underlying surface have a more significant contribution to the meteorological fields, especially on temperature; and (5) for the PBL parameterization schemes, the update of the model version has less impact on the simulation results because each update has small changes and no major changes overall. In general, the configuration of uncertainties needs to be considered comprehensively according to what you need in order to obtain the best simulation results.

The Atmospheric and Oceanic Processes on Thermal Front Variability over the Java Sea

The Java Sea is influenced by various atmospheric and oceanographic factors such as monsoon winds and the inflow of water from adjacent seas, thus leading to the formation of thermal fronts. In this research, the atmospheric process over the thermal front was estimated by using sea surface temperature (SST) gradient and crosswind-SST gradient calculation based on the GLORYS12V1 monthly SST data from Copernicus Marine Service (CMEMS) and the ERA5 monthly surface wind data from European Centre for Medium-Range Weather Forecasts (ECMWF) for 27 years (from 1993 to 2019). Moreover, the oceanic process was determined from the volume transports of three major channels surrounding the Java Sea, which are Karimata Strait, Sunda Strait, and the eastern boundary of the Java Sea. Based on the annual variance analysis, there are four main thermal front areas, such as Northern Java Coast (NJC), Eastern Sumatra Coast (ESC), Western Borneo Coast (WBC), and Eastern Borneo Coast (EBC). The annual variation of the thermal front over the ESC, WBC, and ESC has two peaks in March and October, while the NJC thermal front area has three peaks in March, July, and November. This study results that the thermal front activities over the NJC, WBC, and EBC are significantly correlated to the wind stress curl over those areas, with a coefficient correlation of about 0.62, -0.92, and -0.73, respectively. In addition, the increase in thermal front activity over the NJC area is controlled by negative wind stress curl with no lag time, while the WBC and EBC areas are controlled by positive wind stress curl with a lag time of 1 month (wind stress leads). On the other hand, the thermal front activity over the ESC area is closely related to southward volume transport from the Karimata Strait, with a correlation coefficient of -0.74 and a lag time of 2 months (volume transport leads).

Estimation and Evaluation of Zenith Tropospheric Delay from Single and Multiple GNSS Observations

Multi-Global Navigation Satellite Systems (multi-GNSS) (including GPS, BDS, Galileo, and GLONASS) provide a significant opportunity for high-quality zenith tropospheric delay estimation and its applications in meteorology. However, the performance of zenith total delay (ZTD) retrieval from single- or multi-GNSS observations is not clear, particularly from the new, fully operating BDS-3. In this paper, zenith tropospheric delay is estimated using the single-, dual-, triple-, or four-GNSS Precise Point Positioning (PPP) technique from 55 Multi-GNSS Experiment (MGEX) stations over one year. The performance of GNSS ZTD estimation is evaluated using the International GNSS Service (IGS) standard tropospheric products, radiosonde, and the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA5). The results show that the GPS-derived ZTD time series is more consistent and reliable than those derived from BDS-only, Galileo-only, and GLONASS-only solutions. The performance of the single-GNSS ZTD solution can be enhanced with better accuracy and stability by combining multi-GNSS observations. The accuracy of the ZTD from multi-GNSS observations is improved by 13.8%, 43.8%, 27.6%, and 22.9% with respect to IGS products for the single-system solution (GPS, BDS, Galileo, and GLONASS), respectively. The ZTD from multi-GNSS observations presents higher accuracy and a significant improvement with respect to radiosonde and ERA5 data when compared to the single-system solution.

FuXi: a cascade machine learning forecasting system for 15-day global weather forecast

Over the past few years, the rapid development of machine learning (ML) models for weather forecasting has led to state-of-the-art ML models that have superior performance compared to the European Centre for Medium-Range Weather Forecasts (ECMWF)’s high-resolution forecast (HRES), which is widely considered as the world’s best physics-based weather forecasting system. Specifically, ML models have outperformed HRES in 10-day forecasts with a spatial resolution of 0.25∘. However, the challenge remains in mitigating the accumulation of forecast errors for longer effective forecasts, such as achieving comparable performance to the ECMWF ensemble in 15-day forecasts. Despite various efforts to reduce accumulation errors, such as implementing autoregressive multi-time step loss, relying on a single model has been found to be insufficient for achieving optimal performance in both short and long lead times. Therefore, we present FuXi, a cascaded ML weather forecasting system that provides 15-day global forecasts at a temporal resolution of 6 hours and a spatial resolution of 0.25∘. FuXi is developed using 39 years of the ECMWF ERA5 reanalysis dataset. The performance evaluation demonstrates that FuXi has forecast performance comparable to ECMWF ensemble mean (EM) in 15-day forecasts. FuXi surpasses the skillful forecast lead time achieved by ECMWF HRES by extending the lead time for Z500 from 9.25 to 10.5 days and for T2M from 10 to 14.5 days. Moreover, the FuXi ensemble is created by perturbing initial conditions and model parameters, enabling it to provide forecast uncertainty and demonstrating promising results when compared to the ECMWF ensemble.

Drought Assessment in Greece Using SPI and ERA5 Climate Reanalysis Data

The present work aims to assess the spatial variability and the trends of the annual rainfall and meteorological drought in the entire territory of Greece utilising the ERA5 reanalysis precipitation dataset of the European Centre for Medium-Range Weather Forecasts (ECMWF), which spans from January 1940 to December 2022 (an 83-year period). Drought assessment took place based on the Standardized Precipitation Index (SPI) for timescales ranging from 1 month to 12 months. Evaluation was carried out by calculating SPI using observed rainfall data from five meteorological stations. The annual rainfall and drought severity trends for timescales of 1 (SPI-1), 3 (SPI-3), 6 (SPI-6) and 12 (SPI-12) months were analysed using the Theil–Sen slope method and the Mann–Kendall trend test. The results indicate significant, both increasing and decreasing, annual precipitation trends at the 95% significance level for the Aegean Islands, western Crete and western mainland of Greece. The results also indicate significant drought trends for SPI-12 for the Aegean Islands and western Peloponnese. Trend analysis for SPI-1, SPI-3 and SPI-6 indicate a mixture of non-significantly increasing wetting trends and increasing drought trends at the national scale. In conclusion, the ERA5 dataset seems to be a valuable tool for drought monitoring at the spatial scale.

Characterizing the tropospheric water vapor spatial variation and trend using 2007–2018 COSMIC radio occultation and ECMWF reanalysis data

Abstract. Atmospheric water vapor plays a crucial role in the global energy balance, hydrological cycle, and climate system. High-quality and consistent water vapor data from different sources are vital for weather prediction and climate research. This study assesses the consistency between the Formosa Satellite Mission 3–Constellation Observing System for Meteorology, Ionosphere, and Climate (FORMOSAT-3/COSMIC) radio occultation (RO) and European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis Model 5 (ERA5) water vapor datasets. Comparisons are made across different atmospheric pressure levels (300, 500, and 850 hPa) from 2007 to 2018. Generally, the two datasets show good spatial and temporal agreement. COSMIC's global water vapor retrieval is slightly lower than ERA5's at 500 and 850 hPa, with distinct latitudinal differences between hemispheres. COSMIC exhibits global water vapor increasing trends of 3.47 ± 1.77 % per decade, 3.25 ± 1.25 % per decade, and 2.03 ± 0.65 % per decade at 300, 500, and 850 hPa, respectively. Significant regional variability in water vapor trends, encompassing notable increasing and decreasing patterns, is observable in tropical and subtropical regions. At 500 and 850 hPa, strong water vapor increasing trends are noted in the equatorial Pacific Ocean and the Laccadive Sea, while decreasing trends are evident in the Indo-Pacific Ocean region and the Arabian Sea. Over land, substantial increasing trends at 850 hPa are observed in the southern United States, contrasting with decreasing trends in southern Africa and Australia. The differences between the water vapor trends of COSMIC and ERA5 are primarily negative in the tropical regions at 850 hPa. However, the water vapor increasing trends at 850 hPa estimated from COSMIC are significantly higher than the ones derived from ERA5 data for two low-height stratocumulus-cloud-rich ocean regions west of Africa and South America. These regions with notable water vapor trend differences are located in the Intertropical Convergence Zone (ITCZ) area with frequent occurrences of convection, such as deep clouds. The difference in characterizing water vapor distribution between RO and ERA5 in deep cloud regions may cause such trend differences. The assessment of spatiotemporal variability in RO-derived water vapor and reanalysis of atmospheric water vapor data helps ensure the quality of these datasets for climate studies.

Development and assessment of artificial neural network models for direct normal solar irradiance forecasting using operational numerical weather prediction data

Accurate operational solar irradiance forecasts are crucial for better decision making by solar energy system operators due to the variability of resource and energy demand. Although numerical weather prediction (NWP) models can forecast solar radiation variables, they often have significant errors, particularly in the direct normal irradiance (DNI), which is especially affected by the type and concentration of aerosols and clouds. This paper presents a method based on artificial neural networks (ANN) for generating operational DNI forecasts using weather and aerosol forecasts from the European Center for Medium-range Weather Forecasts (ECMWF) and the Copernicus Atmospheric Monitoring Service (CAMS), respectively. Two ANN models were designed: one uses as input the predicted weather and aerosol variables for a given instant, while the other uses a period of the improved DNI forecasts before the forecasted instant. The models were developed using observations for the location of Évora, Portugal, resulting in 10 min DNI forecasts that for day 1 of forecast horizon showed an improvement over the downscaled original forecasts regarding R2, MAE and RMSE of 0.0646, 21.1 W/m2 and 27.9 W/m2, respectively. The model was also evaluated for different timesteps and locations in southern Portugal, providing good agreement with experimental data.

Subseasonal Variation Characteristics of Low-Cloud Fraction in Southeastern and Northwestern North Pacific

The subseasonal variability of the low-cloud fraction (LCF) over the southeastern North Pacific (SENP) and northwestern North Pacific (NWNP) was studied using satellite observations and the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis. It is found that subseasonal variability of the LCF was closely related to variations in the estimated inversion strength (EIS), sea surface wind speed (SSW), sensible heat flux (SHF), sea surface temperature (SST), surface temperature advection (Tadv), relative humidity (RH), surface level pressure (SLP) and surface air temperature (SAT). An increase in the LCF over the SENP is associated with the development of an anomalous anticyclonic circulation, which is located on the west coast of America. The cold advection, together with the subsidence warming associated with the anticyclonic circulation, strengthens the temperature inversion, favoring the development of the LCF. In the NWNP, the maximum LCF anomaly was also correlated with the stable boundary layer. The southerly wind blows airflow over the Kuroshio Extension from the subtropics, which brings warm and moist air. When air flows to the colder sea surface, it is cooled and condensed by the intensified heat exchange. A lead-lag composite analysis indicates that the mechanisms are different between the SENP and the NWNP, possibly due to the different types of low-level clouds over these two regions. In the SENP, the trade cumulus dominates under a strong capping inversion over the subtropics, whereas fog and stratus often occur under a shallow capping inversion in the NWNP. The effects of atmospheric circulation are also discussed.

Brief communication: On the definition of the low-level jet

Abstract. Low-level jets (LLJs) are examples of non-logarithmic wind speed profiles affecting wind turbine power production, wake recovery, and structural/aerodynamic loading. However, there is no consensus regarding which definition should be applied for jet identification. In this study we argue that a shear definition is more relevant to wind energy than a falloff definition. The shear definition is demonstrated and validated through the development of a European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation reanalysis (ERA5) LLJ climatology for six sites. Identification of LLJs and their morphology, frequency, and intensity is critically dependent on the (i) vertical window of data from which LLJs are extracted and (ii) the definition employed.

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.

GNSS radio occultation excess-phase processing for climate applications including uncertainty estimation

Abstract. Earth observation from space provides a highly valuable basis for atmospheric and climate science, in particular also through climate benchmark data from suitable remote sensing techniques. Measurements by global navigation satellite system (GNSS) radio occultation (RO) qualify to produce such benchmark data records as they globally provide accurate and long-term stable datasets for essential climate variables (ECVs) such as temperature. This requires a rigorous processing of the raw RO measurements to ECVs, with narrow uncertainties. In order to fully exploit this potential, Wegener Center's Reference Occultation Processing System (rOPS) Level 1a (L1a) processing subsystem includes uncertainty estimation in both precise orbit determination (POD) and excess-phase profile derivation. Here we introduce the new rOPS L1a excess-phase processing, the first step in the RO profiles retrieval down to atmospheric profiles, which extracts the atmospheric excess phase from raw SI-traceable RO measurements. This excess-phase processing, for itself algorithmically concise, includes integrated quality control and uncertainty estimation, requiring a complex framework of various subsystems that we first introduce before describing the implementation of the core algorithms. The quality control and uncertainty estimation, computed per RO event, are supported by reliable forward-modeled excess-phase profiles based on the POD orbit arcs and collocated short-range forecast profiles of the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis (ERA5). The quality control removes or alternatively flags excess-phase profiles of insufficient or degraded quality. The uncertainty estimation accounts both for relevant random- and systematic-uncertainty components, and the resulting (total) uncertainty profiles serve as a starting point for the subsequent uncertainty propagation through the retrieval processing chain down to the atmospheric ECV profiles. We also evaluated the quality and reliability of the resulting excess-phase profiles based on Metop-A/B/C (Meteorological Operational) RO datasets for three 3-month periods in 2008, 2013, and 2020 by way of a sensitivity analysis for three representative atmospheric layers (tropo-, strato-, mesosphere), investigating consistency with ERA5-derived profiles, influences of different orbit and clock inputs, and consistency across the different Metop satellites. These consistencies range from centimeter to submillimeter levels, indicating that the new processing can provide highly accurate and robust excess-phase profiles. Furthermore, cross-evaluation and intercomparison with excess-phase data from the established data providers EUMETSAT (European Organisation for the Exploitation of Meteorological Satellites) and UCAR (University Corporation for Atmospheric Research) revealed subtle discrepancies but overall very close agreement, with larger differences compared to UCAR in the boundary layer. The new rOPS L1a processing can hence be considered capable of producing reliable long-term data records including uncertainty estimation for the benefit of climate applications.

Adaptive selection of members for convective-permitting regional ensemble prediction over the western Maritime Continent

A common issue faced by the downscaled regional ensemble prediction systems is the under-dispersiveness of the ensemble forecasts, often attributed to the lack of spread under the initial conditions from the global ensemble. In this study, a novel method that adopts an adaptive approach to selecting global ensemble members for regional downscaling has been developed. Instead of using a fixed set of pre-selected global ensemble members, the adaptive selection performs a sampling algorithm and selects the global ensemble members, which maximizes a fractions skill score (FSS)-based displacement between ensemble members. The method is applied to a convective-permitting ensemble prediction system over the western Maritime Continent, referred to as SINGV-EPS. SINGV-EPS has a grid spacing of 4.5 km and is a 12-member ensemble that is driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) 51-member global ensemble. Month-long trials were conducted in June 2020 to assess the impact of adaptive selection on the ensemble forecast spread and rainfall verification scores. In both fixed pre-selection and adaptive selection experiments, SINGV-EPS was still under-dispersive. However, adaptive selection improved the ensemble spread and reduced the root-mean-square error (RMSE) of the ensemble mean in wind, temperature, and precipitation fields. Further verification of the rainfall forecasts showed that there was a reduction in the Brier score and a higher hit rate in the relative operating characteristic (ROC) curve for all rainfall thresholds when adaptive selection was applied. Additionally, the ensemble mean forecasts from adaptive selection experiments are more accurate beyond 24 h, with a higher FSS for all rainfall thresholds and neighborhood lengths. These results suggest that the adaptive selection is superior to the fixed pre-selection of global ensemble members for downscaled regional ensemble prediction.

On the use of consistent bias corrections to enhance the impact of AeolusLevel‐2B Rayleigh winds on National Oceanic and Atmospheric Administration global forecast skill

AbstractThe operational Aeolus Level‐2B (L2B) horizontal line‐of‐sight (HLOS) retrieved Rayleigh winds, produced by the European Space Agency (ESA), utilize European Centre for Medium‐Range Weather Forecasts (ECMWF) short‐term forecasts of temperature, pressure, and horizontal winds in the Rayleigh–Brillouin and M1 correction procedures. These model fields or backgrounds can contain ECMWF model‐specific errors, which may propagate to the retrieved Rayleigh winds. This study examines the sensitivity of the retrieved Rayleigh winds to the changes in the model backgrounds, and the potential benefit of using the same system, in this case the National Oceanic and Atmospheric Administration's Finite‐Volume Cubed Sphere Global Forecast System (FV3GFS), for both the corrections and the data assimilation and forecast procedures. It is shown that the differences in the model backgrounds (FV3GFS minus ECMWF) can propagate through the Level‐2B horizontal line‐of‐sight Rayleigh wind retrieval process, mainly the M1 correction, resulting in differences in the retrieved Rayleigh winds with mean and standard deviation of magnitude as large as 0.2 m·s−1. The differences reach up to 0.4, 0.6, and 0.7 m·s−1 for the 95th, 99th, and 99.5th percentiles of the sample distribution with maxima of ∼1.4 m·s−1. The numbers of the large differences for the combined lower and upper 5th, 1st, and 0.5th percentile pairs are ∼6,100, 1,220, and 610 between 2.5 and 25 km height globally per day respectively. The ESA‐disseminated Rayleigh wind product (based on the ECMWF corrections) already shows a significant positive impact on the FV3GFS global forecasts. In the observing system experiments performed, compared with the ESA Rayleigh winds, the use of the FV3GFS‐corrected Rayleigh winds lead to ∼0.5% more Rayleigh winds assimilated in the lower troposphere and show enhanced positive impact on FV3GFS forecasts at the day 1–10 range but limited to the Southern Hemisphere.