Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atmospheric Science”

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 937 records · Page 52

Grassland carbon-water cycling is minimally impacted by a photovoltaic array

Agrivoltaic systems, whereby photovoltaic arrays are co-located with crop or forage production, can alleviate the tension between expanding solar development and loss of agricultural land. However, the ecological ramifications of these arrays are poorly known. We used field measurements and a plant hydraulic model to quantify carbon-water cycling in a semi-arid C 3 perennial grassland growing beneath a single-axis tracking solar array in Colorado, USA. Although the agrivoltaic array reduced light availability by 38%, net photosynthesis and aboveground net primary productivity were reduced by only 6–7% while evapotranspiration decreased by 1.3%. The minimal changes in carbon-water cycling occurred largely because plant photosynthetic traits underneath the panels changed to take advantage of the dynamic shading environment. Our results indicate that agrivoltaic systems can serve as a scalable way to expand solar energy production while maintaining ecosystem function in managed grasslands, especially in climates where water is more limiting than light.

14 SOLAR ENERGY↗

Linking repeat lidar with Landsat products for large scale quantification of fire-induced permafrost thaw settlement in interior Alaska

The permafrost–fire–climate system has been a hotspot in research for decades under a warming climate scenario. Surface vegetation plays a dominant role in protecting permafrost from summer warmth, thus, any alteration of vegetation structure, particularly following severe wildfires, can cause dramatic top–down thaw. A challenge in understanding this is to quantify fire-induced thaw settlement at large scales (>1000 km 2 ). In this study, we explored the potential of using Landsat products for a large-scale estimation of fire-induced thaw settlement across a well-studied area representative of ice-rich lowland permafrost in interior Alaska. Six large fires have affected ~1250 km 2 of the area since 2000. We first identified the linkage of fires, burn severity, and land cover response, and then developed an object-based machine learning ensemble approach to estimate fire-induced thaw settlement by relating airborne repeat lidar data to Landsat products. The model delineated thaw settlement patterns across the six fire scars and explained ~65% of the variance in lidar-detected elevation change. Our results indicate a combined application of airborne repeat lidar and Landsat products is a valuable tool for large scale quantification of fire-induced thaw settlement.

54 ENVIRONMENTAL SCIENCES↗

Increased salinity decreases annual gross primary productivity at a Northern California brackish tidal marsh

Tidal marshes sequester 11.4–87.0 Tg C yr –1 globally, but climate change impacts can threaten the carbon capture potential of these ecosystems. Tidal marshes occur across a wide range of salinity, with brackish marshes (0.5–18 ppt (parts per thousand)) dominating global tidal marsh extents. A diverse mix of freshwater- and saltwater-tolerant plant and microbial communities has led researchers to predict that carbon cycling in brackish wetlands may be less sensitive to changes in salinity than fresh- or saltwater wetlands. Rush Ranch, a well-monitored brackish tidal wetland of the San Francisco Bay National Estuarine Research Reserve, experiences highly variable annual salinity regimes. Within a five-year period (2014–2018), Rush Ranch experienced particularly extreme drought-induced salinization during the 2014 and 2015 growing seasons. During drought years, tidal channel salinity rose from a 15 year baseline of 4.7 ppt to growing season peaks of 10.3 ppt and 12.5 ppt. Continuous eddy covariance data from 2014 to 2018 demonstrate that during drought summers, gross primary productivity (GPP) decreased by 24%, whereas ecosystem respiration remained similar among all five years. Stepwise linear regression revealed that salinity, not air temperature or tidal height, was the dominant driver of annual GPP. A random forest model trained to predict GPP based on environmental data from low salinity years (i.e. naive to salinization) significantly over predicted GPP in drought years. When growing season salinities were doubled, annual estimates of net ecosystem exchange of CO 2 decreased by up to 30%. These results provide ecosystem-scale evidence that increased salinity influences CO 2 fluxes dominantly through reductions in GPP. This relationship provides a starting point for incorporating the effect of changes in salinity in wetland carbon models, which could improve wetland carbon forecasting and management for climate resilience.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗

Analytical Investigation of the Role of Lateral Mixing in the Evolution of Nonprecipitating Cumulus. Part I: Developing Clouds

