Engineering Papers⌕ Search

Engineering topics

Katul, Gabriel

Publications and source records attributed to Katul, Gabriel.

33 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The root-zone soil moisture spectrum in a mediterranean ecosystem

Storage of water within soil pores of the root zone introduce memory effects in the dynamics of soil moisture that are considerably longer than the integral timescale of many atmospheric processes. Thus, hydro-climatic states can be “sustained” through land-surface heat and water vapor fluxes primarily because they can “feed off” on this long-term soil moisture memory. Root-zone soil moisture memory is only but one feature characterizing the spectrum of soil moisture dynamics, which is analyzed here using a combination of long-term measurements and models. In particular, the spectrum of root-zone soil moisture content in a Mediterranean ecosystem is examined using 14-years of half-hourly measurements. A distinguishing hydro-climatic feature in such ecosystems is that sources (mainly rainfall) and sinks (mainly evapotranspiration) of soil moisture are roughly out of phase with each other. For over 4 decades of time scales and 7 decades of energy, the canonical shape of the measured soil moisture spectrum is shown to be approximately Lorentzian determined by the soil moisture variance and its memory but with two exceptions: the occurrences of a peak at diurnal-to-daily time scales and a weaker peak at near annual time scales. Model calculations and spectral analysis demonstrate that diurnal and seasonal variations in hydroclimate forcing responsible for variability in evapotranspiration had minor impact on the normalized shape of the soil moisture spectrum. However, their impact was captured by adjustments in the temporal variance. These findings indicate that precipitation and not evapotranspiration variability dominates the multi-scaling properties of soil moisture variability consistent with prior climate model simulations. Furthermore, the soil moisture memory inferred by the annual peak of soil moisture (340 d) is consistent with climate model simulations, while the memory evaluated from the loss function of a linearized mass balance approach leads to a smaller value (50 d), highlighting the effect of weak non-stationarity on soil moisture variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Taylor dispersion in osmotically driven laminar flows in phloem

Sucrose is among the main products of photosynthesis that are deemed necessary for plant growth and survival. It is produced in the mesophyll cells of leaves and translocated to different parts of the plant through the phloem. Progress in understanding this transport process remains fraught with experimental difficulties, thereby prompting interest in theoretical approaches and laboratory studies. The Münch pressure and mass flow model is one of the accepted hypotheses describing the physics of sucrose transport in the phloem. It is based on osmosis creating an energy potential difference between the source and the sink. The flow responding to this energy potential is assumed laminar and described by the Hagen–Poiseuille equation. This study revisits such osmotically driven flows in tubes with membrane walls by including the effects of Taylor dispersion on mass transport. This effect has been overlooked in phloem flow studies. Taylor dispersion can increase the effective transport of solutes by increasing the apparent diffusion coefficient. It is shown that, in addition to the conventional diffusive correction derived for impermeable tube walls, a new adjustment to the mean advective terms arises because of osmotic effects. Because the molecular Schmidt number is very large for sucrose in water, the sucrose front speed and travel times have a direct dependence on the Péclet number for different ranges of the Münch number. Fially, this study establishes upper limits on expected Taylor dispersion enhancement of sucrose transport.

42 ENGINEERING↗
DOE OSTI2019

A Dynamic Optimality Principle for Water Use Strategies Explains Isohydric to Anisohydric Plant Responses to Drought

Optimality principles that underlie models of stomatal kinetics require identifying and formulating the gain and the costs involved in opening stomata. While the gain has been linked to larger carbon acquisition, there is still a debate as to the costs that limit stomatal opening. This work presents an Euler-Lagrange framework that accommodates water use strategy and various costs through the formulation of constraints. The reduction in plant hydraulic conductance due to cavitation is added as a new constraint above and beyond the soil hydrological balance and is analyzed for three different types of whole-plant vulnerability curves. Model results show that differences in vulnerability curves alone lead to relatively iso- and aniso-hydric stomatal behavior. Moreover, this framework explains how the presence of competition (biotic or abiotic) for water alters stomatal response to declining soil water content. This contribution corroborates previous research that predicts that a plant's environment (e.g., competition, soil processes) significantly affects its response to drought and supplies the required mathematical machinery to represent this complexity. The method adopted here disentangles cause and effect of the opening and closure of stomata and complements recent mechanistic models of stomatal response to drought.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Vertical characterization of highly oxygenated molecules (HOMs) below and above a boreal forest canopy

