SEARCH · Engineering Papers
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.
Hybrid Environment Resources and Operations (HERO) is a system of common web services for mission driven research projects; the HERO Python software development kit (SDK) is a tool to programmatically interact with the HERO APIs.
Convective clouds play an important role in Earth’s climate system and are a known source of extreme weather. Gaps in our understanding of convective vertical motions, microphysics, and precipitation across a full range of aerosol and meteorological regimes continue to limit our ability to predict the occurrence and intensity of these cloud systems. To improve predictability, the National Science Foundation (NSF) sponsored a large field experiment entitled “Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE).” ESCAPE took place between 30 May and 30 September 2022 in the vicinity of Houston, Texas, because this area frequently experiences isolated deep convection that interacts with the region’s mesoscale circulations and its range of aerosol conditions. ESCAPE focused on collecting observations of isolated deep convection through innovative sampling and developing novel analysis techniques. This included the deployment of two research aircraft, the National Research Council of Canada Convair-580 and the Stratton Park Engineering Company Learjet, which combined conducted 24 research flights from 30 May to 17 June. On the ground, three mobile X-band radars and one mobile Doppler lidar truck equipped with soundings were deployed from 30 May to 28 June. From 1 August to 30 September 2022, a dual-polarization C-band radar was deployed and operated using a novel, multisensor agile adaptive sampling strategy to track the entire life cycle of isolated convective clouds. Analysis of the ESCAPE observations has already yielded preliminary findings on how aerosols and environmental conditions impact the convective life cycle.
Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.
Abstract An idealized large-eddy simulation of a tropical marine cloud population was performed. At any time, it contained hundreds of clouds, and updraft width in shallow convection emerging from a subcloud layer appeared to be an important indicator of whether specific convective elements deepened. In an environment with 80%–90% relative humidity below the 0°C level, updrafts that penetrated the 0°C level were larger at and above cloud base, which occurred at the lifting condensation level near 600 m. Parcels rising in these updrafts appeared to emerge from boundary layer eddies that averaged ∼200 m wider than those in clouds that only reached 1.5–3 km height. The deeply ascending parcels (growers) possessed statistically similar values of effective buoyancy below the level of free convection (LFC) as parcels that began to ascend in a cloud but stopped before reaching 3000 m (nongrowers). The growers also experienced less dilution above the LFC. Nongrowers were characterized by negative effective buoyancy and rapid deceleration above the LFC, while growers continued to accelerate well above the LFC. Growers occurred in areas with a greater magnitude of background convergence (or weaker divergence) in the subcloud layer, especially between 300 m and cloud base, but whether the convergence actually led to eddy widening is unclear. Significance Statement Cumulonimbus clouds are responsible for many extreme weather phenomena and are important contributors to Earth’s energy balance. However, the processes leading to the growth of individual clouds are not completely understood nor well-represented in weather prediction models. We find that the clouds containing updrafts that start out wider at early stages of their life cycles grow taller, possibly because they are protected more from drier air outside the cloud than narrow clouds. In addition, this work shows how the initial width of clouds might be related to convergence in the lowest part of the atmosphere, at heights where clouds initially develop. However, meteorologists must be careful not to overinterpret these results because numerical simulations inherently include assumptions that may not reflect reality. This reinforces the need to also observe processes occurring at the scales of individual clouds.
Convective instability is affected by both surface heat fluxes into near-surface air parcels and the atmospheric profile encountered as these air parcels are lifted aloft. In this study, we investigate the evolution of buoyancy experienced by surface-based parcels, taking interstorm soil moisture drydowns as fundamental time periods. We use profile observations of atmospheric temperature and humidity from ground-based passive remote sensing and in situ observations available at the U.S. Department of Energy Atmospheric Radiation Measurement Southern Great Plains site in Oklahoma, United States, and examine drydown events during the warm seasons from 2017 to 2022. Results reveal new insights on the emergence and growth of convective available potential energy (CAPE) during drydowns. We observe that the peak convective potential appears in the afternoon and grows in the midtroposphere (68% of 66 drydown events), with the largest CAPE values toward the end of the drydown event. We also find that the near-surface parcel contributes more to the evolution of CAPE than the atmospheric environment (86% of drydown events). These findings highlight the importance of surface heat fluxes and moisture availability from surface-based parcels in modulating the evolution of atmospheric conditions favorable for convection, underscoring the interplay between surface processes and atmospheric dynamics in driving convective instability.
During the RELAMPAGO-CACTI field experiments in 2018-19, an unprecedented number of balloon-borne soundings were collected in Argentina. Radiosondes were launched from both fixed and mobile platforms, yielding 2712 soundings during the period 15 October 2018-30 April 2019. Approximately 20% of these soundings were collected by highly mobile platforms, strategically positioned for each intensive observing period, and launching approximately once per hour. The combination of fixed and mobile soundings capture both the overall conditions characterizing the RELAMPAGO-CACTI campaign, as well as the detailed evolution of environments supporting the initiation and upscale growth of deep convective storms, including some that produced hazardous hail and heavy rainfall. Episodes of frequent convection were characterized by sufficient quantities of moisture and instability for deep convection, along with deep-layer vertical wind shear supportive of organized or rotating storms. Eleven soundings showed most-unstable convective available potential energy (MUCAPE) exceeding 6000 J kg-1, comparable to the extreme instability observed in other parts of the world with intense deep convection. Parameters used to diagnose severe-storm potential showed that conditions were often favorable for supercells and severe hail, but not for tornadoes, primarily owing to insufficient low-level wind shear. High-frequency soundings also revealed the structure and evolution of the boundary layer leading up to convection initiation, convectively generated cold pools, the South American Low-Level Jet (SALLJ), and elevated nocturnal convection. This sounding dataset will enable improved understanding and prediction of convective storms and their surroundings in subtropical South America, as well as comparisons with other heavily studied regions such as the central United States that have not previously been possible.
Intense deep convection and large mesoscale convective systems (MCSs) are known to occur downstream of the Andes in subtropical South America. Deep convection is often focused along the Sierras de Córdoba (SDC) in the afternoon and then rapidly grows upscale and moves to the east overnight. However, how the Andes and SDC impact the life cycle of MCSs under varying synoptic conditions is not well understood. Two sets of terrain-modification experiments using WRF are used to investigate the impact of topography in different synoptic regimes. The first set is run on the 13–14 December 2018 MCS case from RELAMPAGO, which featured a deep synoptic trough, strong lee cyclogenesis near the SDC, an enhanced low-level jet, and rapid upscale growth of an MCS. When the Andes are reduced by 50%, the lee cyclone and low-level jet that develop are weaker than with the full Andes, and the resulting MCS is weak and moves faster to the east. When the SDC are removed, few differences between the environment and resulting MCS relative to the control run are seen. A second set of experiments are run on the 25–26 January 2019 case in which a large MCS developed over the SDC and remained tied there for an extended period under weak synoptic forcing. Finally, the experiment that produces the most similar MCS to the control is when the Andes are reduced by 50% while maintaining the height of the SDC, suggesting the SDC may play a more important role in the MCS life cycle under quiescent synoptic conditions.
This study documents the spatial and temporal distribution of the South American low-level jet (SALLJ) and quantifies its impact on the convective environment using a 6.5-month convection-permitting simulation during the Remote Sensing of Electrification, Lightning, And Mesoscale/Microscale Processes with Adaptive Ground Observations and Clouds, Aerosols, and Complex Terrain Interactions (RELAMPAGO-CACTI) campaigns. Overall, the simulation reproduces the observed SALLJ characteristics in central Argentina near the Sierras de Córdoba (SDC), a focal point for terrain-focused upscale growth. SALLJs most frequently occur in the summer with maxima to the northwest and east of the SDC and minima over the higher terrain. The shallower SALLJs (<1750 m) have a strong overnight skew, while the elevated jets are more equally spread throughout the day. SALLJ periods often have higher amounts of low-level moisture and instability compared to non-SALLJ periods, with these impacts increasing over time when the SALLJ is present and decreasing afterward. The SALLJ may enhance low-level wind shear magnitudes (particularly when accounting for the jet height); however, enhancement is somewhat limited due to the presence of speed shear in most situations. SALLJ periods are associated with low-level directional shear favorable for organized convection and an orientation of cloud-layer wind shear parallel to the terrain, which could favor upscale growth. A case study is shown in which the SALLJ influenced both the magnitude and direction of wind shear concurrent with convective upscale growth near the SDC. This study highlights the complex relationship between the SALLJ and its impacts during periods of widespread convection.
This article explores the role of land cover in relation to housing type and tenure in shaping exposure to extreme urban heat, focusing on residents of mobile and manufactured housing (MH). We hypothesize that MH residents will experience greater exposure to extreme heat than those living in other housing types due to lower levels of proximate vegetation. This hypothesis is based on the unique property relations and tenure regimes that characterize MH, which may disincentivize investments in planting and maintaining trees and vegetation in arid environments. To test this hypothesis, we compare the amount of vegetation on properties across housing types, within-type tenure arrangements and between three urbanized areas in Arizona with different levels of exposure to extreme heat. To conduct this comparison, we combine multispectral land cover data for 1.7 million parcels with tax assessor data and building footprints to measure land cover at a high resolution. We find that MH units have significantly less vegetation than single-family residential properties, and that MH units in MH parks have less vegetation than those on individual lots. We conclude that municipalities should promote (e.g. through incentives or other policy interventions) the planting of more vegetation around MH units where there is a risk of extreme heat exposure. Future research can expand this analysis with a closer examination of how municipal ordinances and policies affect land cover by housing type and tenure.
Existing discontinuous palynological records from coastal and river valley exposures, with varying quality of radiocarbon control, suggested that the regional Holocene climate in southern areas of the Russian Far East was characterized by as many as 10 fluctuations in temperature and/or precipitation. In this study, palynological data from Zerkalnoye Lake, located on the western coast of the Sea of Japan, indicate a gradual decline in temperature through the Middle and Late-Holocene with a Holocene thermal maximum between 8800 and 5500 cal yr BP. Three wetter than present intervals, marked by an increase in Pinus koraiensis, occurred c. 3600–3500 cal yr BP, 2340–2050 cal yr BP, and 1830–1800 cal yr BP. The Zerkalnoye record shows the dominance of Quercus-broadleaf forests during the Middle and Early Holocene, although Quercus shows a gradual decrease as climate cooled during this interval. Diatom analysis of the Zerkalnoye sediments documents that the site was a shallow bay or coastal lagoon until c. 3540 cal yr BP. After that time, a freshwater lake was established, which had variable marine influences probably caused by sea level changes following deglaciation. The diatom data indicate a cool water interval between 2900 and 2580 cal yr BP, also noted in other sites in the region. Changes in the basin’s depositional environment do not affect the palynological record, indicating such sites can provide reliable paleovegetational and paleoclimatic records. In conclusion, the discrepancies of the various regional paleoclimatic scenarios indicate the need for the further collection of continuous records from coastal to alpine zones in this region of northeastern Asia.
The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known. Results Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of >600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group.
Background: Humans emit approximately 30 million microbial cells per hour into their immediate vicinity. However, sampling of aerosolized microbial taxa (aerobiome) remains largely uncharacterized due to the complexity and limitations of sampling techniques, which are highly susceptible to low biomass and rapid sample degradation. Recently, there has been an interest in developing technology that collects naturally occurring water from the atmosphere, even within the built environment. Here, we analyze the feasibility of indoor aerosol condensation collection as a method to capture and analyze the aerobiome. Methods: Aerosols were collected via condensation or active impingement in a laboratory setting over the course of 8 h. Microbial DNA was extracted from collected samples and sequenced (16S rRNA) to analyze microbial diversity and community composition. Dimensional reduction and multivariate statistics were employed to identify significant (p < 0.05) differences in relative abundances of specific microbial taxa observed between the two sampling platforms. Results: Aerosol condensation capture is highly efficient with a yield greater than 95% when compared to expected values. Compared to air impingement, aerosol condensation showed no significant difference (ANOVA, p > 0.05) in microbial diversity. Among identified taxa, Streptophyta and Pseudomonadales comprised approximately 70% of the microbial community composition. Conclusion: The results suggest that condensation of atmospheric humidity is a suitable method for the capture of airborne microbial taxa reflected by microbial community similarity between devices. Future investigation of aerosol condensation may provide insight into the efficacy and viability of this new tool to investigate airborne microorganisms.
Background: The Atribacterota are widely distributed in the subsurface biosphere. Recently, the first Atribacterota isolate was described and the number of Atribacterota genome sequences retrieved from environmental samples has increased significantly; however, their diversity, physiology, ecology, and evolution remain poorly understood. Results: We report the isolation of the second member of Atribacterota, Thermatribacter velox gen. nov., sp. nov., within a new family Thermatribacteraceae fam. nov., and the short-term laboratory cultivation of a member of the JS1 lineage, Phoenicimicrobium oleiphilum HX-OS.bin.34 TS , both from a terrestrial oil reservoir. Physiological and metatranscriptomics analyses showed that Thermatribacter velox B11 T and Phoenicimicrobium oleiphilum HX-OS.bin.34 TS ferment sugars and n-alkanes, respectively, producing H 2 , CO 2 , and acetate as common products. Comparative genomics showed that all members of the Atribacterota lack a complete Wood-Ljungdahl Pathway (WLP), but that the Reductive Glycine Pathway (RGP) is widespread, indicating that the RGP, rather than WLP, is a central hub in Atribacterota metabolism. Ancestral character state reconstructions and phylogenetic analyses showed that key genes encoding the RGP (fdhA, fhs, folD, glyA, gcvT, gcvPAB, pdhD) and other central functions were gained independently in the two classes, Atribacteria (OP9) and Phoenicimicrobiia (JS1), after which they were inherited vertically; these genes included fumarate-adding enzymes (faeA; Phoenicimicrobiia only), the CODH/ACS complex (acsABCDE), and diverse hydrogenases (NiFe group 3b, 4b and FeFe group A3, C). Finally, we present genome-resolved community metabolic models showing the central roles of Atribacteria (OP9) and Phoenicimicrobiia (JS1) in acetate- and hydrocarbon-rich environments. Conclusion: Our findings expand the knowledge of the diversity, physiology, ecology, and evolution of the phylum Atribacterota. This study is a starting point for promoting more incisive studies of their syntrophic biology and may guide the rational design of strategies to cultivate them in the laboratory.
This manuscript investigates the potential effect of a nuclear-disturbed atmospheric environment on the signal attenuation of a ground/satellite transmitter/receiver system for both classical optical and quantum communications applications. Attenuation of a signal transmitted through the rising nuclear cloud and the subsequently transported debris is modeled climatologically for surface-level detonations of 10 kt, 100 kt, and 1 Mt. Attenuation statistics were collected as a function of time after detonation. These loss terms were compared to normal loss sources such as clouds, smoke from fires, and clear sky operation. Finally, the loss was related to the degradation of transmitted entanglement derived from Bayesian mean estimation.
The Bering Sea experiences a seasonal sea ice cover, which is important to the biophysical environment found there. A pool of cold bottom water (<2°C) is formed on the shelf each winter as a result of cooling and vertical mixing due to brine rejection during the predominately local sea ice growth. The extent and distribution of this Cold Pool (CP) is largely controlled by the winter extent of sea ice in the Bering Sea, which can vary considerably and recently has been much lower than average. The cold bottom water of the CP is important for food security because it delineates the boundary between arctic and subarctic demersal fish species. A northward retreat of the CP will likely be associated with migration of subarctic species toward the Chukchi Sea. We use the fully-coupled Regional Arctic System Model (RASM) to examine variability of the extent and distribution of the CP and its relation to change in the sea ice cover in the Bering Sea during the period 1980–2018. RASM results confirm the direct correlation between the extent of sea ice and the CP and show a smaller CP as a consequence of realistically simulated recent declines of the sea ice cover in the Bering Sea. In fact, the area of the CP was found to be only 31% of the long-term mean in July of 2018. In addition, we also find that a low ice year is followed by a later diatom bloom, while a heavy ice year is followed by an early diatom bloom. Finally, the RASM probabilistic intra-annual forecast capability is reviewed, based on 31-member ensembles for 2019–2021, for its potential use for prediction of the winter sea ice cover and the subsequent summer CP area in the Bering Sea.
DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), weights needed to be developed for use in the weighted sum of the different favorability index models produced from geoscientific exploration datasets. This GDR submission includes those weights. The weighting was done using two different approaches: one based on expert opinions, and one based on statistical learning. The weights are intended to describe how useful a particular exploration method is for imaging each component of each play type. They may be adjusted based on the characteristics of the resource under investigation, knowledge of the quality of the dataset, or simply to reduce the impact a single dataset has on the resulting outputs. Within the DEEPEN PFA, separate sets of weights are produced for each component of each play type, since exploration methods hold different levels of importance for detecting each play component, within each play type. The weights for conventional hydrothermal systems were based on the average of the normalized weights used in the DOE-funded PFA projects that were focused on magmatic plays. This decision was made because conventional hydrothermal plays are already well-studied and understood, and therefore it is logical to use existing weights where possible. In contrast, a true PFA has never been applied to superhot EGS or supercritical plays, meaning that exploration methods have never been weighted in terms of their utility in imaging the components of these plays. To produce weights for superhot EGS and supercritical plays, two different approaches were used: one based on expert opinion and the analytical hierarchy process (AHP), and another using a statistical approach based on principal component analysis (PCA). The weights are intended to provide standardized sets of weights for each play type in all magmatic geothermal systems. Two different approaches were used to investigate whether a more data-centric approach might allow new insights into the datasets, and also to analyze how different weighting approaches impact the outcomes. The expert/AHP approach involved using an online tool (https://bpmsg.com/ahp/) with built-in forms to make pairwise comparisons which are used to rank exploration methods against one-another. The inputs are then combined in a quantitative way, ultimately producing a set of consensus-based weights. To minimize the burden on each individual participant, the forms were completed in group discussions. While the group setting means that there is potential for some opinions to outweigh others, it also provides a venue for conversation to take place, in theory leading the group to a more robust consensus then what can be achieved on an individual basis. This exercise was done with two separate groups: one consisting of U.S.-based experts, and one consisting of Iceland-based experts in magmatic geothermal systems. The two sets of weights were then averaged to produce what we will from here on refer to as the "expert opinion-based weights," or "expert weights" for short. While expert opinions allow us to include more nuanced information in the weights, expert opinions are subject to human bias. Data-centric or statistical approaches help to overcome these potential human biases by focusing on and drawing conclusions from the data alone. More information on this approach along with the dataset used to produce the statistical weights may be found in the linked dataset below.
This data package is associated with the manuscript “Disturbance Triggers Non-Linear Microbe-Environment Feedbacks, in revision in Biogeosciences (Sengupta et al. and 2021;https://bg.copernicus.org/preprints/bg-2021-51/). The study used hyporheic zone sediments as a model system to provide an integrated view of how disturbance modulates linkages among microbial ecology, biogeochemistry, and organic matter thermodynamics. Laboratory experiments exposed hyporheic sediment to varying wetting/drying dynamics. Data types include dissolved oxygen rates used to derive respiration rates, Fourier-transform ion cyclotron resonance mass spectrometry (FTICR-MS) data used to derive thermodynamic properties of organic matter, and microbial community assembly metrics derived from amplicon-sequence data of putatively active (cDNA) and whole community (gDNA). The outcomes of the study are condensed into a broadly applicable conceptual model linking external forcing, internal dynamics, and history. This data package is comprised of a file-level metadata (FLMD) csv, metadata csv, and seven folders that contain csv files, R scripts, xml files, and associated documentation: (1) Rates, (2) bNTI, and (3) FTICR, (4) Statistics_Analyses, (5) Raw OTU Beta dispersion Analysis, (6) bMNTD Randomizations, and (7) Data Dictionaries. The FLMD file has a description of each file included in the data package.
Objectives: Climatic drying is predicted for many tropical forests, yet models remain poorly parameterized for tropical forests, hampering predictions of forest-climate feedbacks. We applied an integrated model–experiment approach, parameterizing an ecosystem model Carbon Organisms Rhizosphere and Protection in the Soil Environment (CORPSE) with tropical forest observational data, and comparing model predictions with a field drying manipulation. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis. 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 used the field data to parameterize and run tests in the model.Results: Measured 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. We used this data to parameterize the model, which then predicted increased soil CO2 fluxes in wetter and fertile forests with drying, and decreased fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment in the forests initially suppressed soil CO2 fluxes across forests, with sustained suppression 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), as predicted by the model. The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Code files:CORPSE_array.py: Defines the equations of the CORPSE modelCORPSE_solvers: Functions for running the CORPSE model using either iterative or ordinary differential equation (ODE) solversrun_Panama_sims.py: Read in datasets and run the model simulations for this studyInput data:PanamaGradientEcosystemChem_BT_CPools_20152016CO2_DC_20190615.xlsx: Plot characteristics used in running model simulationsLiCor compiled surface flux only to 2020_03 DC_20200825.xlsx: Surface gas exchange fluxes used in model-data comparisonsPARCHED litterfall data for Ben Sulman LD 20200902.xlsx: Litterfall data used to drive model simulationsInitialization data:state_500y_20190823.csv: Initial state of model pools based on previous spinup runsOutput data:Outputs/prev_moisture_response.csv: Simulations of multiple sites using original model moisture response function.Outputs/updated_moisture_response.csv: Simulations of multiple sites using updated model moisture response function.Outputs/dry15_prev_moisture_response.csv: Simulations with soil moisture reduced by 15%, using original moisture response function.Outputs/dry15_updated_moisture_response.csv: Simulations with soil moisture reduced by 15%, using updated moisture response function.Outputs/dry30_prev_moisture_response.csv: Simulations with soil moisture reduced by 30%, using original moisture response function.Outputs/dry30_updated_moisture_response.csv: Simulations with soil moisture reduced by 30%, using updated moisture response function.Outputs/latestart_prev_moisture_response.csv: Simulations with extended dry season, using original moisture response function.Outputs/latestart_updated_moisture_response.csv: Simulations with extended dry season, using updated moisture response function.Outputs/[site name]_oneyear.csv: One-year simulation for each site in expanded site list using original moisture response function.Outputs/[site name]_oneyear_dried.csv: One-year simulation for each site in expanded site list using original moisture response function, with soil moisture reduced by 25%.Outputs/[site name]_oneyear_updated_moisture_response.csv: One-year simulation for each site in expanded site list using updated moisture response function.Outputs/[site name]_oneyear_updated_moisture_response_dried.csv: One-year simulation for each site in expanded site list using updated moisture response function, with soil moisture reduced by 25%.Field plot location data: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).