Validation of satellite-based sea surface temperature products against in situ observations off the western coast of Sumatra

This study validated the sea surface temperature (SST) datasets from the Group for High-Resolution SST Multi Product Ensemble (GMPE), National Oceanic and Atmospheric Administration (NOAA) Optimal Interpolation (OI) SST version 2 and 2.1 (OIv2 and OIv2.1), and Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) in the area off the western coast of Sumatra against in situ observations. Furthermore, the root mean square differences (RMSDs) of OIv2, OIv2.1, and ECCO2 were investigated with respect to GMPE, whose small RMSD < 0.2 K against in situ observations confirmed its suitability as a reference. Although OIv2 showed a large RMSD (1–1.5 K) with a significant negative bias, OIv2.1 (RMSD < 0.4 K) improved remarkably. In the average SST distributions for December 2017, the differences among the 4 datasets were significant in the areas off the western coast of Sumatra, along the southern coast of Java, and in the Indonesian inland sea. These results were consistent with the ensemble spread distribution obtained with GMPE. The large RMSDs of OIv2 corresponded to high clouds, and it was suggested that the change in the satellites used for SST estimation contributed to the improvement in OIv2.1.

A new pan-European dataset for gridded daily average wind speed based on in-situ observations

&lt;p&gt;&lt;strong&gt;Objective&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;The purpose of the European Copernicus Climate Change Service (C3S) is to support society by providing information about the past, present and future climate. For the service related to &lt;em&gt;in-situ&lt;/em&gt; observations, one of the objectives is to provide high-resolution (0.1x0.1 and 0.25x0.25 degrees) gridded wind speed fields. The gridded wind fields are based on ECA&amp;D daily average station observations for the period 1970-2020.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Research question&lt;/strong&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;We address the following research questions: [1] How efficiently can we provide the gridded wind fields as a statistically reliable ensemble, in order to represent the uncertainty of the gridding? [2] How efficiently can we exploit high-resolution geographical auxiliary variables (e.g. digital elevation model, terrain roughness) to augment the station data from a sparse network, in order to provide gridded wind fields with high-resolution local features?&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Approach&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;In our analysis, we apply greedy forward selection linear regression (FSLR) to include the high-resolution effects of the auxiliary variables on monthly-mean data. These data provide a &amp;#8216;background&amp;#8217; for the daily estimates. We apply cross-validation to avoid FSLR over-fitting and use full-cycle bootstrapping to create FSLR ensemble members. Then, we apply Gaussian process regression (GPR) to regress the daily anomalies. We consider the effect of the spatial distribution of station locations on the GPR gridding uncertainty.&lt;/p&gt;&lt;p&gt;The goal of this work is to produce several decades of daily gridded wind fields, hence, computational efficiency is of utmost importance. We alleviate the computational cost of the FSLR and GPR analyses by incorporating greedy algorithms and sparse matrix algebra in the analyses.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Novelty&lt;/strong&gt;&amp;#160;&amp;#160;&amp;#160;&lt;/p&gt;&lt;p&gt;The gridded wind fields are calculated as a statistical ensemble of realizations. In the present analysis, the ensemble spread is based on uncertainties arising from the auxiliary variables as well as from the spatial distribution of stations.&lt;/p&gt;&lt;p&gt;Cross-validation is used to tune the GPR hyper parameters. Where conventional GPR hyperparameter tuning aims at an optimal prediction of the gridded mean, instead, we tune the GPR hyperparameters for optimal prediction of the gridded ensemble spread.&lt;/p&gt;&lt;p&gt;Building on our experience with providing similar gridded climate data sets, this set of gridded wind fields is a novel addition to the E-OBS climate data sets.&lt;/p&gt;

Coupled parameter perturbations in the SPP scheme of the HarmonEPS

