Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neal”

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 1,261 records · Page 70

Precessing Jet and Large Dust Grains in the V380 Ori NE Star-forming Region

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J=1→0 lines, and dense cores in L1641 were observed in the 2.0–0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of ∼35 km s{sup −1} and a precession period of ∼1600 years. The CO outflow length gives a dynamical timescale of ∼6300 years, and the protostar must be extremely young. The inferred binary separation of 6–70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of β = 0.3–2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of β = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices.

79 ASTRONOMY AND ASTROPHYSICS↗

The Preprocessing of Galaxies in the Early Stages of Cluster Formation in Abell 1882 at z = 0.139

A rare opportunity to distinguish between internal and environmental effects on galaxy evolution is afforded by "SuperGroups," systems that are rich and massive, but include several comparably rich substructures, surrounded by filaments. We present here a multiwavelength photometric and spectroscopic study of the galaxy population in the SuperGroup Abell 1882 (A1882) at z = 0.139, combining new data from the MMT and Hectospec with archival results from the Galaxy And Mass Assembly survey, the Sloan Digital Sky Survey, the Nasa/IPAC Extragalactic Database, the Gemini Multi-Object Spectrograph, and the Galaxy Evolution Explorer. These provide spectroscopic classifications for 526 member galaxies, across wide ranges of local density and velocity dispersion. We identify three prominent filaments along which galaxies seem to be entering the SuperGroup (mostly in E–W directions). A1882 has a well-populated red sequence, containing most galaxies with stellar mass >10 10.5 M Sun , and a pronounced color–density relation even within its substructures. Thus, galaxy evolution responds to the external environment as strongly in these unrelaxed systems as we find in rich and relaxed clusters. From these data, local density remains the primary factor, with a secondary role for distance from the inferred center of the entire structure's potential well. The effects on star formation, as traced by optical and near-UV colors, depend on galaxy mass. We see changes in lower-mass galaxies (M < 10 10.5 M Sun ) at four times the virial radius of major substructures, while the more massive near-UV Green Valley galaxies show low levels of star formation within two virial radii. The suppression of star formation ("quenching") occurs in the infall regions of these structures even before the galaxies enter the denser group environment.

79 ASTRONOMY AND ASTROPHYSICS↗

FORMING CHONDRITES IN A SOLAR NEBULA WITH MAGNETICALLY INDUCED TURBULENCE

Chondritic meteorites provide valuable opportunities to investigate the origins of the solar system. We explore impact jetting as a mechanism of chondrule formation and subsequent pebble accretion as a mechanism of accreting chondrules onto parent bodies of chondrites, and investigate how these two processes can account for the currently available meteoritic data. We find that when the solar nebula is ≤5 times more massive than the minimum-mass solar nebula at a ≃ 2–3 au and parent bodies of chondrites are ≤10{sup 24} g (≤500 km in radius) in the solar nebula, impact jetting and subsequent pebble accretion can reproduce a number of properties of the meteoritic data. The properties include the present asteroid belt mass, the formation timescale of chondrules, and the magnetic field strength of the nebula derived from chondrules in Semarkona. Since this scenario requires a first generation of planetesimals that trigger impact jetting and serve as parent bodies to accrete chondrules, the upper limit of parent bodies’ masses leads to the following implications: primordial asteroids that were originally ≥10{sup 24} g in mass were unlikely to contain chondrules, while less massive primordial asteroids likely had a chondrule-rich surface layer. The scenario developed from impact jetting and pebble accretion can therefore provide new insights into the origins of the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

DETECTION OF H{sub 2}O AND EVIDENCE FOR TiO/VO IN AN ULTRA-HOT EXOPLANET ATMOSPHERE

We present a primary transit observation for the ultra-hot ( T {sub eq} ∼ 2400 K) gas giant expolanet WASP-121b, made using the Hubble Space Telescope Wide Field Camera 3 in spectroscopic mode across the 1.12–1.64 μ m wavelength range. The 1.4 μ m water absorption band is detected at high confidence (5.4 σ ) in the planetary atmosphere. We also reanalyze ground-based photometric light curves taken in the B , r ′, and z ′ filters. Significantly deeper transits are measured in these optical bandpasses relative to the near-infrared wavelengths. We conclude that scattering by high-altitude haze alone is unlikely to account for this difference and instead interpret it as evidence for titanium oxide and vanadium oxide absorption. Enhanced opacity is also inferred across the 1.12–1.3 μ m wavelength range, possibly due to iron hydride absorption. If confirmed, WASP-121b will be the first exoplanet with titanium oxide, vanadium oxide, and iron hydride detected in transmission. The latter are important species in M/L dwarfs and their presence is likely to have a significant effect on the overall physics and chemistry of the atmosphere, including the production of a strong thermal inversion.

