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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 901 records · Page 50

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Accuracy of Some Approximate Gaussian Filters for the Navier--Stokes Equation in the Presence of Model Error

Bayesian state estimation of a dynamical system from a stream of noisy measurements is important in many geophysical and engineering applications where high dimensionality of the state space, sparse observations, and model error pose key challenges. In this work, three computationally feasible, approximate Gaussian data assimilation/filtering algorithms are considered in various regimes of turbulent 2D Navier--Stokes dynamics in the presence of model error. The first source of error arises from the necessary use of reduced models for the forward dynamics of the filters, while a particular type of representation error arises from the finite resolution of observations which mix up information about resolved and unresolved dynamics. Two stochastically parameterized filtering algorithms, referred to as cSPEKF and GCF, are compared with 3DVAR---a prototypical time-sequential algorithm known to be accurate for filtering dissipative systems for a suitably inflated “background” covariance. We provide the first evidence that the stochastically parameterized algorithms, which do not rely on detailed knowledge of the underlying dynamics and do not require covariance inflation, can compete with or outperform an optimally tuned 3DVAR algorithm, and they can overcome competing sources of error in a range of dynamical scenarios.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

On the Need to Modify the Sea Surface Roughness Formulation over Shallow Waters

Abstract The wind stress formulation in an atmospheric model over shallow waters is investigated using year-long observations of the wind profile within the first 100 m of the atmosphere and mesoscale simulations. The model experiments use a range of planetary boundary layer parameterizations to quantify the uncertainty related to the turbulent closure assumptions and thus to isolate the dominant influence of the surface roughness formulation. Results indicate that a positive wind speed bias exists when common open-ocean formulations for roughness are adopted for a region with a water depth of 30 m. Imposition of a wind stress formulation that is consistent with previous shallow-water estimates is necessary to reconcile model wind speeds with observations, providing modeling evidence that supports the increase of surface drag over shallow waters. The possibility of including water depth in the parameterization of roughness length is examined.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification, and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of downgradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over 5 m. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and nighttime entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

Physics-Dynamics Coupling in Weather, Climate, and Earth System Models: Challenges and Recent Progress

Geophysical models of the atmosphere and ocean invariably involve parameterizations. These represent two distinct areas: a) Subgrid processes which the model cannot (yet) resolve, due to its discrete resolution, and b) sources in the equation, due to radiation for example. Hence coupling between these physics parameterizations and the resolved fluid dynamics and also between the dynamics of the different fluids in the system (air and water) is necessary. This coupling is an important aspect of geophysical models. However, often model development is strictly segregated into either physics or dynamics. Hence, this area has many more unanswered questions than in-depth understanding. Furthermore, recent developments in the design of dynamical cores (significant increase of resolution, move to non-hydrostatic equation sets etc), extended process physics (prognostic micro physics, 3D turbulence, non-vertical radiation etc) and predicted future changes of the computational infrastructure (Exascale with its need for task parallelism, data locality and asynchronous time stepping for example) is adding even more complexity and new questions.

Gross, Markus↗

DOE OSTI2020

Assessment of Precipitating Marine Stratocumulus Clouds in the E3SMv1 Atmosphere Model: A Case Study from the ARM MAGIC Field Campaign

This paper presents a process-oriented evaluation of precipitating stratocumulus and its transition to cumulus in version 1 of the Energy Exascale Earth System Model (E3SMv1) using comprehensive case-study observations from a field campaign of the Atmospheric Radiation Measurement program (ARM). The E3SMv1 single-column model (SCM) of the marine boundary layer and its low clouds and precipitation are compared to observations including subcloud drizzle retrievals from a combination of Doppler radar and lidar backscatter measurements. Additionally, the SCM is also compared to a large-eddy simulation (LES) of the same case. The combination of advanced remote sensing observations and LES is a powerful framework to evaluate the physical parameterizations of large-scale models. Given the observed large-scale environment, the E3SMv1 SCM realistically represents the evolution of clouds and boundary layer structure during the stratocumulus-to-cumulus transition. The model well simulates the liquid water path and its diurnal cycle in the stratocumulus period as well as the two-layer vertical thermodynamic structure and lower cloud fraction in the transition period. E3SMv1’s success in simulating the cloud in the stratocumulus period permitted examination of its precipitation processes. Here problems were identified with E3SMv1 producing an unrealistically small subcloud precipitation fraction, an unrealistic double peak in the vertical profiles of precipitation mass, and drizzle that evaporates too close to the surface. Further model diagnostics determined that these unrealistic characteristics resulted from an overly long microphysics time step and an unrealistic parameterization of the precipitation fraction. These results imply that careful consideration of these issues is needed in order to better simulate precipitation processes in marine stratocumulus.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Fast Computational Algorithms for Partial Differential Equations and Uncertainty Quantifications

