Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LaI”

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

SPRUCE: Shrub Leaf and Stem Functional Traits from 10-year Shrub Layer Biomass Harvests (August 2025)

This data set contains Leaf Mass Area (LMA), Specific Leaf Area (SLA), Leaf Area Index (LAI), and Huber value (Hv) measurements from the 10-year Shrub-layer Biomass harvests in the Spruce and Peatland Responses Under Environmental Change (SPRUCE) experiment in August 2025. The total ground-level cross-sectional sapwood area of Chamaedaphne calyculata (CHCA), Rhododendron groenlandicum (RHGR), and Vaccinium angustifolium (VAAN) were calculated for each of the three shrub-layer community survey plots located in each experimental enclosure. Representative current-year leaves were removed, scanned, dried, and weighed to derive SLA. Grab-samples for VAAN were used to calculate SLA when no biomass was harvested in the shrub-layer community plots. These SLA values were used to infer total one-sided leave surface area from the total species lead dry mass for calculation of LAI and Hv.

Birkebak, Joshua [ORNL] (ORCID:0009000955611494)↗

AmeriFlux FLUXNET-1F US-Me2 Metolius mature ponderosa pine

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Me2 Metolius mature ponderosa pine. This is the FLUXNET version of the carbon flux data for the site US-Me2 Metolius mature ponderosa pine produced by applying the standard ONEFlux (1F) software. Site Description - Site Description before Fire (January, 2002 - August, 2020): The mean stand age is 71 years old and the stand age of the oldest 10% of trees is about 108 years old. This site is one of the Metolius core cluster sites with different age and disturbance classes and part of the AmeriFlux network. The overstory is almost exclusively composed of ponderosa pine trees (Pinus ponderosa Doug. Ex P. Laws) with a few scattered incense cedars (Calocedrus decurrens (Torr.) Florin) and has a peak leaf area index (LAI) of 2.1 m2 m-2. Tree height is relatively homogeneous at about 18 m, and the mean tree density is approximately 339 trees ha-1 (Irvine et al., 2008). The understory is sparse with an LAI of 0.2 m2 m-2 and primarily composed of bitterbrush (Purshia tridentata (Push) DC.) and greenleaf manzanita (Arctostaphylos patula Greene). Soils at the site are sandy (69%/24%/7% sand/silt/clay at 0–0.2 m depth and 66%/27%/7% at 0.2–0.5 m depth, and 54%/ 35%/11% at 0.5–1.0 m depth), freely draining with a soil depth of approximately 1.5 m (Irvine et al., 2008; Law et al., 2001b; Schwarz et al., 2004). Green Ridge Fire: On August 20, 2020, the Green Ridge Fire burned through Us-Me2. The fire was ignited by lightning on August 16th, and grew rapidly to the east over the first few days driven by strong, downslope (westerly) afternoon winds. Fire behavior and observed fire effects were highly heterogeneous due to the localized wind pattern carrying the flaming head of the fire forward, and the efforts being made by suppression resources to contain the fire. The site experienced the full range of fire effects, from <1 m high surface fire that charred litter and duff and only consumed shrubs and herbaceous material to full tree (>15 m) crown fire that consumed 100% of needles, small limbs, and surface fuels at high intensity, leaving only ash and bare soil post-fire. Salvage Logging: From late March to late April 2021, salvage logging by the landowner occurred at the site. Almost all trees within the flux footprint were logged except a small area with lower burn severity, where sap flow and automatic soil respiration measurements are continued since the fire. In August 2022 the primary flux system was moved from the top of the damaged tall tower to a nearby shorter tower centering on the salvage/regenerating footprint.

Hanson, Chad↗

