Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “terrestrial models”

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 973 records · Page 54

Meteorite Material Model for Structural Properties

To assess the threat posed by an asteroid entering Earth's atmosphere, one must predict if, when, and how it fragments during entry. A comprehensive understanding of the asteroid material properties is needed to achieve this objective. At present, the meteorite material found on earth are the only objects from an entering asteroid that can be used as representative material and be tested inside a laboratory setting. Due to complex petrology, it is technically challenging and expensive to obtain reliable material properties by means of laboratory test for a family of meteorites. In order to circumvent this challenge, meteorite unit models are developed to determine the effective material properties including Youngs modulus, compressive and tensile strengths and Poissons ratio, that in turn would help deduce the properties of asteroids. The meteorite unit is a representative volume that accounts for diverse minerals, porosity, cracks and matrix composition. The Youngs Modulus and Poissons Ratio in the meteorite units are calculated by performing several hundreds of Monte-Carlo simulations by randomly distributing the various phases inside these units. Once these values are obtained, cracks are introduced in these meteorite units. The size, orientation and distribution of cracks are derived by extensive CT-scans and visual scans of various meteorites from the same family. Subsequently, simulations are performed to attain stress-strain relations, strength and effective modulus values in the presence of these cracks. The meteorite unit models are presented for H, L and LL ordinary chondrites, as well as for terrestrial basalt. In the case of the latter, data from the simulations is compared with experimental data to validate the methodology. These material models will be subsequently used in fragmentation modeling of full scale asteroids.

Meteorite↗

Quantifying the effect of late bombardment on terrestrial zircons

The first 500 million years of Earth history is thought to be a period of intense planetary bombardment, but the timing and flux of this meteorite bombardment is poorly understood. In particular, on the basis of an inferred lunar impact history, some workers have hypothesized a ~3.9 Ga terminal cataclysm (TC) in which there was marked increase in the impact flux affecting the Moon, the Earth and possibly other terrestrial planets. Minerals that survived this enigmatic period offer a way to test early planetary bombardment models as they may contain telltale micro- to nanoscale shock features. Here, we present results from a numerical modeling calculation that assesses the probability that a zircon residing in the crust would escape shock melting or shock deformation during a TC bombardment event. Even with conservative pressure estimates for zircon shock deformation and intermediate bombardment intensities, we find that only ~6% of ≥4.0 Ga crust would be expected to survive a 3.9 Ga cataclysm without experiencing either complete melting or zircon shock metamorphism. In this study, we couple this modeling with a search for shock effects in the oldest zircons from the Acasta Gneiss Complex, which would have been present in the Earth's crust during a putative 3.9 Ga TC. Spatially correlated electron and NanoSIMS ion microscopy of 4.02 Ga igneous zircons from Acasta reveals no evidence of ancient shock. These data, together with similar results from other Hadean zircon suites, confirm that a post-Hadean TC is unlikely to have occurred. We suggest that the dearth of pre-3.9 Ga terrestrial crust and zircons is instead best explained by endogenic processes related to the mechanisms of early crust formation. Our modeling allows us to evaluate bombardment scenarios from the terrestrial zircon record by applying probabilistic interpretations to zircon shock deformation data. This approach will be valuable for other planetary bodies, allowing broader conclusions to be drawn from geographically limited datasets.

58 GEOSCIENCES↗

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible.

Arnold, William↗

Next-Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models possible

Arnold, William R., Sr.↗

Next Generation Lightweight Mirror Modeling Software

The advances in manufacturing techniques for lightweight mirrors, such as EXELSIS deep core low temperature fusion, Corning's continued improvements in the Frit bonding process and the ability to cast large complex designs, combined with water-jet and conventional diamond machining of glasses and ceramics has created the need for more efficient means of generating finite element models of these structures. Traditional methods of assembling 400,000 + element models can take weeks of effort, severely limiting the range of possible optimization variables. This paper will introduce model generation software developed under NASA sponsorship for the design of both terrestrial and space based mirrors. The software deals with any current mirror manufacturing technique, single substrates, multiple arrays of substrates, as well as the ability to merge submodels into a single large model. The modeler generates both mirror and suspension system elements, suspensions can be created either for each individual petal or the whole mirror. A typical model generation of 250,000 nodes and 450,000 elements only takes 5-10 minutes, much of that time being variable input time. The program can create input decks for ANSYS, ABAQUS and NASTRAN. An archive/retrieval system permits creation of complete trade studies, varying cell size, depth, and petal size, suspension geometry with the ability to recall a particular set of parameters and make small or large changes with ease. The input decks created by the modeler are text files which can be modified by any editor, all the key shell thickness parameters are accessible and comments in deck identify which groups of elements are associated with these parameters. This again makes optimization easier. With ANSYS decks, the nodes representing support attachments are grouped into components; in ABAQUS these are SETS and in NASTRAN as GRIDPOINT SETS, this make integration of these models into large telescope or satellite models easier.

