Search NASA⌕ Search

SEARCH · Search NASA

Results for “site balance 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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Stability of the Martian atmosphere: Is heterogeneous catalysis essential?

A comprehensive homogeneous gas phase photochemical model is developed to study the problem of stability of the Martian atmosphere. The one-dimensional model extends from the ground up to 220 km, passing through the homopause at 125 km. The model thus couples the lower (neutral) atmosphere to the ionosphere above which provides significant downward flux of carbon monoxide and oxygen atoms. It is concluded on the basic of currently accepted values for globally and seasonally averaged water vapor abundance, dust opacity and the middle atmospheric eddy mixing coefficient, as well as the relevant laboratory data (particularly the temperature dependence of CO2 absorption cross section and the rate constant for CO+OH reaction), that the rate of re-formation of carbon dioxide exceeds its photolytic destruction rate by about 40%. Furthermore, it is found that this result is virtually independent of the choice of eddy mixing coefficient, unless its value in the middle atmosphere exceeds 10(exp 8) sq cm/sec or is far smaller than 10(exp 5)sq cm/sec, or the dust opacity, unless it exceeds unity, or the water vapor mixing ratio at the surface, unless it is far smaller (less than or = 1 ppm) or far greater (greater than or = 500 ppm) than the average value (approximately 150 ppm). Since none of these extremes represent globally and seasonally averaged conditions on Mars, we propose that the present model requires existence of a mechanism to throttle down the recycling rate of carbon dioxide on Mars. Therefore, it is suggested that a heterogeneous process which provides a sink to the species that participate in the recycling of CO2, i.e., H2O, H2O2, OH, CO or O, in particular, may be necessary to bring about the balance between the CO2 recycling rate and its photolytic destruction rate. Aerosols of dust or ice (pure or doped water or carbon dioxide ice present in the atmosphere of Mars) can provide the appropriate adsorption sites for the above heterogeneous process. Despite our conclusion that some heterogeneous process may be needed, it is important to recognize that one-dimensional models can only provide first-order results which, most likely, represent globally and seasonally averaged conditions. However, it is only after actual temporal, latitudinal and longitudinal variations of relevant atmospheric parameters are included in the model that one can determine fuly whether the problem of atmospheric stabiltiy still continues to persist and whether some heterogeneous process is required to correct it.

Atreya, Sushil K.↗
DOE OSTI2025

Data and code for: Combining eddy covariance towers, field measurements, and the MEMS 2 ecosystem model improves confidence in the climate impacts of bioenergy with carbon capture and storage

BECCS demands extensive land use change that may create positive or negative radiative forcing impacts upstream of the BECCS facility through changes to in situ greenhouse gas fluxes and land surface albedo. When quantifying these upstream climate impacts, even at a single site, different methods can give different estimates. Here we show how three common methods for estimating the net ecosystem carbon balance of bioenergy crops established on former grassland or former cropland can differ in their central estimates and uncertainty.

BECCS↗
DOE OSTI2026

Deposition Height Prediction in Directed Energy Deposition

Using 316L stainless steel as a model material, reduced-order models are developed to predict capture efficiency, deposition height, and site-specific hardness in directed energy deposition. Capture efficiency is predicted over a 15 to 55 pct range using a dimensionless number derived from processing conditions and thermophysical properties. Deposition height is predicted over a 0.3 to 1.3 mm range without in situ sensing or prior training data, using two models based on the same mass and energy-balance principles. Predictions are compared with machine learning approaches. A quantitative relationship links deposition height, primary dendrite arm spacing (PDAS), and hardness: heights of 0.3 to 1.1 mm correspond to PDAS values of 2.7 to 5.1 µm and Vickers hardness (HV) of 160 to 219. Thinner layers cool more rapidly, producing finer microstructures and higher hardness. Samples fabricated with in situ variations in deposition height exhibited up to 55 HV differences between thick and thin regions, demonstrating that local control of deposition height enables predictive, site-specific hardness within a single build. These results establish deposition height prediction as a pathway for a priori process design and property control in directed energy deposition for 316L stainless steel.