This project concerned the construction, testing and analysis of computational algorithms for solving parameterized and stochastic partial differential equations. The study and understanding of equations of this type is of fundamental importance in numerous engineering and scientific applications. Examples include simulation of plasma dynamics in models of electric propulsion and nuclear fusion, simulation of multiphase flows, such as the flow of water, gas and oil in reservoirs, and structural analysis of the dependence of structures on materials. Parametrization is used in such settings when properties of the models such as viscosity of fluids or electric resistivity of materials are not precisely understood and instead are treated as random variables. The resulting solutions are themselves random, and having such solutions will enable engineers to use probabilistic methods to assess the likelihood of events, for example, whether a pollutant in a liquid will exceed a limit, and to use such analyses to develop ways to ensure positive outcomes. Construction of accurate (high resolution) computational solutions is expensive, requiring significant computer time and computational resources, and there is need to reduce computational cost to make simulation useful and effective. The aim of the project was to construct algorithms to efficiently compute surrogate solutions to parameterized problems to allow for efficient and accurate simulation. The technical approach used focused on two related strategies, based on rank-reduction methods and reduced-order models. These methods construct surrogate solutions of parameter-dependent models by projection or interpolation into low-dimensional approximation spaces. Cost savings are achieved if the low-dimensional spaces can be identified and constructed efficiently and if the resulting low-dimensional algebraic systems can be solved cheaply. Accomplishments include: Theoretical and empirical demonstration of the effectiveness of fast multigrid solution strategies for computing low-rank representations of parameter-dependent solutions to discrete partial differential equations, including the first proof establishing so-called textbook convergence properties for low-rank methods. Development of efficient solution algorithms for solving nonlinear parameter-dependent partial differential equations used in models of fluid dynamics. Developent of efficient algorithms for low-rank representation of solutions of time-dependent simulations of fluid dynamics using multi-dimensional tensor representations of solutions.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Understanding Processes Controlling the Temporal and Spatial Variations of PBL Structures Over the ARM SGP Site

The surface heat, moisture, and momentum fluxes are transferred to the atmosphere above through the planetary boundary layer (PBL), where vertical mixing due to turbulent eddies of different sizes plays critical roles. Therefore, reliably representing PBL processes in numerical models is critical for weather, climate, and air quality prediction. Currently, there are over ten PBL schemes that are selectable within the advanced research version of the Weather Research and Forecasting (WRF) model, indicative of the challenges in capturing the impacts of turbulence within the PBL in models. Further improvements in PBL parameterizations are needed for both weather and climate models, as emphasized in many recent national reports, but require an advanced understanding of the underlying boundary layer processes from observations. This project takes advantage of Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) investments in the atmospheric boundary layer observations and Large-Eddy Simulation (LES) ARM Symbiotic Simulation and Observation (LASSO) simulations to characterize PBL structures, understand key physical processes controlling the mixed layer development, and to evaluate PBL parameterization.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Demeter calibration and uncertainty analysis

Demeter is a community spatial downscaling model that disaggregates land use and land cover changes projected by integrated human-Earth system models. Demeter has not been intensively calibrated, and we still lack a good knowledge about its sensitivity to key parameters and the parameter uncertainties. We used long-term global satellite-based land cover records to calibrate key Demeter parameters. The results identified the optimal parameter values and showed that the parameterization substantially improved the model’s performance. The parameters of intensification ratio and selection threshold were the most sensitive and needed to be carefully tuned, especially for regional applications. Further, small parameter uncertainties after calibration can be inflated when propagated into future scenarios, suggesting that users should consider the parameterization equifinality to better account for the uncertainties in the Demeter downscaled products. Our study provides a key reference for Demeter users, and ultimately contribute to reducing the uncertainties in Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Turbulent kinetic energy over large offshore wind farms observed and simulated by the mesoscale model WRF (3.8.1)

