Mapping Drought with Satellite Remote Sensing and Land Surface Modeling
No abstract available
SEARCH · Engineering Papers
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.
No abstract available
Not Available
Observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission have a coarse resolution in time (monthly) and space (roughly 150,000 km(sup 2) at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This work proposes a variant of existing ensemble-based GRACE-TWS data assimilation schemes. The new algorithm differs in how the analysis increments are computed and applied. Existing schemes correlate the uncertainty in the modeled monthly TWS estimates with errors in the soil moisture profile state variables at a single instant in the month and then apply the increment either at the end of the month or gradually throughout the month. The proposed new scheme first computes increments for each day of the month and then applies the average of those increments at the beginning of the month. The new scheme therefore better reflects submonthly variations in TWS errors. The new and existing schemes are investigated here using gridded GRACE-TWS observations. The assimilation results are validated at the monthly time scale, using in situ measurements of groundwater depth and soil moisture across the U.S. The new assimilation scheme yields improved (although not in a statistically significant sense) skill metrics for groundwater compared to the open-loop (no assimilation) simulations and compared to the existing assimilation schemes. A smaller impact is seen for surface and root-zone soil moisture, which have a shorter memory and receive smaller increments from TWS assimilation than groundwater. These results motivate future efforts to combine GRACE-TWS observations with observations that are more sensitive to surface soil moisture, such as L-band brightness temperature observations from Soil Moisture Ocean Salinity (SMOS) or Soil Moisture Active Passive (SMAP). Finally, we demonstrate that the scaling parameters that are applied to the GRACE observations prior to assimilation should be consistent with the land surface model that is used within the assimilation system.
Hydrological models have been calibrated and validated using catchment streamflows. However, using a point measurement does not guarantee correct spatial distribution of model computed heat fluxes, soil moisture and surface temperatures. With the advent of satellites in the late 70s, surface temperature is being measured two to four times a day from various satellite sensors and different platforms. The purpose of this paper is to demonstrate use of satellite surface temperature in (a) validation of model computed surface temperatures and (b) assimilation of satellite surface temperatures into a hydrological model in order to improve the prediction accuracy of soil moistures and heat fluxes. The assimilation is carried out by comparing the satellite and the model produced surface temperatures and setting the "true"temperature midway between the two values. Based on this "true" surface temperature, the physical relationships of water and energy balance are used to reset the other variables. This is a case of nudging the water and energy balance variables so that they are consistent with each other and the true" surface temperature. The potential of this assimilation scheme is demonstrated in the form of various experiments that highlight the various aspects. This study is carried over the Red-Arkansas basin in the southern United States (a 5 deg X 10 deg area) over a time period of a year (August 1987 - July 1988). The land surface hydrological model is run on an hourly time step. The results show that satellite surface temperature assimilation improves the accuracy of the computed surface soil moisture remarkably.
A physics-based risk model was developed to assess the risk associated with thermal protection system failures during the entry, descent and landing phase of a manned spacecraft mission. In the model, entry trajectories were computed using a three-degree-of-freedom trajectory tool, the aerothermodynamic heating environment was computed using an engineering-level computational tool and the thermal response of the TPS material was modeled using a one-dimensional thermal response tool. The model was capable of modeling the effect of micrometeoroid and orbital debris impact damage on the TPS thermal response. A Monte Carlo analysis was used to determine the effects of uncertainties in the vehicle state at Entry Interface, aerothermodynamic heating and material properties on the performance of the TPS design. The failure criterion was set as a temperature limit at the bondline between the TPS and the underlying structure. Both direct computation and response surface approaches were used to compute the risk. The model was applied to a generic manned space capsule design. The effect of material property uncertainty and MMOD damage on risk of failure were analyzed. A comparison of the direct computation and response surface approach was undertaken.
Abstract The diffuse radiation fertilization effect—the increase in plant productivity in the presence of higher diffuse radiation ( K ↓,d )—is an important yet understudied aspect of atmosphere‐biosphere interactions and can modify the terrestrial carbon, energy, and water budgets. The K ↓,d fertilization effect links the carbon cycle with clouds and aerosols, all of which are large sources of uncertainties for our current understanding of the Earth system and for future climate projections. Here we establish to what extent observational and modeling uncertainty in sunlight's diffuse fraction ( k d ) affects simulated gross primary productivity (GPP) and terrestrial evapotranspiration ( λE ). We find only 48 eddy covariance sites with simultaneous sufficient measurements of K ↓,d with none in the tropical climate zone, making it difficult to constrain this mechanism globally using observations. Using a land modeling framework based on the latest version of the Community Land Model, we find that global GPP ranges from 114 Pg C year −1 when using k d forcing from the Modern‐Era Retrospective analysis for Research and Applications, version 2 reanalysis to a ∼7% higher value of 122 Pg C year −1 when using the Clouds and the Earth's Radiant Energy System satellite product, with especially strong differences apparent over the tropical region (mean increase ∼9%). The differences in λE , although smaller (−0.4%) due to competing changes in shaded and sunlit leaf transpiration, can be greater than regional impacts of individual forcing agents like aerosols. Our results demonstrate the importance of comprehensively and systematically validating the simulated k d by atmosphere modules as well as the response to differences in k d within land modules across Earth System Models.
We present a generic flux limiter to account for mass limitations from an arbitrary number of substrates in a biogeochemical reaction network. The flux limiter is based on the observation that substrate (e.g., nitrogen, phosphorus) limitation in biogeochemical models can be represented as to ensure mass conservative and non-negative numerical solutions to the governing ordinary differential equations. Application of the flux limiter includes two steps: (1) formulation of the biogeochemical processes with a matrix of stoichiometric coefficients and (2) application of Liebig's law of the minimum using the dynamic stoichiometric relationship of the reactants. This approach contrasts with the ad hoc down-regulation approaches that are implemented in many existing models (such as CLM4.5 and the ACME (Accelerated Climate Modeling for Energy) Land Model (ALM)) of carbon and nutrient interactions, which are error prone when adding new processes, even for experienced modelers. Through an example implementation with a CENTURY-like decomposition model that includes carbon, nitrogen, and phosphorus, we show that our approach (1) produced almost identical results to that from the ad hoc down-regulation approaches under non-limiting nutrient conditions, (2) properly resolved the negative solutions under substrate-limited conditions where the simple clipping approach failed, (3) successfully avoided the potential conceptual ambiguities that are implied by those ad hoc down-regulation approaches. We expect our approach will make future biogeochemical models easier to improve and more robust.
Vegetation plays a fundamental role not only in theenergy and carbon cycles but also in the global water balanceby controlling surface evapotranspiration (ET). Thus, accu-rately estimating vegetation-related variables has the poten-tial to improve our understanding and estimation of the dy-namic interactions between the water, energy, and carbon cy-cles. This study aims to assess the extent to which a land sur-face model (LSM) can be optimized through the assimilationof leaf area index (LAI) observations at the global scale. Twoobserving system simulation experiments (OSSEs) are per-formed to evaluate the efficiency of assimilating LAI into anLSM through an ensemble Kalman filter (EnKF) to estimateLAI, ET, canopy-interception evaporation (CIE), canopy wa-ter storage (CWS), surface soil moisture (SSM), and terres-trial water storage (TWS). Results show that the LAI dataassimilation framework not only effectively reduces errorsin LAI model simulations but also improves all the modeledwater flux and storage variables considered in this study (ET,CIE, CWS, SSM, and TWS), even when the forcing pre-cipitation is strongly positively biased (extremely wet con-ditions). However, it tends to worsen some of the modeledwater-related variables (SSM and TWS) when the forcingprecipitation is affected by a dry bias. This is attributed tothe fact that the amount of water in the LSM is conservative,and the LAI assimilation introduces more vegetation, whichrequires more water than what is available within the soil
The Lower Mekong Basin (LMB) is an economically and ecologically important region that experiences hydrologic hazards such as floods and droughts, which can directly affect human well-being and limit economic growth and development. To effectively develop long-term plans for addressing hydrologic hazards, the regional hydrological response to climate variability and land cover change needs to be evaluated. This research aims to investigate how climate variability, specifically variations in the precipitation regime, and land cover change will affect hydrologic parameters both spatially and temporally within the LMB. The research goal is achieved by (1) modeling land cover change for a baseline land cover change scenario as well as changes in land cover with increases in forest or agriculture and (2) using projected climate variables and modeled land cover data as inputs into the Variable Infiltration Capacity (VIC) hydrologic model to simulate the changes to the hydrologic system. The VIC model outputs were analyzed against historic values to understand the relative contribution of climate variability and land cover to change, where these changes occur, and to what degree these changes affect the hydrology. This study found that the LMB hydrologic system is more sensitive to climate variability than land cover change. On average, climate variability was found to increase discharge and evapotranspiration (ET) while decreasing water storage. The change in land cover show that increasing forest area will slightly decrease discharge and increase ET while increasing agriculture area increases discharge and decreases ET. These findings will help the LMB by supporting individual country policy to plan for future hydrologic changes as well as policy for the basin as a whole.
NASA's GRACE mission has the potential to be extremely valuable for water resources applications and global water cycle research. What makes GRACE unique among Earth Science satellite systems is that it is able to monitor variations in water stored in all forms, from snow and surface water to soil moisture to groundwater in the deepest aquifers. However, the space and time resolutions of GRACE observations are coarse. GRACE typically resolves water storage changes over regions the size of Nebraska on a monthly basis, while city-scale, daily observations would be more useful for water management, agriculture, and weather prediction. High resolution numerical (computer) hydrology models have been developed, which predict the fates of water and energy after they strike the land surface as precipitation and sunlight. These are similar to weather and climate forecast models, which simulate atmospheric processes. We integrated the GRACE observations into a hydrology model using an advanced technique called data assimilation. The results were new estimates of groundwater, soil moisture, and snow variations, which combined the veracity of GRACE with the high resolution of the model. We tested the technique over the Mississippi River basin, but it will be even more valuable in parts of the world which lack reliable data on water availability.
Abstract Nutrient limitation is widespread in terrestrial ecosystems. Accordingly, representations of nitrogen (N) limitation in land models typically dampen rates of terrestrial carbon (C) accrual, compared with C‐only simulations. These previous findings, however, rely on soil biogeochemical models that implicitly represent microbial activity and physiology. Here we present results from a biogeochemical model testbed that allows us to investigate how an explicit versus implicit representation of soil microbial activity, as represented in the MIcrobial‐MIneral Carbon Stabilization (MIMICS) and Carnegie‐Ames‐Stanford Approach (CASA) soil biogeochemical models, respectively, influence plant productivity, and terrestrial C and N fluxes at initialization and over the historical period. When forced with common boundary conditions, larger soil C pools simulated by the MIMICS model reflect longer inferred soil organic matter (SOM) turnover times than those simulated by CASA. At steady state, terrestrial ecosystems experience greater N limitation when using the MIMICS‐CN model, which also increases the inferred SOM turnover time. Over the historical period, however, warming‐induced acceleration of SOM decomposition over high latitude ecosystems increases rates of N mineralization in MIMICS‐CN. This reduces N limitation and results in faster rates of vegetation C accrual. Moreover, as SOM stoichiometry is an emergent property of MIMICS‐CN, we highlight opportunities to deepen understanding of sources of persistent SOM and explore its potential sensitivity to environmental change. Our findings underscore the need to improve understanding and representation of plant and microbial resource allocation and competition in land models that represent coupled biogeochemical cycles under global change scenarios.
Tremendous progress has been achieved in the development of land models and their inclusion in Earth system models (ESMs). However, we still have very limited knowledge on the performance skills of these land models. This chapter introduces benchmark analysis, which is a procedure to measure performance of models against a set of defined standards. The benchmark analysis includes: (1) defining targeted aspects of model performance to be evaluated; (2) testing model performance in comparison with a set of benchmarks; (3) measuring model performance skill through quantitative metrics; and (4) evaluating model performance and offering suggestions for future model improvement.
Radiocarbon ( 14 C) is a powerful tracer of the global carbon cycle that is commonly used to assess carbon cycling rates in various Earth system reservoirs and as a benchmark to assess model performance. Therefore, it has been recommended that Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 report predicted radiocarbon values for relevant carbon pools. However, a detailed representation of radiocarbon dynamics may be an impractical burden on model developers. Here, we present an alternative approach to compute radiocarbon values from the numerical output of an ESM that does not explicitly represent these dynamics. The approach requires computed 12 C stocks and fluxes among all carbon pools for a particular simulation of the model. From this output, a time-dependent linear compartmental system is computed with its respective state-transition matrix. Using transient atmospheric 14 C values as inputs, the state-transition matrix is then applied to compute radiocarbon values for each pool, the average value for the entire system, and component fluxes. We demonstrate the approach with ELMv1-ECA, the land component of an ESM model that explicitly represents 12 C, and 14 C in 7 soil pools and 10 vertical layers. Results from our proposed method are highly accurate (relative error <0.01%) compared with the ELMv1-ECA 12 C and 14 C predictions, demonstrating the potential to use this approach in CMIP6 and other model simulations that do not explicitly represent 14 C.
Large-scale land models assume unidirectional land-to-river hydrological interactions, without considering feedbacks between channels and land. Using a tested, physically based model with explicit multiway interactions between overland, channel, wetland, and groundwater flows, we assessed how the representation and properties of channels influence simulated land surface hydrologic, biogeochemical, and ecosystem dynamics. A zone near the channels where various fluxes and states are significantly influenced by the channels, referred to as the fan of influence (FoI) of channels, has been identified. We elucidated two mechanisms inducing the model-derived FoI: the base flow mechanism, in which incised, gaining streams lower the water table and induce more base flow, and the relatively more efficient conveyance of the channel network compared to overland flow. We systematically varied drainage density and grid resolution to quantify the size of the FoI, which is found to span a large fraction of the watershed (25-50%) for hydrologic variables including depth to water table and recharge, etc. The FoI is more pronounced with low-resolution simulations but remains noticeable in hyperresolution (25 m) subbasin simulations. The FoI and the channel influence on basin-average fluxes are also similar in simulations with alternative parameter sets. We found that high-order, entrenched streams cause larger FoI. In addition, removing the simulated channels has disproportionally large influence on modeled wetland areas and inundation duration, which has implications for coupled biogeochemical or ecological modeling. Our results suggest that explicit channel representation provides important feedbacks to land surface dynamics which should be considered in meso or large-scale simulations. Since grid refinement incurs prohibitive computational cost, subgrid channel parameterization has advantages in efficiency over grid-based representations that do not distinguish between overland flow and channel flow.
This Modeling Archive is in support of an NGEE Arctic publication "Rapidly changing high-latitude seasonality: Implications for the 21st century carbon cycle in Alaska" DOI: 10.1088/1748-9326/ac4362.The dataset contains "ecosys" model outputs reported in Shirley et al. (2022) that explored how climate warming will shift carbon cycle seasonality in Alaska. The study used analysis of seasonal dynamics to support the prediction that high-latitude ecosystems are carbon sinks and will continue to accumulate carbon throughout the century. The study predicts surprisingly large increases in spring net carbon uptake that result in larger net carbon uptake in spring than in summer by year 2100. Included as *.nc files are modeled daily net ecosystem exchange (NEE), net primary productivity (NPP), heterotrophic respiration (Rh), autotrophic respiration (Ra), gross primary productivity (GPP), and plant nitrogen uptake. The output is provided (1) at 25 km resolution across Alaska; (2) for years 2010 through 2100. Scripts to generate publication figures included as text files (*.m) from MATLAB. A user guide included as *.pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
Studies have shown that irrigated agriculture has substantial effects on local to regional hydrologic and atmospheric processes. These effects include increased evapotranspiration and near surface humidity, decreased air and ground temperature, augmented downwind precipitation, enhanced recharge, depletion of source water stocks, and overall changes to the regional water budget. However, developers of global offline and coupled modeling systems have been slow to incorporate irrigation. This presentation will describe progress towards simulating irrigation worldwide in the Land Information System (LIS) and Global Land Data Assimilation System (GLDAS). We have developed and refined the ability to simulate sprinkler irrigation, furrow irrigation, paddy irrigation, and drip irrigation. We have reconciled multiple global and national irrigation, land cover, and crop type databases in an effort to (1) distinguish irrigated and rainfed agriculture, (2) simulate the proper distribution of crops (as sub-grid tiles) in each model grid pixel, and (3) apply the most likely irrigation method in each tile based on the crop type and country or state level information on the relative usage of sprinkler, flood, and drip irrigation. Using agricultural water withdrawal data, as available, for calibration and evaluation, we show that the implemented approach is promising but that much fine tuning remains to be done. In particular, irrigation model parameters that produce good results in one location often produce less satisfying results in others.
Here, we present an improved procedure of generating initial conditions (ICs) for climate model hindcast experiments with specified sea surface temperature and sea ice. The motivation is to minimize errors in the ICs and lead to a better evaluation of atmospheric parameterizations' performance in the hindcast mode. We apply state variables (horizontal velocities, temperature and specific humidity) from the operational analysis/reanalysis for the atmospheric initial states. Without a data assimilation system, we apply a two-step process to obtain other necessary variables to initialize both the atmospheric (e.g., aerosols and clouds) and land models (e.g., soil moisture). First, we nudge only the model horizontal velocities towards operational analysis/reanalysis values, given a 6-hour relaxation time scale, to obtain all necessary variables. Compared to the original strategy in which horizontal velocities, temperature and specific humidity are nudged, the revised approach produces a better representation of initial aerosols and cloud fields which are more consistent and closer to observations and model's preferred climatology. Second, we obtain land ICs from an offline land model simulation forced with observed precipitation, winds, and surface fluxes. This approach produces more realistic soil moisture in the land ICs. With this refined procedure, the simulated precipitation, clouds, radiation, and surface air temperature over land are improved in the Day 2 mean hindcasts. Following this procedure, we propose a “Core” integration suite which provides an easily repeatable test allowing model developers to rapidly assess the impacts of various parameterization changes on the fidelity of modelled cloud-associated processes relative to observations.