Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regional models”

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 145 records · Page 8

Evaluating a Regional Climate Model Simulation of Greenland Ice Sheet Snow and Firn Density for Improved Surface Mass Balance Estimates

Modeling vertical profiles of snow and firn density near the surface of the Greenland ice sheet (GrIS) is key to estimating GrIS mass balance, and by extension, global sea level change. To understand sources of error in simulated GrIS density, we compare GrIS density profiles from a leading regional climate model with coincident in situ measurements. We identify key contributors to model density and mass balance biases, including underestimated simulated fresh snow density (which leads to underestimation of density in the top 1 m of snow by ~10%). In areas undergoing frequent melting, positive density biases (of 7% in the top 1 m, and 10% between 1 and 10 m) are likely associated with errors in representing meltwater production, retention and refreezing. The results highlight the importance of accurately capturing fresh snow density and meltwater processes in models used to estimate GrIS mass balance change.

P M Alexander↗

Seasonal Monitoring of Melt and Accumulation Within the Deep Percolation Zone of the Greenland Ice Sheet and Comparison with Simulations of Regional Climate Modeling

Increasing melt over the Greenland Ice Sheet (GrIS) recorded over the past several years has resulted in significant changes of the percolation regime of the ice sheet. It remains unclear whether Greenland's percolation zone will act as a meltwater buffer in the near future through gradually filling all pore space or if near-surface refreezing causes the formation of impermeable layers, which provoke lateral runoff. Homogeneous ice layers within perennial firn, as well as near-surface ice layers of several meter thickness have been observed in firn cores. Because firn coring is a destructive method, deriving stratigraphic changes in firn and allocation of summer melt events is challenging. To overcome this deficit and provide continuous data for model evaluations on snow and firn density, temporal changes in liquid water content and depths of water infiltration, we installed an upward-looking radar system (upGPR) 3.4 m below the snow surface in May 2016 close to Camp Raven (66.4779 deg N, 46.2856 deg W) at 2120 m a.s.l. The radar is capable of quasi-continuously monitoring changes in snow and firn stratigraphy, which occur above the antennas. For summer 2016, we observed four major melt events, which routed liquid water into various depths beneath the surface. The last event in mid-August resulted in the deepest percolation down to about 2.3 m beneath the surface. Comparisons with simulations from the regional climate model MAR are in very good agreement in terms of seasonal changes in accumulation and timing of onset of melt. However, neither bulk density of near-surface layers nor the amounts of liquid water and percolation depths predicted by MAR correspond with upGPR data. Radar data and records of a nearby thermistor string, in contrast, matched very well for both timing and depth of temperature changes and observed water percolations. All four melt events transferred a cumulative mass of 56 kg/(sq. m) into firn beneath the summer surface of 2015. We find that continuous observations of liquid water content, percolation depths and rates for the seasonal mass fluxes are sufficiently accurate to provide valuable information for validation of model approaches and help to develop a better understanding of liquid water retention and percolation in perennial firn.

Heilig, Achim↗

Using a Genetic Algorithm to Model Broadband Regional Waveforms for Crustal Structure in the Western United States

In this study, we analyze regional seismograms to obtain the crustal structure in the eastern Great Basin and western Colorado plateau. Adopting a for- ward-modeling approach, we develop a genetic algorithm (GA) based parameter search technique to constrain the one-dimensional crustal structure in these regions. The data are broadband three-component seismograms recorded at the 1994-95 IRIS PASSCAL Colorado Plateau to Great Basin experiment (CPGB) stations and supplemented by data from U.S. National Seismic Network (USNSN) stations in Utah and Nevada. We use the southwestern Wyoming mine collapse event (M(sub b) = 5.2) that occurred on 3 February 1995 as the seismic source. We model the regional seismograms using a four-layer crustal model with constant layer parameters. Timing of teleseismic receiver functions at CPGB stations are added as an additional constraint in the modeling. GA allows us to efficiently search the model space. A carefully chosen fitness function and a windowing scheme are added to the algorithm to prevent search stagnation. The technique is tested with synthetic data, both with and without random Gaussian noise added to it. Several separate model searches are carried out to estimate the variability of the model parameters. The average Colorado plateau crustal structure is characterized by a 40-km-thick crust with velocity increases at depths of about 10 and 25 km and a fast lower crust while the Great Basin has approximately 35- km-thick crust and a 2.9-km-thick sedimentary layer.

Bhattacharyya, Joydeep↗

