Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble data assimilation”

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 181 records · Page 10

Assimilation of SMOS Soil Moisture Retrievals in the Land Information System

Soil moisture is a crucial variable for weather prediction because of its influence on evaporation. It is of critical importance for drought and flood monitoring and prediction and for public health applications. The NASA Short-term Prediction Research and Transition Center (SPoRT) has implemented a new module in the NASA Land Information System (LIS) to assimilate observations from the ESA's Soil Moisture and Ocean Salinity (SMOS) satellite. SMOS Level 2 retrievals from the Microwave Imaging Radiometer using Aperture Synthesis (MIRAS) instrument are assimilated into the Noah LSM within LIS via an Ensemble Kalman Filter. The retrievals have a target volumetric accuracy of 4% at a resolution of 35-50 km. Parallel runs with and without SMOS assimilation are performed with precipitation forcing from intentionally degraded observations, and then validated against a model run using the best available precipitation data, as well as against selected station observations. The goal is to demonstrate how SMOS data assimilation can improve modeled soil states in the absence of dense rain gauge and radar networks.

Blakenship, Clay↗

Assimilation of Lidar Planetary Boundary Layer Height Observations

Lidar backscatter and wind retrievals of the planetary boundary layer height (PBLH) are assimilated into 22 hourly forecasts from the NASA Unified - Weather and Research Forecast (NU-WRF) model during the Plains Elevated Convection Convection at Night (PECAN) campaign on July 11, 2015 in Greensburg, Kansas, using error statistics collected from the model profiles to compute the necessary covariance matrices. Two separate forecast runs using different PBL physics schemes were employed, and comparisons with 6 independent radiosonde profiles were made for each run. Both of the forecast runs accurately predicted the PBLH and the state variable profiles within the planetary boundary layer during the early morning, and the assimilation had a small impact during this time. In the late afternoon, the forecast runs showed decreased accuracy as the convective boundary layer developed. However, assimilation of the Doppler lidar PBLH observations were found to improve the temperature and V velocity profiles relative to independent radiosonde profiles. Water vapor was overcorrected, leading to increased differences with independent data. Errors in the U velocity were made slightly larger. The computed forecast error covariances between the PBLH and state variables were found to rise in the late afternoon, leading to the larger improvements in the afternoon. This work represents the first effort to assimilate PBLH into forecast states using ensemble methods.

Andrew Tangborn↗

Recent Updates in the SMAP Level-4 Soil Moisture Algorithm

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level-4 Soil Moisture (L4_SM) product. The L4_SM data are published with a mean latency of ~2.5 days from the time of observation and provide global, three-hourly, 9-km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. The L4_SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter (EnKF). In 2018, the L4_SM algorithm was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high-quality, space-borne remote sensing observations. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of surface soil moisture from below under non-equilibrium conditions, which brings the model's surface soil moisture more in line with the SMAP Level-2 and in situ soil moisture. Moreover, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4, and the "catchment deficit" model variable was removed from the EnKF state vector to avoid degrading the model's groundwater estimates.Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation-minus-forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4_SM product and the lessons learned from the effort.

Reichle, Rolf H.↗

Quantifying Errors in TRMM-Based Multi-Sensor QPE Products Over Land in Preparation for GPM

Determining uncertainties in satellite-based multi-sensor quantitative precipitation estimates over land of fundamental importance to both data producers and hydro climatological applications. ,Evaluating TRMM-era products also lays the groundwork and sets the direction for algorithm and applications development for future missions including GPM. QPE uncertainties result mostly from the interplay of systematic errors and random errors. In this work, we will synthesize our recent results quantifying the error characteristics of satellite-based precipitation estimates. Both systematic errors and total uncertainties have been analyzed for six different TRMM-era precipitation products (3B42, 3B42RT, CMORPH, PERSIANN, NRL and GSMap). For systematic errors, we devised an error decomposition scheme to separate errors in precipitation estimates into three independent components, hit biases, missed precipitation and false precipitation. This decomposition scheme reveals hydroclimatologically-relevant error features and provides a better link to the error sources than conventional analysis, because in the latter these error components tend to cancel one another when aggregated or averaged in space or time. For the random errors, we calculated the measurement spread from the ensemble of these six quasi-independent products, and thus produced a global map of measurement uncertainties. The map yields a global view of the error characteristics and their regional and seasonal variations, reveals many undocumented error features over areas with no validation data available, and provides better guidance to global assimilation of satellite-based precipitation data. Insights gained from these results and how they could help with GPM will be highlighted.

