Search NASA⌕ Search

SEARCH · Search NASA

Results for “flux vertical gradients”

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.

143 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Alma Measurements of the HNC and HC3N Distributions in Titan's Atmosphere

We present spectrally and spatially resolved maps of HNC and HC3N emission from Titan's atmosphere, obtained using the Atacama Large Millimeter/submillimeter Array on 2013 November 17. These maps show anisotropic spatial distributions for both molecules, with resolved emission peaks in Titan's northern and southern hemispheres. TheHC3N maps indicate enhanced concentrations of this molecule over the poles, consistent with previous studies of Titan's photochemistry and atmospheric circulation. Differences between the spectrally integrated flux distributions of HNC and HC3N show that these species are not co-spatial. The observed spectral line shapes are consistent with HNC being concentrated predominantly in the mesosphere and above (at altitudes z approx.. greater than 400 km), whereas HC3N is abundant at a broader range of altitudes (z approx. equal to 70-600 km). From spatial variations in the HC3N line profile, the locations of the HC3N emission peaks are shown to be variable as a function of altitude. The peaks in the integrated emission from HNC and the line core (upper atmosphere) component of HC3N (at z approx. greater than 300 km) are found to be asymmetric with respect to Titan's polar axis, indicating that the mesosphere may be more longitudinally variable than previously thought. The spatially integrated HNC and HC3N spectra are modeled using the NEMESIS planetary atmosphere code and the resulting best-fitting disk-averaged vertical mixing ratio profiles are found to be in reasonable agreement with previous measurements for these species. Vertical column densities of the best-fitting gradient models for HNC and HC3N are 1.9 × 10(exp 13) per sq.cm and 2.3 × 10(exp 14) per sq.cm, respectively.

ALMA Telescope↗
NASA NTRSDate not supplied

The Response of the Quasi-Biennial Oscillation to Increased CO and its Modulation by Composition Feedbacks

The Quasi-Biennial Oscillation (QBO) impacts the large-scale circulation by altering lower stratospheric meridional temperature gradients which alter the propagation of upward propagating planetary waves. Proposed teleconnections of the QBO include impacts on the stratospheric polar vortices, extratropical surface winter climate, and the Madden-Julian Oscillation. However, long-term projections of the QBO remain highly uncertain. While recent multi-model investigations (Richter et al. 2019, Butchart et al. 2020) show that the amplitude of the QBO weakens robustly among models, changes in QBO period and stability remain highly uncertain. Here we examine the QBO response to increased greenhouse gases using the NASA Goddard Institute for Space Studies Middle Atmosphere Model E2.2 (Rind et al., 2020; Orbe et al. 2020). Compared to lower vertical resolution versions of ModelE, E2.2 also has a higher model top (0.002 hPa) and employs additional interactive non-orographic gravity wave drag sources from convection and shear, which produce a sufficiently realistic QBO, thus rendering it suitable for use in climate change studies. Overall, we find that both the QBO period and amplitude decrease in response to increased CO , the former related to increased lower stratospheric momentum fluxes (associated with convection) and the latter associated partly with a stronger residual mean circulation. Experiments integrated using fully interactive chemistry also reveal that ozone feedbacks significantly impact the magnitude of the QBO amplitude response. Finally, integrations using fixed (pre-industrial) sea surface temperatures (SST) show that the QBO amplitude responds differently to rapid adjustments versus SST feedbacks, as compared to the QBO period.

Quasi-Biennial Oscillation↗
DOE OSTI2023

What controls the interannual variation of Hadley cell extent in the Northern Hemisphere: physical mechanism and empirical model for edge variation

Abstract The Hadley circulation is the most prominent atmospheric meridional circulation, reducing the radiatively driven equator-to-pole temperature gradient. While the Hadley cell extent varies by several degrees from year to year, the detailed dynamical mechanisms behind such variations have not been well elucidated. During the expanded phase of the Hadley cell, many regions on the periphery of the subtropics experience unfavorable climatic conditions. In this study, using ERA5 reanalysis data, we examine the physical chain of events responsible for the interannual variation of the Hadley cell edge (HCE) latitude in the Northern Hemisphere. This variation is mainly caused by changing eddy activity and wave breaking from both stationary and transient waves. In particular, we show that transient waves cause the HCE to shift poleward by increasing the eddy momentum flux divergence (EMFD) and reducing the baroclinicity over 20°–40°N, shifting the region of peak baroclinicity poleward. El Niño/La Niña and the Arctic Oscillation (AO) account for a significant portion (60%) of the interannual fluctuation of the HCE latitude. Through the poleward displacement of eddy activity, La Niña and a positive AO state are associated with the poleward shift of the HCE. The analysis of 28 CMIP5 models reveals statistical relationships between EMFD, vertical shear, and HCE latitude similar to those observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of Transformative Preparation Methods to Push up High Q&G Performance of FRIB Spare HWR Cryomodule Cavities