79 ASTRONOMY AND ASTROPHYSICS↗

The Very Low Albedo of WASP-12b from Spectral Eclipse Observations with Hubble

We present an optical eclipse observation of the hot Jupiter WASP-12b using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. These spectra allow us to place an upper limit of A{sub g}<0.064 (97.5% confidence level) on the planet’s white light geometric albedo across 290–570 nm. Using six wavelength bins across the same wavelength range also produces stringent limits on the geometric albedo for all bins. However, our uncertainties in eclipse depth are ∼40% greater than the Poisson limit and may be limited by the intrinsic variability of the Sun-like host star—the solar luminosity is known to vary at the 10{sup −4} level on a timescale of minutes. We use our eclipse depth limits to test two previously suggested atmospheric models for this planet: Mie scattering from an aluminum-oxide haze or cloud-free Rayleigh scattering. Our stringent nondetection rules out both models and is consistent with thermal emission plus weak Rayleigh scattering from atomic hydrogen and helium. Our results are in stark contrast with those for the much cooler HD 189733b, the only other hot Jupiter with spectrally resolved reflected light observations; those data showed an increase in albedo with decreasing wavelength. The fact that the first two exoplanets with optical albedo spectra exhibit significant differences demonstrates the importance of spectrally resolved reflected light observations and highlights the great diversity among hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Hammermill Performance Through Screen Selection and Hammer Design

Background: Mechanical preprocessing, which includes particle size reduction and mechanical separation, is one of the primary operations in the feedstock supply system for a lignocellulosic biorefinery. It is the means by which raw biomass from the field or forest is mechanically transformed into an on-spec feedstock with characteristics better suited for the fuel conversion process. Results: This work provides a general overview of the objectives and methodologies of mechanical preprocessing and then presents experimental results illustrating (1) improved size reduction via optimization of hammer mill configuration, (2) improved size reduction via pneumatic-assisted hammer milling, and (3) improved control of particle size and particle size distribution through proper selection of grinder process parameters. Conclusion: Optimal grinder configuration for maximal process throughput and efficiency is strongly dependent on feedstock type and properties, such moisture content. Tests conducted using a HG200 hammer grinder indicate that increasing the tip speed, optimizing hammer geometry, and adding pneumatic assist can increase grinder throughput as much as 400%.

09 BIOMASS FUELS↗

Degradation Dynamics of Perovskite Solar Cells Under Fixed Reverse Current Injection

Previous studies of reverse-bias stability in perovskite solar cells have focused primarily on voltage controlled reverse-bias tests. Here we instead present an investigation of perovskite solar cell degradation under well-defined, constant reverse-current stress. We show that the choice of hole-transport layer dictates the dominant degradation pathway: cells using thick poly(triphenylamine) (PTAA) layers with better indium-doped tin oxide (ITO) coverage can tolerate high reverse bias but quickly undergo catastrophic breakdown under fixed reverse current near their one-sun maximum power-point. In contrast, cells modified with the phosphonic-acid interface layer MeO-2PACz, with poorer ITO coverage compared to PTAA, exhibit soft, gradual, and largely recoverable degradation, regardless of the shading conditions. For MeO-2PACz devices, degradation increases with both current magnitude and duration. Importantly, when normalized by injected charge (current times duration), lower currents applied over longer times cause more severe degradation than higher currents over shorter periods. Combining electrical measurements with spatially resolved photoluminescence imaging, we argue against shunt formation and instead support an ion- and charge-mediated interfacial electrochemical degradation mode.

Applied Physics (physics.app-ph)↗

Optimizing cloud motion estimation on the edge with phase correlation and optical flow

Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.

54 ENVIRONMENTAL SCIENCES↗

Practice and philosophy of climate model tuning across six US modeling centers