&lt;p&gt;In an ensemble prediction system (EPS) the uncertainty in the initial atmospheric conditions is usually represented via perturbation of the initial atmospheric state and different boundary conditions at the beginning and throughout the duration of the forecast. These approaches exclude the uncertainty due to the representation of physical processes within the parameterization schemes of a numerical weather prediction model (NWP). Much of the uncertainty in the presentation of physical process arises from uncertain parameter values regulating key physical processes in the boundary-layer and microphysics schemes. This uncertainty can be represented with a Stochastically Perturbed Parameterization (SPP) scheme, where parameter values for the different model grid points are randomly selected from a defined probability density function. The SPP scheme can improve model performance and increase ensemble spread, but may lead to unrealistic parameter values, which can introduce additional model bias. A potential solution is to use coupled/correlated perturbations for relevant SPP parameters to increase the model performance and ensemble spread, while maintaining physically realistic ranges for the parameters. In this study, we investigate the impact of coupled perturbations in key parameters within the boundary-layer and microphysics schemes of the HarmonEPS model using the new SPP scheme. The performance of the coupled perturbations experiment is evaluated against HarmonEPS experiments using independent parameter perturbations, and perturbations in the initial atmospheric state and boundary conditions for both a winter and a summer period. &amp;#160;We find that coupled perturbations in the SPP scheme can decrease model bias and increase the ensemble spread for the 2m temperature and relative humidity, 10m-wind speed and total cloud cover.&lt;/p&gt;

Accounting for land model error in numerical weather prediction ensemble systems: toward ensemble-based coupled land/atmosphere data assimilation

The ensembles used in the NOAA National Centers for Environmental Prediction (NCEP) global data assimilation and numerical weather prediction (NWP) system are under-dispersed at and near the land surface, preventing their use in ensemble-based land data assimilation. Comparison to offline (land-only) data assimilation ensemble systems suggests that while the relevant atmospheric fields are under-dispersed in NCEP’s system, this alone cannot explain the under-dispersed land component, and an additional scheme is required to explicitly account for land model error. This study then investigates several schemes for perturbing the soil (moisture and temperature) states in NCEP’s system, qualitatively comparing the induced ensemble spread to independent estimates of the forecast error standard deviation in soil moisture, soil temperature, 2m temperature, and 2m humidity. Directly adding perturbations to the soil states, as is commonly done in offline systems, generated unrealistic spatial patterns in the soil moisture ensemble spread. Application of a Stochastically Perturbed Physics Tendencies scheme to the soil states is inherently limited in the amount of soil moisture spread that it can induce. Perturbing the land model parameters, in this case vegetation fraction, generated a realistic distribution in the ensemble spread, while also inducing perturbations in the land (soil states) and atmosphere (2m states) that are consistent with errors in the land/atmosphere fluxes. The parameter perturbation method is then recommended for NCEP’s ensemble system, and it is currently being refined within the development of an ensemble-based coupled land/atmosphere data assimilation for NCEP’s NWP system.

Validation of Satellite-Based Sea Surface Temperature Products Against in Situ Observations Off the Western Coast of Sumatra

This study validated the sea surface temperature (SST) datasets from the Group for High-Resolution SST Multi Product Ensemble (GMPE), National Oceanic and Atmospheric Administration (NOAA) Optimal Interpolation (OI) SST version 2 and 2.1 (OIv2 and OIv2.1), and Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) in the area off the western coast of Sumatra against in situ observations. Furthermore, the root mean square differences (RMSDs) of OIv2, OIv2.1, and ECCO2 were investigated with respect to GMPE, whose small RMSD < 0.2 K against in situ observations confirmed its suitability as a reference. Although OIv2 showed a large RMSD (1-1.5 K) with a significant negative bias, OIv2.1 (RMSD < 0.4 K) improved remarkably. In the average SST distributions for December 2017, the differences among the 4 datasets were significant in the areas off the western coast of Sumatra, along the southern coast of Java, and in the Indonesian inland sea. These results were consistent with the ensemble spread distribution obtained with GMPE. The large RMSDs of OIv2 corresponded to high clouds, and it was suggested that the change in the satellites used for SST estimation contributed to the improvement in OIv2.1.

Downscaling of a seasonal ensemble forecast at the convection-permitting resolution over the Horn of Africa using the WRF model

