SEARCH · Engineering Papers

Results for “vegetation model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Applying NIR and MIR spectroscopy for C and soil property prediction in northern cold-region ecosystems. Which approach works better?

Here, developing reliable predictions of soil attributes is necessary to understand northern cold-region climate-soil feedback. Calibration models using near-infrared (NIR) and mid-infrared (MIR) spectroscopy were developed to predict eight commonly measured soil properties for 119 soil samples representing a range of vegetation types, parent materials, and soil types spanning >23° of latitude from southeast Alaska to the Canadian high Arctic. In order to obtain a more accurate prediction, this study compared the performance of linear and non-linear calibration techniques, including lasso regression (Lasso), support vector machine (SVM), random forest (RF) and classic partial least squares (PLS) to predict different soil properties of these soils. Comparing the four models, we noticed that their performance was quite similar for MIR overall, while NIR achieved better results with a PLS model for our dataset. PLS coupled with MIR showed a better performance for soil parameters, such as total organic carbon (TOC), total nitrogen (TN), cation exchange capacity (CEC) and clay (R-squared of 0.9, 0.81, 0.80, and 0.84) when compared with NIR (R-squared of 0.85, 0.72, 0.81 and 0.68). However, using either MIR or NIR spectroscopy, PLS predictions for bulk density (BD) and sand content were not accurate. The variable importance analysis based on the PLS model successfully estimated the relative contribution of wavelengths influencing soil property predictions most. Overall, TOC, TN, CEC and clay mineral predictions are closely related to the occurrence of specific spectral bands in the MIR region. For example, wavelengths at 2978 and 1761 cm -1 for TOC and TN, as well as at 3064 cm -1 for CEC, were selected as the most influential predictor variables. We demonstrated that MIR spectroscopy is a powerful tool for more extensive monitoring in soils of the northern cold climate region; however, NIR could be utilized for rapid estimates when the highest accuracy is not essential.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Haiti Agriculture II: Evaluating the Success of Reforestation Practices in Haiti

The Caribbean country of Haiti has an extensive history of deforestation and environmental degradation stemming from French colonization. Over the past 33 years, the Haiti Reforestation Partnership (HRP) and their partners, Comprehensive Development Program (CODEP), have planted approximately 15.52 million trees. However, these efforts lacked scientific guidance to ensure successful forest stand survival. The NASA DEVELOP team partnered with the HRP to create a habitat suitability model (HSM) by using PlanetScope and Sentinel-2 Multispectral Instrument (MSI) imagery. The team also incorporated Landsat 8 and 9 OLI surface temperature, Centre National de L’Information Geo-Spatiale (CNIGS) Airborne Lidar, and ancillary datasets to analyze areas suitable for future reforestation efforts. The habitat model suggested locations with higher forest stand survival based on topography, soil health, climate, and feasibility to access suggested locations. We confirmed the HSM through a cross-analyzation of high enhanced vegetation index (EVI) values. Areas with lower EVI values and higher suitability based on the HSM were suggested for future planting as they lack well-established forest stands but have optimal conditions for growth. Through the creation of the HSM, the team provided the HRP with static maps of high suitability, a 3D printed elevation model, and a guidebook for animators. Additionally, the team provided a structured video highlighting the HRP’s efforts. Effective reforestation and better forest stand survival would help to achieve the goal of securing community food security. This would also serve as a guide to expand planting efforts into other locations and communities.

Kelli Roberts↗

NASA NTRS2024

Applications of ArcticDEM for measuring volcanic dynamics, landslides, retrogressive thaw slumps, snowdrifts, and vegetation heights