While the role of highly oxygenated molecules (HOMs) in new particle formation (NPF) and secondary organic aerosol (SOA) formation is not in dispute, the interplay between HOM chemistry and atmospheric conditions continues to draw significant research attention. During the Influence of Biosphere-Atmosphere Interactions on the Reactive Nitrogen budget (IBAIRN) campaign in September 2016, profile measurements of neutral HOMs below and above the forest canopy were performed for the first time at the boreal forest SMEAR II station. The HOM concentrations and composition distributions below and above the canopy were similar during daytime, supporting a well-mixed boundary layer approximation. However, much lower nighttime HOM concentrations were frequently observed at ground level, which was likely due to the formation of a shallow decoupled layer below the canopy. Near the ground HOMs were influenced by the changes in the precursors and oxidants and enhancement of the loss on surfaces in this layer, while the HOMs above the canopy top were not significantly affected. Our findings clearly illustrate that near-ground HOM measurements conducted under stably stratified conditions at this site might only be representative of a small fraction of the entire nocturnal boundary layer. This could, in turn, influence the growth of newly formed particles and SOA formation below the canopy where the large majority of measurements are typically conducted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A Structure Function Model Recovers the Many Formulations for Air-Water Gas Transfer Velocity: Air-Water Gas Transfer

Two ideas regarding the structure of turbulence near a clear air-water interface are used to derive a waterside gas transfer velocity kL for sparingly and slightly soluble gases. The first is that k L is proportional to the turnover velocity described by the vertical velocity structure function D ww (r), where r is separation distance between two points. The second is that the scalar exchange between the air-water interface and the waterside turbulence can be suitably described by a length scale proportional to the Batchelor scale lB=ηSc -1/2 , where Sc is the molecular Schmidt number and η is the Kolmogorov microscale defining the smallest scale of turbulent eddies impacted by fluid viscosity. Using an approximate solution to the von Kármán-Howarth equation predicting D ww (r) in the inertial and viscous regimes, prior formulations for k L are recovered including (i) k L =$\sqrt{2/15}$Sc -1/2 $v$, v K is the Kolmogorov velocity defined by the Reynolds number v K η/ν = 1 and ν is the kinematic viscosity of water; (ii) surface divergence formulations; (iii) k L ∝ SC -1/2 u * , where u * is the waterside friction velocity; (iv) k L ∝ Sc -1/2 $\sqrt{gv/u_*}$ for Keulegan numbers exceeding a threshold needed for long-wave generation, where the proportionality constant varies with wave age, g is the gravitational acceleration; and (v) k L = $\sqrt{2/15}$Sc -1/2 (vgβ o q o ) 1/4 in free convection, where qo is the surface heat flux and β o is the thermal expansion of water. The work demonstrates that the aforementioned k L formulations can be recovered from a single structure function model derived for locally homogeneous and isotropic turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Bridging Land-Surface Fluxes and Aerosol Concentrations to Triggering Convective Rainfall

The goal of this project is to investigate yields and properties of aerosol resulting from the photooxidation of plant-emitted hydrocarbons in natural environments of the Amazonia and their potential influences on cloud formation. A combination of intensive and comprehensive field observations and numerical models at multiple scales is being used to investigate the physical and chemical processes relevant to the formation of cloud condensation nuclei (CCN). The final report resulted in 40 peer reviewed manuscripts (from the Duke University co-PI) partly completed in collaboration with the PI and other co-PIs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Competition for light and water in a coupled soil-plant system