Peters-Lidard, Christa D.↗

Assessing the Performance of the South American Land Data Assimilation System Version 2 (SALDAS-2) Energy Balance across Diverse Biomes

Understanding the exchange of energy between the surface and the atmosphere is important in view of the climate scenario. However, it becomes a challenging task due to a sparse network of observations. This study aims to improve the energy balance estimates for the Amazon, Cerrado, and Pampa biomes located in South America using the radiation and precipitation forcing obtained from the Clouds and the Earth’s Radiant Energy System (CERES) and the precipitation CPTEC/MERGE datasets. We employed three surface models—Noah-MP, Community Land Model (CLSM), and Integrated Biosphere Simulator (IBIS)—and conducted modeling experiments, termed South America Land Data Assimilation System (SALDAS-2). The results showed that SALDAS-2 radiation estimates had the smallest errors. Moreover, SALDAS-2 precipitation estimates were better than the Global Land Data Assimilation System (GLDAS) in the Cerrado (MBE = −0.16) and Pampa (MBE = −0.19). Noah-MP presented improvements compared with CLSM and IBIS in 100% of towers located in the Amazon. CLSM tends to overestimate the latent heat flux and underestimate the sensible heat flux in the Amazon. Noah-MP and Ensemble outperformed GLDAS in terms latent and sensible heat fluxes. The potential of SALDAS-2 should be emphasized to provide more accurate estimates of surface energy balance.

modeling↗

Advancing Drought Understanding, Monitoring and Prediction