Evolution of nonprecipitating cumulus clouds (Cu) at the developing stage under the influence of lateral entrainment and mixing is studied analytically using a minimalistic analytical model. We present a model of an ascending cloud volume (a model of developing Cu) whose structure is determined by the processes of droplet diffusion growth/evaporation and entrainment mixing in the horizontal direction. Spatial and time changes of liquid water content, the adiabatic fraction, droplet concentration, and the mean volume droplet radius are calculated. It is shown that the existence of a nondiluted core in a growing cumulus cloud significantly depends on the cloud width and vertical velocity. While at the updraft velocity of 2 m s −1 the core of a 400-m-wide cloud becomes diluted at distances of a few hundred meters above cloud base, the core of a cloud of 1000-m width remains nondiluted at distances up to 1500 m above cloud base. The explanation of this result is simple: the increase in cloud width and the decrease in the updraft velocity increase the time during which the cloud is diluted due to mixing. Since lateral mixing synchronously decreases both the cloud water content and droplet concentration, the variation of the mean volume droplet radius is low inside the cloud. The approximate quantitative condition for cloud formation in updraft is derived. It is shown that a cloud can arise when its vertical velocity exceeds a critical value. To produce clouds, narrow turbulent plumes should ascend at higher velocity as compared to wider plumes. High humidity of the environment air is favorable for formation of clouds from plumes. The comparison of the obtained results with previously published observational data indicates a reasonable agreement. The results can be useful for parameterization purposes.

54 ENVIRONMENTAL SCIENCES↗

Toward a Theory of the Evolution of Drop Morphology and Splintering by Freezing

Abstract The drop freezing process is described by a phase-field model. Two cases are considered: when the freezing is triggered by central nucleation and when nucleation occurs on the drop surface. Depending on the environmental temperature and drop size, different morphological structures develop. Detailed dendritic growth was simulated at the first stage of drop freezing. Independent of the nucleation location, a decrease in temperature within the range from ~−5° to −25°C led to an increase in the number of dendrites and a decrease in their width and the interdendritic space. At temperatures lower than about −25°C, a planar front developed following surface nucleation, while dendrites formed a granular-like structure with small interdendritic distances following bulk nucleation. An ice shell grew in at the same time (but slower) as dendrites following surface nucleation, while it started forming once the dendrites have reached the drop surface in the case of central nucleation. The formed ice morphology at the first freezing stage predefined the splintering probability. We assume that stresses needed to break the ice shell arose from freezing of the water in the interdendritic spaces. Under this assumption, the number of possible splinters/fragments was proportional to the number of dendrites, and the maximum rate of splintering/fragmentation occurred within a temperature range of about −10° to −20°C, is in agreement with available laboratory and in situ measurements. At temperatures < −25°C, freezing did not lead to the formation of significant stresses, making splintering unlikely. The number of dendrites increased with drop size, causing a corresponding increase in the number of splinters. Examples of morphology that favors drop cracking are presented, and the duration of the freezing stages is evaluated. Sensitivity of the freezing process to the surface fluxes is discussed.

54 ENVIRONMENTAL SCIENCES↗

Effects of Vegetation and Topography on the Boundary Layer Structure above the Amazon Forest

Observational data from two field campaigns in the Amazon forest were used to study the vertical structure of turbulence above the forest. The analysis was performed using the reduced turbulent kinetic energy (TKE) budget and its associated two-dimensional phase space. Results revealed the existence of two regions within the roughness sublayer in which the TKE budget cannot be explained by the canonical flat-terrain TKE budgets in the canopy roughness sublayer or in the lower portion of the convective ABL. Data analysis also suggested that deviations from horizontal homogeneity have a large contribution to the TKE budget. Results from LES of a model canopy over idealized topography presented similar features, leading to the conclusion that flow distortions caused by topography are responsible for the observed features in the TKE budget. These results support the conclusion that the boundary layer above the Amazon forest is strongly impacted by the gentle topography underneath.

54 ENVIRONMENTAL SCIENCES↗

A Formula for the Maximum Vertical Velocity in Supercell Updrafts

