Engineering topics
Reichle, R. H.
Publications and source records attributed to Reichle, R. H..
Effective Assimilation of SMAP Observations Using Statistical Techniques
Statistical techniques permit the retrieval of soil moisture estimates in a model climatology while retaining the spatial and temporal signatures of the satellite observations. As a consequence, the need for bias correction prior to an assimilation of these estimates is reduced, which could result in a more effective use of the independent information provided by the satellite observations. In this study, a statistical neural network (NN) retrieval algorithm is calibrated using SMAP brightness temperature observations and modeled soil moisture estimates (similar to those used to calibrate the SMAP Level 4 DA system). Daily values of surface soil moisture are estimated using the NN and then assimilated into the NASA Catchment model. The skill of the assimilation estimates is assessed based on a comprehensive comparison to in situ measurements from the SMAP core and sparse network sites as well as the International Soil Moisture Network. The NN retrieval assimilation is found to significantly improve the model skill, particularly in areas where the model does not represent processes related to agricultural practices. Additionally, the NN method is compared to assimilation experiments using traditional bias correction techniques. The NN retrieval assimilation is found to more effectively use the independent information provided by SMAP resulting in larger model skill improvements than assimilation experiments using traditional bias correction techniques.
Improving Water Level and Soil Moisture Over Peatlands in a Global Land Modeling System
New model structure for peatlands results in improved skill metrics (without any parameter calibration) Simulated surface soil moisture strongly affected by new model, but reliable soil moisture data lacking for validation.
The SMAP Level-4 ECO Project: Linking the Terrestrial Water and Carbon Cycles
The SMAP (Soil Moisture Active Passive) Level-4 projects aims to develop a fully coupled hydrology-vegetation data assimilation algorithm to generate improved estimates of modeled hydrological fields and carbon fluxes. This includes using the new NASA Catchment-CN (Catchment-Carbon-Nitrogen) model, which combines the Catchment land surface hydrology model with dynamic vegetation components from the Community Land Model version 4 (CLM4). As such, Catchment-CN allows a more realistic, fully coupled feedback between the land hydrology and the biosphere. The L4 ECO project further aims to inform the model through the assimilation of Soil Moisture Active Passive (SMAP) brightness temperature observations as well as observations of Moderate Resolution Imaging Spectroradiometer (MODIS) fraction of absorbed photosynthetically active radiation (FPAR). Preliminary results show that the assimilation of SMAP observations leads to consistent improvements in the model soil moisture skill. An evaluation of the Catchment-CN modeled vegetation characteristics showed that a calibration of the model's vegetation parameters is required before an assimilation of MODIS FPAR observations is feasible.
Data Assimilation to Extract Soil Moisture Information From SMAP Observations
Statistical techniques permit the retrieval of soil moisture estimates in a model climatology while retaining the spatial and temporal signatures of the satellite observations. As a consequence, they can be used to reduce the need for localized bias correction techniques typically implemented in data assimilation (DA) systems that tend to remove some of the independent information provided by satellite observations. Here, we use a statistical neural network (NN) algorithm to retrieve SMAP (Soil Moisture Active Passive) surface soil moisture estimates in the climatology of the NASA Catchment land surface model. Assimilating these estimates without additional bias correction is found to significantly reduce the model error and increase the temporal correlation against SMAP CalVal in situ observations over the contiguous United States. A comparison with assimilation experiments using traditional bias correction techniques shows that the NN approach better retains the independent information provided by the SMAP observations and thus leads to larger model skill improvements during the assimilation. A comparison with the SMAP Level 4 product shows that the NN approach is able to provide comparable skill improvements and thus represents a viable assimilation approach.
Assimilation of Global Radar Backscatter and Radiometer Brightness Temperature Observations to Improve Soil Moisture and Land Evaporation Estimates
Active radar backscatter (s◦) observations from the Advanced Scatterometer (ASCAT) and passive radiometer brightness temperature (TB) observations from the Soil Moisture Ocean Salinity (SMOS) mission are assimilated either individually or jointly into the Global Land Evaporation Amsterdam Model (GLEAM) to improve its simulations of soil moisture and land evaporation. To enable s◦ and TB assimilation, GLEAM is coupled to the Water Cloud Model and the L-band Microwave Emission from the Biosphere (L-MEB) model. The innovations, i.e. differences between observations and simulations, are mapped onto the model soil moisture states through an Ensemble Kalman Filter. The validation of surface (0-10 cm) soil moisture simulations over the period 2010-2014 against in situ measurements from the International Soil Moisture Network (ISMN) shows that assimilating s◦ or TB alone improves the average correlation of seasonal anomalies (Ran) from 0.514 to 0.547 and 0.548, respectively. The joint assimilation further improves Ran to 0.559. Associated enhancements in daily evaporative flux simulations by GLEAM are validated based on measurements from 22 FLUXNET stations. Again, the singular assimilation improves Ran from 0.502 to 0.536 and 0.533, respectively for s◦ and TB, whereas the best performance is observed for the joint assimilation (Ran = 0.546). These results demonstrate the complementary value of assimilating radar backscatter observations together with brightness temperatures for improving estimates of hydrological variables, as their joint assimilation outperforms the assimilation of each observation type separately.
The SMAP Level 4 Carbon PRODUCT for Monitoring Terrestrial Ecosystem-Atmosphere CO2 Exchange
The NASA Soil Moisture Active Passive (SMAP) mission Level 4 Carbon (L4_C) product provides model estimates of Net Ecosystem CO2 exchange (NEE) incorporating SMAP soil moisture information as a primary driver. The L4_C product provides NEE, computed as total respiration less gross photosynthesis, at a daily time step and approximate 14-day latency posted to a 9-km global grid summarized by plant functional type. The L4_C product includes component carbon fluxes, surface soil organic carbon stocks, underlying environmental constraints, and detailed uncertainty metrics. The L4_C model is driven by the SMAP Level 4 Soil Moisture (L4_SM) data assimilation product, with additional inputs from the Goddard Earth Observing System, Version 5 (GEOS-5) weather analysis and Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. The L4_C data record extends from March 2015 to present with ongoing production. Initial comparisons against global CO2 eddy flux tower measurements, satellite Solar Induced Canopy Florescence (SIF) and other independent observation benchmarks show favorable L4_C performance and accuracy, capturing the dynamic biosphere response to recent weather anomalies and demonstrating the value of SMAP observations for monitoring of global terrestrial water and carbon cycle linkages.
Evaluating the Utility of Satellite Soil Moisture Retrievals over Irrigated Areas and the Ability of Land Data Assimilation Methods to Correct for Unmodeled Processes
Earth's land surface is characterized by tremendous natural heterogeneity and human-engineered modifications, both of which are challenging to represent in land surface models. Satellite remote sensing is often the most practical and effective method to observe the land surface over large geographical areas. Agricultural irrigation is an important human-induced modification to natural land surface processes, as it is pervasive across the world and because of its significant influence on the regional and global water budgets. In this article, irrigation is used as an example of a human-engineered, often unmodeled land surface process, and the utility of satellite soil moisture retrievals over irrigated areas in the continental US is examined. Such retrievals are based on passive or active microwave observations from the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E), the Advanced Microwave Scanning Radiometer 2 (AMSR2), the Soil Moisture Ocean Salinity (SMOS) mission, WindSat and the Advanced Scatterometer (ASCAT). The analysis suggests that the skill of these retrievals for representing irrigation effects is mixed, with ASCAT-based products somewhat more skillful than SMOS and AMSR2 products. The article then examines the suitability of typical bias correction strategies in current land data assimilation systems when unmodeled processes dominate the bias between the model and the observations. Using a suite of synthetic experiments that includes bias correction strategies such as quantile mapping and trained forward modeling, it is demonstrated that the bias correction practices lead to the exclusion of the signals from unmodeled processes, if these processes are the major source of the biases. It is further shown that new methods are needed to preserve the observational information about unmodeled processes during data assimilation.
Consistency of Estimated Global Water Cycle Variations Over the Satellite Era
Motivated by the question of whether recent indications of decadal climate variability and a possible "climate shift" may have affected the global water balance, we examine evaporation minus precipitation (E-P) variability integrated over the global oceans and global land from three points of view-remotely sensed retrievals / objective analyses over the oceans, reanalysis vertically-integrated moisture convergence (MFC) over land, and land surface models forced with observations-based precipitation, radiation and near-surface meteorology. Because monthly variations in area-averaged atmospheric moisture storage are small and the global integral of moisture convergence must approach zero, area-integrated E-P over ocean should essentially equal precipitation minus evapotranspiration (P-ET) over land (after adjusting for ocean and land areas). Our analysis reveals considerable uncertainty in the decadal variations of ocean evaporation when integrated to global scales. This is due to differences among datasets in 10m wind speed and near-surface atmospheric specific humidity (2m qa) used in bulk aerodynamic retrievals. Precipitation variations, all relying substantially on passive microwave retrievals over ocean, still have uncertainties in decadal variability, but not to the degree present with ocean evaporation estimates. Reanalysis MFC and P-ET over land from several observationally forced diagnostic and land surface models agree best on interannual variations. However, upward MFC (i.e. P-ET) reanalysis trends are likely related in part to observing system changes affecting atmospheric assimilation models. While some evidence for a low-frequency E-P maximum near 2000 is found, consistent with a recent apparent pause in sea-surface temperature (SST) rise, uncertainties in the datasets used here remain significant. Prospects for further reducing uncertainties are discussed. The results are interpreted in the context of recent climate variability (Pacific Decadal Oscillation, Atlantic Meridional Overturning), and efforts to distinguish these modes from longer-term trends.
Consistency of Estimated Global Water Cycle Variations Over the Satellite Era
Motivated by the question of whether recent indications of decadal climate variability and a possible "climate shift" may have affected the global water balance, we examine evaporation minus precipitation (E-P) variability integrated over the global oceans and global land from three points of view-remotely sensed retrievals / objective analyses over the oceans, reanalysis vertically-integrated moisture convergence (MFC) over land, and land surface models forced with observations-based precipitation, radiation and near-surface meteorology. Because monthly variations in area-averaged atmospheric moisture storage are small and the global integral of moisture convergence must approach zero, area-integrated E-P over ocean should essentially equal precipitation minus evapotranspiration (P-ET) over land (after adjusting for ocean and land areas). Our analysis reveals considerable uncertainty in the decadal variations of ocean evaporation when integrated to global scales. This is due to differences among datasets in 10m wind speed and near-surface atmospheric specific humidity (2m qa) used in bulk aerodynamic retrievals. Precipitation variations, all relying substantially on passive microwave retrievals over ocean, still have uncertainties in decadal variability, but not to the degree present with ocean evaporation estimates. Reanalysis MFC and P-ET over land from several observationally forced diagnostic and land surface models agree best on interannual variations. However, upward MFC (i.e. P-ET) reanalysis trends are likely related in part to observing system changes affecting atmospheric assimilation models. While some evidence for a low-frequency E-P maximum near 2000 is found, consistent with a recent apparent pause in sea-surface temperature (SST) rise, uncertainties in the datasets used here remain significant. Prospects for further reducing uncertainties are discussed. The results are interpreted in the context of recent climate variability (Pacific Decadal Oscillation, Atlantic Meridional Overturning), and efforts to distinguish these modes from longer-term trends.
The Impact of Model and Rainfall Forcing Errors on Characterizing Soil Moisture Uncertainty in Land Surface Modeling
The contribution of rainfall forcing errors relative to model (structural and parameter) uncertainty in the prediction of soil moisture is investigated by integrating the NASA Catchment Land Surface Model (CLSM), forced with hydro-meteorological data, in the Oklahoma region. Rainfall-forcing uncertainty is introduced using a stochastic error model that generates ensemble rainfall fields from satellite rainfall products. The ensemble satellite rain fields are propagated through CLSM to produce soil moisture ensembles. Errors in CLSM are modeled with two different approaches: either by perturbing model parameters (representing model parameter uncertainty) or by adding randomly generated noise (representing model structure and parameter uncertainty) to the model prognostic variables. Our findings highlight that the method currently used in the NASA GEOS-5 Land Data Assimilation System to perturb CLSM variables poorly describes the uncertainty in the predicted soil moisture, even when combined with rainfall model perturbations. On the other hand, by adding model parameter perturbations to rainfall forcing perturbations, a better characterization of uncertainty in soil moisture simulations is observed. Specifically, an analysis of the rank histograms shows that the most consistent ensemble of soil moisture is obtained by combining rainfall and model parameter perturbations. When rainfall forcing and model prognostic perturbations are added, the rank histogram shows a U-shape at the domain average scale, which corresponds to a lack of variability in the forecast ensemble. The more accurate estimation of the soil moisture prediction uncertainty obtained by combining rainfall and parameter perturbations is encouraging for the application of this approach in ensemble data assimilation systems.
Assimilation of Passive and Active Microwave Soil Moisture Retrievals
Root-zone soil moisture is an important control over the partition of land surface energy and moisture, and the assimilation of remotely sensed near-surface soil moisture has been shown to improve model profile soil moisture [1]. To date, efforts to assimilate remotely sensed near-surface soil moisture at large scales have focused on soil moisture derived from the passive microwave Advanced Microwave Scanning Radiometer (AMSR-E) and the active Advanced Scatterometer (ASCAT; together with its predecessor on the European Remote Sensing satellites (ERS. The assimilation of passive and active microwave soil moisture observations has not yet been directly compared, and so this study compares the impact of assimilating ASCAT and AMSR-E soil moisture data, both separately and together. Since the soil moisture retrieval skill from active and passive microwave data is thought to differ according to surface characteristics [2], the impact of each assimilation on the model soil moisture skill is assessed according to land cover type, by comparison to in situ soil moisture observations.
Global Intercomparison of 12 Land Surface Heat Flux Estimates
A global intercomparison of 12 monthly mean land surface heat flux products for the period 1993-1995 is presented. The intercomparison includes some of the first emerging global satellite-based products (developed at Paris Observatory, Max Planck Institute for Biogeochemistry, University of California Berkeley, University of Maryland, and Princeton University) and examples of fluxes produced by reanalyses (ERA-Interim, MERRA, NCEP-DOE) and off-line land surface models (GSWP-2, GLDAS CLM/ Mosaic/Noah). An intercomparison of the global latent heat flux (Q(sub le)) annual means shows a spread of approx 20 W/sq m (all-product global average of approx 45 W/sq m). A similar spread is observed for the sensible (Q(sub h)) and net radiative (R(sub n)) fluxes. In general, the products correlate well with each other, helped by the large seasonal variability and common forcing data for some of the products. Expected spatial distributions related to the major climatic regimes and geographical features are reproduced by all products. Nevertheless, large Q(sub le)and Q(sub h) absolute differences are also observed. The fluxes were spatially averaged for 10 vegetation classes. The larger Q(sub le) differences were observed for the rain forest but, when normalized by mean fluxes, the differences were comparable to other classes. In general, the correlations between Q(sub le) and R(sub n) were higher for the satellite-based products compared with the reanalyses and off-line models. The fluxes were also averaged for 10 selected basins. The seasonality was generally well captured by all products, but large differences in the flux partitioning were observed for some products and basins.
Assimilation of Terrestrial Water Storage from GRACE in a Snow-Dominated Basin
Terrestrial water storage (TWS) information derived from Gravity Recovery and Climate Experiment (GRACE) measurements is assimilated into a land surface model over the Mackenzie River basin located in northwest Canada. Assimilation is conducted using an ensemble Kalman smoother (EnKS). Model estimates with and without assimilation are compared against independent observational data sets of snow water equivalent (SWE) and runoff. For SWE, modest improvements in mean difference (MD) and root mean squared difference (RMSD) are achieved as a result of the assimilation. No significant differences in temporal correlations of SWE resulted. Runoff statistics of MD remain relatively unchanged while RMSD statistics, in general, are improved in most of the sub-basins. Temporal correlations are degraded within the most upstream sub-basin, but are, in general, improved at the downstream locations, which are more representative of an integrated basin response. GRACE assimilation using an EnKS offers improvements in hydrologic state/flux estimation, though comparisons with observed runoff would be enhanced by the use of river routing and lake storage routines within the prognostic land surface model. Further, GRACE hydrology products would benefit from the inclusion of better constrained models of post-glacial rebound, which significantly affects GRACE estimates of interannual hydrologic variability in the Mackenzie River basin.
Role of Subsurface Physics in the Assimilation of Surface Soil Moisture Observations
Root zone soil moisture controls the land-atmosphere exchange of water and energy and exhibits memory that may be useful for climate prediction at monthly scales. Assimilation of satellite-based surface soil moisture observations into a land surface model is an effective way to estimate large-scale root zone soil moisture. The propagation of surface information into deeper soil layers depends on the model-specific representation of subsurface physics that is used in the assimilation system. In a suite of experiments we assimilate synthetic surface soil moisture observations into four different models (Catchment, Mosaic, Noah and CLM) using the Ensemble Kalman Filter. We demonstrate that identical twin experiments significantly overestimate the information that can be obtained from the assimilation of surface soil moisture observations. The second key result indicates that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Our experiments also suggest that (faced with unknown true subsurface physics) overestimating surface to root zone coupling in the assimilation system provides more robust skill improvements in the root zone compared with underestimating the coupling. When CLM is excluded from the analysis, the skill improvements from using models with different vertical coupling strengths are comparable for different subsurface truths. Finally, the skill improvements through assimilation were found to be sensitive to the regional climate and soil types.