Having the capacity to monitor droughts in near-real time and providing accurate drought prediction from weeks to seasons in advance can greatly reduce the severity of social and economic damage caused by drought, a leading natural hazard for North America. The congressional mandate to establish the National Integrated Drought Information System (NIDIS; Public Law 109-430) in 2006 was a major impulse to develop, integrate, and provide drought information to meet the challenges posed by this hazard. Significant progress has been made on many fronts. On the research front, efforts by the broad scientific community have resulted in improved understanding of North American droughts and improved monitoring and forecasting tools. We now have a better understanding of the droughts of the twentieth century including the 1930s "Dust Bowl"; we have developed a broader array of tools and datasets that enhance the official North American Drought Monitor based on different methodologies such as state-of-the-art land surface modeling (e.g., the North American Land Data Assimilation System) and remote sensing (e.g., the evaporative stress index) to better characterize the occurrence and severity of drought in its multiple manifestations. In addition, we have new tools for drought prediction [including the new National Centers for Environmental Prediction (NCEP) Climate Forecast System, version 2, for operational prediction and an experimental National Multimodel Ensemble] and have explored diverse methodologies including ensemble hydrologic prediction approaches. Broad NIDIS-inspired progress is influencing the development of a Global Drought Information System (GDIS) under the auspices of the World Climate Research Program. Despite these advances, current drought monitoring and forecasting capabilities still fall short of users' needs, especially the need for skillful and reliable drought forecasts at regional and local scales. To tackle this outstanding challenging problem, focused and coordinated research efforts are needed, drawing from excellence across the broad drought research community. To meet this challenge, National Oceanic and Atmospheric Administration (NOAA)'s Drought Task Force was established in October 2011 with the ambitious goal of achieving significant new advances in the ability to understand, monitor, and predict drought over North America. The Task Force (duration of October 2011-September 2014) is an initiative of NOAA's Climate Program Office Modeling, Analysis, Predictions, and Projections (MAPP) program in partnership with NIDIS. It brings together over 30 leading MAPP-funded drought scientists from multiple academic and federal institutions [involves scientists from NOAA's research laboratories and centers, the National Aeronautics and Space Administration (NASA), U.S. Department of Agriculture, National Center for Atmospheric Research (NCAR), and many universities] in a concerted research effort that builds on individual MAPP research projects. These projects span the wide spectrum of drought research needed to make fundamental advances, from those aimed at the basic understanding of drought mechanisms to those aimed at testing new drought monitoring and prediction tools for operational and service purposes (as part of NCEP's Climate Test Bed). The Drought Task Force provides focus and coordination to MAPP drought research activities and also facilitates synergies with other national and international drought research efforts, including those by the GDIS.

GDIS↗

El Nino-Induced Tropical Ocean/Land Energy Exchange in MERRA-2 and M2AMIP

Studies have shown the correlation and connection of surface temperatures across the globe, ocean and land, related to Tropical SSTs especially El Nino. This climate variability greatly influences regional weather and hydroclimate extremes (e.g. drought and flood). In this paper, we evaluate the relationship of temperatures across the tropical oceans and continents in MERRA-2, and also in a newly developed MERRA-2 AMIP ensemble simulation (M2AMIP). M2AMIP uses the same model and spatial resolution as MERRA-2, producing the same output diagnostics over 10 ensemble members. Composite El Nino temperature data are compared with observations to evaluate the land/sea contrast, variations and phase relationship. The temperature variations are related to surface heat fluxes and the atmospheric temperatures and transport, to identify the processes that lead to the lagged redistribution of heat in the tropics and beyond. Discernable cloud, radiation and data assimilation changes accompany the onset of El Nino affecting continental regions through the progression to and following the peak values. While the model represents these variations in general, regional strengths and weaknesses can be identified.

Bosilovich, Michael G.↗

Improving the Representation of Land Surface Processes Using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here, we emphasize the implementation of tools and approaches to overcome challenges related to land DA to constrain carbon and water cycling. In particular, we discuss the implementation of adaptive inflation to modify ensemble spread in response to time-varying networks of gridded observations. We also discuss methods to generate ensemble spread through boundary condition (meteorology) forcing that can be applied to site-level applications. Next, we describe the application of vertical localization upon surface soil moisture observations, and forward operators specifically designed for the assimilation of snow and solar-induced fluorescence observations. Finally, we discuss the potential benefit of a quantile conserving filter used to update bounded quantities (state or parameter values).

Brett Raczka↗

Version 4 of the SMAP Level-4 Soil Moisture Algorithm and Data Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 3-hourly, 9-km resolution estimates of surface (0-5 cm) and root-zone (0-100 cm) soil moisture with a mean latency of ~2.5 days. The underlying L4_SM algorithm assimilates SMAP radiometer brightness temperature (Tb) observations into the NASA Catchment land surface model using a spatially-distributed ensemble Kalman filter. Version 4 of the L4_SM modeling system includes a reduction in the upward recharge of surface soil moisture from below under non-equilibrium conditions, resulting in reduced bias and improved dynamic range of L4_SM surface soil moisture compared to earlier versions. This change and additional technical modifications to the system reduce the mean and standard deviation of the observation-minus-forecast Tb residuals and overall soil moisture analysis increments while maintaining the skill of the L4_SM soil moisture estimates versus independent in situ measurements; the average, bias-adjusted RMSE in Version 4 is 0.039 m(exp 3) m(exp -3) for surface and 0.026 m(exp 3) m(exp -3) for root-zone soil moisture. Moreover, the coverage of assimilated SMAP observations in Version 4 is near-global owing to the use of additional satellite Tb records for algorithm calibration. L4_SM soil moisture uncertainty estimates are biased low (by 0.01-0.02 m(exp 3) m(exp -3)) against actual errors (computed versus in situ measurements). L4_SM runoff estimates, an additional product of the L4_SM algorithm, are biased low (by 35 mm year (exp -1)) against streamflow measurements. Compared to Version 3, bias in Version 4 is reduced by 46% for surface soil moisture uncertainty estimates and by 33% for runoff estimates.

RMSE↗

The Impact of AMSR-E Soil Moisture Assimilation on Evapotranspiration Estimation

An assessment ofETestimates for current LDAS systems is provided along with current research that demonstrates improvement in LSM ET estimates due to assimilating satellite-based soil moisture products. Using the Ensemble Kalman Filter in the Land Information System, we assimilate both NASA and Land Parameter Retrieval Model (LPRM) soil moisture products into the Noah LSM Version 3.2 with the North American LDAS phase 2 CNLDAS-2) forcing to mimic the NLDAS-2 configuration. Through comparisons with two global reference ET products, one based on interpolated flux tower data and one from a new satellite ET algorithm, over the NLDAS2 domain, we demonstrate improvement in ET estimates only when assimilating the LPRM soil moisture product.

