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,513 records · Page 84

A laboratory study to estimate pore geometric parameters of sandstones using complex conductivity and nuclear magnetic resonance for permeability prediction

Here, we estimate parameters from the Katz and Thompson permeability model using laboratory complex electrical conductivity (CC) and nuclear magnetic resonance (NMR) data to build permeability models parameterized with geophysical measurements. We use the Katz and Thompson model based on the characteristic hydraulic length scale, determined from mercury injection capillary pressure estimates of pore throat size, and the intrinsic formation factor, determined from multisalinity conductivity measurements, for this purpose. Two new permeability models are tested, one based on CC data and another that incorporates CC and NMR data. From measurements made on forty-five sandstone cores collected from fifteen different formations, we evaluate how well the CC relaxation time and the NMR transverse relaxation times compare to the characteristic hydraulic length scale and how well the formation factor estimated from CC parameters compares to the intrinsic formation factor. We find: (1) the NMR transverse relaxation time models the characteristic hydraulic length scale more accurately than the CC relaxation time (R 2 of 0.69 and 0.33 and normalized root mean square errors (NRMSE) of 0.16 and 0.21, respectively); (2) the CC estimated formation factor is well correlated with the intrinsic formation factor (NRMSE50.23). We demonstrate that that permeability estimates from the joint-NMR-CC model (NRMSE50.13) compare favorably to estimates from the Katz and Thompson model (NRMSE50.074). Lastly, this model advances the capability of the Katz and Thompson model by employing parameters measureable in the field giving it the potential to more accurately estimate permeability using geophysical measurements than are currently possible.

36 MATERIALS SCIENCE↗

Contributions of wildland fire to terrestrial ecosystem carbon dynamics in North America from 1990 to 2012

Abstract Burn area and the frequency of extreme fire events have been increasing during recent decades in North America, and this trend is expected to continue over the 21st century. While many aspects of the North American carbon budget have been intensively studied, the net contribution of fire disturbance to the overall net carbon flux at the continental scale remains uncertain. Based on national scale, spatially explicit and long‐term fire data, along with the improved model parameterization in a process‐based ecosystem model, we simulated the impact of fire disturbance on both direct carbon emissions and net terrestrial ecosystem carbon balance in North America. Fire‐caused direct carbon emissions were 106.55 ± 15.98 Tg C/yr during 1990–2012; however, the net ecosystem carbon balance associated with fire was −26.09 ± 5.22 Tg C/yr, indicating that most of the emitted carbon was resequestered by the terrestrial ecosystem. Direct carbon emissions showed an increase in Alaska and Canada during 1990–2012 as compared to prior periods due to more extreme fire events, resulting in a large carbon source from these two regions. Among biomes, the largest carbon source was found to be from the boreal forest, primarily due to large reductions in soil organic matter during, and with slower recovery after, fire events. The interactions between fire and environmental factors reduced the fire‐caused ecosystem carbon source. Fire disturbance only caused a weak carbon source as compared to the best estimate terrestrial carbon sink in North America owing to the long‐term legacy effects of historical burn area coupled with fast ecosystem recovery during 1990–2012.

Chen, Guangsheng↗

Comparison of aerosol optical properties above clouds between POLDER and AeroCom models over the South East Atlantic Ocean during the fire season

Aerosol properties above clouds have been retrieved over the South East Atlantic Ocean during the fire season 2006 using satellite observations from POLDER (Polarization and Directionality of Earth Reflectances). From June to October, POLDER has observed a mean Above-Cloud Aerosol Optical Thickness (ACAOT) of 0.28 and a mean Above-Clouds Single Scattering Albedo (ACSSA) of 0.87 at 550nm. These results have been used to evaluate the simulation of aerosols above clouds in 5 AeroCom (Aerosol Comparisons between Observations and Models) models (GOCART, HadGEM3, ECHAM5-HAM2, OsloCTM2 and SPRINTARS). Most models do not reproduce the observed large aerosol load episodes. The comparison highlights the importance of the injection height and the vertical transport parameterizations to simulate the large ACAOT observed by POLDER. Moreover, some models overestimate the ACSSA. In accordance with recent recommendations of the black carbon refractive index, a higher prescription of the imaginary part allows a better comparison with POLDER’s ACSSA.

