Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stratification”

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 199 records · Page 11

Sr-Nd-Pb isotope systematics of Australasian tektites: Implications for the nature and composition of target materials and possible volatile loss of Pb

The Australasian tektite (AAT) strewn field is the largest strewn field on the Earth with about ~10–30% coverage, both land and ocean, but a clearly identified source impact crater is absent despite the young age of AAT of ca. 790 ka. A genetic link between the Australasian tektites and their unequivocal parental materials is therefore largely impossible to establish. Nevertheless, the nature of the parental materials and the extent of volatilization can be constrained using the splash form tektites, carrying the chemical signatures of high-temperature processes, and the layered (so-called Muong Nong-type) tektites, which are less chemically homogenized and exceptionally abundant in the AAT field. New high-precision Sr, Nd and Pb isotopic measurements were obtained for a chemically and petrographically well-characterized suite of AAT, which included the Muong Nong-type (MN-AAT) with precisely known field locations in Laos and splash forms (SF-AAT) from different parts of the strewn field. Additionally, optically dark and light zones of the MN-AAT were also separately analyzed. Homogeneous ε Nd values from -11.8 to -11.2, combined with a narrow range of two-stage Nd model ages from 1.67 to 1.72 Ga for the entire AAT suite, point to a well-mixed source, in terms of REE, of the crustal segment from which the sedimentary material for tektites was ultimately derived. The Sr isotopic data largely overlap for SF-AAT and MN-AAT ( 87 Sr/ 86 Sr = 0.71636–0.72021) and indicate Paleozoic to Mesozoic sedimentary parentage. However, late Neogene to early Quaternary re-deposition and formation of a thick silt-sized sedimentary section with vertical stratification is required to comply with 10Be data. Lead isotope systematics documents at least three different components which can perhaps be represented by different mineral phases, such as feldspar, zircon, organic matter adsorbed on young sediments etc., sorted during fluvial transport and final deposition. In addition, the SF-AAT have systematically lower Pb contents than the MN-AAT, and generally show isotopically heavier Pb isotopic ratios. This is theoretically consistent with a preferential volatilization of lighter Pb isotopes during evaporation and considerably larger Pb loss from SF-AAT when compared to MN-AAT. Nevertheless, further experimental work would be necessary to unambiguously distinguish kinetic fractionation from source mixing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-site evaluation of stratified and balanced sampling of soil organic carbon stocks in agricultural fields

Estimating soil organic carbon (SOC) stocks in agricultural fields is essential for environmental and agronomic research, management, and policy. Stratified sampling is a classic strategy for estimating mean soil properties, and has recently been codified in SOC monitoring protocols. However, for the specific task of estimating the SOC stock of an agricultural field, concrete guidance is needed for which covariates to stratify on and how much stratification can improve estimation efficiency. It is also unknown how stratified sampling of SOC stocks compares to modern alternatives, notably doubly balanced sampling. To address these gaps, we collected high-density (average of 7 samples ha -1 ) and deep (average of 75 cm) measurements of SOC stocks at eight commercial fields under maize-soybean production in two US Midwestern states. We combined these measurements with a Bayesian geostatistical model to evaluate stratified and balanced sampling strategies that use a set of readily-available geographic, topographic, spectroscopic, and soil survey data. We examined the number of samples needed to achieve a given level of SOC stock estimation accuracy. While stratified sampling using these variables enables an average sample size reduction of 17% (95% CI, 11% to 23%) compared to simple random sampling, doubly balanced sampling is consistently more efficient, reducing sample sizes by 32% (95% CI, 25% to 37%). The data most important to these efficiency gains are a remotely-sensed SOC index, SSURGO estimates of SOC stocks, and the topographic wetness index. We conclude that in order to meet the urgent challenge of climate change, SOC stocks in agricultural fields could be more efficiently estimated by taking advantage of this readily-available data, especially with doubly balanced sampling.

54 ENVIRONMENTAL SCIENCES↗

Observations of multi-ion physics and kinetic effects in a surrogate to the solar CNO reactions

