Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained regions”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tundra Underlain By Thawing Permafrost Persistently Emits Carbon to the Atmosphere Over 15 Years of Measurements

Warming of the Arctic can stimulate microbial decomposition and release of permafrost soil carbon (C) as greenhouse gases, and thus has the potential to influence climate change. At the same time, plant growth can be stimulated and offset C release. This study presents a 15-year time series comprising chamber and eddy covariance measurements of net ecosystem C exchange in a tundra ecosystem in Alaska where permafrost has been degrading due to regional warming. The site was a carbon dioxide source to the atmosphere with a cumulative total loss of 781.6 g C m -2 over the study period. Both gross primary productivity (GPP) and ecosystem respiration (R eco ) were already likely higher than historical levels such that increases in R eco losses overwhelmed GPP gains in most years. This shift to a net C source to the atmosphere likely started in the early 1990s when permafrost was observed to warm and thaw at the site. Shifts in the plant community occur more slowly and are likely to constrain future GPP increases as compared to more rapid shifts in the microbial community that contribute to increased Reco. Observed rates suggest that cumulative net soil C loss of 4.18–10.00 kg C m -2 —8%–20% of the current active layer soil C pool—could occur from 2020 to the end of the century. This amount of permafrost C loss to the atmosphere represents a significant accelerating feedback to climate change if it were to occur at a similar magnitude across the permafrost region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Inclusive Search for Anomalous Single-Photon Production in MicroBooNE

We present an inclusive search for anomalous production of single-photon events from neutrino interactions in the MicroBooNE experiment. The search and its signal definition are motivated by the previous observation of a low-energy excess of electromagnetic shower events from the MiniBooNE experiment. We use the Wire-Cell reconstruction framework to select a sample of inclusive single-photon final-state interactions with a final efficiency and purity of 7.0% and 40.2%, respectively. We leverage simultaneous measurements of sidebands of charged current 𝜈 𝜇 interactions and neutral current interactions producing 𝜋 0 mesons to constrain signal and background predictions and reduce uncertainties. We perform a blind analysis using a dataset collected from February 2016 to July 2018, corresponding to an exposure of 6.34 × 10 20 protons on target from the Booster Neutrino Beam at Fermilab. In the full signal region, we observe agreement between the data and the prediction, with a goodness-of-fit 𝑝 value of 0.11. We then isolate a subsample of these events containing no visible protons, and observe 93 ± 22⁢(stat) ±35⁢(syst) data events above prediction, corresponding to just above 2⁢𝜎 local significance, concentrated at shower energies below 600 MeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Spectroscopy of states in 136 Ba using the 138 Ba(p,t) reaction

Background: The 136 Ba isotope is the daughter nucleus in 136 Xe $ββ$ decay. It also lies in a shape transitional region of the nuclear chart, making it a suitable candidate to test a variety of nuclear models. Purpose: To obtain spectroscopic information on states in 136 Ba, which will allow a better understanding of its low-lying structure. These data may prove useful to constrain future 136 Xe → 136 Ba neutrinoless $ββ$ decay matrix element calculations. Methods: Here, a 138 Ba(p, t) reaction was used to populate states in 136 Ba up to approximately 4.6 MeV in excitation energy. The tritons were detected using a high-resolution Q3D magnetic spectrograph. A distorted wave Born approximation analysis was performed for the measured triton angular distributions. Results: 102 excited states in 136 Ba were observed, out of which 52 are reported for the first time. Definite spin-parity assignments are made for 26 newly observed states, while previously ambiguous assignments for ten other states are resolved. Together with other available data, the results are used to determine level densities in 136 Ba. These were compared with theory predictions, obtained using shell model calculations with Hamiltonians previously used for 136 Xe neutrinoless $ββ$ decay matrix element evaluations. Conclusions: The shell model predicted level densities agree reasonably well for the two Hamiltonians. However the results for theory and experiment are found to agree only at lower energies, diverging from one another for the higher lying states, with the discrepancy increasing with energy. This is presumably because of lower production cross sections for a majority of the higher-lying predicted states and the experimental limitations in resolving a large number of nearly degenerate states predicted by the theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