A high-resolution pseudo-polygon discrete element model for regional sea ice

Here, this work presents a pseudo-polygon discrete element model for high-resolution sea ice simulations. A scale-invariant bonded particle contact model is proposed to model joints between sea ice floes based on the smeared fracture model and a lattice spring beam model, and the Mohr–Coulomb failure criterion is implemented to represent the shearing failure mechanism of sea ice packings under complex loadings. All mechanical parameters of the bond model can be directly determined from laboratory tests. Validations of the proposed model are made by investigations of mechanical response and failure criteria of field sea ice sheets. Compared with the field observations of sea ice from satellite radar and in situ stress sensors, the proposed model is capable of reproducing the typical constitutive behavior and the Coulomb friction envelope of field sea ice. Finally, the proposed discrete element sea ice model is used to study the effect of loading rates on mechanical behavior including failure strength of regional sea ice.

54 ENVIRONMENTAL SCIENCES↗

The utilization of Seasat-A data in numerical weather prediction

The present paper provides the results from several numerical experiments which have been conducted in cases in which SASS (Seasat-A Satellite Scatterometer) winds were either included or excluded. Global scale experiments are discussed along with a regional area model. The regional model is essentially that described by Duffy (1981). A description of the initial conditions used in the regional model experiment is also provided. It is found that Seasat-A data can be important in defining the initial conditions which results in rapid maritime cyclogenesis. Currently conducted work is concerned with the development of a higher horizontal and vertical resolution global model, with an improved planetary boundary layer parameterization.

Duffy, D. G.↗

Runoff evaluation in an Earth System Land Model for permafrost regions in Alaska

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce. In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-based integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM's parameterized representation of total runoff. Results show that (1) ELM's total runoff was the same order of magnitude as the ATS simulations, and both models were similarly variable over time; (2) minor adjustments to coefficients in ELM's runoff parameterization improved the match between the ATS simulation and ELM's parameterized representation of annual and seasonal total runoff; (3) overall, runoff responses in ATS and ELM are more similar in flat hillslope environments compared to steep hillslopes; and (4) shallower active layer thicknesses and higher precipitation simulations resulted in lower correlations between the two models due to greater total runoff. By incorporating the optimized runoff coefficients from the Sagavanirktok River basin into ELM, the simulated total runoff better matched the streamflow observations at a small watershed located on the Seward Peninsula of Alaska. Our findings revealed important insights into the effectiveness of runoff parameterizations in land surface models and pathways for improving runoff coefficients in typical Arctic regions.

54 ENVIRONMENTAL SCIENCES↗

Simulation and Scaling of the Turbulent Vertical Heat Transport and Deep-Cycle Turbulence Across the Equatorial Pacific Cold Tongue

Microstructure observations in the Pacific cold tongue reveal that strong turbulence often penetrates into the thermocline producing hundreds of W/m2 of downward heat transport during nighttime and early morning. However, virtually all observations of this deep-cycle turbulence (DCT) are from 0◦N,140◦W. Here, a hierarchy of ocean process simulations including submesoscale-permitting regional models and turbulence-permitting large eddy simulations (LES) embedded in a regional model provide insight into DCT at and beyond 0◦N,140◦W. A regional hindcast quantifies the spatio-temporal variability of subsurface turbulent heat fluxes throughout the cold tongue from 1999-2016. Climatologically, the mean heat flux is strongest (∼100W/m2) within2◦ of the equator. The heat flux is slightly (∼10 W/m2) stronger in the northern than southern hemisphere throughout the cold tongue. The seasonal cycle of the heat flux, which is not everywhere aligned with the surface flux, ranges from 150 W/m2 near the equator to 30 W/m2 and 10 W/m2 at 4◦N and S respectively. Aseasonal variability of turbulent heat fluxes is logarithmically distributed and temporally uncorrelated with surface heat fluxes. The aseasonal variability of turbulence is highlighted in 34-day LES of Boreal autumn at 0◦N and 3◦N,140◦W. Intense DCT occurs frequently above the undercurrent at 0◦N and intermittently at 3◦N. Daily-mean heat fluxes scale with the bulk vertical shear and the wind stress, which together explain ∼90% of the daily variance across both LES. Observational validation of the scaling at 0◦N,140◦W is encouraging, but observations beyond 0◦N,140◦W can facilitate needed refinement of mixing parameterization in ocean models.

Simulation↗

A Spatiotemporal Prediction Model for Regional Scheduling of Shared Bicycles Based on the INLA Method