Arnold, William R., Sr.↗

The combination of satellite and topographic/isostatic potential models for mean anomaly determinations

A method is presented for the estimation of a global gravity anomaly field using the combination of satellite-derived potential coefficient models and the coefficients implied by the Airy-Heiskanen topographic/isostatic potential (Rummel et al., 1988) from topographic models with a 30-km depth of compensation. Gravity anomalies calculated with this method are compared with a terrestrial 1 x 1 degree anomaly file where the anomaly standard deviations were less than 10 mgals. Using the GEM T1 model (Marsh et al., 1988) to degree 36, the rms anomaly discrepency was + or - 19 mgals, while the rms values for the terrestrial anomalies was + or - 28 mgals.

Rapp, Richard H.↗

Effects of Gravity on the Double-Diffusive Convection during Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe

General 2-D and 3-D finite element model of non-dilute alloy solidification was used to simulate growth of HgCdTe in terrestrial and microgravity conditions. Parametric research was undertaken to investigate effects of gravity level, gravity vector orientation and growth velocity on the pattern of melt convection, shape of crystal/melt interface and radial thermal gradient. Verification of the model was undertaken by comparison with previously published results. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures obtained from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. For gravity levels higher then 10(exp -7) of terrestrial one it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by 50% by choosing proper orientation of the ampoule. The predicted interface shape is in agreement with one obtained experimentally by quenching.

Bune, Andris↗

Detailed gravity anomalies from Geos 3 satellite altimetry data

Detailed gravity anomalies are computed from a combination of Geos 3 satellite altimeter and terrestrial gravity data using least-squares principles. The mathematical model used is based on the Stokes' equation modified for a nonglobal solution. Using Geos 3 data in the calibration area, the effects of several anomaly parameter configurations and data densities/distributions on the anomalies and their accuracy estimates are studied. The accuracy estimates for 1 deg x 1 deg mean anomalies from low density altimetry data are of the order of 4 mgal. Comparison of these anomalies with the terrestrial data and also with Rapp's data derived using collocation techniques show rms differences of 7.2 and 4.9 mgal, respectively. Indications are that the anomaly accuracies can be improved to about 2 mgal with high density data. Estimation of 30 in. x 30 in. mean anomalies indicates accuracies of the order of 5 mgal. Proper verification of these results will be possible only when accurate ground truth data become available.

Gopalapillai, G. S.↗

E3SM‐GCAM: A Synchronously Coupled Human Component in the E3SM Earth System Model Enables Novel Human‐Earth Feedback Research

Abstract Modeling human‐environment feedbacks is critical for assessing the effectiveness of climate change mitigation and adaptation strategies under a changing climate. The Energy Exascale Earth System Model (E3SM) now includes a human component, with the Global Change Analysis Model (GCAM) at its core, that is synchronously coupled with the land and atmosphere components through the E3SM coupling software. Terrestrial productivity is passed from E3SM to GCAM to make climate‐responsive land use and CO 2 emission projections for the next 5‐year period, which are interpolated and passed to E3SM annually. Key variables affected by the incorporation of these feedbacks include land use/cover change, crop prices, terrestrial carbon, local surface temperature, and climate extremes. Regional differences are more pronounced than global differences because the effects are driven primarily by differences in land use. This novel system enables a new type of scenario development and provides a powerful modeling framework that facilitates the addition of other feedbacks between these models. This system has the potential to explore how human responses to climate change impacts in a variety of sectors, including heating/cooling energy demand, water management, and energy production, may alter emissions trajectories and Earth system changes. Plain Language Summary Earth system models help us understand how humans are changing the climate. Currently, these models do not include human systems, so predetermined greenhouse gas, aerosol, and land use change data are input to these models. These data do not reflect human responses to changes projected by Earth system models. We have added a human component to an Earth system model to represent human responses to environmental change and calculate corresponding greenhouse gas and land use change data instead of using predetermined data. Including this human component changes projections of land use, land carbon storage, and regional climate. Key Points We have incorporated a novel, advanceable human component in an Earth system model to simulate human‐Earth feedbacks Including terrestrial productivity feedbacks from the Earth to the human systems affects land change, crop prices, carbon, and climate Regional effects of including terrestrial productivity feedbacks are greater than global effects because land change is the main driver