Model calibration (or tuning ) is a necessary part of developing and testing coupled ocean–atmosphere climate models regardless of their main scientific purpose. There is an increasing recognition that this process needs to become more transparent for both users of climate model output and other developers. Knowing how and why climate models are tuned and which targets are used is essential to avoiding possible misattributions of skillful predictions to data accommodation and vice versa. This paper describes the approach and practice of model tuning for the six major US climate modeling centers. While details differ among groups in terms of scientific missions, tuning targets, and tunable parameters, there is a core commonality of approaches. Furthermore, practices differ significantly on some key aspects, in particular, in the use of initialized forecast analyses as a tool, the explicit use of the historical transient record, and the use of the present-day radiative imbalance vs. the implied balance in the preindustrial era as a target.

58 GEOSCIENCES↗

The path to CAM6: coupled simulations with CAM5.4 and CAM5.5

This paper documents coupled simulations of two developmental versions of the Community Atmosphere Model (CAM) towards CAM6. The configuration called CAM5.4 introduces new microphysics, aerosol, and ice nucleation changes, among others to CAM. The CAM5.5 configuration represents a more radical departure, as it uses an assumed probability density function (PDF)-based unified cloud parameterization to replace the turbulence, shallow convection, and warm cloud macrophysics in CAM. This assumed PDF method has been widely used in the last decade in atmosphere-only climate simulations but has never been documented in coupled mode. Here, we compare the simulated coupled climates of CAM5.4 and CAM5.5 and compare them to the control coupled simulation produced by CAM5.3. We find that CAM5.5 has lower cloud forcing biases when compared to the control simulations. Improvements are also seen in the simulated amplitude of the Niño-3.4 index, an improved representation of the diurnal cycle of precipitation, subtropical surface wind stresses, and double Intertropical Convergence Zone biases. Degradations are seen in Amazon precipitation as well as slightly colder sea surface temperatures and thinner Arctic sea ice. Simulation of the 20th century results in a credible simulation that ends slightly colder than the control coupled simulation. The authors find this is due to aerosol indirect effects that are slightly stronger in the new version of the model and propose a solution to ameliorate this. Altogether, in these early coupled simulations, CAM5.5 produces a credible climate that is appropriate for science applications and is ready for integration into the National Center for Atmospheric Research's (NCAR's) next-generation climate model.

54 ENVIRONMENTAL SCIENCES↗

Regionally refined test bed in E3SM atmosphere model version 1 (EAMv1) and applications for high-resolution modelling

Climate simulation with more accurate process-level representation at finer resolutions is a pressing necessity in order to provide actionable information to policy-makers regarding extreme events in a changing climate. Computational limitation is a major obstacle for building, and running high-resolution (HR, here 0.25° average grid spacing at the equator) models (HRM). A more affordable path to HRM is to use a global regionally refined model (RRM), which only simulates a portion of the globe at HR while the remaining is at low-resolution (LR, 1°). In this study, we compare the Energy Exascale Earth System Model (E3SM) atmosphere model version 1 (EAMv1) RRM with the HR mesh over the contiguous United States (CONUS) to its corresponding globally uniform LR and HR configurations, as well as to observations and reanalysis data. The RRM has a significantly reduced computational cost (roughly proportional to the HR mesh size) relative to the globally uniform HRM. Over the CONUS, we evaluate the simulation of important dynamical and physical quantities as well as various precipitation measures. Differences between the RRM and HRM over the HR region are predominantly small, demonstrating that the RRM reproduces both well- and poorly simulated behaviours of the HRM over the CONUS. Further analysis based on RRM simulations with the LR vs. HR model parameters reveals that RRM performance is greatly influenced by the different parameter choices used in the LR and HR EAMv1. This is a result of the poor scale-aware behaviour of physical parameterizations, especially for variables influencing sub-grid scale physical processes. RRM can serve as a useful framework to test physics schemes across a range of scales, leading to improved consistency in future E3SM versions. Applying nudging-to-observations techniques within the RRM framework also demonstrates significant advantages over a free-running configuration for use as a testbed, and as such represents an efficient and more robust physics testbed capability. Our conclusions provide additional confirmatory evidence that the RRM is an efficient and effective approach for HRM development and hydrologic research.

54 ENVIRONMENTAL SCIENCES↗

GCAM-USA v5.3_water_dispatch: integrated modeling of subnational US energy, water, and land systems within a global framework