3D Continuous Forcing Dataset from 3D Constrained Variational Analysis at SGP

The continuous 3D large-scale forcing (VARANAL3D) data set derived from 3D constrained variational analysis (3DCVA) extends the conventional constrained variational analysis method by incorporating multiple sub-columns within the analysis domain. This advancement introduces spatial variability into the large-scale forcing fields, thereby enriching the data set’s applicability. The VARANAL3D data set spans from 2004 to 2018 and covers a region of 5˚×4.5˚ domain around the ARM SGP site. The analysis domain is divided into 10×9 sub-columns with 0.5˚ resolution. The 3D large-scale forcing data provides necessary variables to drive and evaluate single-column models (SCM), cloud-resolving models (CRM) ,and large-eddy simulations (LES), as well as information for testing model sensitivity to spatial variability of the large-scale forcing data, facilitating more rigorous testing and refinement of physical processes in SCM/CRM/LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Glacial to periglacial transition at the end of the last ice age in the subtropical semiarid Andes

Atmospheric warming and circulation reorganization at the end of the last ice age represent the most important climate change of the last 100,000 years and provide an opportunity to uncover how the southern subtropics cryosphere responded to strong changes in the global climate system. Extensive mapping and chronologic records on cryogenic landforms to better understand the association and interactions between glaciers and viscous creep of ice-rich permafrost landforms (rock glaciers) are widely missing in the region. In this paper, we reconstruct the geomorphic imprint of the Last Glacial Maximum (LGM) and the Termination I in the high Andes of the Río Limarí Basin (30–31°S) in the subtropical semiarid Andes of Chile. 74 new 10 Be surface exposure dating ages constrain the timing of glaciation, deglaciation, and glacial to periglacial transition. Glacial advances occurred first by 41.2 ± 0.6 – 35.0 ± 0.5 ka during Marine Isotope Stage 3, but probably earlier also; then, a second advance occurred during the global LGM between 24.2 ± 0.4 and 18.6 ± 0.2 ka. Deglaciation by 17.6 ± 0.2 ka left extensive hummocky moraines on the main valleys. Characteristic patterns of furrows and ridges typical of rock glaciers and solifluction superimposed on the LGM hummocky moraine indicate ice-rich permafrost in glacial deposits likely between 15.5 ± 0.3 and 13.6 ± 0.3 ka. We propose that moraines deposited by LGM debris-covered glaciers served as a niche for strong seasonal frost and permafrost creep, which substantially modified the original landforms. Finally, our results contribute to a better understanding of major transformations in an ice-rich high mountain area of the southern hemisphere where the interplay of temperature and precipitation changes drove glacial to periglacial transitions.

10Be surface exposure dating↗
DOE OSTI2021

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Herein this paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are benchmarked using a range of test cases and demonstrated using large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian (ILE) method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A simple Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. In the coupling strategy, a separate discretization of the fluid-structure interface is tethered to the volumetric solid mesh via stiff spring-like penalty forces. Our fluid-structure coupling scheme relies on an immersed interface method (IIM) for discrete geometries, which enables the accurate determination of both velocities and stresses along complex internal interfaces. Numerical methods for FSI can suffer from instabilities related to the added mass effect, but computational tests indicate that the methodology introduced here remains stable for selected test cases across a broad range of solid-fluid mass density ratios, including extremely small, nearly equal, equal, and large density ratios. Biomedical FSI demonstration cases include results obtained using this method to simulate the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and to model transport of blood clots in a patient-averaged anatomical model of the inferior vena cava.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on ultralight axion dark matter through galaxy cluster number counts

Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range 10 −32 eV to 10 −24 eV. We put upper bounds on the ultralight axion relic density Ω a in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around m a = 10 −27 eV with Ω a < 0.0035 and m a = 10 −26 eV with Ω a < 0.0079; both are at a 95% confidence level. When combined with cosmic microwave background probes, these bounds are tightened to Ω a < 0.0030 in the m a = 10 −27 eV mass bin and Ω a < 0.0058 in the m a = 10 −26 eV mass bin, with both at a 95% confidence level. This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.

