Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Observed Boundary Layer Controls on Shallow Cumulus at the ARM Southern Great Plains Site

Abstract The boundary layer controls on shallow cumulus (ShCu) convection are examined using a suite of remote and in situ sensors at ARM Southern Great Plains (SGP). A key instrument in the study is a Doppler lidar that measures vertical velocity in the CBL and along cloud base. Using a sample of 138 ShCu days, the composite structure of the ShCu CBL is examined, revealing increased vertical velocity (VV) variance during periods of medium cloud cover and higher VV skewness on ShCu days than on clear-sky days. The subcloud circulations of 1791 individual cumuli are also examined. From these data, we show that cloud-base updrafts, normalized by convective velocity, vary as a function of updraft width normalized by CBL depth. It is also found that 63% of clouds have positive cloud-base mass flux and are linked to coherent updrafts extending over the depth of the CBL. In contrast, negative mass flux clouds lack coherent subcloud updrafts. Both sets of clouds possess narrow downdrafts extending from the cloud edges into the subcloud layer. These downdrafts are also present adjacent to cloud-free updrafts, suggesting they are mechanical in origin. The cloud-base updraft data are subsequently combined with observations of convective inhibition to form dimensionless “cloud inhibition” (CI) parameters. Updraft fraction and liquid water path are shown to vary inversely with CI, a finding consistent with CIN-based closures used in convective parameterizations. However, we also demonstrate a limited link between CBL vertical velocity variance and cloud-base updrafts, suggesting that additional factors, including updraft width, are necessary predictors for cloud-base updrafts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Marine Boundary Layer Clouds Associated with Coastally Trapped Disturbances: Observations and Model Simulations

Modeling marine low clouds and fog in coastal environments remains an outstanding challenge due to the inherently complex ocean–land–atmosphere system. This is especially important in the context of global circulation models due to the profound radiative impact of these clouds. This study utilizes aircraft and satellite measurements, in addition to numerical simulations using the Weather Research and Forecasting (WRF) Model, to examine three well-observed coastally trapped disturbance (CTD) events from June 2006, July 2011, and July 2015. Cloud water-soluble ionic and elemental composition analyses conducted for two of the CTD cases indicate that anthropogenic aerosol sources may impact CTD cloud decks due to synoptic-scale patterns associated with CTD initiation. In general, the dynamics and thermodynamics of the CTD systems are well represented and are relatively insensitive to the choice of physics parameterizations; however, a set of WRF simulations suggests that the treatment of model physics strongly influences CTD cloud field evolution. Specifically, cloud liquid water path (LWP) is highly sensitive to the choice of the planetary boundary layer (PBL) scheme; in many instances, the PBL scheme affects cloud extent and LWP values as much as or more than the microphysics scheme. Results suggest that differences in the treatment of entrainment and vertical mixing in the Yonsei University (nonlocal) and Mellor–Yamada–Janjić (local) PBL schemes may play a significant role. The impact of using different driving models—namely, the North American Mesoscale Forecast System (NAM) 12-km analysis and the NCEP North American Regional Reanalysis (NARR) 32-km products—is also investigated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Convective Transition Statistics over Tropical Oceans for Climate Model Diagnostics: GCM Evaluation

To assess deep convective parameterizations in a variety of GCMs and examine the fast-time-scale convective transition, a set of statistics characterizing the pickup of precipitation as a function of column water vapor (CWV), PDFs and joint PDFs of CWV and precipitation, and the dependence of the moisture–precipitation relation on tropospheric temperature is evaluated using the hourly output of two versions of the GFDL Atmospheric Model, version 4 (AM4), NCAR CAM5 and superparameterized CAM (SPCAM). The 6-hourly output from the MJO Task Force (MJOTF)/GEWEX Atmospheric System Study (GASS) project is also analyzed. Contrasting statistics produced from individual models that primarily differ in representations of moist convection suggest that convective transition statistics can substantially distinguish differences in convective representation and its interaction with the large-scale flow, while models that differ only in spatial–temporal resolution, microphysics, or ocean–atmosphere coupling result in similar statistics. Most of the models simulate some version of the observed sharp increase in precipitation as CWV exceeds a critical value, as well as that convective onset occurs at higher CWV but at lower column RH as temperature increases. While some models quantitatively capture these observed features and associated probability distributions, considerable intermodel spread and departures from observations in various aspects of the precipitation–CWV relationship are noted. For instance, in many of the models, the transition from the low-CWV, nonprecipitating regime to the moist regime for CWV around and above critical is less abrupt than in observations. Additionally, some models overproduce drizzle at low CWV, and some require CWV higher than observed for strong precipitation. For many of the models, it is particularly challenging to simulate the probability distributions of CWV at high temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Triple-Moment Representation of Ice in the Predicted Particle Properties (P3) Microphysics Scheme