Peters-Lidard, Christa D.↗

A Method for Snow Reanalysis: The Sierra Nevada (USA) Example

This work presents a state-of-the art methodology for constructing snow water equivalent (SWE) reanalysis. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data from the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlation greater than 0.95 compared with in situ SWE observations. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (31 years), and accuracy provide a unique dataset for investigating snow processes. This presentation illustrates how the reanalysis dataset was used to provide a basic accounting of the stored snowpack water in the Sierra Nevada over the last 31 years and ultimately improve real-time streamflow predictions.

Bayesian↗

Estimation of Seasonal Snow Water Equivalent Using Landsat Observations

This work presents a methodology for estimating seasonal snow water equivalent (SWE) from the use of remotely sensed Visible and Near Infrared observations from the Landsat mission. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data over the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlations with in situ SWE observations of greater than 0.95. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (32 years), and accuracy provide a unique dataset for investigating snow processes. In particular, this presentation illustrates how the reanalysis dataset was used to provide climatology of the seasonal snowfall accumulation rates, distributions, and variability over the last three decades.

SWE↗

The Representation of Aerosols in GMAO’s Newest Reanalyses

Over the past few years, NASA’s Global Modelling and Assimilation Office has been working on the configuration and production of three new reanalysis products, GEOS-IT, GiOcean, and MERRA-21C. GEOS-IT, or the Goddard Earth Observing System for Instrument Teams, is a 3D variational data assimilation system that runs in a near real time framework however retrospectively provides data back through 1998 to deliver a consistent view of the Earth-atmosphere system for the production of observational NASA products. Retrospective production for GEOS-IT is complete and the meteorology has since been used to produce the one way weakly coupled GiOcean reanalysis. Due to differences in the atmospheric model, particularly related to scavenging, aerosols are not identical in GEOS-IT and GiOcean. MERRA-21C, or the Modern Era Restrospective analysis for Research and Applications in the 21st century, is a hybrid 4D ensemble variational system at a finer horizontal resolution of 0.25 degrees. Although different in their intended use, and therefore configuration, these systems prominently feature coupling between meteorology and aerosols. The differences and similarities in the aerosol configuration between the three systems will be discussed, covering biomass burning and anthropogenic emissions as well as observations used for the assimilation of aerosol optical depth. A large emphasis will be placed on the version of the underlying aerosol module, GOCART, which underwent a complete refactoring and the addition of radiatively active brown carbon between GEOS-IT and MERRA-21C. Independent observations will be used to evaluate the performance of aerosols in both reanalyses, focusing on aerosol optical depth, surface particulate matter, and vertical profiles of aerosol backscatter.

Allison Collow↗