Topographical changes are of fundamental interest to a wide range of Arctic science disciplines faced with the need to anticipate, monitor, and respond to the effects of climate change, including geohazard management, glaciology, hydrology, permafrost, and ecology. This study demonstrates several geomorphological, cryospheric, and biophysical applications of ArcticDEM – a large collection of publicly available, time-dependent digital elevation models (DEMs) of the Arctic. Our study illustrates ArcticDEM's applicability across different disciplines and five orders of magnitude of elevation derivatives, including measuring volcanic lava flows, ice cauldrons, post-failure landslides, retrogressive thaw slumps, snowdrifts, and tundra vegetation heights. We quantified surface elevation changes in different geological settings and conditions using the time series of ArcticDEM. Following the 2014–2015 Bárðarbunga eruption in Iceland, ArcticDEM analysis mapped the lava flow field, and revealed the post-eruptive ice flows and ice cauldron dynamics. The total dense-rock equivalent (DRE) volume of lava flows is estimated to be (1431 ± 2) million m 3 . Then, we present the aftermath of a landslide in Kinnikinnick, Alaska, yielding a total landslide volume of (400 ± 8) × 103 m 3 and a total area of 0.025 km 2 . ArcticDEM is further proven useful for studying retrogressive thaw slumps (RTS). The ArcticDEM-mapped RTS profile is validated by ICESat-2 and drone photogrammetry resulting in a standard deviation of 0.5 m. Volume estimates for lake-side and hillslope RTSs range between 40,000 ± 9000 m 3 and 1,160,000 ± 85,000 m 3 , highlighting applicability across a range of RTS magnitudes. A case study for mapping tundra snow demonstrates ArcticDEM's potential for identifying high-accumulation, late-lying snow areas. The approach proves effective in quantifying relative snow accumulation rather than absolute values (standard deviation of 0.25 m, bias of −0.41 m, and a correlation coefficient of 0.69 with snow depth estimated by unmanned aerial systems photogrammetry). Furthermore, ArcticDEM data show its feasibility for estimating tundra vegetation heights with a standard deviation of 0.3 m (no bias) and a correlation up to 0.8 compared to the light detection and ranging (LiDAR). The demonstrated capabilities of ArcticDEM will pave the way for the broad and pan-Arctic use of this new data source for many disciplines, especially when combined with other imagery products. The wide range of signals embedded in ArcticDEM underscores the potential challenges in deciphering signals in regions affected by various geological processes and environmental influences.

Chunli Dai↗

DOE OSTI2026

Daily evapotranspiration changes during heatwaves at 32 NEON sites, 2019-2021

This dataset provides partitioned evapotranspiration (ET, the combined loss of water from soil and plant surfaces) anomalies during heatwave events—soil evaporation (E) and transpiration (T)—for 268 heatwave events across 32 National Ecological Observatory Network (NEON) flux sites in the contiguous United States from 2019–2021. Using an ensemble of four high-frequency turbulence methods (Flux-variance Similarity, Conditional Eddy Covariance [CEC], CEC with Water-Use Efficiency, and Conditional Eddy Accumulation; see Zahn and Bou-Zeid 2024), half-hourly transpiration-to-evapotranspiration (T/ET) ratios were derived from 20 hertz (Hz, cycles per second) eddy covariance measurements of carbon dioxide (CO₂) and water vapor (H₂O) concentrations. The dataset spans six vegetation types including evergreen and deciduous forests, grasslands, cultivated crops, shrublands, and emergent herbaceous wetlands. Data Package Contents: The dataset includes a single CSV (comma-separated values) file containing daily anomalies (deviations from baseline conditions) for transpiration (Delta_T), evaporation (Delta_E), total evapotranspiration (Delta_ET), and T/ET ratio (Delta_T_ET) during each day of identified heatwave events. The file also includes site codes, dates, heatwave event identifiers, and day-of-heatwave indicators. The CSV file can be opened with spreadsheet software (Microsoft Excel, Google Sheets) or programming environments (Python, R, MATLAB). This resource enables researchers to investigate ecosystem-specific responses to thermal extremes, validate land surface model partitioning of ET fluxes, and examine feedbacks between water cycling and surface energy balance during heatwaves. The dataset is particularly valuable for studies linking vegetation hydraulic strategies to climate resilience, as it captures the divergent responses of shallow-rooted versus deep-rooted ecosystems. Potential applications include improving drought early warning systems, informing irrigation management strategies, and advancing our mechanistic understanding of land-atmosphere interactions under extreme heat conditions.

Day of Heatwave↗

NASA NTRS2004

A Simulation Model of Carbon Cycling and Methane Emissions in Amazon Wetlands