Abstract. Wind farms affect local weather and microclimates;hence, parameterizations of their effects have been developed for numerical weather prediction models.While most wind farm parameterizations (WFPs) include drag effects of wind farms, models differ on whether or not an additionalturbulent kinetic energy (TKE) source should be included in theseparameterizations to simulate the impact of wind farms on the boundary layer.Therefore, we use aircraft measurements above large offshore wind farms in stable conditions to evaluateWFP choices. Of the three case studies we examine, we find the simulated ambient backgroundflow to agree with observations of temperature stratification and winds. This agreement allows us toexplore the sensitivity of simulated wind farm effects with respect to modeling choicessuch as whether or not to include a TKE source, horizontal resolution, vertical resolutionand advection of TKE. For a stably stratified marine atmospheric boundary layer (MABL), a TKE source and a horizontal resolutionon the order of 5 km or finer are necessary to represent the impact of offshore wind farms on the MABL.Additionally, TKE advection results in excessively reduced TKE over the wind farms, which in turn causesan underestimation of the wind speed deficit above the wind farm. Furthermore, using fine verticalresolution increases the agreement of the simulated wind speed with satellite observations of surface wind speed.

17 WIND ENERGY↗

DOE OSTI2014

Perturbed Physics Ensemble Simulations of Cirrus on the Cloud System-resolving Scale

In this study, the effect of uncertainties in the parameterization of ice microphysical processes and initial conditions on the variability of cirrus microphysical and radiative properties are investigated in a series of cloud system-resolving perturbed physics ensemble (PPE) and initial condition ensemble (ICE) simulations. Three cirrus cases representative of mid-latitude, subtropical and tropical cirrus are examined. It is found that the variability in cirrus properties induced by perturbing uncertain parameters in ice microphysics parameterizations outweighs the variability induced by perturbing the initial conditions in midlatitude and subtropical cirrus. However, in tropical anvil cirrus the variability in the PPE and ICE simulations is about the same order of magnitude. The cirrus properties showing the largest sensitivity are ice water content (IWC) and cloud thickness whereas the averaged high cloud cover is only marginally affected. Changes in cirrus ice water path and outgoing longwave radiation are controlled primarily by changes in IWC and cloud thickness but not by changes is the averaged high cloud cover. The change in the vertical distribution of cloud fraction and cloud thickness is caused by changes in cirrus cloud base whereas cloud top is not sensitive to either perturbed physics or perturbed initial conditions. In all cirrus cases, the top three parameters controlling the microphysical variability and radiative impact of cirrus clouds are ice fall speeds, ice autoconversion size thresholds and heterogeneous ice nucleation. Changes in the ice deposition coefficient do not affect the ice water path and outgoing longwave radiation. Similarly, changes in the number concentration of aerosols available for homogeneous freezing have virtually no effect on the microphysical and radiative properties of midlatitude and subtropical cirrus but only little impact on tropical anvil cirrus. Overall, the sensitivity of cirrus microphysical and radiative properties to uncertainties in ice microphysics is largest for midlatitude cirrus and smallest for tropical anvil cirrus.

cirrus↗

DOE OSTI2014

Intercomparison of the Cloud Water Phase among Global Climate Models

Mixed-phase clouds (clouds that consist of both cloud droplets and ice crystals) are frequently present in the Earth’s atmosphere and influence the Earth’s energy budget through their radiative properties, which are highly dependent on the cloud water phase. In this study, the phase partitioning of cloud water is compared among six global climate models (GCMs) and with Cloud and Aerosol Lidar with Orthogonal Polarization retrievals. It is found that the GCMs predict vastly different distributions of cloud phase for a given temperature, and none of them are capable of reproducing the spatial distribution or magnitude of the observed phase partitioning. While some GCMs produced liquid water paths comparable to satellite observations, they all failed to preserve sufficient liquid water at mixed-phase cloud temperatures. Our results suggest that validating GCMs using only the vertically integrated water contents could lead to amplified differences in cloud radiative feedback. The sensitivity of the simulated cloud phase in GCMs to the choice of heterogeneous ice nucleation parameterization is also investigated. The response to a change in ice nucleation is quite different for each GCM, and the implementation of the same ice nucleation parameterization in all models does not reduce the spread in simulated phase among GCMs. The results suggest that processes subsequent to ice nucleation are at least as important in determining phase and should be the focus of future studies aimed at understanding and reducing differences among the models.

cloud water phase↗

DOE OSTI2015

Impacts of a mushy‐layer thermodynamic approach in global sea‐ice simulations using the CICE sea‐ice model