The ‘CNO process’ occurs in heavier stars with finite metallicity in which hydrogen burning is catalyzed in the presence of 12 C. These reactions are more strongly dependent on temperature than the pp cycle reactions, and thus the CNO cycle dominates only in massive stars. For these types of reactions to be studied at ICF facilities such as OMEGA, an implosion platform using heavier nuclei in the fuel and capable of creating ion temperatures on the order of at least 20 keV is required. A potential route to reach these conditions is to take advantage of kinetic effects in low-convergence shock-driven ‘exploding pusher’ implosions. In this experiment, shots were conducted at the OMEGA laser facility using the surrogate reaction 13 C + D. Its cross section is substantially higher than the actual astrophysical CNO reactions. The yield of this reaction in these implosions was much lower than expected. Finally, physical explanations are discussed, with significant species stratification the likely explanation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Logarithmic scaling of higher-order temperature moments in the atmospheric surface layer

A generalized logarithmic law for high-order moments of homogeneous passive scalars is proposed for turbulent boundary layers. This law is analogous to the generalized log law that has been proposed for high-order moments of the turbulent longitudinal velocity and is derived by combining the random sweeping decorrelation hypothesis with a spectral model informed by the attached eddy hypothesis. The proposed theory predicts that the high-order moments of passive scalar fluctuations within the inertial sublayer will vary logarithmically with wall-normal distance (), and is evaluated using high frequency time-series measurements of temperature and streamwise velocity fluctuations obtained in the first meter of the atmospheric surface layer (ASL) under near-neutral thermal stratification. The logarithmic dependence with within the inertial sublayer is observed in both the air temperature and velocity moments, with good agreement to the predictions from the proposed theory. Surprisingly, the proposed theory appears to be as, if not more, valid for transported passive scalars than for the longitudinal velocity.

54 ENVIRONMENTAL SCIENCES↗

Developing Science-based fueling protocols for 250-bar hydrogen tanks onboard hydrogen ferries: Experiments and modeling

Combined modeling and experimental studies are reported of the fueling of a large (28 kg capacity) 250-bar Type IV hydrogen tank of the type being deployed on early hydrogen ferries, such as the MV Sea Change. The primary goal was to determine how such tanks can be successfully fueled with hydrogen (state of charge greater than 97%) within 45 minutes without exceeding the 82 °C temperature limit for such tanks. The modeling studies show that a gas injector is needed to avoid thermal stratification during hydrogen fueling which can result in potential hot spots. Empirically, precooling of the hydrogen to 0 °C was found to be needed in some of the cases examined, as ambient conditions greatly affected the need for a precooling to achieve the 45-minute fill time desired by end users. The experimental results afforded a calibration of the engineering model SOFIL for these large 250-bar tanks, which now enables using SOFIL to predict volume-averaged hydrogen fueling temperatures to an accuracy of ±2.7°C for these tanks. The model can therefore be used to evaluate potential scenarios for development of a standardized fueling methodology for ferries utilizing large Type-IV tanks.

08 HYDROGEN↗

Interfacial antiferromagnetic phase induced two-step magnetization reversal in PbZr 0.52 Ti 0.48 O3/La 0.67 Sr 0.33 MnO 3 superlattices

Artificial multiferroic heterostructures have recently attracted much interests due to the demonstrated magnetoelectric coupling (MEC) and unique functionalities, promising a tantalizing perspective of novel applications in next-generation electronic, memory, sensor, and energy harvesting technologies. Herein, we report a two-step magnetization reversal in PbZr 0.52 Ti 0.48 O 3 /La 0.67 Sr 0.33 MnO 3 (PZT/LSMO) superlattices, which originates from the strongly entangled strain-, ferroelectric (FE)-polarization-, and exchange-dependent effects. Specifically, the preferential occupancy of the in-plane Mn d x 2 -y 2 orbitals is triggered via the collective effects of the large tensile strain and FE polarization, giving rise to an interfacial antiferromagnetic (AFM) layer with strong AFM anisotropy. The strong spin exchange coupling between the AFM layer and the adjacent ferromagnetic (FM) layer facilitates the magnetic stratification of the FM layer, leading to two coercivities, i.e., two-step magnetization reversal. Meanwhile, a sizeable exchange bias (EB) field is induced. The emerged two-step magnetization reversal concomitant with the pronounced EB phenomenon should be a signature of an enhanced MEC in PZT/LSMO superlattices. Finally, our results will stimulate further interests in multiferroic superlattices in applications of multiferroic-based devices.