The FRIB accelerator project construction, a top priority of US nuclear science, was completed in January 2022, and is now moving to user operation. The stable and reliable operation of the accelerating cryomodules is essential in achieving/fulfilling DOE and user expectations. So far, FRIB cryomodules meet all FRIB specifications for cavity performance. However, during the lifetime of machine operation, degradation of cryomodule performance is possible, as reported in similar operating facilities (CEBAF, SNS). If cryomodule degradation is observed at FRIB, the under-performing cryomodule will require replacement/maintenance. In effort to manage operational reliability, FRIB plans to construct a 0.53 half-wave cryomodule to serve as an active spare. In a parallel effort, FRIB will also work toward increasing operational Q and gradient of spare cryomodule cavities to gain an overall performance margin to support future operational reliability. The current FRIB cavity designs have a potential to operate at gradients higher than 8 MV/m, but are currently limited by field emission (FE) and/or high field Q slope (HFQS); known issue in buffered chemical polished (BCP) treated cavities. The proposal looks to develop transformative surface preparation treatments to improve the operational gradient of spare cryomodules higher than 10 MV/m while maintaining high Q. Thus, increasing operational margin by 30 - 50%. With the goal to improve operational reliability set, the proposal will investigate multiple objectives as possible paths forward to achieve an overall increase in cavity performance and gain a better understanding of SRF limiting mechanisms. The proposal will study the application of different chemical surface treatments to 0.53 half-wave cavities, with the addition of low temperature bakes (LTB), and measure their effects on accelerating performance. Proposed chemical treatments to be explored in this proposal include conventional EP acid mixtures, as well as innovated EP and BCP acid mixtures designed to simplify processing paths in migrating FE and HFQS. The proposed transformative treatment wet N-doping also has the potential to replicate recent advancements in SRF technology relating to nitrogen doping and high Q operation without the requirement for an ultra-high vacuum annealing furnace; currently being developed at FNAL and JLAB. In parallel, high Q performance relating to flux trapping will be investigated with the installation of a second layer of magnetic shielding in the vertical test Dewar. The research objectives presented in the proposal, and their corresponding effects on cavity performance, will provide essential knowledge and future guidance to the SRF community and provide possible paths for future SRF based projects and applications.

43 PARTICLE ACCELERATORS↗
NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗
DOE OSTI2026

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Electric Plasma Arc-Lamp Combustor Liner Durability Test System Developed

Silicon carbide matrix composites are candidate materials for high-temperature combustor liners. Because through-the-thickness thermal gradients are the primary cause of stress on combustor liners, a unique test facility was developed at the NASA Glenn Research Center at Lewis Field to simulate in-service pure thermal stress distributions in fiber reinforced silicon carbide cylinders. It was developed initially under Phase II of the High- Speed Research Program. This test stand can accommodate 8-in.-long test cylinders that have outer diameters of 4 in. and a wall thickness of about 0.08 to 0.12 in. One cylinder at a time is loaded vertically into the test stand. Water-cooled plates enclose the open ends of the cylinder and provide cooling. Load plates on the exterior side of the water-cooled plates provide support and compression loads. To evaluate a combustor liner material s potential performance, researchers induce thermal gradients with an axisymmetric, direct-current, electric arc within the cylinder while refrigerated air at a rate of 1.5 lb/sec impinges on the outside surface of the cylinder. The achievable through-the-thickness thermal gradient is predicted to be in excess of 200 C. The 8-in. long, 0.5-in.-diameter plasma arc emits full spectrum visible light; radiant intensity exceeds 300 Watts per square centimeters to produce temperatures in excess of 1500 C on materials with emissivity near unity. Because the system does not rely upon the combustion of fuels to achieve the related thermal conditions, ancillary environmental reactions with the sample are eliminated. The system incorporates a standard mechanical test frame, which can impose constant as well cyclical axial stresses up to 2200 lb upon the test piece. Silicon-carbide-fiber reinforced silicon carbide matrix composite cylinders were instrumented with thin-film thermocouples to obtain through-the-thickness thermal flux measurements. Inside wall temperatures reached 1200 C with only 250 A of current. One of the special features of this configuration is the creation of hoop stress states within the cylinder, which up this point have not been obtainable in planar coupon tests. This facility will allow various operational modes, including accelerated tests of thermal transients simulating the effects of repeated engine ignition as well as prescribed thermal and mechanical histories to simulate various duty cycle profiles. Tests can now be performed on thermal-barrier-coated metallic liners and ceramic composite liners that require a combination of high heat flux and controlled mechanical stresses.

