Search NASA⌕ Search

SEARCH · Search NASA

Results for “canopy conductance”

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 163 records · Page 9

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗

Shortwave infrared detection of vegetation

Shortwave infrared sensors were included on the Thematic Mapper (TM) to observe vegetation reflected radiance patterns related to water leaf content. Analysis of field measurements for corn and soybeans throughout the growing season showed that shortwave infrared measurements enhance discrimination between the species, particularly in midseason. A numerical model of the canopy reflectance showed that differential leaf absorptance can produce the observed patterns. Analysis of coincident studies of leaf optical properties were conducted to generalize the results to other types of vegetation.

Goward, S. N.↗

Monitoring the defoliation of hardwood forests in Pennsylvania using LANDSAT

An automated system for conducting annual gypsy moth defoliation surveys using LANDSAT MSS data and digital processing techniques is described. A two-step preprocessing procedure was developed that uses multitemporal data sets representing forest canopy conditions before and after defoliation to create a digital image in which all nonforest cover types are eliminated or masked out of a LANDSAT image that exhibits insect defoliation. A temporal window for defoliation assessment was identified and a statewide data base was established. A data management system to interface image analysis software with the statewide data base was developed and a cost benefit analysis of this operational system was conducted.

Dottavio, C. L.↗

Canopy cover and microtopography control precipitation-enhanced thaw of ecosystem-protected permafrost

Abstract Northern high-latitudes are projected to get warmer and wetter, which will affect rates of permafrost thaw and mechanisms by which thaw occurs. To better understand the impact of rain, as well as other factors such as snow depth, canopy cover, and microtopography, we instrumented a degrading permafrost plateau in south-central Alaska with high-resolution soil temperature sensors. The site contains ecosystem-protected permafrost, which persists in unfavorable climates due to favorable ecologic conditions. Our study (2020–2022) captured three of the snowiest years and three of the four wettest years since the site was first studied in 2015. Average thaw rates along an across-site transect increased nine-fold from 6 ± 5 cm yr −1 (2015–2020) to 56 ± 12 cm yr −1 (2020–2022). This thaw was not uniform. Hummock locations, residing on topographic high points with relatively dense canopy, experienced only 8 ± 9 cm yr −1 of thaw, on average. Hollows, topographic low points with low canopy cover, and transition locations, which had canopy cover and elevation between hummocks and hollows, thawed 44 ± 6 cm yr −1 and 39 ± 13 cm yr −1 , respectively. Mechanisms of thaw differed between these locations. Hollows had high warm-season soil moisture, which increased thermal conductivity, and deep cold-season snow coverage, which insulated soil. Transition locations thawed primarily due to thermal energy transported through subsurface taliks during individual rain events. Most increases in depth to permafrost occurred below the ∼45 cm thickness seasonally frozen layer, and therefore, expanded existing site taliks. Results highlight the importance of canopy cover and microtopography in controlling soil thermal inputs, the ability of subsurface runoff from individual rain events to trigger warming and thaw, and the acceleration of thaw caused by consecutive wet and snowy years. As northern high-latitudes become warmer and wetter, and weather events become more extreme, the importance of these controls on soil warming and thaw is likely to increase.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of the photochemical reflectance index in AVIRIS imagery

In this paper, we evaluate the potential for extracting the 'photochemical reflectance index' (PRI; previously called the 'physiological reflectance index') from AVIRIS data. This index, which is derived from narrow-band reflectance at 531 and 570 nm, has proven to be a useful indicator of photosynthetic function at the leaf and canopy scales. At the leaf level, PRI varies with photosynthetic capacity, radiation-use efficiency, and vegetation type (unpublished data). This finding is consistent with the hypothesis that vegetation types exhibiting chronically reduced photosynthesis during periods of stress (e.g. drought-tolerant evergreens) invest proportionally more in photoprotective processes than vegetation with high photosynthetic capacity (e.g. crops or deciduous perennials). Vertical transects in tropical and boreal forest canopies have indicated declines in PRI associated with downregulation of photosynthesis at the canopy tops under sunny, dry midday conditions (unpublished data). This reduced PRI in upper canopy levels provides a further basis for examining this signal with the 'view from above' afforded by aircraft overflights. Although many factors could confound interpretation of a subtle physiological signal at the landscape scale, we conducted a preliminary examination of PRI extracted from existing, AVIRIS imagery of Stanford University's Jasper Ridge Biological Preserve obtained on the June 2nd, 1992, overflight. The goal was to use the hyperspectral capabilities of AVIRIS to evaluate the potential of this index for obtaining useful physiological data at the landscape scale. The expectation based on leaf- and canopy-level studies was that regions containing vegetation of reduced photosynthetic capacity (e.g. chaparral or evergreen woodland) would exhibit lower PRI values than regions of high capacity (e.g. deciduous woodland).