36 MATERIALS SCIENCE↗

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the joint cumulative distribution function (CDF) of a vector-valued quantity of interest in problems with random input parameters and initial conditions. Our approach combines MLMC with stratified sampling of the input sample space by replacing standard Monte Carlo at each level with stratified Monte Carlo initialized with proportionally allocated samples. We show that the resulting stratified MLMC (sMLMC) algorithm is more efficient than its standard MLMC counterpart due to the additional variance reduction provided by the stratification of the random parameter's domain, especially at the coarsest levels. Additional computational cost savings are obtained by smoothing the indicator function with a Gaussian kernel, which proves to be an efficient and robust alternative to recently developed polynomial-based techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A critical review of operational strategies for the management of harmful algal blooms (HABs) in inland reservoirs

Occurrences of freshwater harmful algal blooms (HABs) are increasing on a global scale, largely in part due to increased nutrient input and changing climate patterns. While reservoir management strategies that can influence phytoplankton are known, there is no published guideline or protocol for the management of harmful algal blooms. There is a need to establish what factors are the predominant drivers of blooms, and how common reservoir management strategies specifically influence each factor. The following literature review seeks to establish the benefits and drawbacks of operational management strategies that currently exist. The main focus is altering hydrodynamic conditions (hypolimnetic withdrawals, surface flushing, pulsed inflow, artificial mixing), in order to induce environmental changes within the reservoir itself. Since excess nutrients are one of the biggest contributors to worsening bloom conditions, internal nutrient dynamics and reduction are also discussed. Additionally, we review the predominant seasonal factors (stratification, light, temperature, and wind) that affect likelihood of bloom occurrence and duration. In conclusion, the ultimate objective of this review is to increase understanding of the relationships between HAB drivers and reservoir operations in order to inform the development of data, modeling, and management strategies for the prevention and mitigation of blooms.

54 ENVIRONMENTAL SCIENCES↗

Identifying recharge processes into a vast “fossil” aquifer based on dynamic groundwater 81 Kr age evolution

Water in deep aquifers in arid regions is often considered to be "fossil" when modern recharge rates are negligible relative to the reservoir capacity. Over the past five decades, the Nubian Sandstone Aquifer (NSA) in the arid region of the Sinai Peninsula (Egypt) and the Negev Desert (Israel) has been considered to contain fossil water based on 14 C dating, which revealed 14 C ages of about 30 kyr over most of the aquifer. However, this relatively homogeneous age distribution contradicts the expected increase in groundwater age in the direction of decreasing piezometric head along the flow trajectories. Here, dating results with the longer-lived 81 Kr radioisotope (t 1/2 = 229 ± 11 kyr) are presented, highlighting a wide age range of 40 kyr to 630 kyr in the confined sections of the aquifer, all with very low 14 C activity (<1 pmC). Elevated 81Kr and 14 C activities were only observed within or close to the system's recharge areas. These findings support a new perception of groundwater replenishment during different epochs from the early mid- Pleistocene to the Holocene. By tracking the downstream age evolution, rejuvenation was identified in places where the confinement had been breached. At other locations, the existence of an older groundwater body contributing to the aquifer was detected by means of strongly depleted 81 Kr activity. High spatial heterogeneity in groundwater ages close to the discharge zone of the system is attributed to pronounced age stratification with depth. Calculated ages in the more isolated sections of the system were used to assess regional flow velocity, hydraulic conductivity, and their agreement with present recharge rates. We conclude that regional aquifers in which low 14 C activity suggests homogeneous ancient groundwater age should be reevaluated with 81Kr. With an effective age range beyond one million years, this may enable the restoration of recharge history and provide crucial information for the management of groundwater resources.

54 ENVIRONMENTAL SCIENCES↗

Creating an Advanced Sensor Network to calculate real-time, mass-weighted flue gas composition and air heater leakage of a coal-fired utility boiler under dynamic operating conditions

Utilization of renewable energy sources to minimize the environmental impact of energy production has changed the way utility boilers operate, requiring frequent load cycling between full load and partial loads as low as 30%. Dynamic operation of coal-fired utility boilers significantly reduces boiler efficiency when compared to steady state at full load. Data-driven plant optimization has shown success with coal-fired utility boilers under dynamic operating conditions. The purpose of this work was to create an Advanced Sensor Network (ASN) to provide more extensive real-time data to inform dynamic plant optimization of Net Unit Heat Rate (NUHR). The ASN consists of gas sampling grids in the convective pass of the boiler and downstream of the air heater. These sampling grids allow for quantification of spatial variation of flue gas within the boiler and calculation of mass-weighted composition of flue gas through the combination of composition, velocity, and temperature measurements. The comparison of O 2 between the inlet and outlet of the air heater is used to calculate air leakage in real time. Flue gas composition and air heater leakage are both important factors in boiler efficiency and NUHR. Further, the results of this work support the value of mass-weighted averages for determining flue gas composition accurately. The measurements from the ASN show increased composition stratification during dynamic operation, with an average standard deviation 38% higher than observed during steady-state operation. Air heater leakage was also observed to increase from 2.8% to 5.1% following a load change. Prior to the installation of the ASN, these data would not have been available for dynamic control. These real-time data will be leveraged to calculate and optimize for NUHR during dynamic operation in future work.