Here, it is generally accepted that resource availability shapes the structure and function of many ecosystems. Within the soil-plant-atmosphere (SPA) system, resource availability fluctuates in space and time whereas access to resources by individuals is further impacted by plant-to-plant competition. Likewise, transport and transformation of resources within an individual plant is governed by numerous interacting biotic and abiotic processes. The work here explores the co-limitations on water losses and carbon uptake within the SPA arising from fluctuating resource availability and competition. In particular, the goal is to unfold the interplay between plant access and competition for water and light, as well as the impact of transport/redistribution processes on leaf-level carbon assimilation and water fluxes within forest stands. A framework is proposed that couples a three-dimensional representation of soil-root exchanges with a one-dimensional description of stem water flow and storage, canopy photosynthesis, and transpiration. The model links soil moisture redistribution, root water uptake, xylem water flow and storage, leaf potential and stomatal conductance as driven by supply and demand for water and carbon. The model is then used to investigate plant drought resilience of overstory-understory trees simultaneously competing for water and light. Simulation results reveal that understory-overstory interactions increase ecosystem resilience to drought (i.e. stand-level carbon assimilation rates and water fluxes can be sustained at lower root-zone soil water potentials). This resilience enhancement originates from reduced transpiration (due to shading) and hydraulic redistribution in soil supporting photosynthesis over prolonged periods of drought. In particular, the presence of different rooting systems generates localized hydraulic redistribution fluxes that sustain understory transpiration through overstory-understory interactions. Such complex SPA dynamics cannot be properly summarized by equivalent ecosystem-scale Resistor-Capacitor (RC) rep- resentation. However our results show that, with proper averaging across water flow paths, RC models can provide reasonable estimates of stand-level water and carbon fluxes during inter-storm periods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A Kolmogorov-Brutsaert structure function model for evaporation into a turbulent atmosphere

In 1965, Brutsaert proposed a model that predicted mean evaporation rate E¯ from rough surfaces to scale with the 3/4 power law of the friction velocity ( u* ) and the square–root of molecular diffusivity (D m ) for water vapor. Furthermore, in arriving at these results, a number of assumptions were made regarding the surface renewal rate describing the contact durations between eddies and the evaporating surface, the diffusional mass process from the surface into eddies, and the cascade of turbulent kinetic energy sustaining the eddy renewal process itself.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A reduced order model to analytically infer atmospheric CO 2 concentration from stomatal and climate data

To address questions related to the acceleration or deceleration of the global hydrological cycle or links between the carbon and water cycles over land, reliable data for past climatic conditions based on proxies are required. In particular, the reconstruction of palaeoatmospheric CO 2 content (C a ) is needed to assist the separation of natural from anthropogenic C a variability and to explore phase relations between C a and air temperature T a time series. Both T a and C a are needed to fingerprint anthropogenic signatures in vapor pressure deficit, a major driver used to explain acceleration or deceleration phases in the global hydrological cycle. Current approaches to C a reconstruction rely on a robust inverse correlation between measured stomatal density in leaves (ν) of many plant taxa and C a . There are two methods that exploit this correlation: The first uses calibration curves obtained from extant species assumed to represent the fossil taxa, thereby restricting the suitable taxa to those existing today. The second is a hybrid eco-hydrological/physiological approach that determines C a with the aid of systems of equations based on quasi-instantaneous leaf-gas exchange theories and fossil stomatal data collected along with other measured leaf anatomical traits and parameters. In this contribution, a reduced order model (ROM) is proposed that derives C a from a single equation incorporating the aforementioned stomatal data, basic climate (e.g. temperature), estimated biochemical parameters of assimilation and isotope data. The usage of the ROM is then illustrated by applying it to isotopic and anatomical measurements from three extant species. The ROM derivation is based on a balance between the biochemical demand and atmospheric supply of CO 2 that leads to an explicit expression linking stomatal conductance to internal CO 2 concentration (C i ) and C a . The resulting expression of stomatal conductance from the carbon economy of the leaf is then equated to another expression derived from water vapor gas diffusion that includes anatomical traits. When combined with isotopic measurements for long-term C i /C a , C a can be analytically determined and is interpreted as the time-averaged C a that existed over the life-span of the leaf. Key advantages of the proposed ROM are: 1) the usage of isotopic data provides constraints on the reconstructed atmospheric CO 2 concentration from ν, 2) the analytical form of this approach permits direct links between parameter uncertainties and reconstructed C a , and 3) the time-scale mismatch between the application of instantaneous leaf-gas exchange expressions constrained with longer-term isotopic data is reconciled through averaging rules and sensitivity analysis. The latter point was rarely considered in prior reconstruction studies that combined models of leaf-gas exchange and isotopic data to reconstruct C a from ν. The proposed ROM is not without its limitations given the need to a priori assume a parameter related to the control on photosynthetic rate. Our work here further explores immanent constraints for the aforementioned photosynthetic parameter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Multiple mechanisms generate a universal scaling with dissipation for the air-water gas transfer velocity