In the original Predicted Particle Properties (P3) bulk microphysics scheme, all ice-phase hydrometeors are represented by one or more “free” ice categories, where the physical properties evolve smoothly through changes to the four prognostic variables (per category), and with a two-moment representation of the particle size distribution. As such, the spectral dispersion cannot evolve independently, which thus results in limitations in representation of ice—in particular, hail—due to necessary constraints in the scheme to prevent excessive gravitational size sorting. To overcome this, P3 has been modified to include a three-moment representation of the size distribution of each ice category through the addition of a fifth prognostic variable, the sixth moment of the size distribution. The details of the three-moment ice parameterization in P3 are provided. The behavior of the modified scheme, with the single-ice-category configuration, is illustrated through simulations in a simple 1D kinematic model framework as well as with near large-eddy-resolving (250-m grid spacing) 3D simulations of a hail-producing supercell. It is shown that the three-moment ice configuration controls size sorting in a physically based way and leads to an improved capacity to simulate large, heavily rimed ice (hail), including mean and maximum sizes and reflectivity, and thus an overall improvement in the representation of ice-phase particles in the P3 scheme.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

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

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

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Cloud Phase and Relative Humidity Distributions over the Southern Ocean in Austral Summer Based on In Situ Observations and CAM5 Simulations

Cloud phase and relative humidity (RH) distributions at -67° to 0°C over the Southern Ocean during austral summer are compared between in situ airborne observations and global climate simulations. A scale-aware comparison is conducted using horizontally averaged observations from 0.1 to 50 km. Cloud phase frequencies, RH distributions, and liquid mass fraction are found to be less affected by horizontal resolutions than liquid and ice water content (LWC and IWC, respectively), liquid and ice number concentrations (Nc liq and Nc ice , respectively), and ice supersaturation (ISS) frequency. At -10° to 0°C, observations show 27%–34% and 17%–37% of liquid and mixed phases, while simulations show 60%–70% and 3%–4%, respectively. Simulations overestimate (underestimate) LWC and Nc liq in liquid (mixed) phase, overestimate Nc ice in mixed phase, underestimate IWC in ice and mixed phases, and underestimate (overestimate) liquid mass fraction below (above) -5°C, indicating that observational constraints are needed for different cloud phases. RH frequently occurs at liquid saturation in liquid and mixed phases for all datasets, yet the observed RH in ice phase can deviate from liquid saturation by up to 20%–40% at -20° to 0°C, indicating that the model assumption of liquid saturation for coexisting ice and liquid is inaccurate for low liquid mass fractions (<0.1). Simulations lack RH variability for partial cloud fractions (0.1–0.9) and underestimate (overestimate) ISS frequency for cloud fraction <0.1 (≥0.6), implying that improving RH subgrid-scale parameterizations may be a viable path to account for small-scale processes that affect RH and cloud phase heterogeneities. Lastly, two sets of simulations (nudged and free-running) show very similar results (except for ISS frequency) regardless of sample sizes, corroborating the statistical robustness of the model–observation comparisons.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Effect of Tropical Nonconvective Condensation on Uncertainty in Modeled Projections of Rainfall

Abstract We find that part of the uncertainty in the amplitude and pattern of the modeled precipitation response to CO2 forcing traces to tropical condensation not directly involved with parameterized convection. The fraction of tropical rainfall associated with large-scale condensation can vary from a few percent to well over half depending on model details and parameter settings. In turn, because of the coupling between condensation and tropical circulation, the different ways model assumptions affect the large-scale rainfall fraction also affect the patterns of the response within individual models. In two single-model ensembles based on the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM), versions 3.1 and 5.3, we find strong correlations between the fraction of tropical large-scale rain and both climatological rainfall and circulation and the response to CO2 forcing. While the effects of an increasing tropical large-scale rain fraction are opposite in some ways in the two ensembles—for example, the Hadley circulation weakens with the large-scale rainfall fraction in the CAM3.1 ensemble while strengthening in the CAM5.3 ensemble—we can nonetheless understand these different effects in terms of the relationship between latent heating and circulation, and we propose explanations for each ensemble. We compare these results with data from phase 5 of the Coupled Model Intercomparison Project (CMIP5), for which some of the same patterns hold. Given the importance of this partitioning, there is a need for constraining this source of uncertainty using observations. However, since a “large-scale rainfall fraction” is a modeling construct, it is not clear how observations may be used to test various modeling assumptions determining this fraction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating Diurnal and Semi-Diurnal Cycle of Precipitation in CMIP6 Models Using Satellite- and Ground-Based Observations