AmeriFlux FLUXNET-1F US-Lin Lindcove Orange Orchard

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Lin Lindcove Orange Orchard. This is the FLUXNET version of the carbon flux data for the site US-Lin Lindcove Orange Orchard produced by applying the standard ONEFlux (1F) software. Site Description - The experimental site was a citrus orchard about three km west of the UC Lindcove Research and Experiment Station. The site is characterized by a Mediterranean climate typical of Central California, with warm dry summers and cool wet winters. The soil texture for the upper horizon (0-33 cm) was loam with a particle size distribution of 42% sand, 38% silt, and 20% clay (Kearney Ag Center). Soil pH was 7.4 (Valley Tech Soil Agricultural Laboratory Services). The orchard was drip-irrigated with water applied twice per week in the warm season to ensure water avail- ability close to field capacity. The site was kept clean from understory vegetation and mechanical pruning operations took place twice a year to limit the size ofthe orange trees for harvesting and site maintenance. The block of trees in which the tower and instruments were located was ‘Valencia’ orange on trifoliate rootstock, with a planting date in the 1960’s. The square 4 ha block had dimensions of 200 m NeS and EeW. At the end of the study, we harvested a ‘Valencia’ citrus tree from within the study block to measure citrus leaf area index (LAI) and biomass density. The tree height was 3.7 m as measured with a telescoping pole. We measured canopy radii in the four cardinal directions to calculate a planar area (as seen from above) of 12.2 m2. The mean specific leaf area (SLA) of citrus leaves was 85.4 cm2 g-1 obtained by measuring leaf area of five groups of leaf samples with a LiCor Leaf area meter (mod. LI-3100C) followed by drying and weighing. LAI for the orchard was 3.00, derived using the SLA value and the total dry mass of leaves from the harvested tree, multiplied by the plant population of 237 plants per ha as obtained from spacing measurements and Google Maps imagery.

Fares, Silvano↗

Evaluation of the Large-Scale and Regional Climatic Response Across North Africa to Natural Variability in Oceanic Modes and Terrestrial Vegetation Among the CMIP5 Models

Hydrologic variability is a serious threat to the poverty-stricken regions of North Africa. Meanwhile, the scientific community struggles to attribute these extreme climatic episodes to specific oceanic and land drivers. Prior modeling studies have not assessed simulated feedbacks over North Africa against an established observational benchmark, so their results are considered as untested, model-specific findings. The Coupled Model Intercomparison Project Phase Five (CMIP5) archive represents the state-of-the-art in climate projections, with most CMIP5 models now containing interactive vegetation phenology; however, there have been no studies to date of their simulated vegetation feedbacks. The capability of CMIP5 models at accurately simulating the forcing of both sea-surface temperature (SST) and regional leaf area index (LAI) anomalies on North African climate needs to be a key consideration in assessing the models’ overall credibility and determining appropriate model weighting for developing climate projections. The Generalized Equilibrium Feedback Assessment (GEFA) is a promising statistical method, which can be applied to either model output or observations, for isolating the local and remote impacts of individual oceanic or terrestrial forcings on regional climate. We performed a combined observational and modeling assessment of land-ocean-atmosphere interactions across the distinct ecological and moisture gradients of North Africa. We evaluated and demonstrated the reliability of the GEFA statistical method over North Africa using the Community Earth System Model (CESM), applied GEFA to observational data to quantify the observed forcing of ocean basin SST anomalies and regional LAI anomalies on North African climate, evaluated the CMIP5 models’ performance in terms of representing these key observed feedbacks, and formulated CMIP5 feedback performance metrics for weighting North African climate projections. Our study represented the first attempt to separate the observed roles of oceanic and vegetation feedbacks across North Africa, the first systematic assessment and intercomparison of land-ocean-atmosphere feedbacks in CMIP5, and the first exploration of vegetation feedbacks among CMIP5 models. The following work was accomplished. (1) The team generated 14 scientific publications and 30 presentations based on the DOE-funded research. (2) A stepwise version of the GEFA statistical method was developed in which unimportant forcings were dropped to increase the reliability of the results. (3) SGEFA was successfully validated through experiments with CESM in terms of its ability to isolate the atmospheric responses to individual oceanic or land forcings. (4) Observational evidence was revealed for the Sahel’s positive vegetation-rainfall feedback on the seasonal to interannual time scale, and it was attributed to a moisture recycling mechanism rather than an albedo mechanism. (5) An approach with SGEFA was developed in which the individual contributions of soil moisture versus vegetation forcings could be separated, revealing that the former forcing outweighed the latter for sub-Saharan Africa. (6) Tropical ocean temperatures were found to be key regulators of pan-tropical vegetation variability, especially for arid and semi-arid regions, including sub-Saharan Africa. (7) The CMIP5 models largely underestimated the importance of land feedbacks across the Sahel. (8) The general consensus among CMIP5 models indicates, for the late 21st century, a diminished seasonal predictability of sub-Saharan African regional climate and an elevated role of the land surface compared to oceanic drivers in regulating regional climate variability. (9) Land-ocean-atmosphere interactions were demonstrated to be key contributors to the seasonal predictability of African wildfire activity. (10) The El Djouf was determined to be of greater importance than the Bodélé depression in terms of providing trans-Atlantic dust transport to the Americas.

