Search NASA⌕ Search

SEARCH · Search NASA

Results for “ensemble data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

The Seasonal Cycle of Satellite Chlorophyll Fluorescence Observations and its Relationship to Vegetation Phenology and Ecosystem Atmosphere Carbon Exchange

Mapping of terrestrial chlorophyll uorescence from space has shown potentialfor providing global measurements related to gross primary productivity(GPP). In particular, space-based fluorescence may provide information onthe length of the carbon uptake period that can be of use for global carboncycle modeling. Here, we examine the seasonal cycle of photosynthesis asestimated from satellite fluorescence retrievals at wavelengths surroundingthe 740nm emission feature. These retrievals are from the Global OzoneMonitoring Experiment 2 (GOME-2) flying on the MetOp A satellite. Wecompare the fluorescence seasonal cycle with that of GPP as estimated froma diverse set of North American tower gas exchange measurements. Because the GOME-2 has a large ground footprint (40 x 80km2) as compared with that of the flux towers and requires averaging to reduce random errors, we additionally compare with seasonal cycles of upscaled GPP in the satellite averaging area surrounding the tower locations estimated from the Max Planck Institute for Biogeochemistry (MPI-BGC) machine learning algorithm. We also examine the seasonality of absorbed photosynthetically-active radiation(APAR) derived with reflectances from the MODerate-resolution Imaging Spectroradiometer (MODIS). Finally, we examine seasonal cycles of GPP as produced from an ensemble of vegetation models. Several of the data-driven models rely on satellite reflectance-based vegetation parameters to derive estimates of APAR that are used to compute GPP. For forested sites(particularly deciduous broadleaf and mixed forests), the GOME-2 fluorescence captures the spring onset and autumn shutoff of photosynthesis as delineated by the tower-based GPP estimates. In contrast, the reflectance-based indicators and many of the models tend to overestimate the length of the photosynthetically-active period for these and other biomes as has been noted previously in the literature. Satellite fluorescence measurements therefore show potential for improving model GPP estimates.

Joiner, J.↗

Probabilistic Modeling of Aircraft Trajectories for Dynamic Separation Volumes

With a proliferation of new and unconventional vehicles and operations expected in the future, the ab initio airspace design will require new approaches to trajectory prediction for separation assurance and other air traffic management functions. This paper presents an approach to probabilistic modeling of the trajectory of an aircraft when its intent is unknown. The approach uses a set of feature functions to constrain a maximum entropy probability distribution based on a set of observed aircraft trajectories. This model can be used to sample new aircraft trajectories to form an ensemble reflecting the variability in an aircraft's intent. The model learning process ensures that the variability in this ensemble reflects the behavior observed in the original data set. Computational examples are presented.

Lewis, Timothy A.↗

Modelling the Greenland Ice Sheet's Committed Contribution to Sea Level During the 21st Century

Mass loss from the Greenland Ice Sheet can be partitioned between surface mass balance and discharge due to ice dynamics through its marine-terminating outlet glaciers. A perturbation to a glacier terminus (e.g., a calving event) results in both an instantaneous response in velocity and mass loss and a diffusive response due to the evolution of ice thickness over time. This diffusive response means the total impact of a retreat event can take decades to be fully realized. Here we model the committed response of the Greenland Ice Sheet by applying perturbations to the marine-terminating glacier termini that represent recent observed changes, and simulating the response over the 21st century, while holding the climate forcing constant. The sensitivity of the ice sheet response to model parameter uncertainty is explored within an ensemble framework, and Gravity Recovery and Climate Experiment data is used to constrain the results using a Bayesian calibration approach. We find that the Greenland Ice Sheet's committed contribution to 21st century sea level rise is at least 33.5 [17.5 52.4] mm (25th and 75th percentiles in brackets), with at least 6 mm being attributable directly to terminus retreat that occurred between 2007 and 2015. The spread in our projections is driven by uncertainty in the basal friction coefficient. Our results complement the ISMIP6 Greenland projections, which report the ice sheet response to future forcing, excluding the background response. In this way, we can obtain estimates of Greenland's total contribution to sea level rise in 2100.

Isabel J. Nias↗

Development and Evaluation of Ensemble Consensus Precipitation Estimates over High Mountain Asia