Assimilation of Sentinel 1 and SMAP–based satellite soil moisture retrievals into SWAT hydrological model: the impact of satellite revisit time and product spatial resolution on flood simulations in small basins

In runoff generation process, soil moisture plays an important role as it controls the magnitude of the flood events in response to the rainfall inputs. In this study, we investigated the ability of a new era of satellite soil moisture retrievals to improve the Soil & Water Assessment Tool (SWAT) daily discharge simulations via soil moisture data assimilation for two small (< 500 km2) and hydrologically different catchments located in Central Italy. We ingested 1) the Soil Moisture Active and Passive (SMAP) Enhanced L3 Radiometer Global Daily 9 km EASE-Grid soil moisture, 2) the Advanced SCATterometer (ASCAT) H113 soil moisture product released within the EUMETSAT Satellite Application Facility on Support to Operational Hydrology and Water Management (HSAF) which has a nearly daily temporal resolution and sampling of 12.5 km, and 3) a fused ASCAT/Sentinel-1 (S1) satellite soil moisture product named SCATSAR-SWI with temporal and spatial sampling of 1 day and 1 km, respectively into SWAT hydrological model via the Ensemble Kalman Filter (EnKF). Different configurations were tested with the aim of exploring the effect of the hydrological regime, the land use conditions, the spatial sampling and the revisit time of the products (which controls the amount of available data to be potentially ingested). Results show a general improvement of SWAT discharge simulations for all products in terms of error and Nash Sutcliffe efficiency index. In particular, we found a relatively good behavior of both the active and the passive products in terms of low flows improvement especially for the catchment characterized by a higher baseflow component. The benefit of the higher spatial resolution of SCATSAR-SWI obtained via S1 over ASCAT was small, likely due to very challenging areas for the S1 retrieval. Eventually, better performances were obtained for the passive product in the more forested catchment. With the aim of exploring the benefit of having more frequent satellite soil moisture observations to be ingested, we tested the performance of the ASCAT product with a reduced temporal sampling obtained by temporally matching ASCAT observations to that of SMAP. The results show a significant reduction of the performance of ASCAT, suggesting that the correction frequency (due to the higher number of observations available) for small catchments is an important aspect for improving flood forecasting as it helps to adjust more frequently the pre-storm soil moisture conditions.

Shima Azimi↗

Assimilation of Sentinel 1 and SMAP – Based Satellite Soil Moisture Retrievals into SWAT Hydrological Model: the Impact of Satellite Revisit Time and Product Spatial Resolution on Flood Simulations in Small Basins

In runoff generation process, soil moisture plays an important role as it controls the magnitude of the flood events to the rainfall inputs. In this study, we investigated the ability of a new era of satellite soil moisture retrievals to improve the Soil & Water Assessment Tool (SWAT) daily discharge simulations via soil moisture data assimilation for two small (< 500 km2) and hydrologically different catchments located in Central Italy. We ingested 1) the Soil Moisture Active and Passive (SMAP) Enhanced L3 Radiometer Global Daily 9 km EASE-Grid soil moisture, 2) the Advanced SCATterometer (ASCAT) H113 soil moisture product released within the EUMETSAT Satellite Application Facility on Support to Operational Hydrology and Water Management (H-SAF) which has a nearly daily temporal resolution and sampling of 12.5 km, and 3) a fused ASCAT/Sentinel-1 (S1) satellite soil moisture product named SCATSAR-SWI with temporal and spatial sampling of 1 day and 1 km, respectively into SWAT hydrological model via Ensemble Kalman Filter (EnKF). Different configurations were tested with the aim of exploring the effect of the hydrological regime, the land use conditions, the spatial sampling and the revisit time of the products (which controls the amount of available data to be potentially ingested). Results show a general improvement of SWAT discharge simulations for all products in terms of error and Nash Sutcliffe efficiency index. In particular, we found a relatively good behavior of both the active and the passive products in terms of low flows improvement especially for the catchment characterized by a higher baseflow component. The benefit of the higher spatial resolution of SCATSAR obtained via S1 over ASCAT was small likely due to very challenging areas for S1 retrieval. Eventually, better performances were obtained for the passive product in the more forested catchment. With the aim of exploring the benefit of having more frequent satellite soil moisture observations to be ingested, we tested the performance of the ASCAT product with a reduced temporal sampling obtained by temporally matching ASCAT observations to that of SMAP. The results show a significant reduction of the performance of ASCAT, suggesting that the correction frequency (due to the higher number of observations available) for small catchments is an important aspect for improving flood forecasting as it helps to adjust more frequently the soil moisture pre-storm conditions.