We report the diurnal and semi-diurnal cycle of precipitation simulated from CMIP6 models during 1996-2005 are evaluated globally between 60°S and 60°N, as well as at ten selected locations representing three categories of diurnal cycle of precipitation: (1) afternoon precipitation over land, (2) early morning precipitation over ocean, and (3) nocturnal precipitation over land. Three satellite-based and two ground-based rainfall products are used to evaluate the climate models. Globally, the ensemble mean of CMIP6 models shows a diurnal phase of 3 to 4 hours earlier over land and 1 to 2 hours earlier over ocean, when compared with the latest satellite products. These biases are in line with what were found in previous versions of climate models but reduced compared to the CMIP5 ensemble mean. Analysis at the selected locations complimented with in-situ measurements further reinforces these results. Several CMIP6 models have shown a significant improvement in the diurnal cycle of precipitation compared to their CMIP5 counterparts, notably on delaying afternoon precipitation over land. This can be attributed to the use of more sophisticated convective parameterizations. Most models are still unable to capture the nocturnal peak associated with elevated convection and propagating mesoscale convective systems, with a few exceptions that allow convection to be initiated above the boundary layer to capture nocturnal elevated convection. We also quantify an encouraging consistency between the satellite- and ground-based precipitation measurements despite differing spatiotemporal resolutions and sampling periods, which provides confidence in using them to evaluate the diurnal and semi-diurnal cycle of precipitation in climate models.

58 GEOSCIENCES↗
DOE OSTI2022

Linking Large-Scale Double-ITCZ Bias to Local-Scale Drizzling Bias in Climate Models

Tropical precipitation in climate models presents significant biases in both the large-scale pattern (i.e., double intertropical convergence zone bias) and local-scale characteristics (i.e., drizzling bias with too frequent drizzle/convection and reduced occurrences of no and heavy precipitation). By untangling the coupled system and analyzing the biases in precipitation, cloud, and radiation, this study shows that local-scale drizzling bias in atmospheric models can lead to large-scale double-ITCZ bias in coupled models by inducing convective-regime-dependent biases in precipitation and cloud radiative effects (CRE). The double-ITCZ bias consists of a hemispherically asymmetric component that arises from the asymmetric SST bias and a nearly symmetric component that exists in atmospheric models without the SST bias. By increasing light rain but reducing heavy rain, local-scale drizzling bias induces positive (negative) precipitation bias in the moderate (strong) convective regime, leading to the nearly symmetric wet bias in atmospheric models. By affecting the cloud profile, local-scale drizzling bias induces positive (negative) CRE bias in the stratocumulus (convective) regime in atmospheric models. Because the stratocumulus (convective) region is climatologically more pronounced in the southern (northern) tropics, the CRE bias is deemed to be hemispherically asymmetric and drives warm and wet (cold and dry) biases in the southern (northern) tropics when coupled to ocean. Our results suggest that correcting local-scale drizzling bias is critical for fixing large-scale double-ITCZ bias. Finally, the drizzling and double-ITCZ biases are not alleviated in models with mesoscale (0.25°–0.5°) or even storm-resolving (~3 km) resolution, implying that either large-eddy simulation or fundamental improvement in small-scale subgrid parameterizations is needed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

A New Criterion to Improve Operational Drizzle Detection with Ground-Based Remote Sensing

Light shallow precipitation in the form of drizzle is one of the mechanisms for liquid water removal, affecting cloud lifetime and boundary layer dynamics and thermodynamics. The early formation of drizzle drops is of particular interest for quantifying aerosol–cloud–precipitation interactions. In models, drizzle initiation is represented by the autoconversion, that is, the conversion of liquid water from a cloud liquid water category (where particle sedimentation is ignored) into a precipitating liquid water category. Various autoconversion parameterizations have been proposed in recent years, but their evaluation is challenging due to the lack of proper observations of drizzle development in the cloud. This work presents a new algorithm for Classification of Drizzle Stages (CLADS). CLADS is based on the skewness of the Ka-band radar Doppler spectrum. Skewness is sensitive to the drizzle growth in the cloud: the observed Gaussian Doppler spectrum has skewness zero when only cloud droplets are present without any significant fall velocity. Defining downward velocities positive, skewness turns positive when embryonic drizzle forms and becomes negative when drizzle starts to dominate the spectrum. CLADS identifies spatially coherent structures of positive, zero, and negative skewness in space and time corresponding to drizzle seeding, drizzle growth/nondrizzle, and drizzle mature, respectively. We test CLADS on case studies from the Jülich Observatory for Cloud Evolution Core Facility (JOYCE-CF) and the Barbados Cloud Observatory (BCO) to quantitatively estimate the benefits of CLADS compared to the standard Cloudnet target categorization algorithm. We suggest that CLADS can provide additional observational constraints for understanding the processes related to drizzle formation better.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Simulated Polarimetric Fields of Ice Vapor Growth Using the Adaptive Habit Model. Part II: A Case Study from the FROST Experiment