cosmological parameters↗
DOE OSTI2024

Characterization and Optimization of Skipper CCDs for the SOAR Integral Field Spectrograph

We present results from the characterization and optimization of Skipper charge-coupled devices (CCDs) for use in a focal plane prototype for the Southern Astrophysical Research Integral Field Spectrograph (SIFS). We tested eight Skipper CCDs and selected six for SIFS based on performance results. The Skipper CCDs are 6k × 1k, 15 μm pixels, thick, fully depleted, p-channel devices that have been thinned to ~250 μm, backside processed, and treated with an anti-reflective coating. We demonstrate a single-sample readout noise of <4.3 e - rms pixel -1 in all amplifiers. We optimize the readout sequence timing to achieve a readout noise of 0.5 e - rms pixel -1 after 74 non-destructive measurements, which can be accomplished in a region covering 5% of the detector area in a readout time of <4 minutes. We demonstrate single-photon-counting in all 24 amplifiers (four amplifiers per detector) with a readnoise of σ N ~ 0.18 e - rms pixel -1 after N samp = 400 samples, and we constrain the degree of nonlinearity to be ≲1% at low signal levels (0 e - to 50 e - ). Clock-induced charge (CIC) remains an important issue when the Skipper CCD is configured to provide a large full-well capacity. We achieve a CIC rate of <1.45 × 10 -3 e - pixel -1 frame -1 for a full-well capacity of ~900 e - , which increases to a CIC rate of ~3 e - pixel -1 frame -1 for full-well capacities ~40,000-65,000 e - . We also perform conventional CCD characterization measurements such as charge transfer inefficiency (3.44 × 10 -7 on average), dark current (~2 × 10 -4 e - pixel -1 s -1 ), photon transfer curves, cosmetic defects (<0.45% “bad” pixels), and charge diffusion (point-spread function < 7.5 μm) to verify that these properties are consistent with expectations from conventional p-channel CCDs used for astronomy. Furthermore, we provide the first measurements of the brighter-fatter effect and absolute quantum efficiency (≳80% between 450 and 980 nm; ≳90% between 600 and 900 nm) using Skipper CCDs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Search for stealth supersymmetry in final states with two photons, jets, and low missing transverse momentum in proton-proton collisions at s = 13 TeV

The results of a search for stealth supersymmetry in final states with two photons and jets, targeting a phase space region with low missing transverse momentum ( p T miss ), are reported. The study is based on a sample of proton-proton collisions at s = 13 TeV collected by the CMS experiment, corresponding to an integrated luminosity of 138 fb - 1 . As LHC results continue to constrain the parameter space of the minimal supersymmetric standard model, the low p T miss regime is increasingly valuable to explore. To estimate the backgrounds due to standard model processes in such events, we apply corrections derived from simulation to an estimate based on a control selection in data. The results are interpreted in the context of simplified stealth supersymmetry models with gluino and squark pair production. The observed data are consistent with the standard model predictions, and gluino (squark) masses of up to 2150 (1850) GeV are excluded at the 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Proton deflectometry analysis in magnetized plasmas: Magnetic field reconstruction in one dimension

Proton deflectometry is used in magnetized high-energy-density plasmas to observe electromagnetic fields. Here we describe a reconstruction algorithm to recover the electromagnetic fields from proton fluence data in 1-D. The algorithm is verified against analytic solutions and applied to example data. Next, we study the role of source fluence uncertainty for 1-D reconstructions. We show that reconstruction boundary conditions can be used to constrain the source fluence profile and use this to develop a reconstruction using a specified pair of boundary conditions on the magnetic field. From these considerations, we experimentally demonstrate a hybrid mesh-fluence reconstruction technique where fields are reconstructed from fluence data in an interior region with boundary conditions supplied by direct mesh measurements at the boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Environmental controls on isolated convection during the Amazonian wet season