Gamon, John A.↗

Bridging gaps in permafrost-shrub understanding

Permafrost, permanently frozen ground which underlies much of the Arctic and sub-Arctic, is influenced by deciduous shrubs, yet our understanding of permafrost-shrub interactions is limited. This is largely due to a lack of widespread, long-term, high-quality observations across Arctic and sub-Arctic systems, which are difficult to study due to their remote locations. Shrubs are rapidly expanding in many areas of Arctic tundra and can either amplify or inhibit permafrost thaw, making it crucial that we can understand and predict permafrost-shrub interactions. As we have limited time and resources to make observations, we must design our field campaigns to be as impactful as possible. Below we outline the current state of knowledge and give suggestions about how we can maximize our fieldwork’s impact. From field studies conducted around the Arctic, we know that shrubs generally have a cooling effect on permafrost in the summer, largely caused by the shrub canopy shading the soil. However, summertime cooling can be counteracted by a warming effect in winter, caused by the trapping of blowing snow by shrubs, which leads to increased snow depths. In addition to these commonly observed effects, both observations and models show that other local factors such as climate, soil, snow, and disturbances can interact to cause contrasting dominant effects on permafrost temperatures. For example, shrubs that protrude through the snowpack have been observed to cool permafrost in winter via thermal bridging between the air above the snow and the soil below. However, at other research sites, protruding shrubs lower the surface albedo and cause earlier snowmelt in spring, amplifying permafrost thaw by exposing the soil to thawing earlier. Protruding shrubs have also been observed to induce snow melt in autumn, creating ice layers at the snow surface that prevent further snow accumulation and eliminate the ability of shrubs to trap blowing snow, thereby reducing the winter warming effect. While we understand the physical processes that lead to these effects, we cannot explain why we observe different effects between sites.

54 ENVIRONMENTAL SCIENCES↗

Mangrove Carbon Stocks in Pongara National Park, Gabon

Mangroves are recognized for their valued ecosystem services to coastal areas, and the functional linkages between those services and ecosystem carbon stocks have been established. However, spatially explicit inventories are necessary to facilitate management and protection of mangroves, as well as providing a foundation for payment for ecosystem service programs such as REDD+. We conducted an inventory of carbon stocks in mangroves within Pongara National Park (PNP), Gabon using a stratified random sampling design based on forest canopy height derived from TanDEM-X remote sensing data. Ecosystem carbon pools, including aboveground and belowground biomass and necromass, and soil carbon to a depth of 2 m were assessed using measurements and samples from plots distributed among three canopy height classes within the park. There were two mangrove species within the inventory area in PNP, Rhizophora racemosa and R. harrisonii. R. harrisonii was predominant in the sparse, low-stature stands that dominated the west side of the park. In the east side of the park, both species occurred in tall-stature stands, with tree height often exceeding 30 m. Canopy height was an effective means to stratify the inventory area, as biomass was significantly different among the height classes. Despite those differences in aboveground biomass, the soil carbon density was not significantly different among height classes. Soils were the main component of the ecosystem carbon stock, accounting for over 84% of the total. The ecosystem carbon density ranged from 644 to 943 Mg C ha−1 among the three height classes. The ecosystem carbon stock within PNP is estimated to be 40,588 Gg C. The combination of pre-inventory information about stand conditions and their spatial distribution within the assessment area obtained from remote sensing data and a spatial decision support system were fundamental to implementing this relatively large-scale field inventory. This work exemplifies how mangrove carbon stocks can be quantified to augment national C reporting statistics, provide a baseline for projects involving monitoring, reporting and verification (i.e., MRV), and provide data on the forest composition and structure for sustainable management and conservation practices.

Carl C Trettin↗

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗

Drag characteristics of a disk-gap-band parachute with a nominal diameter of 1.65 meters at Mach numbers from 2.0 to 3.0

Supersonic wind-tunnel tests were conducted with disk-gap-band parachute models having a nominal diameter of 1.65 meters and geometric porosities of 10.0, 12.5, and 15.0 percent. Canopy inflation characteristics, angles of attack, and drag performance are presented for deployment behind forebody base extensions which were free to oscillate in pitch and yaw. The effect of increasing suspension-line length on canopy motions and drag performance is included, and the drag performance of a model with 12.5 percent geometric porosity is compared with results from flight tests of a parachute with a nominal diameter of 12.19 meters.

Mayhue, R. J.↗

Diurnal Reflectance Changes in Vegetation Observed with AVIRIS

