Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modelled”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Evaluating Integrated Surface/Subsurface Permafrost Thermal Hydrology Models in ATS (v0.88) Against Observations from a Polygonal Tundra Site: Modeling Archive

This Modeling Archive supports the NGEE Arctic publication "Evaluating integrated surface/subsurface permafrost thermal hydrology models in ATS (v0.88) against observations from a polygonal tundra site" available in the Geoscientific Model Development journal. We evaluated the integrated surface/subsurface permafrost thermal hydrology models in the Advanced Terrestrial Simulator (ATS) against field observations in polygonal tundra at the Barrow Environmental Observatory. We simulated thermal hydrology of three-dimensional ice-wedge polygons with generic but broadly representative surface microtopography. We drove the simulations with meteorological data and observed water table elevations in ice-wedge polygon troughs. With limited calibration of the soil properties and evaporation model parameters, the simulations were found to be consistent with observed water tables, active layer thickness, snowpack depth, evaporation, and observed soil temperatures at several depths in trough, rim, and center. The study builds confidence in the emerging class of integrated surface/subsurface permafrost simulators and provides an optimized set of model parameters for use in projections of permafrost dynamics in a warming climate. This dataset includes python, html, pdf, EXO, H5, xlsx, and csv files. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Model Deformation Measurements of Sonic Boom Models in the NASA Ames 9- by 7-Ft Supersonic Wind Tunnel

The deformations of two sonic-boom models were measured by stereo photogrammetry during tests in the 9- by 7-Ft Supersonic Wind Tunnel at NASA Ames Research Center. The models were geometrically similar but one was 2.75 times as large as the other. Deformation measurements were made by simultaneously imaging the upper surfaces of the models from two directions by calibrated cameras that were mounted behind windows of the test section. Bending and twist were measured at discrete points using conventional circular targets that had been marked along the leading and trailing edges of the wings and tails. In addition, continuous distributions of bending and twist were measured from ink speckles that had been applied to the upper surfaces of the model. Measurements were made at wind-on (M = 1.6) and wind-off conditions over a range of angles of attack between 2.5 deg. and 5.0 deg. At each condition, model deformation was determined by comparing the wind-off and wind-on coordinates of each measurement point after transforming the coordinates to reference coordinates tied to the model. The necessary transformations were determined by measuring the positions of a set of targets on the rigid center-body of the models whose model-axes coordinates were known. Smoothly varying bending and twist measurements were obtained at all conditions. Bending displacements increased in proportion to the square of the distance to the centerline. Maximum deflection of the wingtip of the larger model was about 5 mm (2% of the semispan) and that of the smaller model was 0.9 mm (1% of the semispan). The change in wing twist due to bending increased in direct proportion to distance from the centerline and reached a (absolute) maximum of about -1 at the highest angle of attack for both models. The measurements easily resolved bending displacements as small as 0.05 mm and bending-induced changes in twist as small as 0.05 deg.

model deformation↗
DOE OSTI2016

Hierarchical mixture of experts and diagnostic modeling approach to reduce hydrologic model structural uncertainty

In most water resources applications, a single model structure might be inadequate to capture the dynamic multi-scale interactions among different hydrological processes. Calibrating single models for dynamic catchments, where multiple dominant processes exist, can result in displacement of errors from structure to parameters, which in turn leads to over-correction and biased predictions. An alternative to a single model structure is to develop local expert structures that are effective in representing the dominant components of the hydrologic process and adaptively integrate them based on an indicator variable. In this study, the Hierarchical Mixture of Experts (HME) framework is applied to integrate expert model structures representing the different components of the hydrologic process. Various signature diagnostic analyses are used to assess the presence of multiple dominant processes and the adequacy of a single model, as well as to identify the structures of the expert models. The approaches are applied for two distinct catchments, the Guadalupe River (Texas) and the French Broad River (North Carolina) from the Model Parameter Estimation Experiment (MOPEX), using different structures of the HBV model. Furthermore, the results show that the HME approach has a better performance over the single model for the Guadalupe catchment, where multiple dominant processes are witnessed through diagnostic measures. Whereas, the diagnostics and aggregated performance measures prove that French Broad has a homogeneous catchment response, making the single model adequate to capture the response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Conditional Karhunen–Loève regression model with Basis Adaptation for high-dimensional problems: Uncertainty quantification and inverse modeling