The Amazon rainforest is a vital component of the global climate system, influencing the hydrological cycle and tropical circulation. However, understanding and modeling the evolution of convection in this region remain a scientific challenge. Here, we assess the environmental conditions associated with shallow, congestus, and isolated deep convection days during the wet season (December to April), employing measurements from the Green Ocean Amazon 2014–2015 (GoAmazon2014/5) experiment and large-scale wind fields from the constrained variational analysis. Composites of deep days show moister than average conditions below 3 km early in the morning. Analyzing the water budget at the surface through observations only, we estimated the water vapor convergence term as a residual of the water balance closure. Convergence remains nearly zero during the deep days until early afternoon (13:00 LST), when it becomes a dominant factor in the water budget. At 14:00 LST, the deep days experience a robust upward large-scale vertical velocity, especially above 4 km, which supports the shallow-to-deep convective transition occurring around 16:00–17:00 LST. In contrast, shallow and congestus days exhibit drier pre-convective conditions, along with diurnal water vapor divergence and large-scale subsidence that extend from the surface to the lower free troposphere. Moreover, afternoon precipitation exhibits the strongest linear correlation (0.6) with large-scale vertical velocity, nearly double the magnitude observed for other environmental factors, even moisture, at different levels and periods of the day. Precipitation also exhibits a moderate increase with low-level wind shear, while upper-level shear has a relatively minor negative impact on convection.

Environmental Sciences & Ecology↗
DOE OSTI2023

Estimating Vadose Zone Flow Properties at the 100 K-East Soil Flushing Site Using ERT Monitoring Data: 2023 Interim Report - 100 KE Soil Property Estimation

In situ soil flushing is being using at the Hanford 100 K-East (100 KE) area to transport mobile chromium contamination in the vadose zone to the water table, where it can be collected and treat through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. In practice, it is infeasible to comprehensively monitor which regions of the vadose zone are being infiltrated through direct sampling of pore water. Consequently, there can be significant uncertainty about which regions of the subsurface have been treated, especially if hydrogeologic conditions are favorable for the development of unstable flows and preferred flow pathways through the vadose zone (Jarvis, Koestel, and Larsbo 2016). Current approaches for quantitative monitoring of soil flushing performance rely on contaminant concentration measurements collected from extractions wells. There is no quantitative information on the volume of flush water delivered to targeted regions of the vadose zone at the Hanford Site, leading to significant uncertainty regarding source term removal and long-term impacts to groundwater. If the subsurface hydrogeologic properties at the 100 KE Area were adequately known, qualitative metrics of soil flushing performance could be simulated, thereby negating expenses required to obtain quantitative performance information through borehole drilling/sampling. However, estimating in situ hydrogeologic properties has long proven elusive, due primarily to a lack of sufficient information to constrain heterogeneous property estimates to a useful degree of certainty. Estimating vadose zone hydrogeologic properties is particularly challenging due the dependence of hydraulic conductivity on saturation. This report describes progress toward a first-of-its-kind demonstration using surface time-lapse 3D electrical resistivity tomography (ERT) monitoring data to estimate the hydrogeologic properties that control flush water transport at the 100 KE soil flushing site. The ultimate objective is (1) to verify sufficient information exists in the ERT monitoring data to adequately resolve vadose zone hydraulic properties, and (2) generate a “digital twin” (i.e., a numerical simulator) that can be used to simulate the amount of flush water that has been delivered to each targeted region of the vadose zone, and thereby assess the efficacy of flush water delivery. Resulting performance estimates can be used in leu of comprehensive borehole drilling and direct sampling (or wellbore logging) that would otherwise be required to obtain the same information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Contrasting Regional Carbon Cycle Responses to Seasonal Climate Anomalies Across the East-West Divide of Temperate North America