An integrative carbon study is investigating the hypothesis that measured fluxes of methane from wetlands in the Amazon region can be predicted accurately using a combination of process modeling of ecosystem carbon cycles and remote sensing of regional floodplain dynamics. A new simulation model has been build using the NASA- CASA concept for predicting methane production and emission fluxes in Amazon river and floodplain ecosystems. Numerous innovations area being made to model Amazon wetland ecosystems, including: (1) prediction of wetland net primary production (NPP) as the source for plant litter decomposition and accumulation of sediment organic matter in two major vegetation classes - flooded forests (varzea or igapo) and floating macrophytes, (2) representation of controls on carbon processing and methane evasion at the diffusive boundary layer, through the lake water column, and in wetland sediments as a function of changes in floodplain water level, (3) inclusion of surface emissions controls on wetland methane fluxes, including variations in daily surface temperature and of hydrostatic pressure linked to water level fluctuations. A model design overview and early simulation results are presented.

Potter, Christopher↗

NASA NTRS2012

Forest Canopy Cover and Height from MISR in Topographically Complex Southwestern US Landscape Assessed with High Quality Reference Data

This study addresses the retrieval of spatially contiguous canopy cover and height estimates in southwestern USforests via inversion of a geometric-optical (GO) model against surface bidirectional reflectance factor (BRF) estimates from the Multi-angle Imaging SpectroRadiometer (MISR). Model inversion can provide such maps if good estimates of the background bidirectional reflectance distribution function (BRDF) are available. The study area is in the Sierra National Forest in the Sierra Nevada of California. Tree number density, mean crown radius, and fractional cover reference estimates were obtained via analysis of QuickBird 0.6 m spatial resolution panchromatic imagery usingthe CANopy Analysis with Panchromatic Imagery (CANAPI) algorithm, while RH50, RH75 and RH100 (50, 75, and 100 energy return) height data were obtained from the NASA Laser Vegetation Imaging Sensor (LVIS), a full waveform light detection and ranging (lidar) instrument. These canopy parameters were used to drive a modified version of the simple GO model (SGM), accurately reproducing patterns ofMISR 672 nm band surface reflectance (mean RMSE 0.011, mean R2 0.82, N 1048). Cover and height maps were obtained through model inversion against MISR 672 nm reflectance estimates on a 250 m grid.The free parameters were tree number density and mean crown radius. RMSE values with respect to reference data for the cover and height retrievals were 0.05 and 6.65 m, respectively, with of 0.54 and 0.49. MISR can thus provide maps of forest cover and height in areas of topographic variation although refinements are required to improve retrieval precision.

Chopping, Mark↗

DOE OSTI2024

The Largely Ignored but Important Role of Non‐Rooted, Submersed Vegetation and Small Waterbodies on Greenhouse Gas Emissions

Abstract Small, shallow waterbodies are abundant globally and are associated with greenhouse gas (GHG) emissions disproportionate to their size. One reason for this is the prevalence of non‐rooted, submersed (NRS) macrophytes that often dominate these systems. Despite these plants' significant role in the carbon ecology of their ecosystems, they are often excluded from climate change models due, in part, to the lack of clarity of their effects. Theus et al. (2023, https://doi.org/10.1029/2023JG007758 ) provides new insights into the role that NRS macrophytes have on GHG emissions, demonstrating that increasing densities of these macrophytes lead to increases in both CO 2 and CH 4 gas fluxes. Here, we highlight the importance of both shallow waterbodies and NRS vegetation within the context of the field of aquatic ecosystem biogeochemistry. We identify a need for the inclusion of these drivers into climate change models based partly on the findings of Theus et al. (2023, https://doi.org/10.1029/2023JG007758 ).

Hassett, Erin↗

NASA NTRS2012

Net Primary Production of Terrestrial Ecosystems from 2000 to 2009

The CASA (Carnegie-Ames-Stanford) ecosystem model has been used to estimate monthly carbon fluxes in terrestrial ecosystems from 2000 to 2009, with global data inputs from NASA's Terra Moderate Resolution Imaging Spectroradiometer (MODIS) vegetation cover mapping. Net primary production (NPP) flux for atmospheric carbon dioxide has varied slightly from year-to-year, but was predicted to have increased over short multi-year periods in the regions of the high-latitude Northern Hemisphere, South Asia, Central Africa, and the western Amazon since the year 2000. These CASA results for global NPP were found to be in contrast to other recently published modeling trends for terrestrial NPP with high sensitivity to regional drying patterns. Nonetheless, periodic declines in regional NPP were predicted by CASA for the southern and western Untied States, the southern Amazon, and southern and eastern Africa. NPP in tropical forest zones was examined in greater detail to discover lower annual production values than previously reported in many global models across the tropical rainforest zones, likely due to the enhanced detection of lower production ecosystems replacing primary rainforest.