A large corpus of field and laboratory experiments support the finding that the water side transfer velocity kL of sparingly soluble gases near air-water interfaces scales as k L ~(νε) 1/4 , where ν is the kinematic water viscosity and ε is the mean turbulent kinetic energy dissipation rate. Originally predicted from surface renewal theory, this scaling appears to hold for marine and coastal systems and across many environmental conditions. It is shown that multiple approaches to representing the effects of turbulence on kL lead to this expression when the Kolmogorov microscale is assumed to be the most efficient transporting eddy near the interface. The approaches considered range from simplified surface renewal schemes with distinct models for renewal durations, scaling and dimensional considerations, and a new structure function approach derived using analogies between scalar and momentum transfer. The work offers a new perspective as to why the aforementioned 1/4 scaling is robust.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Linking Meteorology, Turbulence, and Air Chemistry in the Amazon Rain Forest

We describe the salient features of a field study whose goals are to quantify the vertical distribution of plant-emitted hydrocarbons and their contribution to aerosol and cloud condensation nuclei production above a central Amazonian rain forest. Using observing systems deployed on a 50-m meteorological tower, complemented with tethered balloon deployments, the vertical distribution of hydrocarbons and aerosols was determined under different boundary layer thermodynamic states. The rain forest emits sufficient reactive hydrocarbons, such as isoprene and monoterpenes, to provide precursors of secondary organic aerosols and cloud condensation nuclei. Mesoscale convective systems transport ozone from the middle troposphere, enriching the atmospheric boundary layer as well as the forest canopy and surface layer. Through multiple chemical transformations, the ozone-enriched atmospheric surface layer can oxidize rain forest–emitted hydrocarbons. One conclusion derived from the field studies is that the rain forest produces the necessary chemical species and in sufficient amounts to undergo oxidation and generate aerosols that subsequently activate into cloud condensation nuclei.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Soil–plant–atmosphere conditions regulating convective cloud formation above southeastern US pine plantations

Abstract Loblolly pine trees ( Pinus taeda L.) occupy more than 20% of the forested area in the southern United States, represent more than 50% of the standing pine volume in this region, and remove from the atmosphere about 500 g C m per year through net ecosystem exchange. Hence, their significance as a major regional carbon sink can hardly be disputed. What is disputed is whether the proliferation of young plantations replacing old forest in the southern United States will alter key aspects of the hydrologic cycle, including convective rainfall, which is the focus of the present work. Ecosystem fluxes of sensible ( ) and latent heat (LE) and large‐scale, slowly evolving free atmospheric temperature and water vapor content are known to be first‐order controls on the formation of convective clouds in the atmospheric boundary layer. These controlling processes are here described by a zero‐order analytical model aimed at assessing how plantations of different ages may regulate the persistence and transition of the atmospheric system between cloudy and cloudless conditions. Using the analytical model together with field observations, the roles of ecosystem and LE on convective cloud formation are explored relative to the entrainment of heat and moisture from the free atmosphere. Our results demonstrate that cloudy–cloudless regimes at the land surface are regulated by a nonlinear relation between the Bowen ratio and root‐zone soil water content, suggesting that young/mature pines ecosystems have the ability to recirculate available water (through rainfall predisposition mechanisms). Such nonlinearity was not detected in a much older pine stand, suggesting a higher tolerance to drought but a limited control on boundary layer dynamics. These results enable the generation of hypotheses about the impacts on convective cloud formation driven by afforestation/deforestation and groundwater depletion projected to increase following increased human population in the southeastern United States.

Manoli, Gabriele↗