Abstract We perform global simulations of the Los Alamos sea‐ice model, CICE, with a new thermodynamics component that has a fully prognostic, variable bulk salinity vertical profile based on mushy layer physics. The processes of gravity drainage, melt‐water flushing and snow‐ice formation are parameterized to allow the bulk salinity to evolve with time. We analyze the seasonal and spatial variation of sea‐ice bulk salinity, area, volume and thickness and compare these quantities to simulations using the previous thermodynamic component. Adjusting one of the gravity drainage parameters, we find good agreement between simulation results and fieldwork ice‐core observations of sea‐ice bulk salinity. As expected, bulk salinity is highest during periods of ice growth and lowest after periods of ice melt. In the northern hemisphere the new thermodynamics component produces thicker ice than the previous thermodynamics component. Of the nine major differences between the two models, differences in how salinities are calculated and how melt‐pond flushing is parameterized are the principal causes of this thickness difference. Thickness differences are smaller in the southern hemisphere than in the northern hemisphere since a greater fraction of ice melts, and differences cannot accumulate year‐on‐year. Model differences in how ice thickness changes and snow‐ice formation are calculated are the most important causes of the different thickness between the two thermodynamic components in the southern hemisphere. The melt‐pond area and volume are found to be highly sensitive to a parameter choice controlling drainage through macroscopic holes in the ice, in both hemispheres.

Turner, Adrian K.↗

DOE OSTI2014

Evaluation of convection‐permitting model simulations of cloud populations associated with the Madden‐Julian Oscillation using data collected during the AMIE/DYNAMO field campaign

Abstract Regional convection‐permitting model simulations of cloud populations observed during the 2011 Atmospheric Radiation Measurement (ARM) Madden‐Julian Oscillation Investigation Experiment/Dynamics of the Madden‐Julian Oscillation Experiment (AMIE/DYNAMO) field campaign are evaluated against ground‐based radar and ship‐based observations. Sensitivity of model simulated reflectivity, surface rain rate, and cold pool statistics to variations of raindrop breakup/self‐collection parameters in four state‐of‐the‐art two‐moment bulk microphysics schemes in the Weather Research and Forecasting (WRF) model is examined. The model simulations generally overestimate reflectivity from large and deep convective cells, and underestimate stratiform rain and the frequency of cold pools. In the sensitivity experiments, introduction of more aggressive raindrop breakup or decreasing the self‐collection efficiency increases the cold pool occurrence frequency in all of the simulations, and slightly reduces the reflectivity and precipitation statistics bias in some schemes but has little effect on the overall mean surface precipitation. Both the radar observations and model simulations of cloud populations show an approximate power law relationship between convective echo‐top height and equivalent convective cell radius.

Hagos, Samson↗

DOE OSTI2015

Evaluation of cloud-resolving model simulations of midlatitude cirrus with ARM and A-train observations

Cirrus clouds are ubiquitous in the upper troposphere and still constitute one of the largest uncertainties in climate predictions. Our paper evaluates cloud-resolving model (CRM) and cloud system-resolving model (CSRM) simulations of a midlatitude cirrus case with comprehensive observations collected under the auspices of the Atmospheric Radiation Measurements (ARM) program and with spaceborne observations from the National Aeronautics and Space Administration A-train satellites. The CRM simulations are driven with periodic boundary conditions and ARM forcing data, whereas the CSRM simulations are driven by the ERA-Interim product. Vertical profiles of temperature, relative humidity, and wind speeds are reasonably well simulated by the CSRM and CRM, but there are remaining biases in the temperature, wind speeds, and relative humidity, which can be mitigated through nudging the model simulations toward the observed radiosonde profiles. Simulated vertical velocities are underestimated in all simulations except in the CRM simulations with grid spacings of 500 m or finer, which suggests that turbulent vertical air motions in cirrus clouds need to be parameterized in general circulation models and in CSRM simulations with horizontal grid spacings on the order of 1 km. The simulated ice water content and ice number concentrations agree with the observations in the CSRM but are underestimated in the CRM simulations. The underestimation of ice number concentrations is consistent with the overestimation of radar reflectivity in the CRM simulations and suggests that the model produces too many large ice particles especially toward the cloud base. Simulated cloud profiles are rather insensitive to perturbations in the initial conditions or the dimensionality of the model domain, but the treatment of the forcing data has a considerable effect on the outcome of the model simulations. Despite considerable progress in observations and microphysical parameterizations, simulating the microphysical, macrophysical, and radiative properties of cirrus remains challenging. Comparing model simulations with observations from multiple instruments and observational platforms is important for revealing model deficiencies and for providing rigorous benchmarks. But, there still is considerable need for reducing observational uncertainties and providing better observations especially for relative humidity and for the size distribution and chemical composition of aerosols in the upper troposphere.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2014

Precipitation characteristics of CAM5 physics at mesoscale resolution during MC3E and the impact of convective timescale choice