54 ENVIRONMENTAL SCIENCES↗

Thermally driven advection for radioxenon transport from an underground nuclear explosion

Abstract Barometric pumping is a ubiquitous process resulting in migration of gases in the subsurface that has been studied as the primary mechanism for noble gas transport from an underground nuclear explosion (UNE). However, at early times following a UNE, advection driven by explosion residual heat is relevant to noble gas transport. A rigorous measure is needed for demonstrating how, when, and where advection is important. In this paper three physical processes of uncertain magnitude (oscillatory advection, matrix diffusion, and thermally driven advection) are parameterized by using boundary conditions, system properties, and source term strength. Sobol' sensitivity analysis is conducted to evaluate the importance of all physical processes influencing the xenon signals. This study indicates that thermally driven advection plays a more important role in producing xenon signals than oscillatory advection and matrix diffusion at early times following a UNE, and xenon isotopic ratios are observed to have both time and spatial dependence.

Sun, Yunwei↗

Evaporation estimates using weather station data and boundary layer theory

Global estimates of evapotranspiration remain a challenge. In this study, we show that the daily course of air temperature and specific humidity available at routine weather stations can be used to estimate evapotranspiration and the evaporative fraction, the ratio of latent heat flux to available energy at the surface. Indeed, the diurnal increase in air temperature reflects the magnitude of the sensible heat flux and the increase of specific humidity after sunrise reflects the amplitude of evapotranspiration. Here, the method is physically constrained and based on the budget of heat and moisture in the boundary layer. Unlike land surface-based estimates, the proposed boundary layer estimate does not rely on ad hoc surface resistance parameterizations (e.g., Penman-Monteith). The proposed methodology can be applied to data collected at weather stations to estimate evapotranspiration and evaporative fraction under cloudy conditions and in the pre–remote sensing era.

54 ENVIRONMENTAL SCIENCES↗

Quantifying cloud base updraft speeds of marine stratocumulus from cloud top radiative cooling: QUANTIFYING CLOUD BASE UPDRAFTS OF MSC

Marine stratocumulus clouds play a significant role in the Earth's radiation budget. The updrafts at cloud base (Wb) govern the supersaturation and therefore the activation of cloud condensation nuclei, which modifies the cloud and precipitation properties. A statistically significant relationship between Wb and cloud top radiative cooling rate (CTRC) is found from the measurements of the Department of Energy's Atmospheric Radiation Measurement Mobile Facility on board a ship sailing between Honolulu and Los Angeles. A similar relation was found on Graciosa Island but with greater scatter and weaker correlation presumably due to the island effect. Based on the relation, we are able to estimate the cloud base updrafts using a simple formula: W b =-0.44 × CTRC + 22.30 ± 13, where the W b and CTRC have units of cm/s and W/m 2 , respectively. This quantification can be utilized in satellite remote sensing and parameterizations of W b in general circulation models.

54 ENVIRONMENTAL SCIENCES↗

Dependence of EMIC wave parameters during quiet, geomagnetic storm, and geomagnetic storm phase times

As electromagnetic ion cyclotron (EMIC) waves may play an important role in radiation belt dynamics, there has been a push to better include them into global simulations. How to best include EMIC wave effects is still an open question. Recently many studies have attempted to parameterize EMIC waves and their characteristics by geomagnetic indices. However, this does not fully take into account important physics related to the phase of a geomagnetic storm. In this paper we first consider how EMIC wave occurrence varies with the phase of a geomagnetic storm and the SYM-H, AE, and Kp indices. Here we show that the storm phase plays an important role in the occurrence probability of EMIC waves. The occurrence rates for a given value of a geomagnetic index change based on the geomagnetic condition. Then in this study we also describe the typical plasma and wave parameters observed in L and magnetic local time for quiet, storm, and storm phase. These results are given in a tabular format in the supporting information so that more accurate statistics of EMIC wave parameters can be incorporated into modeling efforts.