Across temperate North America, interannual variability (IAV) in gross primary production (GPP) and net ecosystem exchange (NEE) and their relationship with environmental drivers are poorly understood. Here, we examine IAV in GPP and NEE and their relationship to environmental drivers using two state-of-the-science flux products: NEE constrained by surface and space-based atmospheric CO 2 measurements over 2010-2015 and satellite up-scaled GPP from FluxSat over 2001-2017. We show that the arid western half of temperate North America provides a larger contribution to IAV in GPP (104% of east) and NEE (127% of east) than the eastern half, in spite of smaller magnitude of annual mean GPP and NEE. This occurs because anomalies in western ecosystems are temporally coherent across the growing season leading to an amplification of GPP and NEE. In contrast, IAV in GPP and NEE in eastern ecosystems is dominated by seasonal compensation effects, associated with opposite responses to temperature anomalies in spring and summer. Terrestrial biosphere models in the MsTMIP ensemble generally capture these differences between eastern and western temperate North America, although there is considerable spread between models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Groundwater Volume Loss in Mexico City Constrained by InSAR and GRACE Observations and Mechanical Models

Abstract Groundwater withdrawal can cause localized and rapid poroelastic subsidence, spatially broad elastic uplift of low amplitude, and changes in the gravity field. Constraining groundwater loss in Mexico City, we analyze data from the Gravity Recovery and Climate Experiment and its follow‐on mission (GRACE/FO) and Synthetic Aperture Radar (SAR) Sentinel‐1A/B images between 2014 and 2021. GRACE/FO observations yield a groundwater loss of 0.85–3.87 km 3 /yr for a region of ∼300 × 600 km surrounding Mexico City. Using the high‐resolution interferometric SAR data set, we measure >35 cm/yr subsidence within the city and up to 2 cm/yr of uplift in nearby areas. Attributing the long‐term subsidence to poroelastic aquifer compaction and the long‐term uplift to elastic unloading, we apply respective models informed by local geology, yielding groundwater loss of 0.86–12.57 km 3 /yr. Our results suggest Mexico City aquifers have been depleting at faster rates since 2015, exacerbating the socioeconomic and health impacts of long‐term groundwater overdrafts.

58 GEOSCIENCES↗
DOE OSTI2025

Sensitivity of the Shallow‐To‐Deep Convective Transition to Moisture and Wind Shear in the Amazon

Deep convection is the primary influence on weather and climate in tropical regions. However, understanding and simulating the shallow‐to‐deep (STD) convective transition has long been challenging. Here, we conduct high‐resolution numerical simulations to assess the environmental controls on the evolution of isolated convection in the Amazon during the wet season. The large‐scale forcing derived through a constrained variational analysis approach for the GoAmazon2014/5 Experiment is used in the simulations. Through sensitivity experiments, we examine the relative importance of moisture and wind shear in controlling the shallow‐to‐deep convective transition for isolated convective events. Convection exhibits the greatest sensitivity to humidity within the lowest 1.5 km, where a 4 mm reduction in column water vapor nearly suppresses ice water formation on deep convective days. In contrast, a reduction in column water vapor in the free troposphere by a factor of two or more is necessary to produce a comparable impact on convection. Increasing low‐level wind speed from 6 to 9 m s -1 enhances afternoon deep convection, raising the cloud ice mixing ratio by approximately 25%. Conversely, upper‐level wind shear reveals the weakest correlation with daytime convection in our simulations. Our results help characterize the role of moisture and wind shear on the STD transition and our understanding of the underlying mechanisms.

environmental controls↗
DOE OSTI2023

New 26 P( p, γ ) 27 S Thermonuclear Reaction Rate and Its Astrophysical Implications in the rp -process