In supercell environments, previous authors have shown strong connections between the vertical wind shear magnitude, updraft width, and entrainment. Based on these results, it is hypothesized that the influences of entrainment-driven dilution on buoyancy and maximum updraft vertical velocity w in supercell environments are a predictable function of the vertical wind shear profile. It is also hypothesized that the influences of pressure perturbation forces on maximum updraft w are small because of a nearly complete offset between upward dynamic pressure forces and downward buoyant pressure forces. To address these hypotheses, we derive a formula for the maximum updraft w that incorporates the effects of entrainment-driven dilution on buoyancy but neglects pressure gradient forces. Solutions to this formula are compared with output from previous numerical simulations. This formula substantially improves predictions of maximum updraft w over past CAPE-derived formulas for maximum updraft w, which supports the first hypothesis. Furthermore, integrated vertical accelerations along trajectories show substantial offsets between dynamic and buoyant pressure forces, supporting the second hypothesis. It is argued that the new formula should be used in addition to CAPE-derived measures for w in forecast and research applications when accurate diagnosis of updraft speed is required.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evaluating the Conservation of Energy Variables in Simulations of Deep Moist Convection

Abstract It is often assumed in parcel theory calculations, numerical models, and cumulus parameterizations that moist static energy (MSE) is adiabatically conserved. However, the adiabatic conservation of MSE is only approximate because of the assumption of hydrostatic balance. Two alternative variables are evaluated here: MSE − IB and MSE + KE, wherein IB is the path integral of buoyancy ( B ) and KE is kinetic energy. Both of these variables relax the hydrostatic assumption and are more precisely conserved than MSE. This article quantifies the errors that result from assuming that the aforementioned variables are conserved in large-eddy simulations (LES) of both disorganized and organized deep convection. Results show that both MSE − IB and MSE + KE better predict quantities along trajectories than MSE alone. MSE − IB is better conserved in isolated deep convection, whereas MSE − IB and MSE + KE perform comparably in squall-line simulations. These results are explained by differences between the pressure perturbation behavior of squall lines and isolated convection. Errors in updraft B diagnoses are universally minimized when MSE − IB is assumed to be adiabatically conserved, but only when moisture dependencies of heat capacity and temperature dependency of latent heating are accounted for. When less accurate latent heat and heat capacity formulae were used, MSE − IB yielded poorer B predictions than MSE due to compensating errors. Our results suggest that various applications would benefit from using either MSE − IB or MSE + KE instead of MSE with properly formulated heat capacities and latent heats.

54 ENVIRONMENTAL SCIENCES↗

Revisiting the Hail Radar Reflectivity–Kinetic Energy Flux Relation by Combining T-Matrix and Discrete Dipole Approximation Calculations to Size Distribution Observations

The retrieval of hail kinetic energy with weather radars or its simulation in numerical models is challenging because of the shape complexity and variable density of hailstones. We combine 3D scans of individual hailstones with measurements of the particle size distributions (PSD) and T-matrix calculations to understand how hail reflectivity Z changes when approximating hailstones as spheroids, as compared to the realistic shapes obtained by 3D scanning technology. Additionally, recent terminal velocity relations are used to compare Z to the hail kinetic energy flux E ˙ . We parameterize the hail backscattering cross sections at L, S, C, and X bands as a function of size between 0.5 and 5.0 cm, matching the range of the observed PSDs. The scattering calculations use the T-matrix method for size parameters below 1.0 and the discrete dipole approximation (DDA) method otherwise. The DDA calculations are done for 48 digital models of realistic hailstones of sizes between 1 and 5 cm. The DDA cross sections are calculated for multiple orientations and averaged assuming a fully random orientation distribution to provide a single value per hailstone. The T-matrix reflectivity assuming solid ice spheres presents negligible differences to DDA results for size parameters below 1.0. Therefore, T matrix was used to fill in the gaps left by the DDA calculations. The results are mapped to the same size bins of the observed PSDs, allowing the calculation of the radar reflectivity. This is then correlated to E ˙ , allowing a potential improvement of past retrieval methods of E ˙ from Z in multiple wavelengths.

54 ENVIRONMENTAL SCIENCES↗

Total metals & anion concentration data; Slate River floodplain, Crested Butte, CO; May 2020-September 2020

This data package includes processed and undiluted measurements for metal and anion concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and September of 2020. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS:Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA.The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solution (Sigma-Aldrich Multielement standard solution for ICP, St. Louis, MO) using matrix matched to sample solutions (either 2% HCl or HNO3 from AriStar Plus,grade acids from VWR Scientific). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace Elements in Water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. An internal standard (Rh) is added via on-line addition into the sample line using a mixing tee.Analysis by Ion Chromatography (Anions):The protocol follows Method 4110 in Standard Methods for Examination of Water and Wastewater.The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with isocratic method using sodium carbonate eluent. Detection is by chemical suppression of eluent conductivity. Quality control solutions and mixed analyte standards purchased from Inorganic Ventures, Christiansburg, VA.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