Goldsby, Jon C.↗
NASA NTRS2005

Simulations of Dynamics and Transport during the September 2002 Antarctic Major Warming

A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

winds↗
DOE OSTI2023

Throughfall-Reduction Drying Effects on Soil Moisture and Soil Temperature at 30-Minute Intervals Four Lowland Panamanian Forests

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (Cusack et al. 2023). Datasets included: Datasets here include .csv and .xls files for volumetric soil moisture (volume/volume) and soil temperature (°C). Data were taken at 30-minute intervals at 5cm and 20cm depths using permanent probes and data loggers (see methods). Data presented are from 2018-2019. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Satellite Remote Sensing of Snow/Ice Albedo over the Himalayas

The Himalayan glaciers and snowpacks play an important role in the hydrological cycle over Asia. The seasonal snow melt from the Himalayan glaciers and snowpacks is one of the key elements to the livelihood of the downstream densely populated regions of South Asia. During the pre-monsoon season (April-May-June), South Asia not only experiences the reversal of the regional meridional tropospheric temperature gradient (i.e., the onset of the summer monsoon), but also is being bombarded by dry westerly airmass that transports mineral dust from various Southwest Asian desert and arid regions into the Indo-Gangetic Plains in northern India. Mixed with heavy anthropogenic pollution, mineral dust constitutes the bulk of regional aerosol loading and forms an extensive and vertically extended brown haze lapping against the southern slopes of the Himalayas. Episodic dust plumes are advected over the Himalayas, and are discernible in satellite imagery, resulting in dust-capped snow surface. Motivated by the potential implications of accelerated snowmelt, we examine the changes in radiative energetics induced by aerosol transport over the Himalayan snow cover by utilizing space borne observations. Our objective lies in the investigation of potential impacts of aerosol solar absorption on the Top-of-Atmosphere (TOA) spectral reflectivity and the broadband albedo, and hence the accelerated snowmelt, particularly in the western Himalayas. Lambertian Equivalent Reflectivity (LER) in the visible and near-infrared wavelengths, derived from Moderate Resolution Imaging Spectroradiometer radiances, is used to generate statistics for determining perturbation caused due to dust layer over snow surface in over ten years of continuous observations. Case studies indicate significant reduction of LER ranging from 5 to 8% in the 412-860nm spectra. Broadband flux observations, from the Clouds and the Earth's Radiant Energy System, are also used to investigate changes in shortwave TOA flux over dust-laden and dust-free snow covered regions. Additionally, spatio-temporal and intra-seasonal variations of LER, along with snow cover information, are used to characterize the seasonal melt pattern and thus to distinguish the outstanding aerosol-induced snowmelt signal. Results from this observational work are expected to provide better understanding of the radiative impact of aerosols over snow surface, especially its role in the Himalayan hydro-glacialogical variability.

Hsu, N. Christina↗
DOE OSTI2022

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tropical Tree Species Identity and Diameter at Breast Height in a Throughfall-Reduction Drying Experiment in Four Lowland Panamanian Forests

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests to 20cm depths, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Datasets included here include .csv and .xls files for tree species identity and diameter at breast height (dbh) in the study plots. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Throughfall-Reduction Drying Effects on Gravimetric Soil Moisture at Two Depths in Four Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests to 20cm depths, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Datasets included here include .csv and .xls files for gravimetric soil moisture (weight/weight). Soil moisture was collected on a quarterly basis from 0-10 cm and 10-20 cm depths using hand-held constant-volume soil corers. Data are on ~3 month timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Throughfall-Reduction Drying Effects on Forest Floor Plant Tissue Biomass (Oi+Oe soil horizons) and Moisture Content in Four Lowland Panamanian Forests

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We also surveyed all tree species and identified to species in 2018 and 2019.Results: We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Forest floor biomass peaked in the dry season, in contrast. Chronic throughfall exclusion also suppressed soil moisture across the four forests to 20cm depths, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Datasets included here include .csv and .xls files for forest floor biomass (equivalent to the Oi+Oe soil horizons), and a subset include forest floor biomass moisture content (weight/weight) in the study plots. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Throughfall-Reduction Drying Effects on Soil Respiration, Moisture, and Temperature in Four Distinct Lowland Panamanian Forests from 2015-2022