54 ENVIRONMENTAL SCIENCES↗

Estimating Leaf Area Index in Row Crops Using Wheel-Based and Airborne Discrete Return Light Detection and Ranging Data

Leaf area index (LAI) is an important variable for characterizing plant canopy in crop models. It is traditionally defined as the total one-sided leaf area per unit ground area and is estimated by both direct and indirect methods. This paper explores the effectiveness of using light detection and ranging (LiDAR) data to estimate LAI for sorghum and maize with different treatments at multiple times during the growing season from both a wheeled vehicle and Unmanned Aerial Vehicles. Linear and nonlinear regression models are investigated for prediction utilizing statistical and plant structure-based features extracted from the LiDAR point cloud data with ground reference obtained from an in-field plant canopy analyzer (indirect method). Results based on the value of the coefficient of determination ( R 2 ) and root mean squared error for predictive models ranged from ∼0.4 in the early season to ∼0.6 for sorghum and ∼0.5 to 0.80 for maize from 40 Days after Sowing to harvest.

59 BASIC BIOLOGICAL SCIENCES↗

Tropical Tree Crop Simulation with a Process-Based, Daily Timestep Simulation Model (ALMANAC): Description of Model Adaptation and Examples with Coffee and Cocoa Simulations

Coffee (Coffea species) and Cocoa (Theobroma cacao) are important cash crops grown in the tropics but traded globally. This study was conducted to apply the ALMANAC model to these crops for the first time, and to test its ability to simulate them under agroforestry management schemes and varying precipitation amounts. To create this simulation, coffee was grown on a site in Kaua’i, Hawai’i, USA, and cocoa was grown on a site in Sefwi Bekwai, Ghana. A stand-in for a tropical overstory tree was created for agroforestry simulations using altered parameters for carob, a common taller tropical tree for these regions. For both crops, ALMANAC was able to realistically simulate yields when compared to the collected total yield data. On Kaua’i, the mean simulated yield was 2% different from the mean measured yield, and in all three years, the simulated values were within 10% of the measured values. For cocoa, the mean simulated yield was 3% different from the mean measured yield and the simulated yield was within 10% of measured yields for all four available years. When precipitation patterns were altered, in Ghana, the wetter site showed lower percent changes in yield than the drier site in Hawai’i. When agroforestry-style management was simulated, a low Leaf Area Index (LAI) of the overstory showed positive or no effect on yields, but when LAI climbed too high, the simulation was able to show the detrimental effect this competition had on crop yields. These simulation results are supported by other literature documenting the effects of agroforestry on tropical crops. This research has applied ALMANAC to new crops and demonstrated its simulation of different management and environmental conditions. The results show promise for ALMANAC’s applicability to these scenarios as well as its potential to be further tested and utilized in new circumstances.

60 APPLIED LIFE SCIENCES↗

The Effect of Rapid Development on Soil CO2 Efflux in a Cellulosic Biofuel Stand

As awareness of climate change increases, the need for carbon neutral fuel sources is growing. Lignocellulosic biofuel derived from pine trees has been suggested as one potential energy source; however, it requires more research before its efficacy for climate change mitigation can be determined. Due to the large share of forest carbon held in soils and the extensive area of pine plantations in the southeast U.S., a better understanding of plantation soil carbon dynamics is critical for biofuel carbon accounting. This study evaluated the effects of canopy development and productivity on soil CO2 efflux, a proxy for soil respiration (Rs), in an intensively managed loblolly pine (Pinus taeda) stand over a period from May 2015 to December 2019. We found that leaf area index (LAI) and gross ecosystem production (GEP), as well as meteorological variables, had significant effects on Rs, but that both overall Rs and soil carbon pools did not increase over the course of the study. We thus hypothesize that GEP and LAI had intra-annual effects on Rs, and that the lack of change in Rs is the result of an increase in autotrophic respiration (Ra) that offset a decrease in decomposition of the previous stand’s organic matter.