Abstract A new adaptive habit model (AHM) grows ice crystals through vapor deposition while evolving ice particle properties, including shape and effective density. The AHM provides an opportunity to investigate observed microphysical processes through the computation of polarimetric variables and corroboration with microphysical model output. This study is unique because the polarimetric scattering calculations are computed using predicted microphysical parameters rather than a priori assumptions that are imposed within the scattering calculations in the forward simulator, allowing for a more effective comparison to radar observations. Through the simulation of a case in the Front Range of the Rocky Mountains in Colorado using the Advanced Research version of the Weather Research and Forecasting Model, it is found that the AHM approximates ice mass, shape, cloud vertical structure, and temporal evolution as reflected through polarimetric quantities compared to observations. AHM reflectivity magnitudes are similar to those observed with radar and are an improvement over spherical ice crystal assumptions. Further analyses are completed to examine the effect of microphysical processes on the evolution of the differential reflectivity and specific differential phase, both of which are simulated using the AHM. Simulations reveal a polarimetric response to ice crystal mass, number, size, density, and aspect ratio. While results reveal the need for model improvements (e.g., parameterizations for aggregation rate), testing forward-simulated radar fields against observations is a first step in the validation of model microphysical and precipitation processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Anatomy of a Summertime Convective Event over the Arabian Region

This study investigates the structure and evolution of a summertime convective event that occurred on 14 July 2015 over the Arabian region. We use the WRF Model with 1-km horizontal grid spacing and test three PBL parameterizations: the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme; the Asymmetrical Convective Model, version 2, (ACM2) scheme; and the quasi-normal scale-elimination (QNSE) scheme. Convection initiates near the Al Hajar Mountains of northern Oman at around 1100 local time (LT; 0700 UTC) and propagates northwestward. A nonorographic convective band along the west coast of the United Arab Emirates (UAE) develops after 1500 LT as a result of the convergence of cold pools with the sea breeze from the Arabian Gulf. The model simulation employing the QNSE scheme simulates the convection initiation and propagation well. Although the MYNN and ACM2 simulations show convective initiation near the Al Hajar Mountains, they fail to simulate the development of the convective band along the UAE west coast. The MYNN run simulates colder near-surface temperatures and a weaker sea breeze, whereas the ACM2 run simulates a stronger sea breeze but a drier lower troposphere. Sensitivity simulations using horizontal grid spacings of 9 and 3 km show that lower-resolution runs develop broader convective structures and weaker cold pools and horizontal wind divergence, affecting the development of convection along the west coast of the UAE. The 1-km run using the QNSE PBL scheme realistically captures the sequence of events that leads to the moist convection over the UAE and adjacent mountains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Evaluating and Improving Cloud Processes in the Multi-Scale Modeling Framework

The research performed under this grant was intended to improve the embedded cloud model in the Multi-scale Modeling Framework (MMF) for convective clouds by using a 2-moment microphysics scheme rather than the single moment scheme used in all the MMF runs to date. The technical report and associated documents describe the results of testing the cloud resolving model with fixed boundary conditions and evaluation of model results with data. The overarching conclusion is that such model evaluations are problematic because errors in the forcing fields control the results so strongly that variations in parameterization values cannot be usefully constrained

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Toward the Development of a Cold Regions Regional-Scale Hydrologic Model, Final Project Report