Abstract. This paper describes GCAM-USA v5.3_water_dispatch, an open-source model that represents key interactions across economic, energy, water, and land systems in a consistent global framework with subnational detail in the United States. GCAM-USA divides the world into 31 geopolitical regions outside the United States (US) and represents the US economy and energy systems in 51 state-level regions (50 states plus the District of Columbia). The model also includes 235 water basins and 384 land use regions, and 23 of each fall at least partially within the United States. GCAM-USA offers a level of process and temporal resolution rare for models of its class and scope, including detailed subnational representation of US water demands and supplies and sub-annual operations (day and night for each month) in the US electric power sector. GCAM-USA can be used to explore how changes in socioeconomic drivers, technological progress, or policy impact demands for (and production of) energy, water, and crops at a subnational level in the United States while maintaining consistency with broader national and international conditions. This paper describes GCAM-USA's structure, inputs, and outputs, with emphasis on new model features. Four illustrative scenarios encompassing varying socioeconomic and energy system futures are used to explore subnational changes in energy, water, and land use outcomes. We conclude with information about how public users can access the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

54 ENVIRONMENTAL SCIENCES↗

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

54 ENVIRONMENTAL SCIENCES↗

Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1

Abstract. Groundwater plays a key role in meeting water demands, supplying over 40 % of irrigation water globally, with this role likely to grow as water demands and surface water variability increase. A better understanding of the future role of groundwater in meeting sectoral demands requires an integrated hydro-economic evaluation of its cost and availability. Yet substantial gaps remain in our knowledge and modeling capabilities related to groundwater availability, recharge, feasible locations for extraction, extractable volumes, and associated extraction costs, which are essential for large-scale analyses of integrated human–water system scenarios, particularly at the global scale. To address these needs, we developed Superwell, a physics-based groundwater extraction and cost accounting model that operates at sub-annual temporal and at the coarsest 0.5° (≈50 km × 50 km) gridded spatial resolution with global coverage. The model produces location-specific groundwater supply–cost curves that provide the levelized cost to access different quantities of available groundwater. The inputs to Superwell include recent high-resolution hydrogeologic datasets of permeability, porosity, aquifer thickness, depth to water table, recharge, and hydrogeological complexity zones. It also accounts for well capital and maintenance costs, as well as the energy costs required to lift water to the surface. The model employs a Theis-based scheme coupled with an amortization-based cost accounting formulation to simulate groundwater extraction and quantify the cost of groundwater pumping. The result is a spatiotemporally flexible, physically realistic, economics-based model that produces groundwater supply–cost curves. We show examples of these supply–cost curves and the insights that can be derived from them across a set of scenarios designed to explore model outcomes. The supply–cost curves produced by the model show that most (90 %) nonrenewable groundwater in storage globally is extractable at costs lower than USD 0.57 m−3, while half of the volume remains extractable at under USD 0.108 m−3. The global unit cost is estimated to range from a minimum of USD 0.004 m−3 to a maximum of USD 3.971 m−3. We also demonstrate and discuss examples of how these cost curves could be used by linking Superwell's outputs with other models to explore coupled human–environmental system challenges, such as water resources planning and management, or broader analyses of multisectoral feedbacks.

Global Change Analysis Model (GCAM)↗

Moving beyond post hoc explainable artificial intelligence: a perspective paper on lessons learned from dynamical climate modeling

AI models are criticized as being black boxes, potentially subjecting climate science to greater uncertainty. Explainable artificial intelligence (XAI) has been proposed to probe AI models and increase trust. In this review and perspective paper, we suggest that, in addition to using XAI methods, AI researchers in climate science can learn from past successes in the development of physics-based dynamical climate models. Dynamical models are complex but have gained trust because their successes and failures can sometimes be attributed to specific components or sub-models, such as when model bias is explained by pointing to a particular parameterization. We propose three types of understanding as a basis to evaluate trust in dynamical and AI models alike: (1) instrumental understanding, which is obtained when a model has passed a functional test; (2) statistical understanding, obtained when researchers can make sense of the modeling results using statistical techniques to identify input–output relationships; and (3) component-level understanding, which refers to modelers' ability to point to specific model components or parts in the model architecture as the culprit for erratic model behaviors or as the crucial reason why the model functions well. We demonstrate how component-level understanding has been sought and achieved via climate model intercomparison projects over the past several decades. Such component-level understanding routinely leads to model improvements and may also serve as a template for thinking about AI-driven climate science. Currently, XAI methods can help explain the behaviors of AI models by focusing on the mapping between input and output, thereby increasing the statistical understanding of AI models. Yet, to further increase our understanding of AI models, we will have to build AI models that have interpretable components amenable to component-level understanding. We give recent examples from the AI climate science literature to highlight some recent, albeit limited, successes in achieving component-level understanding and thereby explaining model behavior. The merit of such interpretable AI models is that they serve as a stronger basis for trust in climate modeling and, by extension, downstream uses of climate model data.