58 GEOSCIENCES↗

Simultaneous event-specific estimates of transport, loss, and source rates for relativistic outer radiation belt electrons: Event-Specific 1-D Modeling

The most significant unknown regarding relativistic electrons in Earth’s outer Van Allen radiation belt is the relative contribution of loss, transport, and acceleration processes within the inner magnetosphere. Detangling each individual process is critical to improve the understanding of radiation belt dynamics, but determining a single component is challenging due to sparse measurements in diverse spatial and temporal regimes. However, there are currently an unprecedented number of spacecraft taking measurements that sample different regions of the inner magnetosphere. With the increasing number of varied observational platforms, system dynamics can begin to be unraveled. In this work, we employ in-situ measurements during the 13-14 January 2013 enhancement event to isolate transport, loss, and source dynamics in a one dimensional radial diffusion model. We then validate the results by comparing them to Van Allen Probes and THEMIS observations, indicating that the three terms have been accurately and individually quantified for the event. Finally, a direct comparison is performed between the model containing event-specific terms and various models containing terms parameterized by geomagnetic index. Models using a simple 3/Kp loss timescale show deviation from the event specific model of nearly two orders of magnitude within 72 hours of the enhancement event. However, models using alternative loss timescales closely resemble the event specific model.

79 ASTRONOMY AND ASTROPHYSICS↗

Fault geometry inversion and slip distribution of the 2010 M w 7.2 El Mayor-Cucapah earthquake from geodetic data

The 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in Baja, California, and Sonora, Mexico, had primarily right-lateral strike-slip motion and a minor normal-slip component. The surface rupture extended about 120 km in a NW-SE direction, west of the Cerro Prieto fault. We use geodetic measurements including near- to far-field GPS, interferometric synthetic aperture radar (InSAR), and subpixel offset measurements of radar and optical images to characterize the fault slip during the EMC event. We use dislocation inversion methods and determine an optimal nine-segment fault geometry, as well as a subfault slip distribution from the geodetic measurements. With systematic perturbation of the fault dip angles, randomly removing one geodetic data constraint, or different data combinations, we are able to explore the robustness of the inferred slip distribution along fault strike and depth. The model fitting residuals imply contributions of early postseismic deformation to the InSAR measurements as well as lateral heterogeneity in the crustal elastic structure between the Peninsular Ranges and the Salton Trough. We also find that with incorporation of near-field geodetic data and finer fault patch size, the shallow slip deficit is reduced in the EMC event by reductions in the level of smoothing. Finally, these results show that the outcomes of coseismic inversions can vary greatly depending on model parameterization and methodology.

58 GEOSCIENCES↗

The microphysical contributions to and evolution of latent heating profiles in two MC3E MCSs

Abstract The shapes and magnitudes of latent heating profiles have been shown to be different within the convective and stratiform regions of mesoscale convective systems (MCSs). Properly representing these distinctions has significant implications for the atmospheric responses to latent heating on various scales. This study details (1) the microphysical process contributions to latent heating profiles within MCS convective, stratiform, and anvil regions and (2) the time evolution of these profiles throughout the MCS lifetime, using cloud‐resolving model simulations. Simulations of two MCS events that occurred during the Midlatitude Continental Convective Clouds Experiment (MC3E) are conducted. Several features of the simulated MCSs are compared to a suite of observations obtained during the MC3E field campaign, and it is concluded that the simulations reasonably reproduce the MCS events. The simulations show that condensation and deposition are the primary contributors to MCS latent warming, as compared to riming and nucleation processes. In terms of MCS latent cooling, sublimation, melting, and evaporation all play significant roles. It is evident that throughout the MCS lifecycle, convective regions demonstrate an approximately linear decrease in the magnitudes of latent heating rates, while latent heating within stratiform regions is associated with transitions between MCS flow regimes. Such information regarding the temporal evolution of latent heating within convective and stratiform MCS regions could be useful in developing parameterizations representing convective organization.

Marinescu, P. J.↗

Vertical overlap of probability density functions of cloud and precipitation hydrometeors: CLOUD AND PRECIPITATION PDF OVERLAP