09 BIOMASS FUELS↗

Challenging a Global Land Surface Model in a Local Socio-Environmental System

Land surface models (LSMs) predict how terrestrial fluxes of carbon, water, and energy change with abiotic drivers to inform the other components of Earth system models. Here, we focus on a single human-dominated watershed in southwestern Michigan, USA. We compare multiple processes in a commonly used LSM, the Community Land Model (CLM), to observational data at the single grid cell scale. For model inputs, we show correlations (Pearson’s R) ranging from 0.46 to 0.81 for annual temperature and precipitation, but a substantial mismatch between land cover distributions and their changes over time, with CLM correctly representing total agricultural area, but assuming large areas of natural grasslands where forests grow in reality. For CLM processes (outputs), seasonal changes in leaf area index (LAI; phenology) do not track satellite estimates well, and peak LAI in CLM is nearly double the satellite record (5.1 versus 2.8). Estimates of greenness and productivity, however, are more similar between CLM and observations. Summer soil moisture tracks in timing but not magnitude. Land surface reflectance (albedo) shows significant positive correlations in the winter, but not in the summer. Looking forward, key areas for model improvement include land cover distribution estimates, phenology algorithms, summertime radiative transfer modelling, and plant stress responses.

54 ENVIRONMENTAL SCIENCES↗

Global evaluation of terrestrial biogeochemistry in the Energy Exascale Earth System Model (E3SM) and the role of the phosphorus cycle in the historical terrestrial carbon balance

Abstract. The importance of carbon (C)–nutrient interactions to the prediction of future C uptake has long been recognized. The Energy Exascale Earth System Model (E3SM) land model (ELM) version 1 is one of the few land surface models that include both N and P cycling and limitation (ELMv1-CNP). Here we provide a global-scale evaluation of ELMv1-CNP using the International Land Model Benchmarking (ILAMB) system. We show that ELMv1-CNP produces realistic estimates of present-day carbon pools and fluxes. Compared to simulations with optimal P availability, simulations with ELMv1-CNP produce better performance, particularly for simulated biomass, leaf area index (LAI), and global net C balance. We also show ELMv1-CNP-simulated N and P cycling is in good agreement with data-driven estimates. We compared the ELMv1-CNP-simulated response to CO2 enrichment with meta-analysis of observations from similar manipulation experiments. We show that ELMv1-CNP is able to capture the field-observed responses for photosynthesis, growth, and LAI. We investigated the role of P limitation in the historical balance and show that global C sources and sinks are significantly affected by P limitation, as the historical CO2 fertilization effect was reduced by 20 % and C emission due to land use and land cover change was 11 % lower when P limitation was considered. Our simulations suggest that the introduction of P cycle dynamics and C–N–P coupling will likely have substantial consequences for projections of future C uptake.

54 ENVIRONMENTAL SCIENCES↗

Modeling the mechanisms of coastal vegetation dynamics and ecosystem responses to changing water levels

Coastal forests are increasingly experiencing mortality due to inundation by fresh- and seawater, leading to their replacement by marshes. These shifts alter vegetation composition, biogeochemical cycling, carbon storage, and hydrology. Using a hydraulically enabled ecosystem demography model (FATES-Hydro), we conducted numerical experiments to investigate the mechanisms behind inundation-driven forest loss and the ecosystem-scale consequences of forest-to-marsh transitions. We compared mortality processes and their effects across broadleaf and conifer trees at two coastal sites – Lake Erie (freshwater) and Chesapeake Bay (saline). Our simulations show that hydraulic failure, driven by root loss under prolonged flooding, is the primary mortality mechanism across both tree types and sites. Forest replacement by marsh reduced ecosystem-scale leaf area index (LAI), gross primary production (GPP), transpiration, and deep soil water uptake in conifer forests, while broadleaf forests experienced smaller changes due to lower initial LAI and greater marsh compensation. Marsh invasion occurred following canopy thinning driven by tree mortality. These findings suggest that, under similar root loss, hydraulic failure dominates coastal tree mortality regardless of species or water type, with denser forests experiencing stronger ecosystem impacts. Our study identifies key mortality mechanisms and offers testable hypotheses for future empirical studies on coastal vegetation change.