Accurate nuclear reaction rates for 26 P(p, γ) 27 S are pivotal for a comprehensive understanding of the rp-process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of 26 P(p, γ) 27 S because of the lack of nuclear mass and energy level structure information for 27 S. We reevaluate this reaction rate using the experimentally constrained 27 S mass, together with the shell model predicted level structure. It is found that the 26 P(p, γ) 27 S reaction rate is dominated by a direct capture reaction mechanism despite the presence of three resonances at E = 1.104, 1.597, and 1.777 MeV above the proton threshold in 27 S. The new rate is overall smaller than the other previous rates from the Hauser–Feshbach statistical model by at least 1 order of magnitude in the temperature range of X-ray burst interest. In addition, we consistently update the photodisintegration rate using the new 27 S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in the rp-process is explored by postprocessing nucleosynthesis calculations. The final abundance ratio of 27 S/ 26 P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show that the reaction path 26 P(p, γ) 27 S(β + ,ν) 27 P is not as important as previously thought for producing 27 P. The adoption of the new reaction rates for 26 P(p, γ) 27 S only reduces the final production of aluminum by 7.1% and has no discernible impact on the yield of other elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Further improvement and evaluation of nudging in the E3SM Atmosphere Model version 1 (EAMv1): simulations of the mean climate, weather events, and anthropogenic aerosol effects

Abstract. A previous study on the use of nudging in E3SM Atmosphere Model version 1 (EAMv1) had an unresolved issue; i.e., a simulation nudged to EAMv1's own meteorology showed non-negligible deviations from the free-running baseline simulation over some of the subtropical marine stratocumulus and trade cumulus regions. Here, we demonstrate that the deviations can be substantially reduced by (1) changing where the nudging tendency is calculated in the time integration loop of a nudged EAM simulation so as to improve consistency with the free-running baseline and (2) increasing the frequency of the constraining data so as to better capture strong sub-diurnal variations. The fact that modification (2) improves the climate representativeness of the nudged simulations has motivated us to investigate whether the use of newer reanalysis products with higher data frequency can help improve nudged hindcast simulations by better capturing the observed weather events. To answer this question, we present simulations conducted at EAMv1's standard horizontal resolution (approximately 1∘) with nudging towards 6-hourly ERA-Interim reanalysis or 6-hourly, 3-hourly, or hourly ERA5 reanalysis. These simulations are evaluated against the climatology of free-running EAMv1 simulations as well as reanalyses, satellite retrievals, and in situ measurements from the Atmospheric Radiation Measurement user facility. For the 1∘ EAMv1 simulations, we recommend using the relocated nudging tendency calculation and the ERA5 reanalysis at 3-hourly or higher frequency. Simulations used for estimating the anthropogenic aerosol effects often use nudging to help discern signal from noise. The sensitivity of such estimates to the configuration of nudging is investigated in EAMv1, again using the standard 1∘ horizontal resolution. We find that, when estimating the global mean effects, the frequency of constraining data has relatively small impacts, while the choice of nudged variables can change the results substantially. The nudging of air temperature (in addition to horizontal winds) has two non-negligible effects. First, when the constraining data come from reanalysis, the nudging-induced mean bias correction can cause significant changes in the simulated clouds and hence substantially different estimates of the aerosol effects. The impact of the mean bias correction on ice cloud formation has been noted in previous studies and is also seen in EAMv1. For applications like ours, where the preferred configurations of nudging are those capable of providing results consistent with the multi-year free-running simulations, the consequence of the mean bias correction is undesirable. The second important impact of temperature nudging is a significant suppression of adjustments to aerosol forcing, which also causes changes in the estimated aerosol effects. This effect can be seen in simulations nudged to either reanalysis or EAM's own meteorology. These results suggest that nudging horizontal winds but not temperature is a better choice for estimating the anthropogenic aerosol effects.

54 ENVIRONMENTAL SCIENCES↗