42 ENGINEERING↗

Detect and correct bias in multi-site neuroimaging datasets

The desire to train complex machine learning algorithms and to increase the statistical power in association studies drives neuroimaging research to use ever-larger datasets. The most obvious way to increase sample size is by pooling scans from independent studies. However, simple pooling is often ill-advised as selection, measurement, and confounding biases may creep in and yield spurious correlations. In this work, we combine 35,320 magnetic resonance images of the brain from 17 studies to examine bias in neuroimaging. In the first experiment, Name That Dataset, we provide empirical evidence for the presence of bias by showing that scans can be correctly assigned to their respective dataset with 71.5% accuracy. Given such evidence, we take a closer look at confounding bias, which is often viewed as the main shortcoming in observational studies. In practice, we neither know all potential confounders nor do we have data on them. Hence, we model confounders as unknown, latent variables. Kolmogorov complexity is then used to decide whether the confounded or the causal model provides the simplest factorization of the graphical model. Finally, we present methods for dataset harmonization and study their ability to remove bias in imaging features. In particular, we propose an extension of the recently introduced ComBat algorithm to control for global variation across image features, inspired by adjusting for unknown population stratification in genetics. Overall, our results demonstrate that harmonization can reduce dataset-specific information in image features. Further, confounding bias can be reduced and even turned into a causal relationship. However, harmonization also requires caution as it can easily remove relevant subject-specific information. Code is available at https://github.com/ai-med/Dataset-Bias.

42 ENGINEERING↗

Arctic ice-ocean interactions in an 8-to-2 kilometer resolution global model

In the last decades, the Arctic climate has changed dramatically, with the loss of multiyear sea ice one of the clearest consequences. These changes have occurred on relatively rapid timescales, and both accurate short-term Arctic prediction (e.g., 10 days to three months) and climate projection of future Arctic scenarios present ongoing challenges. Here we describe a representation of the Arctic ocean and sea ice in a ultrahigh resolution simulation in which the horizontal grid mesh reduces from 8 km at the equator to 2 km at the poles (UH8to2) for the years 2017-2020. We find the simulation reproduces observed distributions of seasonal sea-ice thickness and concentration realistically, although concentration is biased low in the spring and summer and low biases in thickness are found in the central and eastern basins in the fall. Volume, fresh water, and heat transports through key passages are realistic, lying within observationally determined ranges. Climatological comparisons reveal that the UH8to2 Atlantic Water is shallower, warmer, and saltier than the World Ocean Atlas 2018 climatology for 2005-2017 in the eastern basin. Our analysis suggests that these biases, combined with a lack of stratification in the upper 100 m of the simulated ocean, contribute to the winter biases in modeled sea ice thickness. This relationship between biases in the sea ice and ocean points to a potential positive feedback within the model, illuminating challenges for long term model predictive power in a changing Arctic climate.

58 GEOSCIENCES↗

Using electronic health record metadata to predict housing instability amongst veterans

Housing instability is considered a significant life stressor and preemptive screening should be applied to identify those at risk for homelessness as early as possible so that they can be targeted for specialized care. We developed models to classify patient outcomes for an established VA Homelessness Screening Clinical Reminder (HSCR), which identifies housing instability, in the two months prior to its administration. Logistic Regression and Random Forest models were fit to classify responses using the last 18 months of document activity. We measure concentration of risk across stratifications of predicted probability and observe an enriched likelihood of finding confirmed false negative responses from veterans with diagnosed housing instability. Positive responses were 34 times more likely to be detected within the top 1 % of patients predicted at risk than from those randomly selected. There is a 1 in 4 chance of detecting false negatives within the top 1 % of predicted risk. Machine learning methods can classify between episodes of housing instability using a data-driven approach that does not rely on variables curated from domain experts. This method has the potential to improve clinicians’ ability to identify veterans who are experiencing housing instability but are not captured by HSCR.

60 APPLIED LIFE SCIENCES↗

Identification of combustion mode under MILD conditions using Chemical Explosive Mode Analysis

