European Centre For Medium-Range Weather Forecasts Research Articles

Improving dynamical-statistical subseasonal precipitation forecasts using deep learning: A case study in Southwest China

Subseasonal precipitation forecasting is challenging but critical for water management, energy supply, and disaster prevention. To improve regional subseasonal precipitation prediction, previous studies have proposed a dynamical-statistical projection model (DSPM). In this study, we develop a new method that combines the DSPM and deep learning (DL), called the DL-DSPM. The DSPM is developed using the observed relationships between large-scale atmospheric circulations and regional precipitation, and the dynamical forecasted atmospheric circulations from the European Centre for Medium-Range Weather Forecasts (ECMWF) model. The DL-DSPM improves upon the DSPM by correcting biases in atmospheric circulation forecasts from the ECMWF model using two DL models, namely, residual network and U-Net models. In the case of Southwest China (SWC), DL models can improve atmospheric circulation forecasts at lead times beyond 5 pentads, including large-scale drivers of SWC precipitation variability. The DL-DSPM outperforms the ECMWF model and DSPM forecasts in predicting precipitation anomalies beyond 4 and 5 pentads over most SWC regions, respectively. In addition, the DL-DSPM is more skillful than the ECMWF model and DSPM in predicting extreme precipitation events more than 4 pentads in advance. The successful combination of DL and the DSPM provides a new possible direction for DL applications in subseasonal precipitation forecasting.

Wave retrieval for Sentinel-1 synthetic aperture radar under complex sea state

ABSTRACT In this paper, sea surface waves during complex sea state (i.e. tropical cyclones (TCs) and Kuroshio current) are inverted from dual-polarized (vertical – vertical (VV) and vertical – horizontal (VH)) synthetic aperture radar (SAR) images using an existing algorithm, namely, the parameterized first-guess spectrum method (PFSM). In particular, the wind term is included in two modulation transfer functions (MTFs), i.e. tilt modulation and velocity bunching. More than 3000 Sentinel (S-1) images acquired in interferometric wide (IW) mode and extra wide (EW) mode are collected, and the operational CyclObs wind product provided by the French Research Institute for Exploitation of the Oceans (IFREMER) is collected corresponding to the TC images. Note that these images are located within 500 km away from TC eyes. The wave spectra during the TCs are then inverted using the PFSM algorithm and the SAR-derived wind product. Simultaneously, the 0.05° gridded wave spectra collocated with these images are simulated using a third-generation numeric wave model, called WAVEWATCH-III (WW3), in which the reconstructed wind based on the European Centre for Medium-Range Weather Forecasts (ECMWF), the Copernicus Marine Environment Monitoring Service (CMEMS) current, and the CMEMS sea level are treated as the forcing field. Comparison of the SAR wave retrievals and simulations reveals that the root mean square errors (RMSEs) of the significant wave height (SWH) and mean wave period (MWP) are 0.55 m and 0.43 s, their correlation coefficients (COR) are 0.94 and 0.92, and their scatter index (SI) values are 0.17 and 0.07. These results demonstrate the better performance of the PFSM algorithm using theoretical MTFs, i.e. an RMSE of 0.68 m and COR of 0.90 for the SWH and an RMSE of 0.60 s and a COR of 0.87 for the MWP. Although the accuracy of the wave retrieval from the SAR image under complex sea state is improved when implementing the tilt and velocity bunching MTFs with the strong wind term, SWH retrievals by algorithm PFSM suffer saturation problem at SWH > 8 m.

Performance Evaluation of a National Seven-Day Ensemble Streamflow Forecast Service for Australia