Here, we propose a methodology for improving the accuracy of surrogate models of the observable response of physical systems as a function of the systems’ spatially heterogeneous parameter fields, with applications to uncertainty quantification and parameter estimation in high-dimensional problems. Practitioners often formulate finite-dimensional representations of spatially heterogeneous parameter fields using truncated unconditional Karhunen–Loève expansions (KLEs) for a certain choice of unconditional covariance kernel and construct surrogate models of the observable response with respect to the KLE coefficients. When direct measurements of the parameter fields are available, we propose improving the accuracy of these surrogate models by representing the parameter fields via conditional Karhunen-Loève expansions (CKLEs). CKLEs are constructed by conditioning the covariance kernel of the unconditional expansion on the direct measurements of the parameter field via Gaussian process regression, and then truncating the corresponding KLE. We apply the proposed methodology to constructing surrogate models via the Basis Adaptation (BA) method of the stationary hydraulic head response, measured at spatially discrete observation locations, of a groundwater flow model of the Hanford Site, as a function of the 1000-dimensional representation of the model’s log-transmissivity field. We find that BA surrogate models of the hydraulic head based on CKLEs are more accurate than BA surrogate models based on unconditional expansions for forward uncertainty quantification tasks. Furthermore, we find that inverse estimates of the hydraulic transmissivity field computed using CKLE-based BA surrogate models are more accurate than those computed using unconditional BA surrogate models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

A guideline for appropriate application of vertically-integrated modeling approaches for geologic carbon storage modeling

Mathematical modeling is an essential tool for answering questions related to geologic carbon storage (GCS). The choice of modeling approach depends on the type of questions being asked. In this paper we discuss a series of approaches with a hierarchical complexity including vertically-integrated single-phase flow approaches, vertically-integrated multi-phase flow approaches (with and without vertical equilibrium assumption), three-dimensional multi-phase flow approaches, and fully-coupled multi-phase flow approaches that couple flow with geochemistry and/or geomechanics. Three spatial scales are used to categorize the questions to be addressed by modeling: regional scale (encompasses CO 2 plume extent and majority of area of pressure impact of one or more injection operations), site scale (includes the CO 2 plume extent and some of the area impacted by the pressure increase of a single injection site), and well scale (the immediate vicinity of an injection well). A set of guidelines is developed to help modelers choose the most appropriate modeling approach, and show when simpler modeling approaches may be the better choice. Vertically-integrated single-phase flow models are the most appropriate choice at both the site and regional scales, if the pressure impact outside of the CO 2 plume is of interest. Vertically-integrated multi-phase flow models should be chosen at the regional scale, if the locations of CO 2 plumes are of interest, and at the site scale if vertical segregation of CO 2 and brine is fast or vertical heterogeneity in properties can be presented by distinct, continuous layers. Three-dimensional multi-phase flow models are the appropriate choice at the well and site scales for cases with significant vertical flow components of CO 2 and brine. Fully-coupled multi-phase flow models should only be chosen if pore-space alteration through geochemistry or geomechanics feeds back to fluid flow.

42 ENGINEERING↗
NASA NTRS2021

Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6

The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 °C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 °C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 °C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 °C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.d-2090s.