This project improves meso-scale hydrologic modeling in the boreal forest by: (1) demonstrating the importance of capturing the heterogeneity of the landscape using small scale datasets for parameterization for both small and large basins; (2) demonstrating that in drier parts of the landscape and as the boreal forest dries with climate change, modeling approaches must consider the sensitivity of simulations to soil hydraulic parameters - such as residual water content - that are usually held constant. Thus, variability / flexibility in residual water content must be considered for accurate simulation of hydrologic processes in the boreal forest; (3) demonstrating that assessing climate change impacts on boreal forest hydrology through multiple model integration must account for direct effects of climate change (temperature and precipitation), and indirect effects from climate impacts on landscape characteristics (permafrost and vegetation distribution). Simulations demonstrated that climate change will increase runoff, but will increase ET to a greater extent and result in a drying of the landscape; and (4) vegetation plays a significant role in boreal hydrologic processes in permafrost free areas that have deciduous trees. This landscape type results in a decoupling of ET and precipitation, a tight coupling of ET and temperature, low runoff, and overall soil drying.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving High Resolution Offshore Wind Resource Assessments and Forecasts using Observations in the MA/RI Lease Areas

The third Wind Forecast Improvement Project (WFIP3) sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the Marine Atmospheric Boundary Layer (MABL). WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high use coastal zone, using a 3-D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the MABL over the ocean was done from an air-sea interaction flux tower and extended deployments of a large barge platform. WFIP3 focused on mesoscale and sub-mesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. Numerous critical forecasting phenomena were observed, however the project was terminated prior to the completion of the field observational period and the analysis period.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

ENSO Diversity Simulated in a Revised Cane-Zebiak Model

The El Niño-Southern Oscillation (ENSO) phenomenon features rich sea surface temperature (SST) spatial pattern variations dominated by the Central Pacific (CP) and Eastern Pacific (EP) patterns during its warm phase. Understanding such ENSO pattern diversity has been a subject under extensive research activity. To provide a framework for unveiling the fundamental dynamics of ENSO diversity, an intermediate coupled model based on the Cane-Zebiak-type framework, named RCZ, is established in this study. Compared with the original Cane-Zebiak model, RCZ consists of revised model formulation and well-tuned parameterization schemes. All model components are carefully validated against the observations via the standalone mode, in which the observed anomalous SST (wind stress) forcing is prescribed to drive the atmospheric (oceanic) component. The superiority of RCZ’s model components over those in the original Cane-Zebiak model is evidenced by their better performance in simulating the observations. Coupled simulation with RCZ satisfactorily reproduces aspects of the observed ENSO characteristics, including the spatial pattern, phase-locking, amplitude asymmetry, and, particularly, ENSO diversity/bi-modality. RCZ serves as a promising tool for studying dynamics of ENSO diversity as it resolves most of the relevant processes proposed in the literature, including atmospheric nonlinear convective heating, oceanic nonlinear dynamical heating, and the ENSO/westerly wind burst interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Process-based modeling of soil nitrous oxide emissions from United States corn fields under different management and climate scenarios coupled with evaluation using regional estimates

Direct emissions of soil nitrous oxide during a growing season (N 2 O gs ) can be quantified with process-based models considering interactions between management, climate, and soil moisture when key data are available. We used an adapted “parameterized CENTURY/DAYCENT-model” ( p CENTURY) calibrated with crop growth and soil organic matter decay coefficients at the county-level for the estimation of N 2 O gs in the United States Corn Belt. Model estimated N 2 O-emissions from corn-based biofuels scenarios considering crop rotation, fertilizer inputs, tillage, and weather were compared against meta-summary of field observations from 55 studies. Both model and meta-summary ranked N 2 O gs -emissions to be corn > wheat > soybean phase while model likely underestimated cover crop N 2 O gs -emissions. The N 2 O gs -emissions and the associated emission factors (EFs) were modeled and summarized to be greater after anhydrous ammonia than urea application and from conventional tilled than non-tilled fields. Modeled and observed N 2 O gs -emissions after organic and inorganic fertilizer amendment did not differ due to high variability associated with the treatments. However, the organic fertilizer associated EFs were greater according to meta-summary data because of N input rates. Regionalized weather scenarios indicate hotspots for N 2 O gs -emissions can occur where crop N uptake is limited during dry years and in eastern states also during normal or wet seasons. The p CENTURY-derived N 2 O gs EFs (0.91 ± 0.19%) for counties investigated were only slightly lower than literature (1.07 ± 0.57%) or Tier-1 (1%) values. Our preliminary evaluation of regional soil moisture estimates showed reasonable agreement between monthly soil moisture estimates and the North American Soil Moisture Dataset during the growing season, but overestimation of soil moisture in winter-spring can influence the estimates of annual N 2 O emissions so future work is needed to calibrate soil moisture-associated model parameters. Our work provided scenario-based estimates of climate and management impacts on soil N 2 O gs -emissions together with valuable spatial insights into EFs that will be improved by more accurate information of fertilizer inputs and more temporally refined model evaluation.

54 ENVIRONMENTAL SCIENCES↗