Terrestiral ecosystems↗

NASA NTRS2017

The Impact of Anthropogenic Land Use and Land Cover Change on Regional Climate Extremes

Recent research highlights the role of land surface processes in heat waves, droughts, and other extreme events. Here we use an earth system model (ESM) from the Geophysical Fluid Dynamics Laboratory (GFDL) to investigate the regional impacts of historical anthropogenic land useland cover change (LULCC) on combined extremes of temperature and humidity. A bivariate assessment allows us to consider aridity and moist enthalpy extremes, quantities central to human experience of near-surface climate conditions. We show that according to this model, conversion of forests to cropland has contributed to much of the upper central US and central Europe experiencing extreme hot, dry summers every 2-3 years instead of every 10 years. In the tropics, historical patterns of wood harvesting, shifting cultivation and regrowth of secondary vegetation have enhanced near surface moist enthalpy, leading to extensive increases in the occurrence of humid conditions throughout the tropics year round. These critical land use processes and practices are not included in many current generation land models, yet these results identify them as critical factors in the energy and water cycles of the midlatitudes and tropics.

Near-surface↗

NASA NTRS2021

Northern Great Plains Disasters: Using Earth Observations to Enhance Flood Monitoring on Tribal Lands in the Northern Great Plains

In 2019, the Great Plains experienced unprecedented catastrophic flooding. Large flood events are predicted to increase in frequency and severity, posing risks to communities in this region, particularly Tribal Nations. We used data from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR), imagery from the Sentinel-2 MultiSpectral Instrument (MSI), and digital elevation models (DEMs) from the Shuttle Radar Topography Mission (SRTM) within Google Earth Engine to map historical floods in the region beginning in 2014 with particular attention to the Rosebud Sioux Reservation and the tribal lands of other Great Plains Tribal Water Alliance members. This historical mapping used C-SAR for a combined method approach with a Z-score algorithm in addition to an index for flooded short vegetation. We also developed a flood risk map by weighting different flood predictor variables according to flood risk literature. These variables included soil drainage from the Soil Survey Geographic Database (SSURGO); elevation, slope, and Topographic Wetness Index (TWI) derived from digital elevation models; precipitation from Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS); land cover from the National Land Cover Database (NLDC); and Normalized Difference Vegetation Index (NDVI) derived from Landsat 8 Operational Land Imager (OLI). From the flood extent and risk maps, we identified widespread flooding in short vegetation (including cropland) and noted flood susceptibility in regions exhibiting high social vulnerability and low community resilience (FEMA indices). We created an ArcGIS Online StoryMap to share project background, results, and data. Additionally, we provided a written tutorial so partners may replicate the flood mapping for future flood events.

Anna Ballasiotes↗

NASA NTRS2013

Combining Satellite Data and Models to Assess the Impacts of Urbanization on the Continental US Surface Climate

Urbanization is one of the most important and long lasting forms of land transformation. Urbanization affects the surface climate in different ways: (1) by reduction of the vegetation fraction causing subsequent reduction in photosynthesis and plant s water transpiration, (2) by alternation of surface runoff and infiltration and their impacts on soil moisture and the water table, (3) by change in the surface albedo and surface energy partitioning, and (4) by transformation of the surface roughness length and modification of surface fluxes. Land cover and land use change maps including urban areas have been developed and will be used in a suite of land surface models of different complexity to assess the impacts of urbanization on the continental US surface climate. These maps and datasets based on a full range of available satellite data and ground observations will be used to characterize distant-past (pre-urban), recent-past (2001), present (2010), and near future (2020) land cover and land use changes. The main objective of the project is to assess the impacts of these land transformation on past, current and near-future climate and the potential feedbacks from these changes on the atmospheric, hydrologic, biological, and socio-economic properties beyond the immediate metropolitan regions of cities and their near suburbs. The WRF modeling system will be used to explore the nature and the magnitude of the two-way interactions between urban lands and the atmosphere and assess the overall regional dynamic effect of urban expansion on the northeastern US weather and climate