Total metals, sulfur and organic carbon data; Slate River floodplain, Crested Butte, CO; March 2021-October 2021

This data package includes processed and undiluted measurements for metal, sulfur and organic carbon concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between March and October of 2021. These measurements were all recorded at the Environmental Measurements Facility (EM-1) at Stanford University in Stanford, CA. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 2C until measured at EM-1.Analysis by ICP-MS:Measurements for total metals were performed on an inductively coupled plasma optical emission spectrometer (ICP-OES; iCAP 6300, Thermo Scientific, Cambridge, U.K.). Calibration standards were prepared from the mulit element stock solution (Sigma-Aldrich Multielement standard for ICP, St. Louis, MO) using matrix matched to sample solutions (2% HNO3). Calibration curves included 5 points with correlation coefficients of >0.99. The QC protocol includes a continuing calibration blank and quality control samples that are analyzed just after calibration and again every 20 samples and at the completion of the run. Acceptable QC responses must be between 90-110% of the certified value. Analysis by Total Organic Carbon:Dissolved organic carbon concentrations were quantified on a total organic carbon (TOC) analyzer (TOC-L, Shimadzu, Kyoto, Japan) running the NPOC method. Standard curves were developed using an Organic Carbon Standard from RICCA Chemical Company (Arlington,TX). A series of 4 to 6 standards were automatically diluted by the instrument in a concentration range that spans that of the samples. A blank sample was run just after the calibration curve and at the end of the run. A QC sample was run every 25 samples. Acceptable QC responses must be between 90-110% of the certified value.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗

Dataset for 'Ombadi et al. (2023). A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature'

This package contains the main codes, sample input data and main result files to reproduce the analysis and results presented in the article: “Ombadi et al. (2023), A warming-induced reduction in snow fraction amplifies rainfall extremes, Nature”. The folder consists of the following: (1) “Raw data”: a folder that contains sample input data which is used in some of the codes for demonstration purposes. It also contains data that was not pre-processed such as Elevation data; (2) “Results”: this folder contains files of the main results presented in the paper including: “Annual-Max-Series”, “Change-rainfall-extremes”, “Change-snow-fraction”, “Warming levels_By scenario_model_year” and “Masks”. Description of these folders is detailed in the "Readme.rtf" file; (3) Python jupyter notebooks (Extract_Annual Max Series (AMS).ipynb, Elevation-dependent amplification of rainfall extremes.ipynb, Sensitivity_to_global_warming.ipynb) demonstrate the main steps of analysis. Further description of those notebooks is provided in the "Readme.rtf" file; (4) R code for extreme value analysis (Extreme_Value_Analysis.R). The sample and pre-processed dataset in "Raw data" is obtained from publicly available repositories of CMIP6 and ERA5 datasets; see Methods for more detail. This research was supported by Office of Science, Office of Biological and Environmental Research of the US Department of Energy under contract no. DE-AC02-05CH11231 for the CASCADE Scientific Focus (funded by the Regional and Global Model Analysis Program area within the Earth and Environmental Systems Modeling Program) and the iNAIADS Early Career Research Project (funded by the Environmental Systems Science program).

54 ENVIRONMENTAL SCIENCES↗

Data for publication: "A fresh take: Seasonal changes in terrestrial freshwater inputs impact salt marsh hydrology and vegetation dynamics"

This data repository contains data associated with the manuscript "A fresh take: Seasonal changes in terrestrial freshwater inputs impact salt marsh hydrology and vegetation dynamics". This study was conducted at the Elkhorn Slough National Estuarine Research Reserve in Watsonville, California from October 2019 - June 2022. We sought to understand the role of shallow freshwater inputs from adjacent uplands on salt marsh hydrologic behavior and vegetation productivity. This dataset contains CSV files of the following: daily salt marsh subsurface water level and pore water conductivity, monthly vegetation survey measurements, soil core data. Estuary surface water level, conductivity, and local precipitation were downloaded from the National Estuarine Research Reserve System (Centralized Data Management Office) at https://cdmo.baruch.sc.edu/.

54 ENVIRONMENTAL SCIENCES↗

Changes in the Regional Water Cycle and Their Impact on Societies

ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.

54 ENVIRONMENTAL SCIENCES↗

Atmospheric Instrument Systems and Technology in the Goddard Earth Sciences Division