&lt;p&gt;Studies have shown the benefits of convection-permitting downscaling at the seasonal scale using limited-area models. To evaluate the performance with real forecasts as boundary conditions, four members of the SEAS5 global ensemble were dynamically downscaled over Ethiopia during June, July, and August 2018 at a 3-km resolution. We used a multi&amp;#8208;physics ensemble based on the WRF model to compare the effects of boundary conditions and physics &lt;span&gt;&lt;span&gt;parametrization&lt;/span&gt;&lt;/span&gt; producing 16 ensemble members. With ECMWF analyses as a reference, SEAS5 averaged to a +0.17&amp;#176;C bias over Ethiopia whereas WRF resulted in +1.14&amp;#176;C. With respect to precipitation, the WRF model simulated 264 mm compared to 248 mm for SEAS5 and 236 mm for GPM-IMERG. The maximum northward extension of the tropical rain belt decreased by about 2&amp;#176; in both models. Downscaling enhanced the ensemble spread in precipitation by 60% on average, correcting the SEAS5 underdispersion. The WRF ensemble spread over Ethiopia was mostly generated by the perturbed boundary conditions, as their effect is often 50% larger than the physics&amp;#8208;induced variability. The results indicate that boundary condition perturbations are necessary, although not always sufficient, to generate the right amount of ensemble spread in a limited-area model with complex topography. The next step is to use specific methods to calculate the added value provided by the downscaling.&lt;/p&gt;

Development of stochastically perturbed parameterization scheme for the Noah Land Surface Model with the optimized random forcing parameters using the micro-genetic algorithm

&lt;p&gt;The ensemble data assimilation system is beneficial to represent the initial uncertainties and flow-dependent background error covariance (BEC). In particular, the inevitable model uncertainties can be expressed by ensemble spread, that is the standard deviation of ensemble BEC. However, the ensemble spread generally suffers from under-estimated problems. To alleviate this problem, recent studies employed stochastic perturbation schemes to increases the ensemble spreads by adding the random forcing in the model tendencies (i.e., physical or dynamical tendencies) or parameterization schemes (i.e., PBL, convective scheme, etc.). In this study, we focus on the near-surface uncertainties which are affected by the interactions between the land and atmosphere process. The land surface model (LSM) provides various fluxes as the lower boundary condition to the atmosphere, influencing the accuracy of hourly-to-seasonal scale weather forecasting, but the surface uncertainties were not much addressed yet. In this study, we developed the stochastically perturbed parameterization (SPP) scheme for the Noah LSM. The Weather Research and Forecasting (WRF) ensemble system is used for regional weather forecasting over East Asia, especially over the Korean Peninsula. As a testbed experiment with the newly-developed Noah LSM-SPP system, we first perturbed the soil temperature &amp;#8212; a crucial variable for the near-surface forecasts by affecting sensible heat fluxes, land surface skin temperature and surface air temperature, and hence lower-tropospheric temperature. Here, the random forcing used in perturbation is made by the tuning parameters for amplitude, length scale, and time scales: they are commonly determined empirically by trial and error. In order to find optimal tuning parameter values, we applied a global optimization algorithm &amp;#8212; the micro-genetic algorithm (micro-GA) &amp;#8212; to achieve the smallest root-mean-squared errors. Our results indicate that optimization of the random forcing parameters contributes to an increase in the ensemble spread and a decrease in the ensemble mean errors in the near-surface and lower-troposphere uncertainties. Further experiments will be conducted by including soil moisture in the testbed.&lt;/p&gt;

Dynamical Structures of Cross-Domain Forecast Error Covariance of a Simulated Tropical Cyclone in a Convection-Permitting Coupled Atmosphere–Ocean Model