Bounoua, L.↗

NASA NTRS2017

Development of High-Resolution Dynamic Dust Source Function - A Case Study with a Strong Dust Storm in a Regional Model

A high-resolution dynamic dust source has been developed in the NASA Unified-Weather Research and Forecasting (NU-WRF) model to improve the existing coarse static dust source. In the new dust source map, topographic depression is in 1-km resolution and surface bareness is derived using the Normalized Difference Vegetation Index (NDVI) data from Moderate Resolution Imaging Spectroradiometer (MODIS). The new dust source better resolves the complex topographic distribution over the Western United States where its magnitude is higher than the existing, coarser resolution static source. A case study is conducted with an extreme dust storm that occurred in Phoenix, Arizona in 0203 UTC July 6, 2011. The NU-WRF model with the new high-resolution dynamic dust source is able to successfully capture the dust storm, which was not achieved with the old source identification. However the case study also reveals several challenges in reproducing the time evolution of the short-lived, extreme dust storm events.

Dynamic dust source↗

NASA NTRS2014

Impact of Dust Radiative Forcing upon Climate

Dust aerosols perturb the atmospheric radiative flux at both solar and thermal wavelengths, altering the energy and water cycles. The climate adjusts by redistributing energy and moisture, so that local temperature perturbations, for example, depend upon the forcing over the entire extent of the perturbed circulation. Within regions frequently mixed by deep convection, including the deep tropics, dust particles perturb the surface air temperature primarily through radiative forcing at the top of the atmosphere (TOA). Many models predict that dust reduces global precipitation. This reduction is typically attributed to the decrease of surface evaporation in response to dimming of the surface. A counterexample is presented, where greater shortwave absorption by dust increases evaporation and precipitation despite greater dimming of the surface. This is attributed to the dependence of surface evaporation upon TOA forcing through its influence upon surface temperature and humidity. Perturbations by dust to the surface wind speed and vegetation (through precipitation anomalies) feed back upon the dust aerosol concentration. The current uncertainty of radiative forcing attributed to dust and the resulting range of climate perturbations calculated by models remain a useful test of our understanding of the mechanisms relating dust radiative forcing to the climate response.

Feedback upon dust mobilization↗

NASA NTRS2002

Sources of Uncertainty in the Prediction of LAI / fPAR from MODIS

To explicate the sources of uncertainty in the prediction of biophysical variables over space, consider the general equation: where z is a variable with values on some nominal, ordinal, interval or ratio scale; y is a vector of input variables; u is the spatial support of y and z ; x and u are the spatial locations of y and z , respectively; f is a model and B is the vector of the parameters of this model. Any y or z has a value and a spatial extent which is called its support. Viewed in this way, categories of uncertainty are from variable (e.g. measurement), parameter, positional. support and model (e.g. structural) sources. The prediction of Leaf Area Index (LAI) and the fraction of absorbed photosynthetically active radiation (fPAR) are examples of z variables predicted using model(s) as a function of y variables and spatially constant parameters. The MOD15 algorithm is an example of f, called f(sub 1), with parameters including those defined by one of six biome types and solar and view angles. The Leaf Canopy Model (LCM)2, a nested model that combines leaf radiative transfer with a full canopy reflectance model through the phase function, is a simpler though similar radiative transfer approach to f(sub 1). In a previous study, MOD15 and LCM2 gave similar results for the broadleaf forest biome. Differences between these two models can be used to consider the structural uncertainty in prediction results. In an effort to quantify each of the five sources of uncertainty and rank their relative importance for the LAI/fPAR prediction problem, we used recent data for an EOS Core Validation Site in the broadleaf biome with coincident surface reflectance, vegetation index, fPAR and LAI products from the Moderate Resolution Imaging Spectrometer (MODIS). Uncertainty due to support on the input reflectance variable was characterized using Landsat ETM+ data. Input uncertainties were propagated through the LCM2 model and compared with published uncertainties from the MOD15 algorithm.