Coarse-resolution climate models increasingly rely on probability density functions (PDFs) to represent subgrid-scale variability of prognostic variables. While PDFs characterize the horizontal variability, a separate treatment is needed to account for the vertical structure of clouds and precipitation. When sub-columns are drawn from these PDFs for microphysics or radiation parameterizations, appropriate vertical correlations must be enforced via PDF overlap specifications. This study evaluates the representation of PDF overlap in the Subgrid Importance Latin Hypercube Sampler (SILHS) employed in the assumed PDF turbulence and cloud scheme called the Cloud Layers Unified By Binormals (CLUBB). PDF overlap in CLUBB-SILHS simulations of continental and tropical oceanic deep convection is compared with overlap of PDF of various microphysics variables in cloud-resolving model (CRM) simulations of the same cases that explicitly predict the 3D structure of cloud and precipitation fields. CRM results show that PDF overlap varies significantly between different hydrometeor types, as well as between PDFs of mass and number mixing ratios for each species, - a distinction that the current SILHS implementation does not make. In CRM simulations that explicitly resolve cloud and precipitation structures, faster falling species, such as rain and graupel, exhibit significantly higher coherence in their vertical distributions than slow falling cloud liquid and ice. These results suggest that to improve the overlap treatment in the sub-column generator, the PDF correlations need to depend on hydrometeor properties, such as fall speeds, in addition to the currently implemented dependency on the turbulent convective length scale.

58 GEOSCIENCES↗

Investigating dominant characteristics of fires across the Amazon during 2005-2014 through satellite data synthesis of combustion signatures

Estimates of fire emissions remain uncertain due to limited constraints on the variations in fire characteristics. Here we demonstrate the utility of space-based observations of smoke constituents in addressing this limitation. We introduce a satellite-derived smoke index (SI) as an indicator of the dominant phase of large-scale fires. This index is calculated as the ratio of the geometric mean of observed fractional enhancements (due to fire) in carbon monoxide and aerosol optical depth to that of nitrogen dioxide. Here we assess the usefulness of this index on fires in the Amazon. We analyze the seasonal, regional, and interannual joint distribution of SI and fire radiative power (FRP) in relation to fire hotspots, land cover, Drought Severity Index, and deforestation rate estimates. We also compare this index with an analogous quantity derived from field data or emission inventories. Our results show that SI changes from low (more flaming) to high (more smoldering) during the course of a fire season, which is consistent with the changes in observed maximum FRPs from high to low. We also find that flaming combustion is more dominant in areas where deforestation fires dominate, while smoldering combustion has a larger influence during drought years when understory fires are more likely enhanced. Lastly, we find that the spatiotemporal variation in SI is inconsistent with current emission inventories. Although we recognize some limitations of this approach, our results point to the utility of SI as a proxy for overall combustion efficiency in the parameterization of fire emission models.

54 ENVIRONMENTAL SCIENCES↗

Land-atmosphere coupling and climate prediction over the U.S. Southern Great Plains

Biases in land-atmosphere coupling in climate models can contribute to climate prediction biases, but land models are rarely evaluated in the context of this coupling. Here, we tested land-atmosphere coupling and explored effects of land surface parameterizations on climate prediction in a single-column version of the National Center for Atmospheric Research Community Earth System Model (CESM1.2.2) and an off-line Community Land Model (CLM4.5). The correlation between leaf area index (LAI) and surface evaporative fraction (ratio of latent to total turbulent heat flux) was substantially underpredicted compared to observations in the U.S. Southern Great Plains, while the correlation between soil moisture and evaporative fraction was overpredicted by CLM4.5. To estimate the impacts of these errors on climate prediction, we modified CLM4.5 by prescribing observed LAI, increasing soil resistance to evaporation, increasing minimum stomatal conductance, and increasing leaf reflectance. The modifications improved the predicted soil moisture-evaporative fraction (EF) and LAI-EF correlations in off-line CLM4.5 and reduced the root-mean-square error in summer 2 m air temperature and precipitation in the coupled model. The modifications had the largest effect on prediction during a drought in summer 2006, when a warm bias in daytime 2 m air temperature was reduced from +6°C to a smaller cold bias of -1.3°C, and a corresponding dry bias in precipitation was reduced from -111 mm to -23 mm. The role of vegetation in droughts and heat waves is underpredicted in CESM1.2.2, and improvements in land surface models can improve prediction of climate extremes.