The Australian Bureau of Meteorology offers a national operational 7-day ensemble streamflow forecast service covering regions of high environmental, economic, and social significance. This semi-automated service generates streamflow forecasts every morning and is seamlessly integrated into the Bureau’s Hydrologic Forecasting System (HyFS). Ensemble rainfall forecasts, European Centre for Medium-Range Weather Forecasts (ECMWF), and Poor Man’s Ensemble (PME), available in the Numerical Weather Prediction (NWP) suite, are used to generate these streamflow forecasts. The NWP rainfall undergoes pre-processing using the Catchment Hydrologic Pre-Processer (CHyPP) before being fed into the GR4H rainfall–runoff model, which is embedded in the Short-term Water Information Forecasting Tools (SWIFT) hydrological modelling package. The simulated streamflow is then post-processed using Error Representation and Reduction In Stages (ERRIS). We evaluated the performance of the operational rainfall and streamflow forecasts for 96 catchments using four years of operational data between January 2020 and December 2023. Performance evaluation metrics included the following: CRPS, relative CRPS, CRPSS, and PIT-Alpha for ensemble forecasts; NSE, PCC, MAE, KGE, PBias, and RMSE; and three categorical metrics, CSI, FAR, and POD, for deterministic forecasts. The skill scores, CRPS, relative CRPS, CRPSS, and PIT-Alpha, gradually decreased for both rainfall and streamflow as the forecast horizon increased from Day 1 to Day 7. A similar pattern emerged for NSE, KGE, PCC, MAE, and RMSE as well as for the categorical metrics. Forecast performance also progressively decreased with higher streamflow volumes. Most catchments showed positive performance skills, meaning the ensemble forecast outperformed climatology. Both streamflow and rainfall forecast skills varied spatially across the country—they were generally better in the high-runoff-generating catchments, and poorer in the drier catchments situated in the western part of the Great Dividing Range, South Australia, and the mid-west of Western Australia. We did not find any association between the model forecast skill and the catchment area. Our findings demonstrate that the 7-day ensemble streamflow forecasting service is robust and draws great confidence from agencies that use these forecasts to support decisions around water resource management.

Accelerating Lagrangian transport simulations on graphics processing units: performance optimizations of Massive-Parallel Trajectory Calculations (MPTRAC) v2.6

Abstract. Lagrangian particle dispersion models are indispensable tools for the study of atmospheric transport processes. However, Lagrangian transport simulations can become numerically expensive when large numbers of air parcels are involved. To accelerate these simulations, we made considerable efforts to port the Massive-Parallel Trajectory Calculations (MPTRAC) model to graphics processing units (GPUs). Here we discuss performance optimizations of the major bottleneck of the GPU code of MPTRAC, the advection kernel. Timeline, roofline, and memory analyses of the baseline GPU code revealed that the application is memory-bound, and performance suffers from near-random memory access patterns. By changing the data structure of the horizontal wind and vertical velocity fields of the global meteorological data driving the simulations from structure of arrays (SoAs) to array of structures (AoSs) and by introducing a sorting method for better memory alignment of the particle data, performance was greatly improved. We evaluated the performance on NVIDIA A100 GPUs of the Jülich Wizard for European Leadership Science (JUWELS) Booster module at the Jülich Supercomputing Center, Germany. For our largest test case, transport simulations with 108 particles driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis, we found that the runtime for the full set of physics computations was reduced by 75 %, including a reduction of 85 % for the advection kernel. In addition to demonstrating the benefits of code optimization for GPUs, we show that the runtime of central processing unit (CPU-)only simulations is also improved. For our largest test case, we found a runtime reduction of 34 % for the physics computations, including a reduction of 65 % for the advection kernel. The code optimizations discussed here bring the MPTRAC model closer to applications on upcoming exascale high-performance computing systems and will also be of interest for optimizing the performance of other models using particle methods.

Sun from Behind the Clouds: the Appeals Board of the European Centre for Medium-Range Weather Forecasts

Abstract One of the lesser-known regional international organizations in Europe is the European Centre for Medium-Range Weather Forecasts (ECMWF). It has a mandate to undertake tasks such as the production of weather forecasts, conduct scientific research to improve meteorological projections, maintain extensive datasets, and to better coordinate European meteorological infrastructure. The Centre has three sites across three European states where it has over 400 staff, and given it is immune from suit, its staff members only have an internal justice mechanism to resort to in the event of employment/pensions disputes. Little is known about the ECMWF Appeals Board given its untransparent nature. As uncovered in this article, it has undergone partial reform, and simultaneously, there has been an uptick in the number of cases. This article analyses these developments and highlights that more must be done to align the ECMWF Appeals Boards with the features of other international administrative tribunals.