E3SM↗

Reduction of structural impacts and distinction of photosynthetic pathways in a global estimation of GPP from space-borne solar-induced chlorophyll fluorescence

Quantifying global photosynthesis remains a challenge due to a lack of accurate remote sensing proxies. Solar-induced chlorophyll fluorescence (SIF) has been shown to be a good indicator of photosynthetic activity across various spatial scales. However, a global and spatially challenging estimate of terrestrial gross primary production (GPP) based on satellite SIF remains unresolved due to the confounding effects of species-specific physical and physiological traits and external factors, such as canopy structure or photosynthetic pathway (C3orC4). Here we analyze an ensemble of far-red SIF data from OCO-2 satellite and ground observations at multiple sites, using the spectral invariant theory to reduce the effects of canopy structure and to retrieve a structure-corrected total canopy SIF emission (SIF total). We find that the relationships between observed canopy-leavingSIF and ecosystem GPP vary significantly among biomes. In contrast, the relationships between SIF total and GPP converge around two unique models, one for C3and one for C4plants. We show that the two single empirical models can be used to globally scale satellite SIF observations to terrestrial GPP. We obtain an independent estimate of global terrestrial GPP of 129.56 ± 6.54 PgC/year for the 2015–2017 period, which is consistent with the state-of-the-art data- and process-oriented models. The new GPP product shows improved sensitivity to previously undetected ‘hotspots’ of productivity, being able to resolve the double-peak in GPP due to rotational cropping systems. We suggest that the direct scheme to estimate GPP presented here, which is based on satelliteSIF, may open up new possibilities to resolve the dynamics of global terrestrial GPP across space and time.

scale satellite SIF observations to terrestrial GP↗

Understanding Subseasonal Moisture Recycling from Extreme MCSs Using a Land–Atmosphere Coupled Water Tracer Model

Extreme precipitation events often produce severe flooding and pose significant threats to our society. However, their subseasonal impacts on the terrestrial and atmospheric systems are not well understood. Here, we use a new land–atmosphere coupled water tracer tool embedded in the Weather Research and Forecasting (WRF) Model to “tag” precipitation produced by a series of extreme mesoscale convective systems (MCSs) in May 2015 over the southern Great Plains to reveal their contributions to different terrestrial water storages and moisture recycling during May–July. During May, we find that precipitation from earlier MCSs can contribute to 20%–50% of total evapotranspiration (ET) in later May and contribute ∼10% of local precipitation for later MCSs. After May, the water “tagged” to preceding May MCS events has the most active contribution to moisture recycling during the first 5 days in June, but this contribution diminishes quickly afterward. Both periods indicate a quick turnover of tagged precipitation, which may be representative of the southern Great Plains region. The quick turnover is closely associated with the direct evaporation from the soil surface due to the predominant shrubland and grassland coverage over this region. Over some forested areas, however, the tracer-ET flux is dominated by transpiration that stably supplies a small fraction (<1%) of moisture to the atmosphere throughout June and July. In conclusion, the tracer-contributed plume (through evaporation and transpiration) in atmospheric moisture can extend 500–700 km downstream between May and July, highlighting the far-reaching contribution to moisture recycling and precipitation from the tagged MCS events in the southern Great Plains.

Atmosphere-land interaction↗

Volcanic Plume Heights on Mars: Limits of Validity for Convective Models