Kunkel, William [Univ. of Wisconsin, Madison, WI (↗
DOE OSTI2020

ORBIT: Offshore Renewables Balance-of-System and Installation Tool

This report describes the Offshore Renewables Balance-of-system Installation Tool ("ORBIT"), a new model developed by the National Renewable Energy Laboratory (NREL) to evaluate the balance-of-system (BOS) costs associated with offshore wind projects. In the context of wind energy projects BOS costs encompass all expenses required to construct a project other than the capital expenditures (CapEx) of the turbines and towers, including the procurement costs for all other components (such as substructures, cables, and electrical infrastructure), offshore and land-based construction costs, port costs, site surveying fees, permitting fees, and leasing fees are all categorized as BOS costs. BOS costs significantly contribute to the levelized cost of energy (LCOE) for offshore wind, typically comprising over 50\% of the CapEx for a fixed-bottom offshore wind project and 60\% for a floating project. In addition, technology solutions and installation methods vary drastically between projects as they are impacted by factors such as vessel availability, geographic considerations and site geotechnical conditions. These effects require a model with appropriate fidelity to understand how these costs scale as turbine rating increases and the offshore wind supply chain, particularly offshore construction vessels, is expanded. For offshore wind, cost savings attributed to increased turbine rating are primarily realized through BOS procurement and project installation as fewer substructures and less cable are required. BOS costs represent both a modeling challenge, as well as an opportunity, for project developers to optimize solutions to reduce costs. It is critical to understand how these costs are affected by novel technologies, innovative installation processes, and operational constraints in order identify meaningful cost reductions for offshore wind energy.

17 WIND ENERGY↗
NASA NTRS2022

Evaluation of remote sensing-based evapotranspiration products at low-latitude eddy covariance sites

Remote sensing-based evapotranspiration (ET) products have been evaluated primarily using data from northern middle latitudes; therefore, little is known about their performance at low latitudes. To address this bias, an evaluation dataset was compiled using eddy covariance data from 40 sites between latitudes 30° S and 30° N. The flux data were obtained from the emerging network in Mexico (MexFlux) and from openly available databases of FLUXNET, AsiaFlux, and OzFlux. This unique reference dataset was then used to evaluate remote sensing-based ET products in environments that have been underrepresented in earlier studies. The evaluated products were: MODIS ET (MOD16, both the discontinued collection 5 (C5) and the latest collection (C6)), Global Land Evaporation Amsterdam Model (GLEAM) ET, and Atmosphere-Land Exchange Inverse (ALEXI) ET. Products were compared with unadjusted fluxes (ETorig) and with fluxes corrected for the lack of energy balance closure (ETebc). Three common statistical metrics were used: coefficient of determination (R2), root mean square error (RMSE), and percent bias (PBIAS). The effect of a vegetation mismatch between pixel and site on product evaluation results was investigated by examining the relationship between the statistical metrics and product-specific vegetation match indexes. Evaluation results of this study and those published in the literature were used to examine the performance of the products across latitudes. Differences between the MOD16 collection 5 and 6 datasets were generally smaller than differences with the other products. Performance and ranking of the evaluated products depended on whether ETorig or ETebc was used. When using ETorig, GLEAM generally had the highest R2, smallest PBIAS, and best RMSE values across the studied land cover types and climate zones. Neither MOD16 nor ALEXI performed consistently better than the other. When using ETebc, none of the products stood out in terms of both low bias and strong correlations. The use of ETebc instead of ETorig affected the biases more than the correlations. The product evaluation results showed no significant relationship with the degree of match between the vegetation at the pixel and site scale. The latitudinal comparison showed tendencies of lower R2 (all products) but better PBIAS and normalized RMSE values (MOD16 and GLEAM) for forests at low latitudes than for forests at northern middle latitudes. For non-forest vegetation, the products showed no clear latitudinal differences in performance.

Diego Salazar-Martínez↗
DOE OSTI2020

Conceptual Model Update of Mercury Sources and Flux at Y-12 and Upper East Fork Poplar Creek, Oak Ridge,Tennessee

The interactions between the key processes that control the fate and transport of mercury (Hg) in the sediments, groundwater, and stream water within the Y-12 National Security Complex (Y-12) in Oak Ridge, Tennessee are complex, and many aspects are not well understood. In contaminated environments such as Y-12, conceptual models are frequently developed to aid in visualizing and understanding the dynamic nature of the hydrologic, geochemical, and physical environment. These conceptual models integrate data in an internally consistent manner to understand processes that control the fate and transport of contaminants. Over the past few decades of environmental investigation at Y-12, a number of conceptual models have been developed to identify and define various technical processes at various scales. The main purpose of the activities described in this report is to update the 2011 conceptual model for Upper East Fork Poplar Creek using the most recent Hg concentration and flux data and informed scientific interpretation. Emphasis was placed on improving the previous model by implementing a more detailed spatial and temporal approach to visualize transport pathways in the watershed and trends in flux and concentration over time, and to compare baseflow and stormflow system dynamics. Detailed descriptions of historical and current Hg sources and transport pathways are also provided for the east and west ends of the facility. This conceptual model will allow the US Department of Energy to evaluate past and present remedial activities and provide a strong technical basis for prioritizing and optimizing remedial responses in a cost-effective and efficient manner. Twenty years ago, Hg flux from Outfall (OF) 200 at the headwaters of East Fork Poplar Creek represented approximately 20% of the overall flux leaving Y-12. By the time of the 2011 conceptual model report, that percentage had increased to 70%–80%. Flux estimates compiled for this report using data from 2009 to 2018 show that OF200 accounted for roughly 60% of the flux leaving Y-12. The relative role of the four storm drain conduits to OF200 appeared similar to 2011, with OF163 being the greatest contributor to downstream flux. An increase in Hg concentrations and flux occurred throughout and downstream of west end storm drains in 2011 in response to a storm drain cleanout. Leading up to the 2011 storm drain cleanout, annual baseflow flux was higher at OF200A6 than at Station 17. From 2011 to 2018, the trend reversed with concentrations and flux being higher at Station 17 than at OF200A6. A spike in concentrations occurred in July 2018, potentially due to a one-time influx of Hg that occurred during COLEX decontamination and decommissioning activities at Y-12. This report includes recommendations for improving our conceptual understanding of Hg sources, transport pathways, and flux at Y-12. Important recommendations include further evaluation of transport pathways, concurrent measurement of Hg concentration and flux under baseflow and stormflow conditions at multiple sites, and research to understand the connections among periphyton and methylmercury concentration, bioavailability, and bioaccumulation. Ongoing facility decontamination and decommissioning efforts, site characterization, remedial actions, and research are facilitating the collection of new Hg data at Y-12. The conceptual model update provided herein takes advantage of the extensive recent Hg sampling to provide an updated assessment of Hg mass balance and support a refined understanding of Hg behavior at or near Y-12. This assessment should assist in future environmental management decisions and in mitigating the impacts of Hg on the surrounding environment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2010

A Microwave Radiance Assimilation Study for a Tundra Snowpack

Recent studies have begun exploring the assimilation of microwave radiances for the modeling and retrieval of snow properties. At a point scale, and for short durations (i week), radiance assimilation (RA) results are encouraging. However, in order to determine how practical RA might be for snow retrievals when applied over longer durations, larger spatial scales, and/or different snow types, we must expand the scope of the tests. In this paper we use coincident microwave radiance measurements and station data from a tundra site on the North Slope of Alaska. The field data are from the 3rd Radio-brightness Energy Balance Experiment (REBEX-3) carried out in 1994-95 by the University of Michigan. This dataset will provide a test of RA over months instead of one week, and for a very different type of snow than previous snow RA studies. We will address the following questions: flow well can a snowpack physical model (SM), forced with local weather, match measured conditions for a tundra snowpack?; How well can a microwave emission model, driven by the snowpack model, match measured microwave brightnesses for a tundra snowpack?; How well does RA increase or decrease the fidelity of estimates of snow depth and temperatures for a tundra snowpack?

Kim, Edward↗
NASA NTRS2003

The Broken Belt: Meteorite Concentrations on Stranded Ice

Since the first Antarctic meteorite concentrations were discovered more than 25 years ago, many theories regarding the role of iceflow in the production of meteorite concentrations have been put forward, and most agree on the basic principles. These models suggest that as the East Antarctic icesheet flows toward the margins of the continent, meteorites randomly located within the volume of ice are transported toward the icesheet margin. Where mountains or subsurface obstructions block glacial flow, diversion of ice around or over an obstruction reduces horizontal ice movement rates adjacent to the barriers and creates a vertical (upward) component of movement. If local mechanisms for ice loss (ablation) exist at such sites, an equilibrium surface will develop according to the balance between ice supply and loss, and the cargo of meteorites is exhumed on a blue ice surface. The result is a conceptual conveyor belt bringing meteorite-bearing volumes of ice from the interior of the continent to stagnant or slowmoving surfaces where ice is then lost and a precious cargo is left as a lag deposit. Cassidy et al. provides an excellent overview of how this model has been adapted to several Antarctic stranding surfaces.

Harvey, R. P.↗
DOE OSTI2021

SU(4) Chiral Spin Liquid, Exciton Supersolid, and Electric Detection in Moiré Bilayers

We propose a moiré bilayer as a platform where exotic quantum phases can be stabilized and electrically detected. Moiré bilayers consist of two separate moiré superlattice layers coupled through the interlayer Coulomb repulsion. In the small distance limit, an SU(4) spin can be formed by combining layer pseudospin and the real spin. As a concrete example, we study an SU(4) spin model on triangular lattice in the fundamental representation. By tuning a three-site ring exchange term K ~ ( t 3 /U 2 ), we find the SU(4) symmetric crystallized phase and an SU (4) 1 chiral spin liquid at the balanced filling. We also predict two different exciton supersolid phases with interlayer coherence at imbalanced filling under displacement field. Especially, the system can simulate an SU(2) Bose-Einstein condensation by injecting interlayer excitons into the magnetically ordered Mott insulator at the layer polarized limit. Smoking gun evidences of these phases can be obtained by measuring the pseudospin transport in the counterflow channel.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Simulation of Boreal Ecosystem Carbon and Water Budgets: Scaling from Local to Regional Extents

A coupled water and energy balance model is developed. This model can predict the partitioning of water and energy between major source, sink and storage elements within the Boreal-Ecosystem-Atmospheric Study (BOREAS) areas. The results of testing the model against data collected at BOREAS tower sites during Intensive Field Campaigns and remotely sensed data collected across the BOREAS region are presented.

Wood, Eric F.↗
DOE OSTI2025

Bringing Inorganic Carbon to Life: Developing Model Metalloenzymes for C1 Conversion Reactions

The global carbon cycle is carefully balanced through the use of specialized enzymes in plants, algae, bacteria, and archaea. A primordial metabolic pathway for the conversion of inorganic carbon into cellular biomass uses a large, nickel-containing enzyme called carbon monoxide dehydrogenase (CODH)/acetyl coenzyme A synthase (ACS). This system fixes carbon dioxide (CO 2 ) into carbon monoxide (CO) at the CODH site. The carbon monoxide is then used to generate acetyl coenzyme A, a biological building block, through a key carbon-carbon bond forming step at the nickel site in ACS. Despite the significance of these processes in the context of energy conversion, the fundamental chemistry underlying these transformations has remained elusive, in part due to the complexity of the natural enzyme. To better understand the principles governing these biological processes, which represent one-carbon (C 1 ) activation reactions, this project will develop functional models of CODH and ACS based on modifications to a small metalloprotein scaffold. Comprehensive characterization of these systems using advanced spectroscopic and biophysical techniques will reveal key elements responsible for conferring high levels of activity to the model systems, with implications for gaining insight into the mechanisms of the natural enzymes. Moreover, the principles learned from this research can be used to guide design of robust catalysts for efficient conversion of CO 2 and CO into liquid fuels.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Ecosystem-Level Energy and Water Budgets Are Resilient to Canopy Mortality in Sparse Semiarid Biomes

Climate-driven woody vegetation mortality is a defining feature of semiarid biomes that drives fundamental changes in ecosystem structure. However, the observed impacts of woody mortality on ecosystem-scale energy and water budgets and the responses of surviving vegetation are highly variable among studies in water-limited environments. A previous girdling manipulation experiment in a piñon-juniper woodland suggested that although ecosystem-scale evapotranspiration was not altered by large-scale piñon mortality, soil water content decreased and the surviving juniper experienced greater water stress than juniper in an undisturbed woodland. Here we experimentally explored to what extent mortality-induced changes in energy balance components can explain these results. We compared energy fluxes measured above two adjacent piñon-juniper woodlands where piñon girdling was implemented at one site and the other subsequently experienced large-scale natural piñon mortality. We found that the mortality-induced decrease in canopy area was not sufficient to alter surface reflectance, roughness, and partitioning between energy budget components at both sites. A radiative transfer model estimated that because of the sparse premortality canopy, surface reflectance is more sensitive to a large increase in understory leaf area than further loss of crown area. Increased water stress in the remaining juniper following both mortality events can be explained by an increase in radiation on the ground that promoted higher soil temperature and evaporation. We found similar responses of ecosystem and tree-level functions to both girdling and natural mortality. This suggests that girdling is an appropriate approach to explore the impact of tree mortality on ecosystem structure, function, and energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Power capture and power take-off load of a self-balanced dual-flap oscillating surge wave energy converter

Wave energy converters are an important part of future renewable energy infrastructure. Predicting their power matrix, capture width ratio, and power take-off loads at a targeted site is required for performance assessment before deployment. Because their testing is very expensive, numerical modeling and simulations play a significant role in those assessments. Linear potential flow theory has limited accuracy under large amplitude wave forcing. More accurate predictions can be obtained by using higher-fidelity models, which are computationally expensive. We present a framework for multi-fidelity numerical simulations to determine the hydrodynamic response, wave capture capability, and power take-off load of a full-scale dual-flap oscillating surge wave energy converter. This design exploits out-of-phase motion by setting the distance between the flaps to half the wavelength of the most occurring wave. The simulations are validated using a 1:10 model experiments in a wave tank. Based on these validations, it was determined that Euler simulations provide an acceptable prediction with 90% reduction in computational time with only 11% error. Utilizing Euler simulations at full-scale, the results demonstrate that the annual electrical energy output is 1.79 GWh under regular wave conditions. Here, one significant improvement over single-flap designs is the capture width ratio which exceeds unity.

16 TIDAL AND WAVE POWER↗
NASA NTRSDate not supplied

TPSAS-NF1676L-33912-DND

The agricultural, renewable energy management, and science communities need global surface and top-of-atmosphere (TOA) radiative fluxes on a near real-time basis for uses such as building energy performance monitoring and seasonal crop yield modeling. The Clouds and Earth's Radiant Energy System (CERES) FLASHFlux (Fast Longwave and SHortwave radiative Flux) data products address this need by enhancing the speed of CERES processing using simplified calibration, averaging techniques and fast radiation parameterizations to produce global fluxes within a week of satellite observations. The CERES FLASHFlux provides two data products: 1) an hourly Level 2 Single Scanner Footprint (SSF) data separately for Aqua and Terra observations, and 2) a daily Level 3 Time Interpolated and Spatially Averaged (TISA) 1degree x 1degree gridded data that combines Aqua and Terra observations. The CERES FLASHFlux Version4A includes an update to process MODIS Collection 6.1 in an improved cloud retrieval algorithm, and an update to the CERES Edition 4 angular distribution models. We inter-compared FLASHFlux Version4A, FLASHFlux Version3C, and CERES EBAF (Energy Balanced and Filled) Top-of-Atmosphere fluxes to evaluate these improvements. We also compare both SSF overpass and TISA radiative flux products to a set of surface measurement sites that are globally distributed and assess agreement. FLASHFlux TISA and CERES EBAF are used in the analysis of TOA anomalies for the annual State of Climate report. This also poster highlights the methodology use to merge the two data products by correcting the bias in FLASHFlux TISA and inter-calibrate to CERES EBAF.

Parnchai K Sawaengphokhai↗
NASA NTRS2024

A Physics-Based Antarctic Melt Detection Technique: Combining Advanced Microwave Scanning Radiometer 2, Radiative-Transfer Modeling, and Firn Modeling

Surface melt on ice shelves has been linked to hydrofracture and subsequent ice shelf breakup. Since the 1990s, scientists have been using microwave radiometers to detect melt on ice shelves and ice sheets by applying various statistical thresholding techniques to identify significant increases in brightness temperature that are associated with melt. In this study, instead of using a fixed threshold, we force the Snow Radiative Transfer model (SMRT) with outputs from the Community Firn Model (CFM) to create a dynamic, physics-based threshold for melt. In the process, we also combine our method with statistical thresholding techniques and produce microwave grain size information in the process. We run this “Hybrid Method” across the Larsen C ice shelf as well as thirteen sites on the Antarctic Ice Sheet. Melt and non-melt days from the Hybrid Method and three statistical thresholding techniques match with surface energy balance within 94±1%; the effect of melt on the passive microwave is mostly binary and thus largely detectable by statistical thresholding techniques as well as physics-based techniques. Rather than always replacing statistical thresholding techniques with the Hybrid Method, we recommend using the Hybrid Method in studies where melt volume or grain size is of interest. In this study, we show that the Hybrid Method can be used to (a) model dry snow brightness temperatures of Antarctic snow and (b) derive a measure of grain size; therefore, it is an important step forwards towards using firn and radiative transfer modeling to quantify melt rather than to simply detect melt days.

Marissa E Dattler↗
DOE OSTI2024

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Statistical framework to assess long-term spatio-temporal climate changes: East River mountainous watershed case study

Abstract Evaluation of long-term temporal and spatial climatic change in mountainous regions is a critical challenge because of the interactive effects of multiple land and climatic factors and processes. Here we present the application of the statistical framework to the assessment of changes of climatic conditions, using data from 17 meteorological stations across the East River watershed near Crested Butte, Colorado, USA, and spanning the period from 1966 to 2021. The framework is developed based on (1) a time-series analysis of daily, monthly, and yearly averaged meteorological parameters (temperature, relative humidity, precipitation, wind speed, etc.), (2) evaluation and time series analysis of potential evapotranspiration (ET o ), actual evapotranspiration (ET), aridity index (AI), standard precipitation index (SPI) and standard precipitation-evapotranspiration index (SPEI), and (3) a temporal-spatial climatic zonation of the studied area based on the hierarchical clustering and PCA analysis of the SPEI, because the SPEI can be considered an integrative characteristic of the changes of climatic conditions. The Budyko model, with the application of the Penman–Monteith equation for the estimation of ET o , was used to determine the ET. The time series analysis of the AI is used to identify the periods with energy limited and water limited conditions. Hierarchical clustering of site locations for the three temporal segments of the SPEI showed a significant temporal-spatial shifts, indicating that dynamic climatic processes drive zonation patterns. Therefore, the watershed climatic zonation requires periodic re-evaluation based on the structural time series analysis of meteorological and water balance data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