54 ENVIRONMENTAL SCIENCES↗

Exploring the effects of a nonhydrostatic dynamical core in high-resolution aquaplanet simulations

This study explores the impact of a nonhydrostatic dynamical core in high-resolution regional climate simulations using an aquaplanet framework. The Weather Research and Forecasting (WRF) model is used to conduct simulations with both hydrostatic (H) and nonhydrostatic (NH) solvers at horizontal grid spacings (?x) of 36, 12, and 4 km. The differences between the H and NH simulated precipitation (?P) are notable even at ?x = 12 km in the inter-tropical convergence zone and the transition region to the drier subtropics. At gray-zone grid spacing (12 km and 4 km) over the tropics, ?P is sensitive to whether a cumulus parameterization scheme is used or not. With an idealized Witch of Agnesi land mountain, differences in the precipitation and circulation (vertical velocity) between the H and NH simulations are significant even at ?x = 36 km in the tropics due largely to the strong feedbacks related to moist processes. The differences increase as the model grid spacing and mountain half width (a) are reduced, accompanied by a shift towards a more nonhydrostatic flow regime at a = 24 km. Latent heat release drastically enhances the differences between the NH and H simulations and extends the effect of nonhydrostatic dynamics to a broader region over the mountain and downstream over the tropics. Overall, differences exist between H and NH simulations even at resolutions between 12 and 36 km, but the differences are sensitive to the representations of moist physics and other features such as horizontal diffusion used in the WRF model.

nonhydrostatic dynamics, latent heating, nonhydros↗

The diurnal cycle of clouds and precipitation at the ARM SGP site: Cloud radar observations and simulations from the multiscale modeling framework

Millimeter Wavelength Cloud Radar (MMCR) data from December 1996 to December 2010, collected at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site, are used to examine the diurnal cycle of hydrometeor occurrence. These data are categorized into clouds (-40 dBZ e ≤ reflectivity < -10 dBZ e ), drizzle and light precipitation (-10 dBZ e ≤ reflectivity < 10 dBZ e ), and heavy precipitation (reflectivity ≥ 10 dBZ e ). The same criteria are implemented for the observation-equivalent reflectivity calculated by feeding outputs from a Multiscale Modeling Framework (MMF) climate model into a radar simulator. The MMF model consists of the National Center for Atmospheric Research Community Atmosphere Model with conventional cloud parameterizations replaced by a cloud-resolving model. We find that a radar simulator combined with the simple reflectivity categories can be an effective approach for evaluating diurnal variations in model hydrometeor occurrence. It is shown that the MMF only marginally captures observed increases in the occurrence of boundary layer clouds after sunrise in spring and autumn and does not capture diurnal changes in boundary layer clouds during the summer. Above the boundary layer, the MMF captures reasonably well diurnal variations in the vertical structure of clouds and light and heavy precipitation in the summer but not in the spring.

54 ENVIRONMENTAL SCIENCES↗

Can nudging be used to quantify model sensitivities in precipitation and cloud forcing?

Efficient simulation strategies are crucial for the development and evaluation of high resolution climate models. This paper evaluates simulations with constrained meteorology for the quantification of parametric sensitivities in the Community Atmosphere Model version 5 (CAM5). Two parameters are perturbed as illustrating examples: the convection relaxation time scale (TAU), and the threshold relative humidity for the formation of low-level stratiform clouds (rhminl). Results suggest that the fidelity and computational efficiency of the constrained simulations depend strongly on 3 factors: the detailed implementation of nudging, the mechanism through which the perturbed parameter affects precipitation and cloud, and the magnitude of the parameter perturbation. In the case of a strong perturbation in convection, temperature and/or wind nudging with a 6-hour relaxation time scale leads to non-negligible side effects due to the distorted interactions between resolved dynamics and parameterized convection, while a 1-year free running simulation can satisfactorily capture the annual mean precipitation sensitivity in terms of both global average and geographical distribution. In the case of a relatively weak perturbation the large-scale condensation scheme, results from 1-year free-running simulations are strongly affected by noise associated with internal variability, while nudging winds effectively reduces the noise, and reasonably reproduces the response of precipitation and cloud forcing to parameter perturbation. These results indicate that caution is needed when using nudged simulations to assess precipitation and cloud forcing sensitivities to parameter changes in general circulation models. We also demonstrate that ensembles of short simulations are useful for understanding the evolution of model sensitivities.

