Programming active cohesive granular matter with mechanically induced phase changes
Interactions of simple robots can be mapped to a discrete algorithm that accurately predicts the ensemble’s emergent behaviors.
Li, Shengkai↗
SEARCH · Engineering Papers
Results for “ensemble prediction”
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.
Exascale Dislocation Simulator
ExaDiS, or Exascale Dislocation Simulator, is a set of software modules written to enable numerical simulations of large groups of moving and interacting dislocations, line defects in crystals responsible for crystal plasticity. By tracking time evolution of sufficiently large dislocation ensembles, ExaDiS predicts plasticity response and plastic strength of crystalline materials. ExaDiS is built around a library of core functions for Discrete Dislocation Dynamics method specifically written to perform efficient computations on GPUs.
Bertin, Nicolas↗
Improving Solution Accuracy and Convergence for Stochastic Physics Parameterizations with Colored Noise
Stochastic parameterizations are used in numerical weather prediction and climate modeling to help capture the uncertainty in the simulations and improve their statistical properties. Convergence issues can arise when time integration methods originally developed for deterministic differential equations are applied naively to stochastic problems. In previous studies, it has been demonstrated that a correction term, known in stochastic analysis as the Itô correction, can help improve solution accuracy for various deterministic numerical schemes and ensure convergence to the physically relevant solution without substantial computational overhead. The usual formulation of the Itô correction is valid only when the stochasticity is represented by white noise. In this study, a generalized formulation of the Itô correction is derived for noises of any color. The formulation is applied to a test problem described by an advection–diffusion equation forced with a spectrum of fast processes. We present numerical results for cases with both constant and spatially varying advection velocities to show that, for the same time step sizes, the introduction of the generalized Itô correction helps to substantially reduce time integration error and significantly improve the convergence rate of the numerical solutions when the forcing term in the governing equation is rough (fast varying); alternatively, for the same target accuracy, the generalized Itô correction allows for the use of significantly longer time steps and, hence, helps to reduce the computational cost of the numerical simulation.
54 ENVIRONMENTAL SCIENCES↗
Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models
Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.
Geochemistry & Geophysics↗
Impacts of Mesh Refinement and Number of Ensemble Members on Climate Variability Predictions in E3SM
Explore the source record for details and available documents.
Houchens, Brent C.↗
Statistical mechanics of 'negative temperature' states
Consideration of the dynamics of a two-dimensional guiding center plasma, recently shown by Taylor and McNamara (1971) to be identical to the dynamics of the discrete vortex model of Onsager (1949). A semirigorous application of the methods of equilibrium statistical mechanics to the guiding center plasma (or equivalently, the line vortex system) is presented. An adaptation of the apparatus of the theory of probability is attempted, in the form given by Khinchin (1949) to obtain ensemble-average predictions for the states of the guiding center plasma. Interest focuses primarily on the regime in which the interaction energy is high enough to be above the Onsager 'negative temperature' threshold.
Montgomery, D.↗
Soil Moisture and Snow Cover: Active or Passive Elements of Climate?
A key question in the study of the hydrologic cycle is the extent to which surface effects such as soil moisture and snow cover are simply passive elements or whether they can affect the evolution of climate on seasonal and longer time scales. We have constructed ensembles of predictability studies using the NCAR CCM3 in which we compared the relative roles of initial surface and atmospheric conditions over the central and western U.S. GAPP region in determining the subsequent evolution of soil moisture and of snow cover. We have also made sensitivity studies with exaggerated soil moisture and snow cover anomalies in order to determine the physical processes that may be important. Results from simulations with realistic soil moisture anomalies indicate that internal climate variability may be the strongest factor, with some indication that the initial atmospheric state is also important. The initial state of soil moisture does not appear important, a result that held whether simulations were started in late winter or late spring. Model runs with exaggerated soil moisture reductions (near-desert conditions) showed a much larger effect, with warmer surface temperatures, reduced precipitation, and lower surface pressures; the latter indicating a response of the atmospheric circulation. These results suggest the possibility of a threshold effect in soil moisture, whereby an anomaly must be of a sufficient size before it can have a significant impact on the atmospheric circulation and hence climate. Results from simulations with realistic snow cover anomalies indicate that the time of year can be crucial. When introduced in late winter, these anomalies strongly affected the subsequent evolution of snow cover. When introduced in early winter, however, little or no effect is seen on the subsequent snow cover. Runs with greatly exaggerated initial snow cover indicate that the high reflectivity of snow is the most important process by which snow cover can impact climate, through lower surface temperatures and increased surface pressures. In early winter, the amount of solar radiation is very small and so this albedo, effect is inconsequential while in late winter, with the sun higher in the sky and period of daylight longer, the effect is much stronger. The results to date were obtained for model runs with present-day conditions. We are currently analyzing runs made with projected forcings for the 21st century to see if these results are modified in any way under likely scenarios of future climate change.
Oglesby, Robert J.↗
Soil Moisture and Snow Cover: Active or Passive Elements of Climate?
A key question in the study of the hydrologic cycle is the extent to which surface effects such as soil moisture and snow cover are simply passive elements or whether they can affect the evolution of climate on seasonal and longer time scales. We have constructed ensembles of predictability studies using the NCAR CCM3 in which we compared the relative roles of initial surface and atmospheric conditions over the central and western U.S. in determining the subsequent evolution of soil moisture and of snow cover. We have also made sensitivity studies with exaggerated soil moisture and snow cover anomalies in order to determine the physical processes that may be important. Results from simulations with realistic soil moisture anomalies indicate that internal climate variability may be the strongest factor, with some indication that the initial atmospheric state is also important. The initial state of soil moisture does not appear important, a result that held whether simulations were started in late winter or late spring. Model runs with exaggerated soil moisture reductions (near-desert conditions) showed a much larger effect, with warmer surface temperatures, reduced precipitation, and lower surface pressures; the latter indicating a response of the atmospheric circulation. These results suggest the possibility of a threshold effect in soil moisture, whereby an anomaly must be of a sufficient size before it can have a significant impact on the atmospheric circulation and hence climate. Results from simulations with realistic snow cover anomalies indicate that the time of year can be crucial. When introduced in late winter, these anomalies strongly affected the subsequent evolution of snow cover. When introduced in early winter, however, little or no effect is seen on the subsequent snow cover. Runs with greatly exaggerated initial snow cover indicate that the high reflectively of snow is the most important process by which snow cover cart impact climate, through lower surface temperatures and increased surface pressures. In early winter, the amount of solar radiation is very small and so this albedo effect is inconsequential while in late winter, with the sun higher in the sky and period of daylight longer, the effect is much stronger.
Oglesby, Robert J.↗
Soil Moisture and Snow Cover: Active or Passive Elements of Climate
A key question is the extent to which surface effects such as soil moisture and snow cover are simply passive elements or whether they can affect the evolution of climate on seasonal and longer time scales. We have constructed ensembles of predictability studies using the NCAR CCM3 in which we compared the relative roles of initial surface and atmospheric conditions over the central and western U.S. in determining the subsequent evolution of soil moisture and of snow cover. Results from simulations with realistic soil moisture anomalies indicate that internal climate variability may be the strongest factor, with some indication that the initial atmospheric state is also important. Model runs with exaggerated soil moisture reductions (near-desert conditions) showed a much larger effect, with warmer surface temperatures, reduced precipitation, and lower surface pressures; the latter indicating a response of the atmospheric circulation. These results suggest the possibility of a threshold effect in soil moisture, whereby an anomaly must be of a sufficient size before it can have a significant impact on the atmospheric circulation and climate. Results from simulations with realistic snow cover anomalies indicate that the time of year can be crucial. When introduced in late winter, these anomalies strongly affected the subsequent evolution of snow cover. When introduced in early winter, however, little or no effect is seen on the subsequent snow cover. Runs with greatly exaggerated initial snow cover indicate that the high reflectivity of snow is the most important process by which snow cover can impact climate, through lower surface temperatures and increased surface pressures. The results to date were obtained for model runs with present-day conditions. We are currently analyzing runs made with projected forcings for the 21st century to see if these results are modified in any way under likely scenarios of future climate change. An intriguing new statistical technique involving 'clustering' is developed to assist in this analysis.
Oglesby, Robert J.↗
Estimating Single-Trial Responses in EEG
Accurate characterization of single-trial field potential responses is critical from a number of perspectives. For example, it allows differentiation of an evoked response from ongoing EEG. We previously developed the multiple component Event Related Potential (mcERP) algorithm to improve resolution of the single-trial evoked response. The mcERP model states that multiple components, each specified by a stereotypic waveform varying in latency and amplitude from trial to trial, comprise the evoked response. Application of the mcERP algorithm to simulated data with three independent, synthetic components has shown that the model is capable of separating these components and estimating their variability. Application of the model to single trial, visual evoked potentials recorded simultaneously from all V1 laminae in an awake, fixating macaque yielded local and far-field components. Certain local components estimated by the model were distributed in both granular and supragranular laminae. This suggests a linear coupling between the responses of thalamo-recipient neuronal ensembles and subsequent responses of supragranular neuronal ensembles, as predicted by the feedforward anatomy of V1. Our results indicate that the mcERP algorithm provides a valid estimation of single-trial responses. This will enable analyses that depend on trial-to-trial variations and those that require separation of the evoked response from background EEG rhythms
Shah, A. S.↗
ENSO Bred Vectors in Coupled Ocean-Atmosphere General Circulation Models
The breeding method has been implemented in the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Coupled General Circulation Model (CGCM) with the goal of improving operational seasonal to interannual climate predictions through ensemble forecasting and data assimilation. The coupled instability as cap'tured by the breeding method is the first attempt to isolate the evolving ENSO instability and its corresponding global atmospheric response in a fully coupled ocean-atmosphere GCM. Our results show that the growth rate of the coupled bred vectors (BV) peaks at about 3 months before a background ENSO event. The dominant growing BV modes are reminiscent of the background ENSO anomalies and show a strong tropical response with wind/SST/thermocline interrelated in a manner similar to the background ENSO mode. They exhibit larger amplitudes in the eastern tropical Pacific, reflecting the natural dynamical sensitivity associated with the presence of the shallow thermocline. Moreover, the extratropical perturbations associated with these coupled BV modes reveal the variations related to the atmospheric teleconnection patterns associated with background ENSO variability, e.g. over the North Pacific and North America. A similar experiment was carried out with the NCEP/CFS03 CGCM. Comparisons between bred vectors from the NSIPP CGCM and NCEP/CFS03 CGCM demonstrate the robustness of the results. Our results strongly suggest that the breeding method can serve as a natural filter to identify the slowly varying, coupled instabilities in a coupled GCM, which can be used to construct ensemble perturbations for ensemble forecasts and to estimate the coupled background error covariance for coupled data assimilation.
Yang, S. C.↗
Symmetry, Statistics and Structure in MHD Turbulence
Here, we examine homogeneous MHD turbulence in terms of truncated Fourier series. The ideal MHD equations and the associated statistical theory of absolute equilibrium ensembles are symmetric under P, C and T. However, the presence of invariant helicities, which are pseudoscalars under P and C, dynamically breaks this symmetry. This occurs because the surface of constant energy in phase space has disjoint parts, called components: while ensemble averages are taken over all components, a dynamical phase trajectory is confined to only one component. As the Birkhoff-Khinchin theorem tells us, ideal MHD turbulence is thus non-ergodic. This non-ergodicity manifests itself in low-wave number Fourier modes that have large mean values (while absolute ensemble theory predicts mean values of zero). Therefore, we have coherent structure in ideal MHD turbulence. The level of non-ergodicity and amount of energy contained in the associated coherent structure depends on the values of the helicities, as well as on the presence, or not, of a mean magnetic field and/or overall rotation. In addition to the well known cross and magnetic helicities, we also present a new invariant, which we call the parallel helicity, since it occurs when mean field and rotation axis are aligned. The question of applicability of these results to real (i.e., dissipative) MHD turbulence is also examined. Several long-time numerical simulations on a 64(exp 3) grid are given as examples. It is seen that coherent structure begins to form before decay dominates over nonlinearity. The connection of these results with inverse spectral cascades, selective decay, and magnetic dynamos is also discussed.
Shebalin, John V.↗
A Method for Snow Reanalysis: The Sierra Nevada (USA) Example
This work presents a state-of-the art methodology for constructing snow water equivalent (SWE) reanalysis. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data from the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlation greater than 0.95 compared with in situ SWE observations. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (31 years), and accuracy provide a unique dataset for investigating snow processes. This presentation illustrates how the reanalysis dataset was used to provide a basic accounting of the stored snowpack water in the Sierra Nevada over the last 31 years and ultimately improve real-time streamflow predictions.
Bayesian↗
Soil Organic Carbon and Nitrogen Feedbacks on Crop Yields Under Climate Change
A critical omission from climate change impact studies on crop yield is the interaction between soil organic carbon (SOC), nitrogen (N) availability, and carbon dioxide (CO2). We used a multimodel ensemble to predict the effects of SOC and N under different scenarios of temperatures and CO2 concentrations on maize (Zea mays L.) and wheat (Triticum aestivum L.) yield in eight sites across the world. We found that including feedbacks from SOC and N losses due to increased temperatures would reduce yields by 13% in wheat and 19% in maize for a 3°C rise temperature with no adaptation practices. These losses correspond to an additional 4.5% (+3°C) when compared to crop yield reductions attributed to temperature increase alone. Future CO2 increase to 540 ppm would partially compensate losses by 80% for both maize and wheat at +3°C, and by 35% for wheat and 20% for maize at +6°C, relative to the baseline CO2 scenario.
Basso, B.↗
Estimation of Seasonal Snow Water Equivalent Using Landsat Observations
This work presents a methodology for estimating seasonal snow water equivalent (SWE) from the use of remotely sensed Visible and Near Infrared observations from the Landsat mission. The method is comprised of two main components: (1) a coupled land surface model and snow depletion curve model, which is used to generate an ensemble of predictions of SWE and snow cover area for a given set of (uncertain) inputs, and (2) a reanalysis step, which updates estimation variables to be consistent with the satellite observed depletion of the fractional snow cover time series. This method was applied over the Sierra Nevada (USA) based on the assimilation of remotely sensed fractional snow covered area data over the Landsat 5-8 record (1985-2016). The verified dataset (based on a comparison with over 9000 station years of in situ data) exhibited mean and root-mean-square errors less than 3 and 13 cm, respectively, and correlations with in situ SWE observations of greater than 0.95. The method (fully Bayesian), resolution (daily, 90-meter), temporal extent (32 years), and accuracy provide a unique dataset for investigating snow processes. In particular, this presentation illustrates how the reanalysis dataset was used to provide climatology of the seasonal snowfall accumulation rates, distributions, and variability over the last three decades.
SWE↗
Deep Learning for Environmental and Ecological Prediction, Evaluation and Insight With Ensembles of Water Quality (Deep-View)
Explore the source record for details and available documents.
Stephanie Schollaert Uz↗
A Multiple-Model Convection-Permitting Ensemble Examination of the Probabilistic Prediction of Tropical Cyclones: Hurricanes Sandy (2012) and Edouard (2014)
Explore the source record for details and available documents.
Christopher Melhauser↗