Improving forecasts of precipitation extremes over northern and central Italy using machine learning

AbstractThe accurate prediction of intense precipitation events is one of the main objectives of operational weather services. This task is even more relevant nowadays, with the rapid progression of global warming which intensifies these events. Numerical weather prediction models have improved continuously over time, providing uncertainty estimation with dynamical ensembles. However, direct precipitation forecasting is still challenging. Greater availability of machine‐learning tools paves the way to a hybrid forecasting approach, with the optimal combination of physical models, event statistics, and user‐oriented postprocessing. Here we describe a specific chain, based on a random‐forest (RF) pipeline, specialised in recognising favourable synoptic conditions leading to precipitation extremes and subsequently classifying extremes into predefined types. The application focuses on northern and central Italy, taken as a testbed region, but is seamlessly extensible to other regions and time‐scales. The system is called MaLCoX (Machine Learning model predicting Conditions for eXtreme precipitation) and is running daily at the Italian regional weather service of ARPAE Emilia‐Romagna. MalCoX has been trained with the ARCIS gridded high‐resolution precipitation dataset as the target truth, using the last 20 years of the European Centre for Medium‐Range Weather Forecasts (ECMWF) reforecast dataset as input predictors. We show that, with a long enough training period, the optimal blend of larger‐scale information with direct model output improves the probabilistic forecast accuracy of extremes in the medium range. In addition, with specific methods, we provide a useful diagnostic to convey to forecasters the underlying physical storyline which makes a meteorological event extreme.

Using CloudSat to Advance the Global Precipitation Climatology Project (GPCP) over Antarctica

Remote sensing-based precipitation products face several challenges in high latitudes and specifically over frozen surfaces (i.e., snow and ice). Consequently, precipitation estimates tend to be lower in quality over these regions, including Antarctica, the coldest continent on Earth. In this study, we developed a method for adjusting precipitation estimates over Antarctica by leveraging CloudSat's ability to capture snowfall compared to other satellite products over snow and ice surfaces. We addressed limitations of CloudSat, such as poor spatiotemporal sampling, noise, and incomplete coverage near the poles. We utilized the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th Generation (ERA5) to guide development, particularly during the period when the TIROS Operational Vertical Sounder (TOVS) was contributing. We adjusted monthly and seasonal TOVS and Atmospheric Infrared Sounder – Infrared (AIRS-IR) precipitation biases over Antarctica at the pixel level, as these two products are the primary inputs for the Global Precipitation Climatology Project (GPCP) in high latitudes. We assessed the adjusted TOVS and AIRS-IR through analyses of geographical maps and time series of monthly and seasonal mean precipitation rates. The results are encouraging, indicating that the proposed approach could replace the current approach used in the GPCP for adjusting precipitation estimates from AIRS-IR and TOVS over Antarctica. Adjusting precipitation estimates from TOVS and AIRS-IR using the proposed approach improves the Kling-Gupta efficiency (KGE) over their entire period by 162% and 147%, respectively. Moreover, the proposed approach can be applied to adjust other precipitation products over Antarctica at the pixel level, including satellites and reanalysis products.

Climate classification for Northeast Brazil using reanalysis data and the Absolute Aridity Index