Previous studies have overestimated volcanic plume heights on Mars. In this work, we demonstrate that volcanic plume rise models, as currently formulated, have only limited validity in any environment. These limits are easily violated in the current Mars environment and may also be violated for terrestrial and early Mars conditions. We indicate some of the shortcomings of the model with emphasis on the limited applicability to current Mars conditions. Specifically, basic model assumptions are violated when (1) vertical velocities exceed the speed of sound, (2) radial expansion rates exceed the speed of sound, (3) radial expansion rates approach or exceed the vertical velocity, or (4) plume radius grossly exceeds plume height. All of these criteria are violated for the typical Mars example given here. Solutions imply that the convective rise, model is only valid to a height of approximately 10 kilometers. The reason for the model breakdown is hat the current Mars atmosphere is not of sufficient density to satisfy the conservation equations. It is likely that diffusion and other effects governed by higher-order differential equations are important within the first few kilometers of rise. When the same criteria are applied to eruptions into a higher-density early Mars atmosphere, we find that eruption rates higher than 1.4 x 10(exp 9) kilograms per second also violate model assumptions. This implies a maximum extent of approximately 65 kilometers for convective plumes on early Mars. The estimated plume heights for both current and early Mars are significantly lower than those previously predicted in the literature. Therefore, global-scale distribution of ash seems implausible.

Glaze, Lori S.↗

Use of Reference Frames for Interplanetary Navigation at JPL

Navigation of interplanetary spacecraft is typically based on range, Doppler, and differential interferometric measurements made by ground-based telescopes. Acquisition and interpretation of these observations requires accurate knowledge of the terrestrial reference frame and its orientation with respect to the celestial frame. Work is underway at JPL to reprocess historical VLBI and GPS data to improve realizations of the terrestrial and celestial frames. Improvements include minimal constraint alignment, improved tropospheric modeling, better orbit determination, and corrections for antenna phase center patterns.

International Celestial Reference Frame (ICRF)↗

Modeling and Analysis of Realistic Fire Scenarios in Spacecraft

An accidental fire inside a spacecraft is an unlikely, but very real emergency situation that can easily have dire consequences. While much has been learned over the past 25+ years of dedicated research on flame behavior in microgravity, a quantitative understanding of the initiation, spread, detection and extinguishment of a realistic fire aboard a spacecraft is lacking. Virtually all combustion experiments in microgravity have been small-scale, by necessity (hardware limitations in ground-based facilities and safety concerns in space-based facilities). Large-scale, realistic fire experiments are unlikely for the foreseeable future (unlike in terrestrial situations). Therefore, NASA will have to rely on scale modeling, extrapolation of small-scale experiments and detailed numerical modeling to provide the data necessary for vehicle and safety system design. This paper presents the results of parallel efforts to better model the initiation, spread, detection and extinguishment of fires aboard spacecraft. The first is a detailed numerical model using the freely available Fire Dynamics Simulator (FDS). FDS is a CFD code that numerically solves a large eddy simulation form of the Navier-Stokes equations. FDS provides a detailed treatment of the smoke and energy transport from a fire. The simulations provide a wealth of information, but are computationally intensive and not suitable for parametric studies where the detailed treatment of the mass and energy transport are unnecessary. The second path extends a model previously documented at ICES meetings that attempted to predict maximum survivable fires aboard space-craft. This one-dimensional model implies the heat and mass transfer as well as toxic species production from a fire. These simplifications result in a code that is faster and more suitable for parametric studies (having already been used to help in the hatch design of the Multi-Purpose Crew Vehicle, MPCV).

life support↗

Formation of Outer Planets: Overview

An overview of current theories of planetary formation, with emphasis on giant planets is presented. The most detailed models are based upon observation of our own Solar System and of young stars and their environments. Terrestrial planets are believe to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. According to the prevailing core instability model, giant planets begin their growth by the accumulation of small solid bodies, as do terrestrial planets. However, unlike terrestrial planets, the growing giant cores become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk disspates. The primary questions regarding the core instability model is whether planets with small cores can accrete gaseous enveloples within the lifetimes of gaseous protoplanetary disks. The main alternative giant planet formation model is the disk instability model, in which gaseous planets form directly via gravitational instabilities within protoplanetary disks. Formation of giant planets via gas instability has never been demonstrated for realistic disk conditions. Moreover, this model has difficulty explaining the supersolar abundances of heavy elements in Jupiter and Saturn, and it does not explain the orgin of planets like Uranus and Neptune.

Lissauer, Jack↗

Revisiting particle dry deposition and its role in radiative effect estimates