Dungan, Jennifer L.↗

DOE OSTI2023

Machine learning models inaccurately predict current and future high-latitude C balances

The high-latitude carbon (C) cycle is a key feedback to the global climate system, yet because of system complexity and data limitations, there is currently disagreement over whether the region is a source or sink of C. Recent advances in big data analytics and computing power have popularized the use of machine learning (ML) algorithms to upscale site measurements of ecosystem processes, and in some cases forecast the response of these processes to climate change. Due to data limitations, however, ML model predictions of these processes are almost never validated with independent datasets. To better understand and characterize the limitations of these methods, we develop an approach to independently evaluate ML upscaling and forecasting. We mimic data-driven upscaling and forecasting efforts by applying ML algorithms to different subsets of regional process-model simulation gridcells, and then test ML performance using the remaining gridcells. In this study, we simulate C fluxes and environmental data across Alaska using ecosys, a process-rich terrestrial ecosystem model, and then apply boosted regression tree ML algorithms to training data configurations that mirror and expand upon existing AmeriFLUX eddy-covariance data availability. We first show that a ML model trained using ecosys outputs from currently-available Alaska AmeriFLUX sites incorrectly predicts that Alaska is presently a modeled net C source. Increased spatial coverage of the training dataset improves ML predictions, halving the bias when 240 modeled sites are used instead of 15. However, even this more accurate ML model incorrectly predicts Alaska C fluxes under 21st century climate change because of changes in atmospheric CO 2 , litter inputs, and vegetation composition that have impacts on C fluxes which cannot be inferred from the training data. Our results provide key insights to future C flux upscaling efforts and expose the potential for inaccurate ML upscaling and forecasting of high-latitude C cycle dynamics.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1985

The role of isoprene oxidation in the tropospheric ozone budget in the tropics

A comprehensive chemical mechanism for the oxidation of isoprene (a hydrocarbon, C5H8 emitted primarily by vegetation) by OH and O3 in the troposphere was developed and incorporated into a one-dimensional steady-state photochemical model of the troposphere. Flux boundary conditions for NOx (NO + NO2), HNO3, O3, and CO were used to investigate the changes produced in the tropospheric concentrations and integrated column of ozone from including isoprene chemistry in the model. Two calculations were performed at 15 deg N latitude for annual conditions using identical flux boundary conditions for NOx, HNO3, O3, and CO; in one calculation, the chemistry describing isoprene oxidation was included while in the other it was not. Both sets of calculations included reactions describing the chemistry of anthropogenic nonmethane hydrocarbons. The calculations showed decreases in concentrations of ozone throughout the troposphere when isoprene chemistry was included. Concentrations of NOx and HNO3 increased in the lower troposphere and decreased in the upper troposphere while concentrations of CO and PAN increased throughout the troposphere when isoprene chemistry was included. Implications of this study to the budgets of these species in the tropics is discussed.

Brewer, D. A.↗

NASA NTRS1987

Spectral characteristics and the extent of paleosols of the Palouse formation

Three spectral models defining the spatial distribution of soil areas by levels of amorphous iron, organic carbon, and the ratio of amorphous iron to organic carbon were developed and field verification studies were conducted. The models used particular Thematic Mapper band ratios selected by statistical correlation with soil chemical data. The ability of the models to indicate erosion severity and to differentiate between iron enriched and carbonate paleosols is discussed. In addition, the effect of vegetation cover on paleosols is addressed.

Frazier, B. E.↗

NASA NTRS1987

Concerning the relationship between evapotranspiration and soil moisture

The relationship between the evapotranspiration and soil moisture during the drying, supply-limited phase is studied. A second scaling parameter, based on the evapotranspirational supply and demand concept of Federer (1982), is defined; the parameter, referred to as the threshold evapotranspiration, occurs in vegetation-covered surfaces just before leaf stomata close and when surface tension restricts moisture release from bare soil pores. A simple model for evapotranspiration is proposed. The effects of natural soil heterogeneities on evapotranspiration computed from the model are investigated. It is observed that the natural variability in soil moisture, caused by the heterogeneities, alters the relationship between regional evapotranspiration and the area average soil moisture.

Wetzel, Peter J.↗