Scenario Model Intercomparison Project (ScenarioMI↗
DOE OSTI2017

Influence of landscape heterogeneity on water available to tropical forests in an Amazonian catchment and implications for modeling drought response: Water Available to Tropical Forest

The Amazon basin experienced periodic droughts in the past, and climate models projected more intense and frequent droughts in the future. How tropical forests respond to drought may depend on water availability, which is modulated by landscape heterogeneity. Using the one-dimensional ACME Land Model (ALM) and the three-dimensional ParFlow variably saturated flow model, a series of numerical experiments were performed for the Asu catchment in central Amazon to elucidate processes that influence water available for plant use and provide insights for improving Earth system models. Results from ParFlow show that topography has a dominant influence on groundwater table and runoff through lateral flow. Without any representations of lateral processes, ALM simulates very different seasonal variations in groundwater table and runoff compared to ParFlow even if it is able to reproduce the long-term spatial average groundwater table of ParFlow through simple parameter calibration. In the ParFlow simulations, the groundwater table is evidently deeper and the soil saturation is lower in the plateau compared to the valley. However, even in the plateau during the dry season in the drought year of 2005, plant transpiration is not water stressed in the ParFlow simulations as the soil saturation is still sufficient to maintain a soil matric potential for the stomata to be fully open. This finding is insensitive to uncertainty in atmospheric forcing and soil parameters, but the empirical wilting formulation used in the models is an important factor that should be addressed using observations and modeling of coupled plant hydraulics-soil hydrology processes in future studies.

Fang, Yilin↗
DOE OSTI2021

Assessment of a tomography-informed polyhedral discrete element modelling approach for complex-shaped granular woody biomass in stress consolidation

The design of handling equipment for granular biomass primarily requires experiments. Discrete element models (DEM) can provide designers with detailed insight into the behaviour of granular materials. However, granular biomass comprising complex-shaped particles is difficult to model with DEM. A tomography-informed DEM approach and an exhaustive assessment of the ability of the approach to predict the bulk behaviour of milled pines using experimental data is presented. Nano-CT scan was conducted to obtain 3D particle surface geometries as the basis for particle shape approximation by a polyhedral model and a sphero-polyhedral model. These models were applied in the simulation of a compressibility test. Our parametric study showed that particle Young's modulus and restitution coefficient are the two main properties influencing the simulated bulk behaviour of the DEM particles and that the level of approximating DEM particle shape for real particles is critical for accurately replicating the bulk behaviour of milled pines. The polyhedral model demonstrated better suitability than its sphero-polyhedral counterpart for modelling pine particles. The polyhedral model calibrated in a compressibility test was then applied in the simulations of a friction test without additional parameter tuning. Interlocking was dominant in the shear of bulk pine particles. Remarkably, the polyhedral model predicted the frictional behaviour of the pine particles when compared to the experimental results. The limitations of the model, as well as possible ways for enhancement, are discussed. This work provided novel insights into the suitability of complex-shaped DEM models for granular woody biomass.

Biomass flow↗
DOE OSTI2018

Role of Troposphere-Convection-Land Coupling in the Southwestern Amazon Precipitation Bias of the Community Earth System Model version 1 (CESM1)

The Amazon forest plays an important role in the tropical circulation and global water and biogeochemical cycles. However, the common Amazon precipitation biases in global climate models undermine our ability to understand future changes in the Amazon and their global consequences. The Community Earth System Model exhibits a persistent wet season dry bias over the southwestern Amazon. The bias is also characterized by inaccurate diurnal timing and intensity of rainfall events. Analysis of a diverse set of simulations suggests that the bias originates from local processes. A known land model bias in the partitioning of surface energy fluxes is found not to contribute significantly to the rainfall bias. Key model processes are identified to be the low-level divergent circulation controlling horizontal moisture flow and the sensitivity of the modeled deep convection to the lower tropospheric moisture. The lower tropospheric moisture divergence dictates the spatiotemporal distributions of convective precipitation through parameterized entrainment mixing, which generally improves the rainfall diurnal cycle. However, it also limits convective precipitation over the areas with frequent and large-scale divergence in the lower troposphere. The southwestern Amazon is such an area, particularly in the late afternoon, although uncertainty exists in the low-level flow in the simulation and reanalyses. One way to sustain rainfall in the presence of low-level divergence is a convective trigger based on both the moisture and temperature changes by large-scale processes, which is shown to help in reducing the dry bias.

Sakaguchi, Koichi↗
DOE OSTI2025

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine Learning Augmented Predictive and Generative Model for Rupture Life in Ferritic and Austenitic Steels

The Larson-Miller parameter (LMP) offers an efficient and fast scheme to estimate the creep rupture life of alloy materials for high temperature applications. However, owing to poor generalizability and dependence on the constant C, which is typically not known a-priori, estimations using the Larson-Miller parameter often result in suboptimal performance for a wide range of materials. At best it is useful in comparing alloys of similar composition. In this work, three machine learning (ML) schemes were developed for rupture life prediction for 9-12% Cr ferritic-martensitic steels and austenitic stainless steels, i.e., a hierarchical model to parameterize LMP using the LMP constant C to compute rupture life, a hierarchical model to parameterize both C and LMP to compute rupture life, and a direct prediction of rupture life. Specifically, we show that the third scheme, using a gradient boosting algorithm, can be used to train ML models for very accurate prediction of rupture life in a variety of alloys (Pear-son Correlation Coefficient > 0.9 for 9-12% Cr and > 0.8 for austenitic stainless steels). In addition, the Shapley value was used to quantify feature importance, making the model interpretable by identifying the effect of various features on the model performance. Furthermore, a variational autoencoder-based generative model was built by conditioning on the experimental dataset to sample hypothetical synthetic candidate alloys from the learnt joint distribution not existing in both 9-12 % Cr ferritic-martensitic steel and austenitic stainless steel datasets. Finally, based on the predictive and generative model, a reinforcement learning strategy has been proposed to guide experimentalists into designing better heat resistant alloys.

Mamun, Md Osman G.↗
DOE OSTI2014

Guidelines for the formulation of Lagrangian stochastic models for particle simulations of single-phase and dispersed two-phase turbulent flows

In this paper, we establish a set of criteria which are applied to discuss various formulations under which Lagrangian stochastic models can be found. These models are used for the simulation of fluid particles in single-phase turbulence as well as for the fluid seen by discrete particles in dispersed turbulent two-phase flows. The purpose of the present work is to provide guidelines, useful for experts and non-experts alike, which are shown to be helpful to clarify issues related to the form of Lagrangian stochastic models. A central issue is to put forward reliable requirements which must be met by Lagrangian stochastic models and a new element brought by the present analysis is to address the single- and two-phase flow situations from a unified point of view. For that purpose, we consider first the single-phase flow case and check whether models are fully consistent with the structure of the Reynolds-stress models. In the two-phase flow situation, coming up with clear-cut criteria is more difficult and the present choice is to require that the single-phase situation be well-retrieved in the fluid-limit case, elementary predictive abilities be respected and that some simple statistical features of homogeneous fluid turbulence be correctly reproduced. This analysis does not address the question of the relative predictive capacities of different models but concentrates on their formulation since advantages and disadvantages of different formulations are not always clear. Indeed, hidden in the changes from one structure to another are some possible pitfalls which can lead to flaws in the construction of practical models and to physically unsound numerical calculations. A first interest of the present approach is illustrated by considering some models proposed in the literature and by showing that these criteria help to assess whether these Lagrangian stochastic models can be regarded as acceptable descriptions. A second interest is to indicate how future developments can be safely built, which is also relevant for stochastic subgrid models for particle-laden flows in the context of Large Eddy Simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Surrogate modelling the Baryonic Universe II: On forward modelling the colours of individual and populations of galaxies

ABSTRACT Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model fsps, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model universemachine, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that universemachine predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with per cent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.

Chaves-Montero, Jonás (ORCID:0000000295534261)↗
DOE OSTI2025

Graph reinforcement learning for exploring model spaces beyond the standard model

We present a methodology for performing scans of beyond the standard model (BSM) parameter spaces with reinforcement learning. We identify a novel procedure using graph neural networks that is capable of exploring spaces of models without the user specifying a fixed particle content, allowing broad classes of BSM models to be explored—in theory, the technique is applicable to nearly any model space with a prespecified gauge group. We provide a generic procedure by which a suitable graph grammar can be developed for any BSM model that features user-specified symmetry groups and a finite number of different possible particle species, the use of which is applicable to a variety of machine learning tasks over the actions of BSM theories beyond our particular reinforcement learning use case. As a proof of concept, we construct the graph grammar for theories with vectorlike leptons that may or may not be charged under a dark U ( 1 ) group, inspired by portal matter extensions of the sub-GeV vector portal/kinetic mixing simplified dark matter models. We then use this graph grammar to create a reinforcement learning environment tasked with creating models with these vectorlike leptons that are consistent with a list of a variety of precision observables. The reinforcement learning agent succeeds in developing models that can address the observed muon anomalous magnetic moment discrepancy while remaining consistent with flavor violation and electroweak precision observables, including both constructions that have previously been studied as well as new models that have not, to our knowledge, previously been identified. By inspecting the resulting ensembles of models that the agent produces and experimenting with different configurations for our reinforcement learning environment and graph grammar, we also infer various lessons about the development of these environments that can be transferable to reinforcement learning scans of more complicated model spaces and comment on future directions for the development of this technique into a more mature tool. Published by the American Physical Society 2025

Wojcik, George N.↗
DOE OSTI2015

A revised partiality model and post-refinement algorithm for X-ray free-electron laser data

An updated partiality model and post-refinement algorithm for XFEL snapshot diffraction data is presented and confirmed by observing anomalous density for S atoms at an X-ray wavelength of 1.3 Å. Research towards using X-ray free-electron laser (XFEL) data to solve structures using experimental phasing methods such as sulfur single-wavelength anomalous dispersion (SAD) has been hampered by shortcomings in the diffraction models for X-ray diffraction from FELs. Owing to errors in the orientation matrix and overly simple partiality models, researchers have required large numbers of images to converge to reliable estimates for the structure-factor amplitudes, which may not be feasible for all biological systems. Here, data for cytoplasmic polyhedrosis virus type 17 (CPV17) collected at 1.3 Å wavelength at the Linac Coherent Light Source (LCLS) are revisited. A previously published definition of a partiality model for reflections illuminated by self-amplified spontaneous emission (SASE) pulses is built upon, which defines a fraction between 0 and 1 based on the intersection of a reflection with a spread of Ewald spheres modelled by a super-Gaussian wavelength distribution in the X-ray beam. A method of post-refinement to refine the parameters of this model is suggested. This has generated a merged data set with an overall discrepancy (by calculating the R{sub split} value) of 3.15% to 1.46 Å resolution from a 7225-image data set. The atomic numbers of C, N and O atoms in the structure are distinguishable in the electron-density map. There are 13 S atoms within the 237 residues of CPV17, excluding the initial disordered methionine. These only possess 0.42 anomalous scattering electrons each at 1.3 Å wavelength, but the 12 that have single predominant positions are easily detectable in the anomalous difference Fourier map. It is hoped that these improvements will lead towards XFEL experimental phase determination and structure determination by sulfur SAD and will generally increase the utility of the method for difficult cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Greenland Surface Mass Balance as Simulated by the Community Earth System Model. Part I: Model Evaluation and 1850–2005 Results

The modeling of the surface mass balance (SMB) of the Greenland Ice Sheet (GIS) requires high-resolution models in order to capture the observed large gradients in the steep marginal areas. Until now, global climate models have not been considered suitable to model ice sheet SMB owing to model biases and insufficient resolution. This study analyzes the GIS SMB simulated for the period 1850–2005 by the Community Earth System Model (CESM), which includes a new ice sheet component with multiple elevation classes for SMB calculations. The model is evaluated against observational data and output from the regional model Regional Atmospheric Climate Model version 2 (RACMO2). Because of a lack of major climate biases, a sophisticated calculation of snow processes (including surface albedo evolution) and an adequate downscaling technique, CESM is able to realistically simulate GIS surface climate and SMB. CESM SMB agrees reasonably well with in situ data from 475 locations ($r$ = 0.80) and output from RACMO2 ($r$ = 0.79). The simulated mean SMB for 1960–2005 is 359 ± 120 Gt yr -1 in the range of estimates from regional climate models. The simulated seasonal mass variability is comparable with mass observations from the Gravity Recovery and Climate Experiment (GRACE), with synchronous annual maximum (May) and minimum (August–September) and similar amplitudes of the seasonal cycle. CESM is able to simulate the bands of precipitation maxima along the southeast and northwest margins, but absolute precipitation rates are underestimated along the southeastern margin and overestimated in the high interior. The model correctly simulates the major ablation areas. Total refreezing represents 35% of the available liquid water (the sum of rain and melt).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

A Physically Based Runoff Routing Model for Land Surface and Earth System Models

A new physically based runoff routing model, called the Model for Scale Adaptive River Transport (MOSART), has been developed to be applicable across local, regional, and global scales. Within each spatial unit, surface runoff is first routed across hillslopes and then discharged along with subsurface runoff into a ‘‘tributary subnetwork’’ before entering the main channel. The spatial units are thus linked via routing through the main channel network, which is constructed in a scale-consistent way across different spatial resolutions. All model parameters are physically based, and only a small subset requires calibration.MOSART has been applied to the Columbia River basin at 1/ 168, 1/ 88, 1/ 48, and 1/ 28 spatial resolutions and was evaluated using naturalized or observed streamflow at a number of gauge stations. MOSART is compared to two other routing models widely used with land surface models, the River Transport Model (RTM) in the Community Land Model (CLM) and the Lohmann routing model, included as a postprocessor in the Variable Infiltration Capacity (VIC) model package, yielding consistent performance at multiple resolutions. MOSART is further evaluated using the channel velocities derived from field measurements or a hydraulic model at various locations and is shown to be capable of producing the seasonal variation and magnitude of channel velocities reasonably well at different resolutions. Moreover, the impacts of spatial resolution on model simulations are systematically examined at local and regional scales. Finally, the limitations ofMOSART and future directions for improvements are discussed.

runoff routing, scale adaptive↗
DOE OSTI2020

Two-Layer Dynamic Recycling Model (2L-DRM): Learning from Moisture Tracking Models of Different Complexity

Atmospheric moisture tracking models are used to identify and quantify sources and sinks of water in the atmospheric branch of the hydrologic cycle. These models are primarily used to investigate the origin of moisture resulting in precipitation for particular regions around the globe. Moisture tracking models vary widely in their level of complexity, depending on the number of physical processes represented. Complex models are comprehensive in their physical representation, but computationally much more expensive than simple models, which only focus on specific physical processes and use simplifying assumptions. We present the mathematical derivation of the new two-layer dynamical recycling model (2L-DRM), a simple analytical moisture tracking model that relaxes the vertically-integrated formulation of the original one-layer DRM (Dominguez et al., 2006). By comparing the simple DRM to a very complex moisture tracking model that uses water vapor tracers embedded within the Weather Research and Forecast regional climate model (WRF-WVT) for the North American Monsoon region, we pinpoint the absence of vertical wind shear as the main deficiency in the simple DRM. When comparing both simple models (DRM and 2L-DRM) to the WRF-WVT (which we treat as ‘truth’), the 2L-DRM better captures the spatial extent, the net amount and the temporal variability of precipitation that originates from oceanic and local terrestrial sources. The 2L-DRM is well suited to study the large-scale climatological sources of moisture, and for these applications, performs on-par with the much more complex and computationally demanding WRF-WVT model.

Dominguez, Francina↗