Wet and dry deposition remove aerosols from the atmosphere, and these processes control aerosol lifetime and thus impact climate and air quality. Dry deposition is a significant source of aerosol uncertainty in global chemical transport and climate models. Dry deposition parameterizations in most global models were developed when few particle deposition measurements were available. However, new measurement techniques have enabled more size-resolved particle flux observations. We combined literature measurements with data that we collected over a grassland in Oklahoma and a pine forest in Colorado to develop a dry deposition parameterization. We find that relative to observations, previous parameterizations overestimated deposition of the accumulation and Aitken mode particles, and underestimated in the coarse mode. These systematic differences in observed and modeled accumulation mode particle deposition velocities are as large as an order of magnitude over terrestrial ecosystems. As accumulation mode particles form most of the cloud condensation nuclei (CCN) that influence the indirect radiative effect, this model-measurement discrepancy in dry deposition alters modeled CCN and radiative forcing. We present a revised observationally driven parameterization for regional and global aerosol models. Using this revised dry deposition scheme in the Goddard Earth Observing System (GEOS)-Chem chemical transport model, we find that global surface accumulation-mode number concentrations increase by 62% and enhance the global combined anthropogenic and natural aerosol indirect effect by −0.63 W m −2 . Our observationally constrained approach should reduce the uncertainty of particle dry deposition in global chemical transport models.

54 ENVIRONMENTAL SCIENCES↗

Geospatial coherence of surface-atmosphere fluxes in the upper Great Lakes region

Surface-atmosphere fluxes are known to vary at multiple time scales, but uncertainty is high as to how fluxes change spatially within regions. With an increase in the number of eddy covariance towers, we are now able to examine the geospatial coherence of ecosystem fluxes, using time-series correlation. Eighteen sites from Michigan and Wisconsin were used in this study, ranging from 100 m to 600 km apart. Surface-atmosphere fluxes from a six-month period were used to quantify spatial coherence on a pair-wise basis. Using geospatial statistics, carbon and sensible heat (H) fluxes were found to be 95% correlated directly outside of their flux footprint and 56% correlated up to a distance of ~35 km. Latent (LE) and momentum (τ) fluxes were less correlated, 83% directly outside of their flux footprint and 40% at a distance of ~130 km albeit, at a much larger spatial distance than for the carbon and sensible heat fluxes. All fluxes showed strong spectral resonance at diel and seasonal timescales, with 1-, 2- and 3-month periods being common modes of variability among H, LE, and τ fluxes. Results based on Empirical Orthogonal Function show distinct transitions of net ecosystem exchange from fall to winter before photosynthesis or respiration while H and τ do not exhibit coherent trends. Here, this work demonstrates the potential of quantifying geospatial coherence of surface-atmosphere fluxes in the Midwestern United States, with the ability to predict fluxes beyond the spatial limit of a single flux tower footprint. Ultimately, expanding the flux measurements to larger scales would allow better spatial scaling of terrestrial surface-atmosphere fluxes between tower footprint and modeling or remote sensing scales.

54 ENVIRONMENTAL SCIENCES↗

Salinity exposure affects lower-canopy specific leaf area of upland trees in a coastal deciduous forest

Sea level rise and increasing storm surges are likely to affect the canopy physiology, ecology, and structure of coastal forests, even well in advance of tree mortality. Laboratory and greenhouse studies have documented that saltwater exposure can trigger changes in leaf-level physiology and morphology, but few in situ studies have examined how tree-specific leaf area (SLA), the ratio of leaf area to mass and a crucial trait and model parameter, is affected by saline soils. We conducted an observational study of SLA in a mid-Atlantic (USA) coastal deciduous forest, taking advantage of a natural gradient in salinity along a tidal creek. Measured SLA of the 239 trees and seven species sampled ranged from Carya glabra (N = 6 trees, mean SLA = 277.9±36.3cm 2 /g) to Fagus grandifolia (N=60, 321.9±62.9cm 2 g); as expected, trees species and canopy position (sun versus shade) significantly affected SLA. For trees (N=100) directly exposed to the tidal creek, salinity was highly significant after accounting for species (P<0.001), with trees in the lower reaches of the creek having lower SLA. Leaf area index (LAI), computed from SLA and litter traps, ranged from 4.8 to 15.8 and was inversely related to salinity exposure; the spatial variability in leaf litter production contributed much more to LAI uncertainty than did SLA variability. These in situ results are correlative but consistent with the hypothesis, based on previous greenhouse studies, that the stress of chronic salinity exposure changes species’ leaf morphology. Our findings are useful for understanding the growing effects of saltwater intrusion into upland forests, as well as parameterizing and testing ecosystem-scale models simulating forest stressors and disturbances at the terrestrial-aquatic interface.

54 ENVIRONMENTAL SCIENCES↗