54 ENVIRONMENTAL SCIENCES↗

Toward a Minimal Representation of Aerosols in Climate Models: Description and Evaluation in the Community Atmosphere Model CAM5

A modal aerosol module (MAM) has been developed for the Community Atmosphere Model version 5 (CAM5), the atmospheric component of the Community Earth System Model version 1 (CESM1). MAM is capable of simulating the aerosol size distribution and both internal and external mixing between aerosol components, treating numerous complicated aerosol processes and aerosol physical, chemical and optical properties in a physically based manner. Two MAM versions were developed: a more complete version with seven-lognormal modes (MAM7), and a three-lognormal mode version (MAM3) for the purpose of long-term (decades to centuries) simulations. Major approximations in MAM3 include assuming immediate mixing of primary organic matter (POM) and black carbon (BC) with other aerosol components, merging of the MAM7 fine dust and fine sea salt modes into the accumulation mode, merging of the MAM7 coarse dust and coarse sea salt modes into the single coarse mode, and neglecting the explicit treatment of ammonia and ammonium cycles. Simulated sulfate and secondary organic aerosol (SOA) mass concentrations are remarkably similar between MAM3 and MAM7 as most ({approx}90%) of these aerosol species are in the accumulation mode. Differences of POM and BC concentrations between MAM3 and MAM7 are also small (mostly within 10%) because of the assumed hygroscopic nature of POM, so that freshly emitted POM and BC are wet-removed before mixing internally with soluble aerosol species. Sensitivity tests with the POM assumed to be hydrophobic and with slower aging process increase the POM and BC concentrations, especially at high latitudes (by several times). The mineral dust global burden differs by 10% and sea salt burden by 30-40% between MAM3 and MAM7 mainly due to the different size ranges for dust and sea salt modes and different standard deviations of log-normal size distribution for sea salt modes between MAM3 and MAM7. The model is able to qualitatively capture the observed geographical and temporal variations of aerosol mass and number concentrations, size distributions, and aerosol optical properties. However, there are noticeable biases, e.g., simulated sulfate and mineral dust concentrations at surface over the oceans are too low. Simulated BC concentrations are significant low in the Arctic. There is a low bias in modeled aerosol optical depth on the global scale, especially in the developing counties. There biases in aerosol simulations clearly indicate the need for improvements of aerosol processes (e.g., emission fluxes of anthropogenic aerosols and precursor gases in developing countries, boundary layer nucleation) and properties (e.g., primary aerosol emission size, POM hygroscopicity). In addition the critical role of cloud properties (e.g., liquid water content, cloud fraction) responsible for the wet scavenging of aerosol is highlighted.

54 ENVIRONMENTAL SCIENCES↗

A new ensemble-based consistency test for the Community Earth System Model (pyCECT v1.0)

Abstract. Climate simulation codes, such as the Community Earth System Model (CESM), are especially complex and continually evolving. Their ongoing state of development requires frequent software verification in the form of quality assurance to both preserve the quality of the code and instill model confidence. To formalize and simplify this previously subjective and computationally expensive aspect of the verification process, we have developed a new tool for evaluating climate consistency. Because an ensemble of simulations allows us to gauge the natural variability of the model's climate, our new tool uses an ensemble approach for consistency testing. In particular, an ensemble of CESM climate runs is created, from which we obtain a statistical distribution that can be used to determine whether a new climate run is statistically distinguishable from the original ensemble. The CESM ensemble consistency test, referred to as CESM-ECT, is objective in nature and accessible to CESM developers and users. The tool has proven its utility in detecting errors in software and hardware environments and providing rapid feedback to model developers.

58 GEOSCIENCES↗