Objectives: Climatic drying is predicted for many tropical forests, yet effects on soil properties across moisture and soil gradients within tropical forests remain poorly characterized, hampering predictions of forest-climate feedbacks. We hypothesized that throughfall reduction (i.e. chronic drying) would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests by further reductions in soil moisture, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis and soil saturation. We conducted a field study in four distinct lowland Panamanian forests that vary in background rainfall and soil fertility. In each forest we constructed throughfall reduction structures that reduce moisture inputs to plots by ~50% (see methods). In drying and paired control plots we measured soil CO2 fluxes, soil moisture (volume/volume), soil temperature (°C), air temperature (°C), and forest floor biomass during wet-dry cycles (2015 – 2022). We also collected litterfall biomass and surveyed all tree species for diameter at breast height (dbh), and identified tree species in 2018 and 2019.Results:We found that soil moisture peaked in the wet season and declined in the dry season. Measured soil CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. Chronic throughfall exclusion also suppressed soil moisture across the four forests, and also initially suppressed soil CO2 fluxes across forests. There was sustained suppression of soil CO2 fluxes after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season). The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils, as the drying effect increased with time. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks. Datasets included: Datasets included here include .csv and .xls files for soil respiration (carbon dioxide flux) collected using the LiCor (see methods), air temperature °C , soil moisture (vl/vl), and soil temperature °C. Soil moisture and soil temperature were collected at 0-10cm depths using hand-held probes (see methods). Data are on ~3 week timescales from 2015-2022 with some gaps. There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site). No special software is needed to open these files.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2004

Determining Sizes of Particles in a Flow from DPIV Data

A proposed method of measuring the size of particles entrained in a flow of a liquid or gas would involve utilization of data from digital particle-image velocimetry (DPIV) of the flow. That is to say, with proper design and operation of a DPIV system, the DPIV data could be processed according to the proposed method to obtain particle sizes in addition to particle velocities. As an additional benefit, one could then compute the mass flux of the entrained particles from the particle sizes and velocities. As in DPIV as practiced heretofore, a pulsed laser beam would be formed into a thin sheet to illuminate a plane of interest in a flow field and the illuminated plane would be observed by means of a charge-coupled device (CCD) camera aimed along a line perpendicular to the illuminated plane. Unlike in DPIV as practiced heretofore, care would be taken to polarize the laser beam so that its electric field would lie in the illuminated plane, for the reason explained in the next paragraph. The proposed method applies, more specifically, to transparent or semitransparent spherical particles that have an index of refraction different from that of the fluid in which they are entrained. The method is based on the established Mie theory, which describes the scattering of light by diffraction, refraction, and specular reflection of light by such particles. In the case of a particle illuminated by polarized light and observed in the arrangement described in the preceding paragraph, the Mie theory shows that the image of the particle on the focal plane of the CCD camera includes two glare spots: one attributable to light reflected toward the camera and one attributable to light refracted toward the camera. The distance between the glare spots is a known function of the size of the particle, the indices of refraction of the particle material, and design parameters of the camera optics. Hence, the size of a particle can be determined from the distance between the glare spots. The proposed method would be implemented in an algorithm that would automatically identify, and measure the distance between, the glare spots for each particle for which a suitable image has been captured in a DPIV image frame. The algorithm (see figure) would begin with thresholding of data from the entire image frame to reduce noise, thereby facilitating discrimination of particle images from the background and aiding in the separation of overlapping particles. It is important not to pick a threshold level so high that the light intensity between a given pair of glare spots does not fall below the threshold value, leaving the glare spots disconnected. The image would then be scanned in a sequence of rows and columns of pixels to identify groups of adjacent pixels that contain nonzero brightnesses and that are surrounded by pixels of zero brightness. Each such group would be assumed to constitute the image of one particle. Each such group would be further analyzed to determine whether the image was saturated; saturated particle images must be rejected because the locations of glare spots in saturated images cannot accurately be determined. Within each unsaturated particle image, the centroids (deemed to be the locations) of the glare spots would be determined by means of gradients of brightness distributions and three-point horizontal and three-point vertical Gaussian estimates based on the brightness values of the brightest pixels and the pixels adjacent to them. If the brightness of a given particle image contained only one peak, then it would be assumed that a second glare spot did not exist and that image would be rejected.

Wernet, M. P.↗