The climate of the Northeast Region of Brazil (NEB) has been intensively studied and analyzed for climate classification. The aridity index of the United Nations Environment Programme (UNEP) (AIUNEP) has been used for this purpose, but without fully satisfactory results. The input variables needed for its calculation are precipitation and reference potential evapotranspiration (ET0). However, although rainfall stations recording routine measurements of precipitation are well distributed in the NEB, they do not provide the necessary variables for estimating ET0. Thus, interpolation is used to calculate ET0, but this can generate errors. Another objective climate classification approach is the Thornthwaite method, based on the determination of the moisture index (Im), whose calculation also requires weather station data. Thus, seeking to circumvent the problem of paucity of stations and improve the spatial distribution of information on meteorological variables in the NEB, the present work had as one of its objectives to validate reanalysis data from ERA5 of the European Center for Medium-range Weather Forecast (ECMWF) and the unified gauge-based analysis of global daily precipitation project of the Climate Prediction Center/National Oceanic and Atmospheric Administration (CPC/NOAA). After validation, climate classifications were developed for the NEB using the AIUNEP and Im. It was observed that the Thornthwaite climate classification overestimated the aridity in the NEB, while the IAUNEP tended to underestimate it. For this reason, a new climate classification index, called the absolute aridity index (Iab), was suggested, yielding satisfactory results.

Evaluation of precipitation extremes in ERA5 reanalysis driven regional climate simulations over the CORDEX-Australasia domain

Reanalysis-driven regional climate simulations using the Weather Research and Forecasting (WRF) model in New South Wales (NSW) and Australian Regional Climate Modelling (NARCliM) Version 2.0 are assessed for capturing precipitation extreme indices. Seven configurations of the WRF model driven by ECMWF (European Centre for Medium-Range Weather Forecasts) Reanalysis v5 (ERA5) for Australia from 1979 to 2020 at 20 km resolution are evaluated. We assess the spatiotemporal patterns of six selected Expert Team on Sector-Specific Climate Indices (ET-SCI) precipitation extremes by comparing regional climate model (RCM) simulations against gridded observations. The RCMs evaluated have varying levels of accuracy in simulating precipitation extremes. While they capture climatology and coefficient of variation of precipitation extremes relatively well, temporal correlation and trend reproduction present challenges. Some RCMs perform more effectively for specific extreme indices, while others encounter challenges in accurately replicating them. No single RCM excels in all aspects, highlighting the need to consider specific strengths when selecting RCMs for global climate model (GCM) driven simulations.

Assimilation of GNSS tropospheric gradients into the Weather Research and Forecasting (WRF) model version 4.4.1

Abstract. In this study, we have incorporated tropospheric gradient observations from a Global Navigation Satellite Systems (GNSS) ground station network into the Weather Research and Forecasting (WRF) model through a newly developed observation operator. The experiments aim at testing the functionality of the developed observation operator and at analyzing the impact of tropospheric gradients on the sophisticated data assimilation (DA) system. The model was configured for a 0.1° mesh over Germany with 50 vertical levels up to 50 hPa. Our initial conditions were obtained from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) data at 0.25° resolution, and conventional observations were obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF), restricted to mainly surface stations and radiosondes. We selected approximately 100 GNSS stations with high data quality and availability covering Germany. We performed DA every 6 h for June and July 2021. Four experiments were conducted: (1) a control run assimilating only conventional observations; (2) an impact run assimilating zenith total delays (ZTDs) on top of the control run; (3) an impact gradient run assimilating ZTDs and gradients on top of the control run; and (4) a gradient run assimilating only gradients on top of the control run. The error for the impact run was reduced by 32 % and 10 % for ZTDs and gradients, whereas the error for the impact gradient run was reduced by 35 % and 18 %, respectively. The gradient errors for the gradient run were nearly equal to those of the impact gradient. Overall, the newly developed operator for the WRFDA system works as intended. In particular, the combined assimilation of gradients and the ZTDs led to a notable improvement in the humidity field at altitudes above 2.5 km. With the operator codes developed and freely available to the WRF users, we aim to trigger further GNSS tropospheric gradient assimilation studies.

The investigation of retrieval method for aerosol and water vapor profiles based on MAX-DOAS technology