Precipitation estimates are highly uncertain in complex regions such as High-Mountain Asia (HMA), where ground measurements are very difficult to obtain, and atmospheric dynamics poorly understood. Though gridded products derived from satellite-based observations and/or reanalysis can provide temporally and spatially distributed estimates of precipitation, there are significant inconsistencies in these products. As such, to date, there is little agreement in the community on the best and most accurate gridded precipitation product in HMA, which is likely area dependent because of HMA’s strong heterogeneities and complex orography. Targeting these gaps, this article presents the development of a consensus ensemble precipitation product using three gridded precipitation datasets (the Integrated Multi-satellitE Retrieals for Global Precipitation Measurement IMERG, the Climate Hazards group Infrared Precipitation with Stations CHIRPS, and the ECMWF Reanalysis ERA5) with a localized probability matched mean (LPM) approach. We evaluate the performance of the LPM estimate along with a simple ensemble mean (EM) estimate to overcome the differences and disparities of the three selected constituent products on long-term averages and trends in HMA. Our analysis demonstrates that LPM reduces the high biases embedded in the ensemble members and provides more realistic spatial patterns compared to EM. LPM is also a good alternative for merging data products with different spatio-temporal resolutions. By filtering disparities among the individual ensemble members, LPM overcomes the problem of a certain product performing well only in a particular area and provides a consensus estimate with plausible temporal trends.

Fadji Z Maina↗

A Scalable Real-Time Data Assimilation Framework for Predicting Turbulent Atmosphere Dynamics

AI-based foundation models like FourCastNet, GraphCast are revolutionizing weather and climate predictions but are not yet ready for operational use. Their limitation lies in the absence of a data assimilation system to incorporate real-time Earth system observations, crucial for accurately forecasting events like tropical cyclones. To overcome these obstacles, we introduce a generic real-time data assimilation framework and demonstrate its end-to-end performance on the Frontier supercomputer. This framework comprises two primary modules: an ensemble score filter (EnSF), which significantly outperforms the state-of-the-art data assimilation method, and a vision transformer-based surrogate capable of real-time adaptation through the integration of observational data. We demonstrate both the strong and weak scaling of our framework up to 1024 GPUs on the Exascale supercomputer, Frontier. Our results not only illustrate the framework's exceptional scalability on high-performance computing systems, but also demonstrate the importance of supercomputers in real-time data assimilation for weather and climate predictions.

Lu, Dan↗

Interannual Variability of OLR as Observed by AIRS and CERES

This paper compares spatial anomaly time series of OLR (Outgoing Longwave Radiation) and OLR(sub CLR) (Clear Sky OLR) as determined using observations from CERES Terra and AIRS over the time period September 2002 through June 2011. Both AIRS and CERES show a significant decrease in global mean and tropical mean OLR over this time period. We find excellent agreement of the anomaly time-series of the two OLR data sets in almost every detail, down to 1 deg X 1 deg spatial grid point level. The extremely close agreement of OLR anomaly time series derived from observations by two different instruments implies that both sets of results must be highly stable. This agreement also validates to some extent the anomaly time series of the AIRS derived products used in the computation of the AIRS OLR product. The paper also examines the correlations of anomaly time series of AIRS and CERES OLR, on different spatial scales, as well as those of other AIRS derived products, with that of the NOAA Sea Surface Temperature (SST) product averaged over the NOAA Nino-4 spatial region. We refer to these SST anomalies as the El Nino Index. Large spatially coherent positive and negative correlations of OLR anomaly time series with that of the El Nino Index are found in different spatial regions. Anomalies of global mean, and especially tropical mean, OLR are highly positively correlated with the El Nino Index. These correlations explain that the recent global and tropical mean decreases in OLR over the period September 2002 through June 2011, as observed by both AIRS and CERES, are primarily the result of a transition from an El Nino condition at the beginning of the data record to La Nina conditions toward the end of the data period. We show that the close correlation of global mean, and especially tropical mean, OLR anomalies with the El Nino Index can be well accounted for by temporal changes of OLR within two spatial regions which lie outside the NOAA Nino-4 region, in which anomalies of cloud cover and mid-tropospheric water vapor are both highly negatively correlated with the El Nino Index. Agreement of the AIRS and CERES OLR(sub CLR) anomaly time series is less good, which may be a result of the large sampling differences in the ensemble of cases included in each OLR(sub CLR) data set.

Susskind, Joel↗

An Inverse Chance-constrained Approach to the Calibration of Robust Models