54 ENVIRONMENTAL SCIENCES↗

Coupling Spectral-bin Cloud Microphysics with the MOSAIC Aerosol Model in WRF-Chem: Methodology and Results for Marine Stratocumulus Clouds

Aerosol-cloud interaction processes can be represented more physically with bin cloud microphysics relative to bulk microphysical parameterizations. However, due to computational power limitations in the past, bin cloud microphysics was often run with very simple aerosol treatments. The purpose of this study is to represent better aerosol-cloud interaction processes in the Chemistry version of Weather Research and Forecast model (WRF-Chem) at convection-permitting scales by coupling spectral-bin cloud microphysics (SBM) with the MOSAIC sectional aerosol model. A flexible interface is built that exchanges cloud and aerosol information between them. The interface contains a new bin aerosol activation approach, which replaces the treatments in the original SBM. It also includes the modified aerosol resuspension and in-cloud wet removal processes with the droplet loss tendencies and precipitation fluxes from SBM. The newly coupled system is evaluated for two marine stratocumulus cases over the Southeast Pacific Ocean with either a simplified aerosol setup or full-chemistry. We compare the aerosol activation process in the newly-coupled SBM-MOSAIC against the SBM simulation without chemistry using a simplified aerosol setup, and the results show consistent activation rates. A longer time simulation reinforces that aerosol resuspension through cloud drop evaporation plays an important role in replenishing aerosols and impacts cloud and precipitation in marine stratocumulus clouds. Evaluation of the coupled SBM-MOSAIC with full-chemistry using aircraft measurements suggests that the new model works realistically for the marine stratocumulus clouds, and improves the simulation of cloud microphysical properties compared to a simulation using MOSAIC coupled with the Morrison two-moment microphysics.

Gao, Wenhua↗

Effects of explicit convection on global land-atmosphere coupling in the superparameterized CAM

Here, conventional global climate models are prone to producing unrealistic land-atmosphere coupling signals. Cumulus and convection parameterizations are natural culprits but the effect of bypassing them with explicitly resolved convection on global land-atmosphere coupling dynamics has not been explored systematically. We apply a suite of modern land-atmosphere coupling diagnostics to isolate the effect of cloud Superparameterization in the Community Atmosphere Model (SPCAM) v3.5, focusing on both the terrestrial segment (i.e., soil moisture and surface turbulent fluxes interaction) and atmospheric segment (i.e., surface turbulent fluxes and precipitation interaction) in the water pathway of the landatmosphere feedback loop. At daily timescales, SPCAM produces stronger uncoupled terrestrial signals (negative sign) over tropical rainforests in wet seasons, reduces the terrestrial coupling strength in the Central Great Plain in American, and reverses the coupling sign (from negative to positive) over India in the boreal summer season—all favorable improvements relative to reanalysis-forced land modeling. Analysis of the triggering feedback strength (TFS) and amplification feedback strength (AFS) shows that SPCAM favorably reproduces the observed geographic patterns of these indices over North America, with the probability of afternoon precipitation enhanced by high evaporative fraction along the eastern United States and Mexico, while conventional CAM does not capture this signal. We introduce a new diagnostic called the Planetary Boundary Layer (PBL) Feedback Strength (PFS), which reveals that SPCAM exhibits a tight connection between the responses of the lifting condensation level, the PBL height, and the rainfall triggering to surface turbulent fluxes; a triggering disconnect is found in CAM.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