Among the most important short-term dynamic biological processes are diurnal changes in canopy water relations. Plant regulation of water transport through stomatal openings affects other gaseous transport processes, often dramatically decreasing photosynthetic fixation of carbon dioxide during periods of water stress. Water stress reduces stomatal conductance of water vapor through the leaf surface and alters the diurnal timing of stomatal opening. Under non-water stressed conditions, stomates typically open soon after dawn and transpire water vapor throughout the daylight period. During stress periods, stomates may close for part of the day, generally near mid-day. Under prolonged stress conditions, stomatal closure shifts to earlier times during the day; stomates may close by mid-morning and remain closed until the following morning - or remain closed entirely. Under these conditions the relationship between canopy greenness (e.g., measured with a vegetation index or by spectral mixture analysis) and photosynthetic fixation of carbon is lost and the remotely sensed vegetation metric is a poor predictor of gas exchange. Prediction of stomatal regulation and exchange of water and trace gases is critical for ecosystem and climate models to correctly estimate budgets of these gases and understand or predict other processes like gross and net ecosystem primary production. Plant gas exchange has been extensively studied by physiologists at the leaf and whole plant level and by biometeorologists at somewhat larger scales. While these energy driven processes follow a predictable if somewhat asymmetric diurnal cycle dependent on soil water availability and the constraints imposed by the solar energy budget, they are nonetheless difficult to measure at the tree and stand levels using conventional methods. Ecologists have long been interested in the potential of remote sensing for monitoring physiological changes using multi-temporal images. Much of this research has focused on day-to-day changes in water use, especially for agricultural applications. Ustin et al. showed seasonal changes in canopy water content in chaparral shrub could be estimated using optical methods. Vanderbilt et al. followed asymmetric diurnal changes in the reflectance of a walnut orchard, but could not attribute specific reflectance changes to specific changes in canopy architecture or physiology. Forests and shrub lands in California experience prolonged periods of drought, sometimes extending six months without precipitation. The conifer and evergreen chaparral communities common to the foothill region around the central valley of California retain their foliage throughout the summer and have low transpiration rates despite high net radiation and temperature conditions. In contrast, grasslands and drought resistant deciduous species in the same habitat are seasonally dormant in summer. Because of differences in the mechanisms of drought tolerance, rooting depth and physiology between different plant communities in the region, it is likely that they display differences in diurnal water relations. The presence of diverse plant communities provides an opportunity to investigate possible diurnal landscape patterns in water relations that could be observed by an airborne hyperspectral scanner. This investigation of AVIRIS data collected over forest and shrub land represents the continuation of a prior investigation involving spectral mixture analysis of diurnal effects in the same AVIRIS data set.

Vanderbilt, V. C.↗

Stem hydraulic conductivity and vulnerability to cavitation for 26 tree species in Panama

Stem hydraulic conductivity and vulnerability to cavitation were measured for 26 tree species located in Panama. The data were generated to better understand the ecology of the focal tree species. Complementary NGEE-Tropics datasets for these species include sap flow, leaf-level gas exchange, and leaf water potential. Stem samples were collected from distal branches of canopy trees, brought to the Smithsonian Tropical Research Institute laboratory in Gamboa, Panama, and allowed to dry to various water potentials before measurements. For each species, a Weibull function was fit to the 90% quantile of the relationship between stem area specific hydraulic conductivity (Ks) and stem water potential. From these functions, maximum Ks and vulnerability parameters were derived. The data files in the package include raw data, derived parameters, and the R script used for analysis.

54 ENVIRONMENTAL SCIENCES↗

Street-level temperature estimation using graph neural networks: Performance, feature embedding and interpretability

Estimating street-level air temperature is a challenging task due to the highly heterogeneous urban surfaces, canyon-like street morphology, and the diverse physical processes in the built environment. Though pioneering studies have embarked on investigations via data-driven approaches, many questions remain to be answered. Here, in this study, we leveraged an innovative framework and redefined the street-level temperature estimation problem using Graph Neural Networks (GNN) with spatial embedding techniques. The results showed that GNN models are more capable and consistent of estimating street-level temperature among tested locations, benefiting from its unique strength in handling extensive data over unstructured graph topology. In addition, we conducted in-depth analysis of feature importance to enhance the model interpretability. Among the urban features analyzed in this study, the time-variant canopy density and meter-level land use data emerge as crucial factors. Our findings highlight GNN 's high potential in capturing the complex dynamics between urban elements and their impacts on microclimate, thus offering valuable insights for comprehensive urban data collection and urban climate modeling in general. Collectively, this study also contributes to urban planning and policy by providing avenues to enhance city resilience against climate change, thereby advancing the agenda for environmental stewardship and urban sustainability.

54 ENVIRONMENTAL SCIENCES↗

C-band microwave scattering from small balsam fir