Small basins↗

Basin-Scale Assessment of the Land Surface Water Budget in the National Centers for Environmental Prediction Operational and Research NLDAS-2 Systems

The purpose of this study is to evaluate the components of the land surface water budget in the four land surface models (Noah, SAC-Sacramento Soil Moisture Accounting Model, (VIC) Variable Infiltration Capacity Model, and Mosaic) applied in the newly implemented National Centers for Environmental Prediction (NCEP) operational and research versions of the North American Land Data Assimilation System version 2 (NLDAS-2). This work focuses on monthly and annual components of the water budget over 12 National Weather Service (NWS) River Forecast Centers (RFCs). Monthly gridded FLUX Network (FLUXNET) evapotranspiration (ET) from the Max-Planck Institute (MPI) of Germany, U.S. Geological Survey (USGS) total runoff (Q), changes in total water storage (dS/dt, derived as a residual by utilizing MPI ET and USGS Q in the water balance equation), and Gravity Recovery and Climate Experiment (GRACE) observed total water storage anomaly (TWSA) and change (TWSC) are used as reference data sets. Compared to these ET and Q benchmarks, Mosaic and SAC (Noah and VIC) in the operational NLDAS-2 overestimate (underestimate) mean annual reference ET and underestimate (overestimate) mean annual reference Q. The multimodel ensemble mean (MME) is closer to the mean annual reference ET and Q. An anomaly correlation (AC) analysis shows good AC values for simulated monthly mean Q and dS/dt but significantly smaller AC values for simulated ET. Upgraded versions of the models utilized in the research side of NLDAS-2 yield largely improved performance in the simulation of these mean annual and monthly water component diagnostics. These results demonstrate that the three intertwined efforts of improving (1) the scientific understanding of parameterization of land surface processes, (2) the spatial and temporal extent of systematic validation of land surface processes, and (3) the engineering-oriented aspects such as parameter calibration and optimization are key to substantially improving product quality in various land data assimilation systems.

NLDAS-2↗

Connecting Satellite Observations with Water Cycle Variables Through Land Data Assimilation: Examples Using the NASA GEOS-5 LDAS

A land data assimilation system (LDAS) can merge satellite observations (or retrievals) of land surface hydrological conditions, including soil moisture, snow, and terrestrial water storage (TWS), into a numerical model of land surface processes. In theory, the output from such a system is superior to estimates based on the observations or the model alone, thereby enhancing our ability to understand, monitor, and predict key elements of the terrestrial water cycle. In practice, however, satellite observations do not correspond directly to the water cycle variables of interest. The present paper addresses various aspects of this seeming mismatch using examples drawn from recent research with the ensemble-based NASA GEOS-5 LDAS. These aspects include (1) the assimilation of coarse-scale observations into higher-resolution land surface models, (2) the partitioning of satellite observations (such as TWS retrievals) into their constituent water cycle components, (3) the forward modeling of microwave brightness temperatures over land for radiance-based soil moisture and snow assimilation, and (4) the selection of the most relevant types of observations for the analysis of a specific water cycle variable that is not observed (such as root zone soil moisture). The solution to these challenges involves the careful construction of an observation operator that maps from the land surface model variables of interest to the space of the assimilated observations.

Land data assimilation↗

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