This paper proposes a strategy to calibrate computational models according to uncertain input-output data. To this end, uncertainty in the data is first quantified by creating adversarial data sets. Samples drawn from such sets are then mapped from the input-output space to the parameter space using an inverse mapping. This mapping minimizes the collective output spread of an ensemble of point predictions while satisfying a set of individual data-matching requirements. The distribution of the resulting parameter points, which often exhibits strong parameter dependencies, is then modeled using sliced-normals. The chance-constrained formulation used to learn this distribution enables the analyst to trade-off a greater likelihood for most of the data against a lower likelihood for some of the data thereby relaxing the conservatism of the calibrated model. This formulation not only neglects the worst-performing quantiles of each adversarial distribution but also eliminates the potentially serious effects that outliers might have on the resulting model. This calibration approach not only has a considerably lower computational cost than the standard forward approach but it also allows for the identification of suitable distribution classes, which in turn yield better calibrated models.

Calibration↗

File Specification for M2AMIP Products: Version Number - 1.0

The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) is an atmospheric reanalysis computed with the Goddard Earth Observing (EOS) System, Version 5.12.4 (GEOS) data assimilation system (Gelaro et al., 2017). To supplement the reanalysis, the GEOS General Circulation Model (GCM) used in MERRA-2 has been used to generate a 10-member ensemble of simulations, configured following the convention of the Atmospheric Model Intercomparison Project (AMIP; Gates et al., 1992). Each ensemble member was initialized using meteorological fields from a different date in November 1979. The AMIP simulations used the sea-surface temperature (SST) and sea-ice boundary conditions that were used in MERRA-2 (Bosilovich et al., 2016). This 10-member ensemble of AMIP simulations, denoted M2AMIP, is available for download in a group of self-describing files, which are documented in this office note. All data collections are provided on the same horizontal grid as MERRA-2. This grid has 576 points in the longitudinal direction and 361 points in the latitudinal direction, corresponding to a resolution of 0.625 degrees by 0.5 degrees. Although data collections are available at this grid, all fields are computed on a cubed-sphere grid with an approximate resolution of 50 km by 50 km and are then spatially interpolated to the latitude-longitude grid. There are no changes in the vertical grids used: variables are provided on either the native vertical grid of 72 model layers, or interpolated to 42 standard pressure levels. Unlike MERRA, no data collections are available at the vertical layer edges. More details on the grid are provided in Section 4. MERRA-2 introduced observation-based precipitation forcing for the land surface parameterization and the corresponding variable PRECTOTCORR in the MERRA-2 FLX (surface turbulent fluxes and related quantities) and LFO (land-surface forcing) collections (see Section 6; Reichle et al., 2017). While this variable is still available for M2AMIP, there was no observation-based forcing, making the value identical to the model derived precipitation, PRECTOT. Similarly, without data assimilation, the values for the analysis increments, D*DTANA, in the tendency and vertically integrated file collections are zero. The M2AMIP data are available for download online through the NASA Center for Climate Simulation (NCCS) DataPortal (https://portal.nccs.nasa.gov/datashare/gmao_m2amip/). Data are arranged in subdirectories based on ensemble member, followed by year and month. Control files that are compatible with the Grid Analysis and Display System (GrADS) are available in the ctl_daily and ctl_monthly directories for the hourly, three hourly, and monthly mean data. Control files for the monthly mean diurnal cycle can be found in the ctl_diurnal subdirectory within the directory for each individual ensemble member.

M2AMIP↗

Crop Model Improvement Reduces the Uncertainty of the Response to Temperature of Multi-Model Ensembles

To improve climate change impact estimates and to quantify their uncertainty, multi-model ensembles (MMEs) have been suggested. Model improvements can improve the accuracy of simulations and reduce the uncertainty of climate change impact assessments. Furthermore, they can reduce the number of models needed in a MME. Herein, 15 wheat growth models of a larger MME were improved through re-parameterization and/or incorporating or modifying heat stress effects on phenology, leaf growth and senescence, biomass growth, and grain number and size using detailed field experimental data from the USDA Hot Serial Cereal experiment (calibration data set). Simulation results from before and after model improvement were then evaluated with independent field experiments from a CIMMYT worldwide field trial network (evaluation data set). Model improvements decreased the variation (10th to 90th model ensemble percentile range) of grain yields simulated by the MME on average by 39% in the calibration data set and by 26% in the independent evaluation data set for crops grown in mean seasonal temperatures greater than 24 C. MME mean squared error in simulating grain yield decreased by 37%. A reduction in MME uncertainty range by 27% increased MME prediction skills by 47%. Results suggest that the mean level of variation observed in field experiments and used as a benchmark can be reached with half the number of models in the MME. Improving crop models is therefore important to increase the certainty of model-based impact assessments and allow more practical, i.e. smaller MMEs to be used effectively.

damage assessment↗

An Ensemble-Based Smoother with Retrospectively Updated Weights for Highly Nonlinear Systems

Monte Carlo computational methods have been introduced into data assimilation for nonlinear systems in order to alleviate the computational burden of updating and propagating the full probability distribution. By propagating an ensemble of representative states, algorithms like the ensemble Kalman filter (EnKF) and the resampled particle filter (RPF) rely on the existing modeling infrastructure to approximate the distribution based on the evolution of this ensemble. This work presents an ensemble-based smoother that is applicable to the Monte Carlo filtering schemes like EnKF and RPF. At the minor cost of retrospectively updating a set of weights for ensemble members, this smoother has demonstrated superior capabilities in state tracking for two highly nonlinear problems: the double-well potential and trivariate Lorenz systems. The algorithm does not require retrospective adaptation of the ensemble members themselves, and it is thus suited to a streaming operational mode. The accuracy of the proposed backward-update scheme in estimating non-Gaussian distributions is evaluated by comparison to the more accurate estimates provided by a Markov chain Monte Carlo algorithm.

Monte Carlo↗

A Generalized, Compactly-Supported Correlation Function for Data Assimilation Applications

Correlation functions play an essential role in modern data assimilation, where they are used to model covariances given a set of tunable parameters or applied as tapering functions to localize covariances in ensemble-based schemes. One of the most widely-used correlation functions in data assimilation is the Gaspari and Cohn (1999) piecewise-rational, compactly-supported parametric correlation function (hereafter referred to as GC99). The GC99 correlation function is useful due to its tunable cut-off parameter c and Gaussian-like shape achieved when the parameter a is set to one-half. These properties are attractive for tapering functions in data assimilation applications. However, the GC99 correlation function is homogeneous over Euclidean 3-space and isotropic when restricted to the sphere, properties that may be less than ideal for some geophysical applications. GC99 is also compactly-supported on a sphere of fixed radius, which requires tuning of the cut-off parameter c that can depend on the specific application. This work presents a generalization of the GC99 correlation function that allows the cut-off parameter c and shape parameter a to vary over space to gain more flexibility in shape while maintaining its compact support property. The function, which we call the Generalized Gaspari Cohn (GenGC) correlation function, introduces inhomogeneity in Euclidean 3-space and anisotropy when restricted to the sphere by allowing both parameters c and a to vary, as functions, over the spatial domain. The GC99 correlation function is a special case of GenGC where the functions c and a are held constant, as fixed parameters rather than functions. The GenGC correlation function also generalizes the follow-on to the work of Gaspari and Cohn (1999) presented in Gaspari et al. (2006), which allowed a to vary while keeping c fixed. We illustrate through simple one- and two-dimensional examples the variety of inhomogeneous and anisotropic correlation functions GenGC can produce by varying c and a over space, and suggest applications where they may be useful in data assimilation, such as covariance modeling or localization. In particular, we describe how the GenGC correlation function can be used to construct covariances using correlation length and variance fields derived from dynamics. For example, the correlation length field for advective dynamics is governed by a partial differential equation (PDE) in N spatial dimensions, where N is the number of space dimensions of the state. Correlation length fields can be determined from this PDE and used with GenGC to construct the corresponding correlations. We can then approximate the full covariance by rescaling by the variance, which also satisfies a PDE in N spatial dimensions for advective dynamics. Thus we can approximate the full covariance without solving the covariance PDE, which is in 2N spatial dimensions, by solving just two PDEs each in only N spatial dimensions. This approach to evolving the correlation length and variance fields, then reconstructing the correlations using GenGC, is suggested as an alternative to current methods of covariance modeling in data assimilation algorithms.

GC99↗

Visualizing 2D Probability Distributions from Satellite Image-Derived Data

Creating maps of biophysical and geophysical variables using Earth Observing System (EOS) satellite image data is an important component of Earth science. These 2D maps have a single value at every location and standard techniques are used to visualize them. Current tools fall short, however, when it is necessary to describe a distribution of values at each location. Distributions may represent a frequency of occurrence over time, frequency of occurrence from multiple runs of an ensemble forecast or possible values from an uncertainty model. 'Distribution data sets' are described, then a case study is presented to visualize such 2D distributions. Distribution data sets are different from multivariate data sets in the sense that the values are for a single variable instead of multiple variables. Our case study data consists of multiple realizations of percent forest cover, generated using a geostatistical technique that combines ground measurements and satellite imagery to model uncertainty about forest cover. We present several approaches for analyzing and visualizing such data sets. The first is a pixel-wise analysis of the probability density functions for the 2D image while the second is an analysis of features identified within the image. Such pixel-wise and feature-wise views will give Earth scientists a more complete understanding of distribution data sets.

Kao, David↗

Impact of Satellite Microwave Radiance Data Assimilation on GEOS Atmospheric Analysis and Forecasts in Tropics

The NASA Global Modeling and Assimilation Office (GMAO) has developed a system to assimilate all-sky microwave radiance data in the Goddard Earth Observing System (GEOS). The system provides additional constraints on the analysis process near the storm regions and adjusts the geophysical parameters such as precipitation, cloud, moisture, surface pressure, and wind by combining information from microwave radiance measurements and model forecasts in an optimal manner. The system proved that assimilating the all-sky microwave radiance data improve the GEOS atmospheric analyses and forecasts. This all-sky data framework has been included in the GEOS Forward Processing (FP) system and currently assimilating all-sky GMI data in real-time for GEOS global analysis and forecast production at the GMAO. This presentation describes the methodologies to assimilate cloud- and precipitation-affected microwave radiances in GEOS data assimilation system based on hybrid four-dimensional ensemble-variational (4D-EnVar) configuration. In the current operational GEOS-FP system, the addition of all-sky GMI radiances has the largest impact in the Tropics. Specific humidity is significantly improved in the short term (0-72 hours) forecasts. Similar improvements are seen in the Tropical and lower tropospheric temperature and winds. More detailed results on all-sky microwave radiance data impact on GEOS analysis and forecasts are discussed in this presentation.

Kim, Min-Jeong↗

Gearbox bearing crack growth prognostics and uncertainty quantification with physics-informed machine learning

This paper introduces the extreme theory of functional connections (X-TFC), a physics-informed machine learning algorithm, and tailors it to estimate the remaining useful life (RUL) of wind turbine gearbox bearings experiencing fatigue crack growth. Unlike purely data-driven methods, X-TFC embeds a physics model, based on Head's theory in this work, into its training objective. The core of X-TFC is a random-projection single-layer neural network trained via an extreme learning machine, which requires only limited damage progression data and solves for output weights with a least-squares optimization algorithm. A composite loss function balances the network's fit to observed degradation data against the residuals of the governing crack growth differential equation, ensuring the learned damage trajectory remains physically plausible. When applied to a vibration-based health-index (HI) dataset measured during the growth of a crack on the inner ring of a high-speed bearing in a wind turbine gearbox (Bechhoefer and Dubé, 2020), X-TFC achieves near-zero prediction bias. Even when trained on only the first 10 %–20 % of the damage progression data, with sufficient physics weighting its predictions remain monotonic and smooth, delivering high prognosability and trendability. To quantify the epistemic uncertainty, we employ a Monte Carlo ensemble of independently initialized X-TFC models trained on noise-perturbed data, which yields confidence intervals around each RUL estimate and captures both model-parameter and epistemic uncertainty. In addition to a vibration-based HI, we demonstrate that the proposed framework can be directly applied to a supervisory control and data acquisition (SCADA) data-based HI (Eftekhari Milani et al., 2026) measured during similar wind turbine gearbox bearing crack faults, preserving its accuracy and interpretability. This extension shows the versatility of our approach, which is applicable to bearings of multiple gearbox manufacturers, models, and ratings using only SCADA data. By integrating domain knowledge with machine learning, X-TFC offers a rapid, reliable tool for crack prognostics. Its adaptability to other bearing failure modes, such as pitch bearing ring cracks, positions X-TFC as a powerful enabler of data-driven, physics-informed asset management in the wind energy sector and beyond.

17 WIND ENERGY↗

Multihierarchy Gaussian Process Models for Probabilistic Aerodynamic Databases using Uncertain Nominal and Off-Nominal Configuration Data

Probabilistic aerodynamic databases are a crucial component of the development lifecycle for aerospace vehicles. A key challenge when building aerodynamic databases is that most data used to construct them represent various simplifications of the real flight vehicle. For example, wind tunnel models often simplify the vehicle geometry and surface roughness characteristics, while CFD computations often make simplifications to the physics being modeled, such as fully laminar or turbulent calculations. Multifidelity data fusion models rely on a user being able to define a hierarchy of fidelity levels anchored to some "truth" data. This approach is unsatisfactory when no data can be considered to accurately reflect real flight conditions. In this work, we provide an alternative approach by presenting a consistent mathematical framework for building probabilistic aerodynamic databases in the form of a conditional probability distribution described by an ensemble of multifidelity Gaussian Processes. Instead of relying on a single hierarchy of data fidelity levels, the presented framework identifies a "nominal" configuration and potential corrections to the nominal which represent specific physical phenomena not represented in the nominal data. The nominal and correction functions themselves are constructed as multifidelity Gaussian Processes and linearly combined to form an ensemble model which fuses the uncertainties associated nominal and correction models. Results obtained using the proposed framework on a simplified Orion Crew Module wind tunnel dataset demonstrate the predictive capability of the multihierarchy framework. We further demonstrate the benefits of such a probabilistic aerodynamic database approach through function sampling and computing the conditional distributions of derived quantities, such as the trim angle of attack and aerodynamic coefficients at trim.

Gaussian Processes↗

Flight effects on fan noise with static and wind tunnel comparisons

A flight test program utilizing a JT15D-1 turbofan engine has been conducted with the objectives of studying flight effects on fan noise and evaluating the simulation effectiveness of both a wind tunnel and a static test configuration that incorporated an inlet control device (ICD). In conjunction with synchronized laser-radar tracking and meteorological profile information, data obtained from a linear array of ground microphones was narrowband-analyzed and ensemble-averaged to yield highly accurate far-field flight acoustic results. Utilizing appropriate corrections, flight, wind tunnel, and static acoustic data were normalized to a static-equivalent, 100-foot radius, lossless reference condition. Data comparisons showed that both the static test with ICD and wind tunnel were generally very effective in simulating flight blade-passage-frequency (BPF) noise levels. However, differences were observed in broadband noise levels and in the details of the multiple-pure-tone harmonics.

Preisser, J. S.↗

The first ensemble of kilometer-scale simulations of a hydrological year over the third pole

An accurate understanding of the current and future water cycle over the Third Pole is of great societal importance, given the role this region plays as a water tower for densely populated areas downstream. An emerging and promising approach for skillful climate assessments over regions of complex terrain is kilometer-scale climate modeling. As a foundational step towards such simulations over the Third Pole, we present a multi-model and multi-physics ensemble of kilometer-scale regional simulations for the hydrological year of October 2019 to September 2020. The ensemble consists of 13 simulations performed by an international consortium of 10 research groups, configured with a horizontal grid spacing ranging from 2.2 to 4 km covering all of the Third Pole region. These simulations are driven by ERA5 and are part of a Coordinated Regional Climate Downscaling EXperiment Flagship Pilot Study on Convection-Permitting Third Pole. The simulations are compared against available gridded and in-situ observations and remote-sensing data, to assess the performance and spread of the model ensemble compared to the driving reanalysis during the cold and warm seasons. Although ensemble evaluation is hindered by large differences between the gridded precipitation datasets used as a reference over this region, we show that the ensemble improves on many warm-season precipitation metrics compared with ERA5, including most wet-day and hour statistics, and also adds value in the representation of wet spells in both seasons. As such, the ensemble will provide an invaluable resource for future improvements in the process understanding of the hydroclimate of this remote but important region.

54 ENVIRONMENTAL SCIENCES↗

Impact of Design Factors for ESA CCI Satellite Soil Moisture Data Assimilation over Europe

In this study, soil moisture retrievals of the combined active–passive ESA Climate Change Initiative (CCI) soil moisture product are assimilated into the Noah-MP land surface model over Europe using a one-dimensional ensemble Kalman filter and an 18-yr study period. The performance of the data assimilation (DA) system is evaluated by comparing it with a model-only experiment (at in situ sites) and by assessing statistics of innovations and increments as DA diagnostics (over the entire domain). For both assessments, we explore the impact of three design choices, resulting in the following insights. 1) The magnitude of the assumed observation errors strongly affects the skill improvements evaluated against in situ stations and internal diagnostics. 2) Choosing between climatological or monthly cumulative distribution function matching as the observation bias correction method only has a marginal effect on the in situ skill of the DA system. However, the internal diagnostics suggest a more robust system parameterization if the observations are rescaled monthly. 3) The choice of atmospheric reanalysis dataset to force the land surface model affects the model-only skill and the DA skill improvements. The model-only skill is higher with input from the MERRA-2 than with input from the ERA5 reanalysis, resulting in larger DA skill improvements for the latter. Additionally, we show that the added value of the DA strongly depends on the quality of the satellite retrievals and land cover, with the most substantial soil moisture skill improvements occurring over croplands and skill degradation occurring over densely forested areas.

Soil moisture↗