Simultaneously obtaining the vertical distribution and profile of aerosol and water vapor is crucial for in-depth understanding of the Earth's radiation balance, regional water cycle, and the causes of frequent haze weather in China. This study developed a system and inversion method based on multi-axis differential optical absorption spectroscopy (MAX-DOAS) to simultaneously retrieve tropospheric aerosol and water vapor profiles. The absorption of dioxy (O4) was obtained by multi-elevation measurement using the developed ground-based MAX-DOAS system. By minimizing a cost function combining an atmospheric radiative transfer model, aerosol extinction profiles were first inverted, followed by water vapor profiles. Optimized parameter selection in the look-up table method reduced reliance on prior profile information, enhancing pollutant distribution retrieval accuracy. In the monitoring period in Huaibei region, aerosols were mainly below 1.5 km, while water vapor concentrations decreased with altitude. Comparisons of water vapor vertical column densities (VCDs) retrieved through the look-up table method with the European Centre for Medium-Range Weather Forecasts(ECMWF) ERA5 model and geometric approximation showed good agreement, with correlation coefficients of 0.93 and 0.98, respectively. To comprehend water vapor sources at various vertical layers, a 24-h backward trajectory clustering analysis was conducted using the HYSPLIT model based on observed wind fields. Findings revealed that at 500 m altitude, water vapor primarily emanated from the southeast, whereas at 1 km and 2 km altitudes, it predominantly originated from the southwest. The study advances simultaneous tropospheric aerosol and water vapor profile retrieval, providing reliable technical support for the investigation of regional pollution.

How well can global ensemble forecasts predict tropical cyclones in the southwest Indian Ocean?

AbstractThe southwest Indian Ocean (SWIO) recently experienced its most active, costliest and deadliest cyclone season on record (2018–2019). The anticipation and forecasting of natural hazards, such as tropical cyclones, are crucial to preparing for their impacts, but it is important to understand how well forecasting systems can predict them. Despite the vulnerability of the SWIO to tropical cyclones, comparatively little research has focused on this region, including understanding the ability of numerical weather prediction systems to predict cyclones and their impacts in southeast Africa. In this study, we evaluate ensemble probabilistic and high‐resolution deterministic forecasts of tropical cyclones in the SWIO from 2010 to 2020, using two state‐of‐the‐art global forecasting systems: one from the European Centre for Medium‐Range Weather Forecasts (ECMWF) and the other from the U.K. Met Office. We evaluate predictions of the track, assessing the location of the centre of each storm and its speed of movement, as well as its intensity, looking at maximum wind speeds and minimum central pressure, and discuss how the forecasts have evolved over the 10‐year period. Overall, ECMWF typically provides more accurate forecasts, but both systems tend to underestimate translation speed and intensity. We also investigate the impact of the Madden‐Julian Oscillation (MJO) on tropical cyclones and their forecasts. The MJO impacts where and when tropical cyclones form, their tracks and intensities, which in turn impacts forecast skill. These results are intended to provide an increased understanding of the ability of global forecasting systems to predict tropical cyclones in the SWIO, for the purpose of decision making and anticipatory action.

Deep Learning of a 200-member Ensemble with a Limited Historical Training to Improve the Prediction of Extreme Precipitation Events

Abstract This study introduces a deep learning (DL) scheme to generate reliable and skillful probabilistic quantitative precipitation forecasts (PQPFs) in a postprocessing framework. Enhanced machine learning model architecture and training mechanisms are proposed to improve reliability and skill of PQPFs while permitting computationally efficient model fitting using a short training dataset. The methodology is applied to postprocessing of 24h accumulated PQPFs from an ensemble forecast system recently introduced by the Center for Western Weather and Water Extremes (CW3E) and for lead times from one to six days. The ensemble system was designed based on a high-resolution version of the Weather Research and Forecasting (WRF) model, named West-WRF, to produce a 200-member ensemble in near real-time (NRT) over the Western U.S. during the boreal cool seasons to support forecastinformed reservoir operations (FIRO) and studies of prediction of heavy-to-extreme events. Postprocessed PQPFs are compared with those from the raw West-WRF ensemble, the operational Global Ensemble Forecast System (GEFSv12) as well as the ensemble from the European Centre for Medium-Range Weather Forecasts (ECMWF). As an additional baseline, we provide PQPFs verification metrics from a recently developed neural network postprocessing scheme proposed by Ghazvinian et al. (2022). Results demonstrate that skill of postprocessed forecasts significantly outperform PQPFs and deterministic forecasts from raw ensembles and the recently developed algorithm. The resulting PQPFs broadly improve upon the reliability and skill of baselines in predicting heavy to extreme precipitation (e.g., > 75 mm) across all lead times while maintaining the spatial structure of the high-resolution raw ensemble.