Abstract The physics suite of the Community Atmosphere Model version 5 (CAM5) has recently been implemented in the Weather Research and Forecasting (WRF) model to explore the behavior of the parameterization suite at high resolution and within the controlled setting of a limited area model. The initial paper documenting this capability characterized the behavior for a northern high‐latitude region. This paper characterizes the precipitation characteristics for continental, midlatitude, springtime conditions during the Midlatitude Continental Convective Clouds Experiment (MC3E) over the central United States. This period exhibited a range of convective conditions from those driven strongly by large‐scale synoptic regimes to more locally driven convection. The study focuses on the precipitation behavior at 32 km grid spacing to better anticipate how the physics will behave in a global model when used at similar grid spacing in the coming years. Importantly, one change to the Zhang‐McFarlane deep convective parameterization when implemented in WRF was to make the convective timescale parameter an explicit function of grid spacing. This study examines the sensitivity of the precipitation to the default value of the convective timescale in WRF, which is 600 s for 32 km grid spacing, to the value of 3600 s used for 2° grid spacing in CAM5. For comparison, a 1200 s and an infinite convective timescale are also used. The results show that the 600 s timescale gives the most accurate precipitation amount over the central United States. However, this setting has the worst precipitation diurnal cycle, with the convection too tightly linked to the daytime surface heating. Longer timescales greatly improve the diurnal cycle but result in less precipitation and produce a low bias. An analysis of rain rates shows the accurate precipitation amount with the shorter timescale is assembled from an over abundance of drizzle combined with too few heavy rain events. With longer timescales, one can improve the frequency distribution, particularly for the extreme rain rates. Ultimately, without changing other aspects of the physics, one must decide between accurate diurnal timing and rain amount when choosing an appropriate convective timescale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Impact of the dynamical core on the direct simulation of tropical cyclones in a high-resolution global model

Our paper examines the impact of the dynamical core on the simulation of tropical cyclone (TC) frequency, distribution, and intensity. The dynamical core, the central fluid flow component of any general circulation model (GCM), is often overlooked in the analysis of a model's ability to simulate TCs compared to the impact of more commonly documented components (e.g., physical parameterizations). The Community Atmosphere Model version 5 is configured with multiple dynamics packages. This analysis demonstrates that the dynamical core has a significant impact on storm intensity and frequency, even in the presence of similar large-scale environments. In particular, the spectral element core produces stronger TCs and more hurricanes than the finite-volume core using very similar parameterization packages despite the latter having a slightly more favorable TC environment. Furthermore, these results suggest that more detailed investigations into the impact of the GCM dynamical core on TC climatology are needed to fully understand these uncertainties. Key Points The impact of the GCM dynamical core is often overlooked in TC assessments The CAM5 dynamical core has a significant impact on TC frequency and intensity A larger effort is needed to better understand this uncertainty

58 GEOSCIENCES↗

DOE OSTI2015

How does increasing horizontal resolution in a global climate model improve the simulation of aerosol-cloud interactions?

The Community Atmosphere Model Version 5 is run at horizontal grid spacing of 2, 1, 0.5, and 0.25 degrees, with the meteorology nudged towards the Year Of Tropical Convection analysis, and cloud simulators and the collocated A-Train satellite observations are used to explore the resolution dependence of aerosol-cloud interactions. The higher-resolution model produces results that agree better with observations, showing an increase of susceptibility of cloud droplet size, indicating a stronger first aerosol indirect forcing (AIF), and a decrease of susceptibility of precipitation probability, suggesting a weaker second AIF. The resolution sensitivities of AIF are attributed to those of droplet nucleation and precipitation parameterizations. The annual average AIF in the northern hemisphere mid-latitudes (where most anthropogenic emissions occur) in the 0.25° model is reduced by about 1 W m⁻² (-30%) compared to the 2° model, leading to a 0.26 W m⁻² reduction (-15%) in the global annual average AIF.

Ma, Po -Lun↗

DOE OSTI2015

Quantifying sources and sinks of reactive gases in the lower atmosphere using airborne flux observations

Atmospheric composition is governed by the interplay of emissions, chemistry, deposition, and transport. Substantial questions surround each of these processes, especially in forested environments with strong biogenic emissions. Utilizing aircraft observations acquired over a forest in the southeast U.S., we calculate eddy covariance fluxes for a suite of reactive gases and apply the synergistic information derived from this analysis to quantify emission and deposition fluxes, oxidant concentrations, aerosol uptake coefficients, and other key parameters. Evaluation of results against state-of-the-science models and parameterizations provides insight into our current understanding of this system and frames future observational priorities. As a near-direct measurement of fundamental process rates, airborne fluxes offer a new tool to improve biogenic and anthropogenic emissions inventories, photochemical mechanisms, and deposition parameterizations.

54 ENVIRONMENTAL SCIENCES↗