Understanding the dynamics of the flow-dependent forecast error covariance across the air–sea interface is beneficial toward revealing the potential influences of strongly coupled data assimilation on tropical cyclone (TC) initialization in coupled models, and the fundamental dynamics associated with TC air–sea interactions. A 200-member ensemble of convection-permitting forecasts from a coupled atmosphere–ocean regional model is used to investigate the forecast error covariance across the oceanic and atmospheric domains during the rapid intensification of Hurricane Florence (2018). Forecast uncertainties in both atmospheric and oceanic domains, from an Eulerian perspective, increase with forecast lead time, mainly from TC displacement errors. In a storm-relative framework, the ensemble forecast uncertainties in both domains are predominantly caused by differences in the simulated storm intensity and structure. The largest ensemble spread in the atmospheric pressure, temperature, and wind fields can be found within the TC inner-core region. Alternatively, the largest ensemble spread in the upper-ocean currents and temperature fields are located along the cold wake behind the storm. Cross-domain ensemble correlations between simulated atmospheric (oceanic) observations and oceanic (atmospheric) state variables in the storm-relative coordinates are highly anisotropic, variable dependent, and ultimately driven by the dynamics of TC air–sea interactions. Meaningful and dynamically consistent cross-domain ensemble correlations suggest that it is possible to use atmospheric and oceanic observations to simultaneously update state variables associated with the coupled ocean–atmosphere prediction of TCs using strongly coupled data assimilation. Sensitivity experiments demonstrate that at least 60–80 ensemble members are required to represent physically consistent cross-domain correlations and minimize sampling errors.

Cross-basin Interactions between the Tropical Atlantic and Pacific in the ECMWF Hindcasts

Using the ensemble hindcasts of the European Centre for Medium-Range Weather Forecasts (ECMWF) coupled model for the period of 1980-2005, spatio-temporal evolution in the covariability of sea surface temperature (SST) and low-level winds in the ensemble mean and spread over the tropical Atlantic is investigated with the month-reliant singular value decomposition (SVD) method, which treats the variables in a given monthly sequence as one time step. The leading mode of the ensemble mean represents a co-evolution of SST and winds over the tropical Atlantic associated with a phase transition of El Niño from the peak to decay phase, while the second mode is related to a phase transition from El Niño to La Niña, indicating a precursory role of the north tropical Atlantic (NTA) SST warming in La Niña development.The leading mode of ensemble spread in SST and winds further illustrates that an NTA SST anomaly acts as a precursor for El Niño-Southern Oscillation (ENSO). A north-tropical pathway for the delayed effect of the NTA SST anomaly on the subsequent ENSO event is identified; the NTA SST warming induces the subtropical Northeast Pacific SST cooling through the modulation of a zonal-vertical circulation, setting off a North Pacific Meridional Mode (NPMM). The coupled SST-wind anomalies migrate southwestward to the central equatorial Pacific and eventually amplify into a La Niña event in the following months due to the equatorial Bjerknes feedback. Ensemble spread greatly increases the sample size and affords insights into the inter-basin interactions between the tropical Atlantic and Pacific, as demonstrated here in the NTA SST impact on ENSO.

What Does a Convection-Allowing Ensemble of Opportunity Buy Us in Forecasting Thunderstorms?

The High Resolution Ensemble Forecast v2.1 (HREFv2.1), an operational convection-allowing model (CAM) ensemble, is an “ensemble of opportunity” wherein forecasts from several independently designed deterministic CAMs are aggregated and postprocessed together. Multiple dimensions of diversity in the HREFv2.1 ensemble membership contribute to ensemble spread, including model core, physics parameterization schemes, initial conditions (ICs), and time lagging. In this study, HREFv2.1 forecasts are compared against the High Resolution Rapid Refresh Ensemble (HRRRE) and the Multiscale data Assimilation and Predictability (MAP) ensemble, two experimental CAM ensembles that ran during the 5-week Spring Forecasting Experiment (SFE) in spring 2018. The HRRRE and MAP are formally designed ensembles with spread achieved primarily through perturbed ICs. Verification in this study focuses on composite radar reflectivity and updraft helicity to assess ensemble performance in forecasting convective storms. The HREFv2.1 shows the highest overall skill for these forecasts, matching subjective real-time impressions from SFE participants. Analysis of the skill and variance of ensemble member forecasts suggests that the HREFv2.1 exhibits greater spread and more effectively samples model uncertainty than the HRRRE or MAP. These results imply that to optimize skill in forecasting convective storms at 1–2-day lead times, future CAM ensembles should employ either diverse membership designs or sophisticated perturbation schemes capable of representing model uncertainty with comparable efficacy.