Climate Variability of Air Temperature and Its Warming Trends in the Arabian Gulf

AbstractOver a 41-year period (1979–2019), the climate of the Arabian Gulf (Gulf) exhibited notable variability. This variation, influenced by mesoscale and globalscale meteorological and climatic determinants, showed temperature fluctuations, predominantly in the Northern Gulf (N Gulf). The N Gulf recorded an average temperature of 25.9 °C, a noticeable 3 °C higher than the broader Gulf average. This period marked extreme temperature records, with a peak at 53.8 °C in July 2016, paralleling the highest known temperatures in Asia. Data reveals a warming trend (0.1–0.6 °C per decade) in the region, with the N Gulf’s rise surpassing global trends. The years 2014–2019 were notably warmer, impacting overall trends, and 1997 was the transitional year shifting from cooler to warmer climates. The N Gulf’s seasonal dynamics are also transforming, with extended summers and shortened winters. In this study, we employed empirical orthogonal function (EOF) decomposition analysis to examine the interannual to multidecadal fluctuations in air temperature within the Gulf from 1979 to 2019, using hourly air temperature data from the Watch Forcing Dataset (WFDE5) generated from the fifth generation (ERA5) European Centre for Medium-Range Weather Fore-casts (ECMWF). Through this method, we discerned three principal EOF modes, which together accounted for 82% of the air temperature variability in the Gulf, highlighting external influences such as North Atlantic Oscillation (NAO), El Nino-Southern Oscillation (ENSO), and Atlantic Multidecadal Oscillation (AMO). The data for these external influences were sourced from the National Oceanic and Atmospheric Administration (NOAA) Physical Sciences Laboratory. These findings stress the pressing need for detailed climate studies and mitigation strategies for a rapidly warming Gulf region, especially the N Gulf.

Use of CAMS near Real-Time Aerosols in the HARMONIE-AROME NWP Model

Near real-time aerosol fields from the Copernicus Atmospheric Monitoring Services (CAMS), operated by the European Centre for Medium-Range Weather Forecasts (ECMWF), are configured for use in the HARMONIE-AROME Numerical Weather Prediction model. Aerosol mass mixing ratios from CAMS are introduced in the model through the first guess and lateral boundary conditions and are advected by the model dynamics. The cloud droplet number concentration is obtained from the aerosol fields and used by the microphysics and radiation schemes in the model. The results show an improvement in radiation, especially during desert dust events (differences of nearly 100 W/m2 are obtained). There is also a change in precipitation patterns, with an increase in precipitation, mainly during heavy precipitation events. A reduction in spurious fog is also found. In addition, the use of the CAMS near real-time aerosols results in an improvement in global shortwave radiation forecasts when the clouds are thick due to an improved estimation of the cloud droplet number concentration.

Opposite effects of aerosols and meteorological parameters on warm clouds in two contrasting regions over eastern China