Dock-less bicycle-sharing programs have been widely accepted as an efficient mode to benefit health and reduce congestions. And modeling and prediction has always been a core proposition in the field of transportation. Most of the existing demand prediction models for shared bikes take regions as research objects; therefore, a POI-based method can be a beneficial complement to existing research, including zone-level, OD-level, and station-level techniques. Point of interest (POI) is the location description of spatial entities, which can reflect the cycling route characteristics for both commuting and noncommuting trips to a certain extent, and is also the main generating point and attraction point of shared-bike travel flow. In this study, we make an effort to model a POI-level cycling demand with a Bayesian hierarchical method. The proposed model combines the integrated nested Laplace approximation (INLA) and random partial differential equation (SPDE) to cope with the huge computation in the modeling process. In particular, we have adopted the dock-less bicycle-sharing rental records of Mobike as a case study to validate our method; the study area was one of the fastest growing urban districts in Shanghai in August 2016. The operation results show that the method can help better understand, measure, and characterize spatiotemporal patterns of bike-share ridership at the POI level and quantify the impact of the spatiotemporal effect on bicycle-sharing use.

Yu, Zhuoran↗

Tidal simulation using regional ocean modeling systems (ROMS)

The purpose of our research is to test the capability of ROMS in simulating tides. The research also serves as a necessary exercise to implement tides in an operational ocean forecasting system. In this paper, we emphasize the validation of the model tide simulation. The characteristics and energetics of tides of the region will be reported in separate publications.

tide modeling↗

Yield Estimates of three Historical Atmospheric Nuclear Explosions at Lop Nor, from Reduced Order Models of Regionally Recorded Rayleigh Waveforms and One Seismic Station: A Brief Communication

We use a reduced order model (ROM) for Rayleigh waveforms sourced by large atmospheric explosions along with data collected from a single seismic station (TLG) to estimate the yields of three atmospheric nuclear tests with order-megaton (MT) yields. These historical tests were conducted by China between 1973 and 1980 (CHIC 15, CHIC 16, and CHIC 26). We use our ROM to estimate yields of 2.34 MT (CHIC 15), 0.51 MT (CHIC 16), and 1.15 (CHIC 26) in the 20.5 s period band that are consistent with publicly accepted values of 3 MT (CHIC 15), 0.50 MT (CHIC 16), and 0.730 MT (CHIC 26).

58 GEOSCIENCES↗

Acoustic intensity calculations for axisymmetrically modeled fluid regions

An algorithm for calculating acoustic intensities from a time harmonic pressure field in an axisymmetric fluid region is presented. Acoustic pressures are computed in a mesh of NASTRAN triangular finite elements of revolution (TRIAAX) using an analogy between the scalar wave equation and elasticity equations. Acoustic intensities are then calculated from pressures and pressure derivatives taken over the mesh of TRIAAX elements. Intensities are displayed as vectors indicating the directions and magnitudes of energy flow at all mesh points in the acoustic field. A prolate spheroidal shell is modeled with axisymmetric shell elements (CONEAX) and submerged in a fluid region of TRIAAX elements. The model is analyzed to illustrate the acoustic intensity method and the usefulness of energy flow paths in the understanding of the response of fluid-structure interaction problems. The structural-acoustic analogy used is summarized for completeness. This study uncovered a NASTRAN limitation involving numerical precision issues in the CONEAX stiffness calculation causing large errors in the system matrices for nearly cylindrical cones.

Hambric, Stephen A.↗

Japanese MAGSAT team

Construction of a model of the regional magnetic field and investigation of the local magnetic anomalies and their origin were approaches used in attempts to study the crustal structure near Japan and its Antarctic bases. Spatial properties of the regional magnetic field and comparison of the regional model with that derived from MAGSAT data are discussed. Possible causes of the magnetic anomalies, and results of aeromagnetic surveys incorporating gravity and seismic data are explored. Ionospheric and magnetospheric contributions to geomagnetic variations, field-aligned currents, magnetic geomagnetic pulsations, and hydromagnetic waves by analysis of MAGSAT data are also examined.

Fukushima, N.↗

Regional Carbon Dioxide and Water Vapor Exchange Over Heterogeneous Terrain