An experiment to examine the C-band backscattering characteristics of conifer trees was conducted using a truck-mounted scatterometer. Small (1 m tall) balsam fir (Abies balsamea) were arranged at various equidistant spacings on a platform to present canopies of varying density to the radar. C-band backscattering measurements of a range of canopy densities were acquired under different polarizations and incidence angles. The measured backscattering coefficient from the tree canopies increased with increasing biomass, but approached a maximum at a LAI of 2.5 and fresh biomass of 3.3 kg/sq m. A backscatter model was implemented using measured canopy attributes and showed close agreement with scatterometer measurements over the range of canopy densities. Model results indicated that branches were the prime scatterers of the radar while needles were found to only slightly attenuate the radar signal.

Ranson, K. J.↗

BOREAS TE-5 Leaf Gas Exchange Data

The BOREAS TE-5 team collected measurements in the NSA and SSA on gas exchange, gas composition, and tree growth. The leaf photosynthetic gas exchange data were collected in the BOREAS NSA and the SSA from 06-Jun- 1994 to 13-Sep- 1994 using a LI-COR 6200 portable photosynthesis system. The data were collected to compare the photosynthetic capacity, stomata] conductance, and leaf intercellular CO, concentrations among the major tree species at the BOREAS sites. The data are average values from diurnal measurements on the upper canopy foliage (sun leaves). The data are available in tabular ASCII files. The data files are available on a CD-ROM (see document number 20010000884), or from the Oak Ridge National Laboratory (ORNL) Distributed Activity Archive Center (DAAC).

Hall, Forrest G.↗

Flow Visualization Techniques in Wind Tunnel Tests of a Full-Scale F/A-18 Aircraft

The proposed paper presents flow visualization performed during experiments conducted on a full-scale F/A-18 aircraft in the 80- by 120-Foot Wind-Tunnel at NASA Ames Research Center. The purpose of the flow-visualization experiments was to document the forebody and leading edge extension (LEX) vortex interaction along with the wing flow patterns at high angles of attack and low speed high Reynolds number conditions. This investigation used surface pressures in addition to both surface and off-surface flow visualization techniques to examine the flow field on the forebody, canopy, LEXS, and wings. The various techniques used to visualize the flow field were fluorescent tufts, flow cones treated with reflective material, smoke in combination with a laser light sheet, and a video imaging system for three-dimension vortex tracking. The flow visualization experiments were conducted over an angle of attack range from 20 deg to 45 deg and over a sideslip range from -10 deg to 10 deg. The various visualization techniques as well as the pressure distributions were used to understand the flow field structure. The results show regions of attached and separated flow on the forebody, canopy, and wings as well as the vortical flow over the leading-edge extensions. This paper will also present flow visualization comparisons with the F-18 HARV flight vehicle and small-scale oil flows on the F-18.

Lanser, Wendy R.↗

Full-Waveform, Wide-Swath Lidar Imaging of Forested and Urban Areas in Leaf-On Conditions: Development, Results and Future Direction

Full-Waveform lidar measurements provide unprecedented views of the vertical and horizontal structure of vegetation and the topography of the Earth s surface. Utilizing a high signal-to-noise ratio lidar system, larger than typical laser footprints (10-20 m), and the recorded time history of interaction between a short-duration (approx. 10 ns) pulse of laser light and the surface of the Earth, full-waveform lidar is able to simultaneously image sub-canopy topography as well as the vertical structure of any overlying vegetation. These data reveal the true 3-D vegetation structure in leaf-on conditions enabling important biophysical parameters such as above-ground biomass to be estimated with unprecedented accuracy. An airborne lidar mission was conducted July-August 2003 in support of the North America Carbon Program. NASA s Laser Vegetation Imaging Sensor (LVIS) was used to image approximately 2,000 km$^2$ in Maine, New Hampshire, Massachusetts and Maryland. Areas with available ground and other data were included (e.g., experimental forests, FLUXNET sites) in order to facilitate as many bio- and geophysical investigations as possible. Data collected included ground elevation and canopy height measurements for each laser footprint, as well as the vertical distribution of intercepted surfaces. Data will be publicly distributed within 6- 12 months of collection. Further details of the mission, including the lidar system technology, the locations of the mapped areas, and examples of the numerous data products that can be derived from the return waveform data products will be presented. Future applications including detection of ground and vegetation canopy changes and a spaceborne implementation of wide-swath, full-waveform imaging lidar will also be discussed.

Blair, James B.↗

Field testing of thermal canopy models in a spruce-fir forest

Recent advances in remote sensing technology allow the use of the thermal infrared region to gain information about vegetative surfaces. Extending existing models to account for thermal radiance transfers within rough forest canopies is of paramount importance. This is so since all processes of interest in the physical climate system and biogeochemical cycles are thermally mediated. Model validation experiments were conducted at a well established boreal forest; northern hardwood forest ecotone research site located in central Maine. Data was collected to allow spatial and temporal validation of thermal models. Emphasis was placed primarily upon enhancing submodels of stomatal behavior, and secondarily upon enhancing boundary layer resistance submodels and accounting for thermal storage in soil and vegetation.

Source record↗