Ding, Junyan [Occidental College, Los Angeles, CA ↗

Development of a plant carbon–nitrogen interface coupling framework in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0)

Plant and microbial nitrogen (N) dynamics and N availability regulate the photosynthetic capacity and capture, allocation, and turnover of carbon (C) in terrestrial ecosystems. Studies have shown that a wide divergence in representations of N dynamics in land surface models leads to large uncertainties in the biogeochemical cycle of terrestrial ecosystems and then in climate simulations as well as the projections of future trajectories. In this study, a plant C–N interface coupling framework is developed and implemented in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0). The main concept and structure of this plant C–N framework and its coupling strategy are presented in this study. This framework takes more plant N-related processes into account. The dynamic ratio (CNR) for each plant functional type (PFT) is introduced to consider plant resistance and adaptation to N availability to better evaluate the plant response to N limitation. Furthermore, when available N is less than plant N demand, plant growth is restricted by a lower maximum carboxylation capacity of RuBisCO (V c,max ), reducing gross primary productivity (GPP). In addition, a module for plant respiration rates is introduced by adjusting the respiration with different rates for different plant components at the same N concentration. Since insufficient N can potentially give rise to lags in plant phenology, the phenological scheme is also adjusted in response to N availability. All these considerations ensure a more comprehensive incorporation of N regulations to plant growth and C cycling. This new approach has been tested systematically to assess the effects of this coupling framework and N limitation on the terrestrial carbon cycle. Long-term measurements from flux tower sites with different PFTs and global satellite-derived products are employed as references to assess these effects. The results show a general improvement with the new plant C–N coupling framework, with more consistent emergent properties, such as GPP and leaf area index (LAI), compared to the observations. The main improvements occur in tropical Africa and boreal regions, accompanied by a decrease in the bias in global GPP and LAI by 16.3 % and 27.1 %, respectively.

54 ENVIRONMENTAL SCIENCES↗

Benchmarking soil moisture and its relationship to ecohydrologic variables in Earth System Models

Soil moisture (SM) is a key regulator of ecosystem biogeophysics, influencing plant water relations and land-atmosphere energy exchanges. We evaluate the representation of SM in 16 Earth System Models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) using the International Land Model Benchmarking (ILAMB) framework, focusing on surface (0–5, 0–10 cm) and rootzone (0–100 cm) depths, as well as key ecohydrological variables like gross primary productivity (GPP), leaf area index (LAI), and evapotranspiration (ET), and their coupling. Models are benchmarked against multiple observational and assimilated datasets to assess both state variables and cross-variable relationships. Surface SM is generally well represented (r > 0.87), while rootzone SM variability is systematically overestimated (normalized standard deviation > 1). ET shows strong agreement with observations (r > 0.9), whereas GPP and LAI exhibit larger inter-model spread. Skill in individual variables does not guarantee realistic SM–ecohydrology coupling, which varies strongly across models and depends on the reference dataset. Köppen-based regional analyses reveal strong regime dependence, with several models performing well in Tropical and Temperate regions but degrading in Continental (high-latitude) zones. Across both global and regional benchmarks, models cluster by land surface framework, indicating that structural choices in soil hydrology and soil–plant coupling exert a first-order control on performance. These results provide process-relevant benchmarks and suggest that improving the representation of vertical soil structure, rooting depth distributions, and soil–plant hydraulic coupling will be central to advancing soil moisture realism in next-generation Earth system models.

CMIP6↗

Further tests of the suits reflectance model