In spite of setbacks due to forest fires, eviction after a change of landowners and unanticipated need to upgrade and replace much of the instrumentation, substantial progress has been made during the past three years, resulting in major new findings. Although most of the results are in manuscript form, three papers have been published and a fourth was recently submitted. The data has been subjected to extensive quality control. Extra attention has been devoted to the influence of tilt rotation and flux-calculation method, particularly with respect to nocturnal fluxes. Previous/standard methods for calculating nocturnal fluxes with moderate and strong stability are inadequate and lead to large random fluxes errors for individual records, due partly to inadvertent inclusion of mesoscale motions that strongly contaminant the estimation of fluxes by weak turbulence. Such large errors are serious for process studies requiring carbon dioxide fluxes for individual records, but are substantially reduced when averaging fluxes over longer periods as in calculation of annual NEE budgets. We have employed a superior method for estimating fluxes in stable conditions with a variable averaging width . Mesoscale fluxes are generally unimportant except for events and are generally not systematic or predictable. Mesoscale or regional models of our region are not able to reproduce important aspects of the diurnally varying wind field

Mahrt, Larry J.↗

Understanding the Effect of Land Cover Classification on Model Estimates of Regional Carbon Cycling in the Boreal Forest Biome

The original objectives of this proposed 3-year project were to: 1) quantify the respective contributions of land cover and disturbance (i.e., wild fire) to uncertainty associated with regional carbon source/sink estimates produced by a variety of boreal ecosystem models; 2) identify the model processes responsible for differences in simulated carbon source/sink patterns for the boreal forest; 3) validate model outputs using tower and field- based estimates of NEP and NPP; and 4) recommend/prioritize improvements to boreal ecosystem carbon models, which will better constrain regional source/sink estimates for atmospheric C02. These original objectives were subsequently distilled to fit within the constraints of a 1 -year study. This revised study involved a regional model intercomparison over the BOREAS study region involving Biome-BGC, and TEM (A.D. McGuire, UAF) ecosystem models. The major focus of these revised activities involved quantifying the sensitivity of regional model predictions associated with land cover classification uncertainties. We also evaluated the individual and combined effects of historical fire activity, historical atmospheric CO2 concentrations, and climate change on carbon and water flux simulations within the BOREAS study region.

Kimball, John↗

Assimilation of GEO and LEO Satellite Retrievals in WRF-Chem (15 km and 4 km) with ‘Top-Down’ Emissions Estimation

We are constraining concentrations and emissions for all criteria pollutants (CO, O3, NO2, SO2, PM10, and PM2.5) with WRF-Chem/DART in applications for FRAPPE (15-km grid) and COLORADO (4-km grid). Our results show that: (i) dynamic emissions estimation at medium resolutions (15 km) improves forecast skill; (ii) At high resolutions, the results look good, but we do not have validation data; (iii) for what may be the first time, we assimilate O3 retrieval profiles in a regional model. The problem is that the averaging kernels generally extend above the upper boundary of regional models. We use O3 upper boundary conditions from the global model to solve that problem; and (iv) In the high resolution experiments, we document the potential benefits of assimilating TEMPO O3 profile and NO2 tropospheric column retrievals. Here, we assimilate proxy TEMPO retrievals from a GEOS-Chem nature run, so there’s a conceptual problem due to the potential bias of the proxy retrievals, but the point is to demonstrate our ability to assimilate TEMPO and identify its potential impacts.

Chemical data assimilation↗

Soil moisture and hydrology projections of the permafrost region – a model intercomparison

Abstract. This study investigates and compares soil moisture and hydrology projections of broadly used land models with permafrost processes and highlights the causes and impacts of permafrost zone soil moisture projections. Climate models project warmer temperatures and increases in precipitation (P) which will intensify evapotranspiration (ET) and runoff in land models. However, this study shows that most models project a long-term drying of the surface soil (0–20 cm) for the permafrost region despite increases in the net air–surface water flux (P-ET). Drying is generally explained by infiltration of moisture to deeper soil layers as the active layer deepens or permafrost thaws completely. Although most models agree on drying, the projections vary strongly in magnitude and spatial pattern. Land models tend to agree with decadal runoff trends but underestimate runoff volume when compared to gauge data across the major Arctic river basins, potentially indicating model structural limitations. Coordinated efforts to address the ongoing challenges presented in this study will help reduce uncertainty in our capability to predict the future Arctic hydrological state and associated land–atmosphere biogeochemical processes across spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗

WMLES of Boundary-Layer Transition

The potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) is investigated through analysis of Direct Numerical Simulation (DNS) and WMLES results. The analysis indicates satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with a WMLES computation that yields reasonable agreement with DNS results; however, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation, boundary-layer↗