Abstract. The sensitivity (S) of cloud parameters to the influence of different aerosol and meteorological parameters has in most previous aerosol–cloud interaction (aci) studies been addressed using traditional statistical methods. In the current study, relationships between cloud droplet effective radius (CER) and aerosol optical depth (AOD, used as a proxy for cloud condensation nuclei, CCN), i.e., the sensitivity (S) of CER to AOD, are investigated with different constraints of AOD and cloud liquid water path (LWP). In addition to traditional statistical methods, the geographical detector method (GDM) is applied in this study to quantify the relative importance of the effects of aerosol and meteorological parameters, as well as their interaction, on S. Moderate Resolution Imaging Spectroradiometer (MODIS) C6 L3 data and European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-5 reanalysis data, for the period from 2008 to 2022, were used to investigate aci over eastern China. Two contrasting areas were selected: the heavily polluted Yangtze River Delta (YRD) and a relatively clean area over the East China Sea (ECS). Linear regression analysis shows that CER decreases with the increase in AOD (negative S) in the moderately polluted atmosphere (0.1<AOD<0.3) over the ECS, whereas, in contrast, CER increases with increasing AOD (positive S) in the polluted atmosphere (AOD>0.3) over the YRD. Evaluation of S as function of the LWP shows that in the moderately polluted atmosphere over the ECS, S is negative in the LWP interval [40 g m−2, 200 g m−2], and the sensitivity of CER to AOD is substantially stronger as LWP is larger. In contrast, in the polluted atmosphere over the YRD, S is positive in the LWP interval [0 g m−2, 120 g m−2] and does not change notably as function of LWP in this interval. The study further shows that over the ECS, the CER is larger for higher low tropospheric stability (LTS) and relative humidity (RH) but lower for higher pressure vertical velocity (PVV). Over the YRD, there is no significant influence of LTS on the relationship between CER and AOD. The GDM has been used as an independent method to analyze the sensitivity of cloud parameters to AOD and meteorological parameters (RH, LTS and PVV). The GDM has also been used to analyze the effects of interactions between two parameters and thus obtain information on confounding meteorological effects on the aci. Over the ECS, cloud parameters are sensitive to almost all parameters considered except for cloud top pressure (CTP), and the sensitivity to AOD is larger than that to any of the meteorological factors. Among the meteorological factors, the cloud parameters are most sensitive to PVV and least sensitive to RH. Over the YRD, the explanatory power of the sensitivity of cloud parameters to AOD and meteorological parameters is much smaller than over the ECS, except for RH, which has a statistically significant influence on CTP and can explain 74 % of the variation of CTP. The results from the GDM analysis show that cloud parameters are more sensitive to the combination of aerosol and a meteorological parameter than to each parameter alone, but confounding effects due to co-variation of both parameters cannot be excluded.

A review of the BuFeng-1 GNSS-R mission: calibration and validation results of sea surface and land surface

ABSTRACT In this paper, we will conclude the results of Bufeng-1 (BF-1) A/B data processing, calibration workflow, and validation of the calibrated sea surface winds, land surface soil moisture, and sea surface height measurements. Since 2019, the BF-1 mission has operated in-orbit for over 4 years. The Earth reflected delay Doppler maps (DDMs) are continuously collected to perform global sea surface and land observations. At the same time, the intermediate frequency (IF) raw data are also obtained for 12 seconds every pass in diagnostic mode. To begin with, a brief description of the spaceborne Global Navigation Satellite System Reflectometry (GNSS-R) technique will be provided in the introduction. Next, we will present the overview of Chinese BF-1 mission and the data specifications used in our research. In the next section, the BF-1 mission-related spaceborne power calibration and validation are presented to show the support to power DDM observable production for sea surface and land surface applications. Then, the status of Chinese Beidou System (BDS) Equivalent Isotropic Radiated Power (EIRP) acquisition programme is then introduced. Furthermore, the latest sea surface height (SSH) measurements results including two modes (group delay and carrier phase) and wind speed derivation based on machine learning (ML) method will be spatial-temporal aligned and validated with auxiliary datasets including Denmark Technology University (DTU) mean sea surface (MSS) products and European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis. The previous published results of sea surface winds retrieval under Hurricane conditions and soil moisture retrieval are also reviewed for the BF-1 mission applications. Finally, the conclusion of BF-1 derived results will be discussed to draw out ongoing/future works.

Spatial-temporal characteristics analysis of solar irradiance forecast errors in Europe and North America