Experiments performed by stacking cotton leaves in the port of a spectroradiometer indicate that single leaf reflectance ceases to vary with more than two leaves in the visible region and eight leaves in the infrared region. Chance and LeMaster have shown that the Suits spectral reflectance model predicts an asymptotic dependence of crop reflectance on leaf area index (LAI) with crop reflectance static for leaf area indices in excess of two in the visible regions and six in the infrared regions of the spectrum. These results are experimentally verified in the field for Milam and Penjamo spring wheat, and a theoretical relationship is discussed that relates crop reflectance at 650 nm to crop canopy LAI. Experimental data are given that relate observer zenith angle to crop reflectance for wheat. The Suits reflectance model calculations for wheat fail to agree with this data.

Lemaster, E. W.↗

Evaluating soil moisture and yield of winter wheat in the Great Plains using Landsat data

Locating areas where soil moisture is limiting to crop growth is important for estimating winter-wheat yields on a regional basis. In the 1975-76 growing season, we evaluated soil-moisture conditions and winter-wheat yields for a five-state region of the Great Plains using Landsat estimates of leaf area index (LAI) and an evapotranspiration (ET) model described by Kanemasu et al (1977). Because LAI was used as an input, the ET model responded to changes in crop growth. Estimated soil-water depletions were high for the Nebraska Panhandle, southwestern Kansas, southeastern Colorado, and the Texas Panhandle. Estimated yields in five-state region ranged from 1.0 to 2.9 metric ton/ha.

Heilman, J. L.↗

A three-part geometric model to predict the radar backscatter from wheat, corn, and sorghum

A model to predict the radar backscattering coefficient from crops must include the geometry of the canopy. Radar and ground-truth data taken on wheat in 1979 indicate that the model must include contributions from the leaves, from the wheat head, and from the soil moisture. For sorghum and corn, radar and ground-truth data obtained in 1979 and 1980 support the necessity of a soil moisture term and a leaf water term. The Leaf Area Index (LAI) is an appropriate input for the leaf contribution to the radar response for wheat and sorghum, however the LAI generates less accurate values for the backscattering coefficient for corn. Also, the data for corn and sorghum illustrate the importance of the water contained in the stalks in estimating the radar response.

Ulaby, F. T.↗

Relating the radar backscattering coefficient to leaf-area index

The relationship between the radar backscattering coefficient of a vegetation canopy, sigma(0) sub can, and the canopy's leaf area index (LAI) is examined. The relationship is established through the development of a model for corn and sorghum and another for wheat. Both models are extensions of the cloud model of Attema and Ulaby (1978). Analysis of experimental data measured at 8.6, 13.0, 17.0, and 35.6 GHz indicates that most of the temporal variations of sigma(0) sub can can be accounted for through variations in green LAI alone, if the latter is greater than 0.5.

Ulaby, F. T.↗

Relating the microwave backscattering coefficient to leaf area index

This paper examines the relationship between the microwave backscattering coefficient of a vegetation canopy, sigma (can, 0) and the canopy's leaf area index (LAI). The relationship is established through the development of one model for corn and sorghum and another for wheat. Both models are extensions of the cloud model of Attema and Ulaby (1978). Analysis of experimental data measured at 8.6, 13.0, 17.0, and 35.6 GHz indicates that most of the temporal variations of sigma (can, 0) can be accounted for through variations in green LAI alone, if the latter is greater than 0.5.

Ulaby, F. T.↗

Spectral estimators of absorbed photosynthetically active radiation in corn canopies

Most models of crop growth and yield require an estimate of canopy leaf area index (LAI) or absorption of radiation. Relationships between photosynthetically active radiation (PAR) absorbed by corn canopies and the spectral reflectance of the canopies were investigated. Reflectance factor data were acquired with a LANDSAT MSS band radiometer. From planting to silking, the three spectrally predicted vegetation indices examined were associated with more than 95% of the variability in absorbed PAR. The relationships developed between absorbed PAR and the three indices were evaluated with reflectance factor data acquired from corn canopies planted in 1979 through 1982. Seasonal cumulations of measured LAI and each of the three indices were associated with greater than 50% of the variation in final grain yields from the test years. Seasonal cumulations of daily absorbed PAR were associated with up to 73% of the variation in final grain yields. Absorbed PAR, cumulated through the growing season, is a better indicator of yield than cumulated leaf area index. Absorbed PAR may be estimated reliably from spectral reflectance data of crop canopies.

Gallo, K. P.↗