Direct Numerical Simulations (DNS) data of Moderate or Intense Low-oxygen Dilution (MILD) combustion are analysed in this paper to identify the contributions of the autoignition and flame modes. This is performed using an extended Chemical Explosive Mode Analysis (CEMA) which accounts for diffusion effects allowing it to discriminate between deflagration and autoignition. This analysis indicates that in premixed MILD combustion conditions, the main combustion mode is ignition for all dilution and turbulence levels and for the two reactant temperature conditions considered. In non-premixed conditions, the preponderance of the ignition mode was observed to depend on the axial location and mixture fraction stratification. With a large mixture fraction lengthscale, ignition is more preponderant in the early part of the domain while the deflagrative mode increases further downstream. On the other hand, when the mixture fraction lengthscale is small, sequential autoignition is observed. Finally, the various combustion modes are observed to correlate strongly with mixture fraction where lean mixtures are more likely to autoignite while stoichiometric and rich mixtures are more likely to react as deflagrative structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Influence of pilot-fuel mixing on the spatio-temporal progression of two-stage autoignition of diesel-sprays in low-reactivity ambient fuel-air mixture

The spatial and temporal locations of autoignition for direct-injection compression-ignition engines depend on fuel chemistry, temperature, pressure, and mixing trajectories in the fuel jets. Dual-fuel systems can provide insight into both fuel-chemistry and physical effects by varying fuel reactivities and engine operating conditions. In this context, the spatial and temporal progression of two-stage autoignition of a diesel-fuel surrogate, n-heptane, in a lean-premixed charge of synthetic natural-gas (NG) and air is imaged in an optically accessible heavy-duty diesel engine. The lean-premixed charge of NG is prepared by fumigation upstream of the engine intake manifold. Optical diagnostics include high-speed (15kfps) cool-flame chemiluminescence-imaging as an indicator of low-temperature heat-release (LTHR) and OH* chemiluminescence-imaging as an indicator high-temperature heat-release (HTHR). NG prolongs the ignition delay of the pilot fuel and increases the combustion duration. Zero-dimensional chemical-kinetics simulations provide further understanding by replicating a Lagrangian perspective for mixtures evolving along streamlines originating either at the fuel nozzle or in the ambient gas, for which the pilot-fuel concentration is either decreasing or increasing, respectively. The zero-dimensional simulations predict that LTHR initiates most likely on the air streamlines before transitioning to HTHR, either on fuel-streamlines or on air-streamlines in regions of near-constant φ. Due to the relatively short pilot-fuel injection-durations, the transient increase in entrainment near the end of injection (entrainment wave) is important for quickly creating auto-ignitable mixtures. To achieve desired combustion characteristics, e.g., multiple ignition-kernels and favorable combustion phasing and location (e.g., for reducing wall heat-transfer or optimizing charge stratification), adjusting injection parameters could tailor mixing trajectories to offset changes in fuel ignition chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LES of HCCI combustion of iso-octane/air in a flat-piston rapid compression machine

Homogeneous Charge Compression Ignition (HCCI) engines promise better efficiency and cleaner emissions than conventional piston engines, but can be challenging to control. Rapid compression machines (RCM) provide a simplified configuration for investigating HCCI combustion behavior, which is necessary for effective control of engine ignition timing and peak pressures. Here, in this study, we assess the utility of large eddy simulations (LES) for predicting HCCI combustion in a 3-D configuration. To this end, LES with finite-rate chemistry employing a 99-species iso-octane/air mechanism of two RCM operating conditions are performed. The RCM configuration under consideration was designed by Strozzi et al. with a flat piston to introduce large amounts of thermal stratification representative of realistic HCCI engine conditions, through the generation of corner vortices. It is shown that the simulation provides reasonable agreement with temperature fluctuations (7% difference), as well as ignition delay in the short ignition case (1 ms difference), while the long ignition case (35 ms difference) highlights more substantial deficiencies that are still within the expected uncertainty from the employed chemical mechanism. Flame propagation modes predicted by LES agree with experimental observations: spontaneous ignition is seen in the short ignition case, while deflagration is more prominent in the long ignition case. Analysis of global and local quantities classify the short ignition case in a mixed ignition regime, and the long ignition case in the mild ignition regime. These results demonstrate the utility of FRC-LES for investigations of multimode combustion regimes of HCCI combustion in a 3-D configuration.

33 ADVANCED PROPULSION SYSTEMS↗

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