Accurate photovoltaic (PV) power forecast is crucial to power systems. Solar irradiance forecast is the fundamentals of PV power forecast. Current researches merely focus on improving the forecast accuracy. There hasn't been comprehensive study on the spatial-temporal characteristics of solar irradiance forecast errors. Thus, this paper analyzed the error characteristics of solar irradiance forecast provided by European Centre for Medium-Range Weather Forecasts (ECMWF), one of the most accurate numerical weather prediction (NWP) products. ECMWF's forecasts were compared with ERA5 solar irradiance reanalysis data to reveal the error distribution characteristics. Four-year (2017–2020) data from the geographical region bounded by 63°N, −126°W, 21°S, and 36°E (covering most parts of Europe and North America) were studied. Experiments show that solar irradiance forecast errors peak at noon of the local time. Furthermore, correlations between solar irradiance forecast errors and other weather variables were also revealed. Experiments show that solar irradiance forecast errors have negative correlation with low cloud cover forecast errors and positive correlation with air temperature forecast errors. The possible reasons for the correlation relationship were also analyzed in detail. This paper systematically reveals the spatial-temporal characteristics of solar irradiance forecast errors and provides a useful guideline for solar PV system operation.

Summer snowline altitude gradients in Western Norway are influenced by maritime climate

ABSTRACT Mountain snowline dynamics are relatively underreported with few studies exploring spatial snowline dynamics. Whilst clear regional-scale relationships between snowline location and temperature exist in European mountains, recent research at higher latitudes reports no response to climate change. In maritime mountains, snowlines occupy complex environmental gradients. Using timeseries of satellite data from Landsat missions 5–8 (151 images between 1984 and 2021), we explored sub-regional summer snowline dynamics across the maritime-continental climate gradient in the Western Norwegian mountains. We characterize spatio-temporal snowline altitude dynamics and investigate the climate factors altering snowline patterns. Summer snowline altitudes were found to increase inland at around double the rate of the 0°C summer isotherm. Data from the European Centre for Medium-Range Weather Forecasts (ECMWF) land component of the fifth generation of European Reanalysis (ERA5-Land), showed a potential ‘maritime-mountain’ effect with coastal orographic snowfall and cloud cover-induced surface solar downwelling radiation amplifying maritime-continental snowline altitude gradients alongside surface atmospheric temperature. This was replicated in the Canadian Rocky Mountains. Between 1984 and 2021, we found spatial summer snowline gradients in Norway decreased and propose multiple climate forcings are responsible, potentially masking links between snowlines and climate change. Although non-significant, the data also suggest regional summer snowline altitudes increased. This study demonstrates the complex spatial heterogeneity in snow-climate relationships and highlights how long-term snow dynamics can be queried using fine-grain (Landsat) resolution satellite data. We share our approach through a Google Earth Engine web-app that rapidly executes spatial snowline analyses for global mountain regions via a graphical user interface.

Improving the statistical downscaling performance of climatic parameters with convolutional neural networks

Abstract This study examines two downscaling techniques, convolutional neural networks (CNNs) and feedforward neural networks for predicting precipitation and temperature, alongside statistical downscaling as a benchmark model. The daily climate predictors were extracted from the European Center for Medium-range Weather Forecast (ECMWF) ERA5 dataset spanning from 1979 to 2010 for Tabriz city, located in the northwest of Iran. The biases in precipitation data of ERA5 predictors were corrected through the empirical quantile mapping method. Also, two nonlinear predictor screening methods, random forest and mutual information, were employed, alongside linear correlation coefficient. While these methods facilitate identification of dominant regional climate change drivers, it is essential to consider their limitations, such as sensitivity to parameter settings, assumptions about data relationships, potential biases in handling redundancy and correlation, challenges in generalizability across datasets, and computational complexity. Evaluation results indicated that CNN, when applied without predictor screening, achieves coefficient of determination of 0.98 for temperature and 0.71 for precipitation. Ultimately, future projections were employed under two shared socioeconomic pathways (SSPs), SSP2-4.5 and SSP5-8.5, and concluded that the most increase in temperature by 2.9 °C and decrease in precipitation by 3.5 mm may occur under SSP5-8.5.