Studies of the Earth’s atmosphere require a comprehensive set of observations that rely on instruments flown on spacecraft, aircraft, and balloons as well as those deployed on the surface. Within NASA’s Goddard Space Flight Center (GSFC) Earth Sciences Division-Atmospheres, laboratories and offices maintain an active program of instrument system development and observational studies that provide: 1) information leading to a basic understanding of atmospheric processes and their relationships with the Earth’s climate system, 2) prototypes for future flight instruments, 3) instruments to serve as calibration references for satellite missions, and 4) instruments for future field validation campaigns that support ongoing space missions. Our scientists participate in all aspects of instrument activity, including component and system design, calibration techniques, retrieval algorithm development, and data processing systems. The Atmospheres Program has well-equipped labs and test equipment to support the development and testing of instrument systems, such as a radiometric calibration and development facility to support the calibration of ultraviolet and visible (UV/VIS), space-borne solar backscatter instruments. This document summarizes the features and characteristics of 46 instrument systems that currently exist or are under development. The report is organized according to active, passive, or in situ remote sensing across the electromagnetic spectrum. Most of the systems are considered operational in that they have demonstrated performance in the field and are capable of being deployed on relatively short notice. Other systems are under study or of low technical readiness level (TRL). The systems described herein are designed mainly for surface or airborne platforms. However, two Cubesat systems also have been developed through collaborative efforts. The Solar Disk Sextant (SDS) is the single balloon-borne instrument. The lidar systems described herein are designed to retrieve clouds, aerosols, methane, water vapor pressure, temperature, and winds. Most of the lasers operate at some wavelength combination of 355, 532, and 1064 nm. The various systems provide high sensitivity measurements based on returns from backscatter or Raman scattering including intensity and polarization. Measurements of the frequency (Doppler) shift of light scattered from various atmospheric constitutes can also be made. Microwave sensors consist of both active (radar) and passive (radiometer) systems. These systems are important for studying processes involving water in various forms. The dielectric properties of water affect microwave brightness temperatures, which are used to retrieve atmospheric parameters such as rainfall rate and other key elements of the hydrological cycle. Atmosphere radar systems operate in the range from 9.6 GHz to 94 GHz and have measurement accuracies from -5 to 1 dBZ; radiometers operate in the 50 GHz to 874 GHz range with accuracies from 0.5 to 1 degree K; conical and cross-track scan modes are used. Our passive optical sensors, consisting of radiometers and spectrometers, collectively operate from the UV into the infrared. These systems measure energy fluxes and atmospheric parameters such as trace gases, aerosols, cloud properties, or altitude profiles of various species. Imager spatial resolution varies from 37 m to 400 m depending on altitude; spectral resolution is as small as 0.5 nm. Many of the airborne systems have been developed to fly on multiple aircraft.

Platnick, Steven E.↗

Predicting spring green-up across diverse North American grasslands

Vegetation phenology influences many ecosystem and climate processes, such as carbon uptake and energy and water cycles. Thus, understanding drivers of vegetation phenology is crucial for predicting current and future impacts of climate change on ecological systems. Existing models can accurately predict the date of spring green-up in temperate forests but tend to perform poorly in grassland systems. We hypothesize this is because most do not incorporate water availability, a primary limiting factor for grassland plants. In this study, we used long-term datasets of digital imagery from the PhenoCam Network of 43 diverse North American grassland sites (195 site-years) to test existing spring phenology models, as well as develop several new models that incorporate precipitation or soil moisture (53 models). As a result, most of the new models performed substantially better, with the best model requiring sufficient accumulated precipitation followed by warm temperatures to trigger spring onset (root mean square error, RMSE, between predicted and observed dates = 16.0 days). Importantly, the best model performed well across all grassland types using a single set of parameters, from temperate to arid grasslands. Since plants are adapted to their local climates, model performance was further improved when parameters were independently optimized for four separate climate regions (RMSE = 10.4 days). Therefore, both sufficient precipitation and temperature are required for grassland green-up, but optimal thresholds vary by region. Running the top model with projected climate data (representative concentration pathway 8.5) suggests that, depending on the climate region, spring onset will occur up to 12 days earlier within 100 years in temperature-limited sites, but the trend is unclear for precipitation-limited sites (3.5 ± 8.0 days later). This new phenology model improves our ability to understand and predict grassland dynamics, with implications for both current and future ecosystem processes related to carbon and water cycling.

54 ENVIRONMENTAL SCIENCES↗

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