Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Long-Term Behavior of the Atlantic Interhemispheric SST Gradient in the CMIP5 Historical Simulations

Multidecadal and longer changes to the Atlantic interhemispheric sea surface temperature gradient (AITG) in phase 5 of the Coupled Model Intercomparison Project (CMIP5) historical simulations are investigated. Observations show a secular trend to this gradient over most of the twentieth century, with the southern lobe warming faster relative to its northern counterpart. A previous study of phase 3 of the CMIP (CMIP3) suggests that this trend is partially forced by anthropogenic sulfate aerosols. This analysis collectively confirms the partially forced trend for the CMIP5 and by anthropogenic aerosols. Like the CMIP3, the CMIP5 also simulates a reversal in the AITG trend in the late 1970s, which was attributed to a leveling offof the anthropogenic aerosol influence and increased influence of greenhouse gases in the late twentieth century. Two (of 25) CMIP5 models, however, systematically simulate a twentieth-century trend opposite to observed, leading to some uncertainty regarding the forced nature of the AITG trend. Here, the observed AITG also exhibits a pronounced multidecadal modulation on top of the trend, associated with the Atlantic multidecadal oscillation (AMO). Motivated by a recent suggestion that the AMO is a forced response to aerosols, the causes of this multidecadal behavior were also examined. A few of the CMIP5 models analyzed do produce multidecadal AITG variations that are correlated to the observed AMO-like variation, but only one, the Hadley Centre Global Environmental Model, version 2 (HadGEM2), systematically simulates AMO-like behavior with both the requisite amplitude and phase. The CMIP5 simulations thus point to a robust aerosol influence on the historical AITG trend but not to the AMO-like multidecadal behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

GRIZZLY MODEL OF FULLY COUPLED HEAT TRANSFER, MOISTURE DIFFUSION, ALKALI-SILICA REACTION AND FRACTURING PROCESSES IN CONCRETE

Idaho National Laboratory (INL) has been developing a fully coupled thermo-hydro-mechanical-chemical (THMC) concrete model, referred as Grizzly, with the ultimate goal of reliably simulating and predicting the long-term performance of concrete structures subjected to various aging and degradation mechanisms. In this paper, we present a fully coupled concrete model implemented in our Grizzly code for simulating coupled moisture/heat transfer, stress-dependent ASR swelling and elastic/inelastic mechanics. Example simulation results and preliminary model validation results against experimental data reported in the literature will be provided also in this paper. The close matches between the simulation results and experiments clearly reveal the potential of the fully coupled solution strategy for solving tightly coupled THMC concrete models. We also present simulation results of using a rate-dependent damaging mechanics model for simulating the ASR swelling induced degradation of concrete structures and its feedback to moisture and heat transport processes.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS2008

Gas Core Reactor Numerical Simulation Using a Coupled MHD-MCNP Model

Analysis is provided in this report of using two head-on magnetohydrodynamic (MHD) shocks to achieve supercritical nuclear fission in an axially elongated cylinder filled with UF4 gas as an energy source for deep space missions. The motivation for each aspect of the design is explained and supported by theory and numerical simulations. A subsequent report will provide detail on relevant experimental work to validate the concept. Here the focus is on the theory of and simulations for the proposed gas core reactor conceptual design from the onset of shock generations to the supercritical state achieved when the shocks collide. The MHD model is coupled to a standard nuclear code (MCNP) to observe the neutron flux and fission power attributed to the supercritical state brought about by the shock collisions. Throughout the modeling, realistic parameters are used for the initial ambient gaseous state and currents to ensure a resulting supercritical state upon shock collisions.

Kazeminezhad, F.↗
DOE OSTI2016

Simulation of Hydraulic and Natural Fracture Interaction Using a Coupled DFN-DEM Model

The presence of natural fractures will usually result in a complex fracture network due to the interactions between hydraulic and natural fracture. The reactivation of natural fractures can generally provide additional flow paths from formation to wellbore which play a crucial role in improving the hydrocarbon recovery in these ultra-low permeability reservoir. Thus, accurate description of the geometry of discrete fractures and bedding is highly desired for accurate flow and production predictions. Compared to conventional continuum models that implicitly represent the discrete feature, Discrete Fracture Network (DFN) models could realistically model the connectivity of discontinuities at both reservoir scale and well scale. In this work, a new hybrid numerical model that couples Discrete Fracture Network (DFN) and Dual-Lattice Discrete Element Method (DL-DEM) is proposed to investigate the interaction between hydraulic fracture and natural fractures. Based on the proposed model, the effects of natural fracture orientation, density and injection properties on hydraulic-natural fractures interaction are investigated.

Discrete Fracture Network↗
NASA NTRSDate not supplied

Dynamical Sensitivity in Response to a Wide Range of Abrupt CO 2 Forcings

An improved understanding of dynamical variability is important for projecting future changes in extratropical weather and the interaction between the extratropical troposphere and the tropics, stratosphere, cryosphere and ocean. Despite their simplicity, the abrupt 2x- and 4xCO 2 forcing simulations from the Coupled Model Intercomparison Project (CMIP) Phase 6 DECK experiments enable a mechanistic look into the forcing and feedback response characteristics of models that can be unambiguously attributed to an increase in carbon dioxide concentrations. Thus, while typically used to evaluate the climate sensitivity in models, with a primary focus on global surface temperature change, here we focus on several measures of extratropical variability, including projected changes in the storm tracks and in stratospheric polar vortex variability. Results are primarily based on simulations produced using low- and high-top versions of the NASA Goddard Institute for Space Studies Model (ModelE) but results are also presented from the larger CMIP6 multi-model ensemble. In addition to the 2x- and 4xCO 2 simulations, we also explore the linearity of the response of extratropical dynamical variability in ModelE to varying levels of CO 2 spanning the range 1/8-8xCO 2 . In particular, we show that the expansion of the jet streams varies nonlinearly with increasing CO 2 , especially in the Northern Hemisphere, and that this can be interpreted in terms of variations in ocean heat transport. The impact of composition feedbacks on changes in variability is also discussed.

dynamical variability↗
NASA NTRS2019

Assessment of Satellite Sea Surface Salinity Products Using a Coupled ENSO Prediction Model

Much work has gone into revising and updating algorithms for converting satellite-measured radiances to useful ocean variables like sea surface salinity (e.g. SMOS - Boutin et al., 2017, SMAP - Fore et al., 2016 and Aquarius - Meissner et al., 2018). As part of our Ocean Salinity Science Team work, we utilize an intermediate-complexity air/sea coupled model as a transfer function to test if more mature satellite SSS model algorithms actually improve ENSO forecast skill. For initialization of the coupled forecast, we demonstrate that the positive impact of SSS assimilation is brought about by surface freshening near the eastern edge of the western Pacific warm pool and density changes that lead to shallower mixed layer between 10S-5N. In addition, salting near the ITCZ leads to a deepening of the mixed layer and thermocline near 8N. These patterns together provide the background state to amplify equatorial Kelvin waves and improve ENSO hindcasts (Hackert et al., 2019). Here we extend this work to compare the impact of various pairs of original and improved satellite SSS algorithms. For instance we compare SMAP V4.1 with the latest, SMAP V4.2, to see what impact algorithm improvements may have on ENSO forecasts. SSS observations are tested on seasonal to interannual variability of tropical Indo-Pacific Ocean dynamics as well as on dynamical ENSO forecasts by initializing twelve-month forecasts for each month of available data. All experiments assimilate satellite sea level (SL), sea surface temperature (SST), and in situ subsurface temperature and salinity observations (Tz, Sz). Additionally various satellite, blended, and in-situ SSS products are assimilated. We find that including satellite SSS significantly improves Niño3.4 sea surface temperature anomaly validation, more mature SSS model algorithms are generally improving ENSO forecasts over time, and more satellite SSS data coverage helps to extend useful forecasts.

Ocean↗
DOE OSTI2016

Changes in future air quality, deposition, and aerosol-cloud interactions under future climate and emission scenarios

The prospect of global climate change will have wide scale impacts, such as ecological stress and human health hazards. One aspect of concern is future changes in air quality that will result from changes in both meteorological forcing and air pollutant emissions. In this study, the GU-WRF/Chem model is employed to simulate the impact of changing climate and emissions following the IPCC AR4 SRES A1B scenario. An average of 4 future years (2020, 2030, 2040, and 2050) is compared against an average of 2 current years (2001 and 2010). Under this scenario, by the Mid-21st century global air quality is projected to degrade with a global average increase of 2.5 ppb in the maximum 8-hr O 3 level and of 0.3 mg m 3 in 24-hr average PM2.5. However, PM2.5 changes are more regional due to regional variations in primary aerosol emissions and emissions of gaseous precursor for secondary PM2.5. Increasing NOx emissions in this scenario combines with a wetter climate elevating levels of OH, HO 2 , H 2 O 2 , and the nitrate radical and increasing the atmosphere’s near surface oxidation state. This differs from findings under the RCP scenarios that experience declines in OH from reduced NOx emissions, stratospheric recovery of O 3 , and increases in CH 4 and VOCs. Increasing NO x and O 3 levels enhances the nitrogen and O 3 deposition, indicating potentially enhanced crop damage and ecosystem stress under this scenario. The enhanced global aerosol level results in enhancements in aerosol optical depth, cloud droplet number concentration, and cloud optical thickness. This leads to dimming at the Earth’s surface with a global average reduction in shortwave radiation of 1.2 W m 2 . This enhanced dimming leads to a more moderate warming trend and different trends in radiation than those found in NCAR’s CCSM simulation, which does not include the advanced chemistry and aerosol treatment of GU-WRF/Chem and cannot simulate the impacts of changing climate and emissions with the same level of detailed treatments. This study indicates that effective climate mitigation and emission control strategies are needed to prevent future health impact and ecosystem stress. Further, studies that are used to develop these strategies should use fully coupled models with sophisticated chemical and aerosol-interaction treatments that can provide a more realistic representation of the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

Warming Overwhelms the Efficacy of Wet Conditions to Moderate Extreme Heat and Atmospheric Aridity Across the Central Plains

While the relationships between dry land surface conditions, heat, and aridity have been well-established, few studies have addressed whether global warming will affect the ability of wet conditions to moderate high temperatures and atmospheric aridity. Using Coupled Model Intercomparison Project Phase 6 models, we demonstrate that absolute changes in the monthly maximum temperature distribution during Central North American summers strongly outweigh the historical cooling effect of high precipitation and soil moisture conditions. Although wet conditions nearly always prevent concurrent extreme temperatures in the baseline period, these conditions are 40%–48% and 96%–98% less effective at 1 and 2° of global warming, respectively. However, high precipitation and soil moisture partially retain the ability to constrain concurrent high vapor pressure deficit conditions below historical thresholds at 1–2° of warming. Our results highlight the growing vulnerability of Central North America to warmer temperatures and drier atmospheric conditions, even during periods of high precipitation and soil moisture.

global warming↗
DOE OSTI2021

Assessing the Representation of Synoptic Variability Associated With California Extreme Precipitation in CMIP6 Models

Abstract Days of extreme precipitation over California are evaluated in Coupled Model Intercomparison Project Phase 6 (CMIP6) models and the ERA‐Interim reanalysis. In the current climate, the model spread in composited precipitation on extreme precipitation days is closely related to the magnitude of composited integrated vapor transport (IVT) across models, a proxy for the intensity of atmospheric rivers. Most models underestimate the magnitude of IVT associated with extreme precipitation, according to ERA‐Interim. This is due mostly to the contribution of moisture, which almost all models overestimate, while the contribution of lower‐tropospheric wind speed is generally closer to the reanalyses. Moreover, most models underestimate the variance in the latitude of maxima of numerous variables among days of extreme California precipitation. That is, in the general circulation models there is a lack of diversity in the latitude of the disturbances bringing winter precipitation to California. In the future climate, most models project a decrease in the frequency of southward‐displaced disturbances among California extreme precipitation days. Hence, the greatest increases in extreme precipitation are over northern California. However, the historical underestimate of the latitudinal variance of disturbances calls into question the reliability of these projections. This bias should be especially considered for dynamical downscaling efforts over the region.

Norris, Jesse↗
DOE OSTI2015

Application of an Online-Coupled Regional Climate Model, WRF-CAM5, over East Asia for Examination of Ice Nucleation Schemes: Part II. Sensitivity to Heterogeneous Ice Nucleation Parameterizations and Dust Emissions

Aerosol particles can affect cloud microphysical properties by serving as ice nuclei (IN). Large uncertainties exist in the ice nucleation parameterizations (INPs) used in current climate models. In this Part II paper, to examine the sensitivity of the model predictions to different heterogeneous INPs, WRF-CAM5 simulation using the INP of Niemand et al. (N12) [1] is conducted over East Asia for two full years, 2006 and 2011, and compared with simulation using the INP of Meyers et al. (M92) [2], which is the original INP used in CAM5. M92 calculates the nucleated ice particle concentration as a function of ice supersaturation, while N12 represents the nucleated ice particle concentration as a function of temperature and the number concentrations and surface areas of dust particles. Compared to M92, the WRF-CAM5 simulation with N12 produces significantly higher nucleated ice crystal number concentrations (ICNCs) in the northern domain where dust sources are located, leading to significantly higher cloud ice number and mass concentrations and ice water path, but the opposite is true in the southern domain where temperatures and moistures play a more important role in ice formation. Overall, the simulation with N12 gives lower downward shortwave radiation but higher downward longwave radiation, cloud liquid water path, cloud droplet number concentrations, and cloud optical depth. The increase in cloud optical depth and the decrease in downward solar flux result in a stronger shortwave and longwave cloud forcing, and decreases temperature at 2-m and precipitation. Changes in temperature and radiation lower surface concentrations of OH, O3, SO42−, and PM2.5, but increase surface concentrations of CO, NO2, and SO2 over most of the domain. By acting as cloud condensation nuclei (CCN) and IN, dust particles have different impacts on cloud water and ice number concentrations, radiation, and temperature at 2-m and precipitation depending on whether the dominant role of dust is CCN or IN. These results indicate the importance of the heterogeneous ice nucleation treatments and dust emissions in accurately simulating regional climate and air quality.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Linking process and structure in the friction stir scribe joining of dissimilar materials: A computational approach with experimental support

We present that friction stir welding (FSW) is a popular technique to join dissimilar materials in numerous applications. The solid state nature of the process enables joining materials with strikingly different physical properties. For welds in lap configuration, an enhancement to this technology is made by introducing a short, hard insert, referred to as a cutting-scribe, at the bottom of the tool pin. The cutting-scribe induces deformation in the bottom plate which leads to the formation of mechanical interlocks or hook like structures at the interface of two materials. A thermo-mechanical computational model employing a coupled Eulerian-Lagrangian approach is developed to quantitatively capture the morphology of these interlocks during the FSW process. Simulations using this model are validated by experimental observations. In conclusion, the identified interface morphology coupled with the predicted temperature field from this process–structure model can be used to estimate the post-weld microstructure and joint strength.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Reproducible, Component-based Modeling with TopoFlow, A Spatial Hydrologic Modeling Toolkit

Modern geoscientists have online access to an abundance of different data sets and models, but these resources differ from each other in myriad ways and this heterogeneity works against interoperability as well as reproducibility. The purpose of this paper is to illustrate the main issues and some best practices for addressing the challenge of reproducible science in the context of a relatively simple hydrologic modeling study for a small Arctic watershed near Fairbanks, Alaska. This study requires several different types of input data in addition to several, coupled model components. All data sets, model components and processing scripts (e.g. for preparation of data and figures, and for analysis of model output) are fully documented and made available online at persistent URLs. Similarly, all source code for the models and scripts is open-source, version controlled and made available online via GitHub. Each model component has a Basic Model Interface (BMI) to simplify coupling and its own HTML help page that includes a list of all equations and variables used. The set of all model components (TopoFlow) has also been made available as a Python package for easy installation. Three different graphical user interfaces for setting up TopoFlow runs are described, including one that allows model components to run and be coupled as web services.

58 GEOSCIENCES↗
NASA NTRS2016

Modeling the Ionosphere-Thermosphere Response to a Geomagnetic Storm Using Physics-based Magnetospheric Energy Input: OpenGGCM-CTIM Results

The magnetosphere is a major source of energy for the Earth's ionosphere and thermosphere (IT) system. Current IT models drive the upper atmosphere using empirically calculated magnetospheric energy input. Thus, they do not sufficiently capture the storm-time dynamics, particularly at high latitudes. To improve the prediction capability of IT models, a physics-based magnetospheric input is necessary. Here, we use the Open Global General Circulation Model (OpenGGCM) coupled with the Coupled Thermosphere Ionosphere Model (CTIM). OpenGGCM calculates a three-dimensional global magnetosphere and a two-dimensional high-latitude ionosphere by solving resistive magnetohydrodynamic (MHD) equations with solar wind input. CTIM calculates a global thermosphere and a high-latitude ionosphere in three dimensions using realistic magnetospheric inputs from the OpenGGCM. We investigate whether the coupled model improves the storm-time IT responses by simulating a geomagnetic storm that is preceded by a strong solar wind pressure front on August 24, 2005. We compare the OpenGGCM-CTIM results with low-earth-orbit satellite observations and with the model results of Coupled Thermosphere-Ionosphere-Plasmasphere electrodynamics (CTIPe). CTIPe is an up-to-date version of CTIM that incorporates more IT dynamics such as a low-latitude ionosphere and a plasmasphere, but uses empirical magnetospheric input. OpenGGCMCTIM reproduces localized neutral density peaks at approx. 400 km altitude in the high-latitude dayside regions in agreement with in situ observations during the pressure shock and the early phase of the storm. Although CTIPe is in some sense a much superior model than CTIM, it misses these localized enhancements. Unlike the CTIPe empirical input models, OpenGGCM-CTIM more faithfully produces localized increases of both auroral precipitation and ionospheric electric fields near the high-latitude dayside region after the pressure shock and after the storm onset, which in turn effectively heats the thermosphere and causes the neutral density increase at 400 km altitude.

Connor, Hyunju K.↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2014

Evidence of resonant mode coupling and the relationship between low and high frequencies in a rapidly rotating a star

In the theory of resonant mode coupling, the parent and child modes are directly related in frequency and phase. The oscillations present in the fast rotating δ Sct star KIC 8054146 allow us to test the most general and generic aspects of such a theory. The only direct way to separate the parent and coupled (child) modes is to examine the correlations in amplitude variability between the different frequencies. For the dominant family of related frequencies, only a single mode and a triplet are the origins of nine dominant frequency peaks ranging from 2.93 to 66.30 cycles day{sup –1} (as well as dozens of small-amplitude combination modes and a predicted and detected third high-frequency triplet). The mode-coupling model correctly predicts the large amplitude variations of the coupled modes as a product of the amplitudes of the parent modes, while the phase changes are also correctly modeled. This differs from the behavior of 'normal' combination frequencies in that the amplitudes are three orders of magnitude larger and may exceed even the amplitudes of the parent modes. We show that two dominant low frequencies at 5.86 and 2.93 cycles day{sup –1} in the gravity-mode region are not harmonics of each other, and their properties follow those of the almost equidistant high-frequency triplet. We note that the previously puzzling situation of finding two strong peaks in the low-frequency region related by nearly a factor of two in frequency has been seen in other δ Sct stars as well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A Case Study of Urbanization Impact on Summer Precipitation in the Greater Beijing Metropolitan Area. Urban Heat Island Versus Aerosol Effects

Convection-resolving ensemble simulations using the WRF-Chem model coupled with a single-layer Urban Canopy Model (UCM) are conducted to investigate the individual and combined impacts of land use and anthropogenic pollutant emissions from urbanization on a heavy rainfall event in the Greater Beijing Metropolitan Area (GBMA) in China. The simulation with the urbanization effect included generally captures the spatial pattern and temporal variation of the rainfall event. An improvement of precipitation is found in the experiment including aerosol effect on both clouds and radiation. The expanded urban land cover and increased aerosols have an opposite effect on precipitation processes, with the latter playing a more dominant role, leading to suppressed convection and rainfall over the upstream (northwest) area, and enhanced convection and more precipitation in the downstream (southeast) region of the GBMA. In addition, the influence of aerosol indirect effect is found to overwhelm that of direct effect on precipitation in this rainfall event. Increased aerosols induce more cloud droplets with smaller size, which favors evaporative cooling and reduce updrafts and suppress convection over the upstream (northwest) region in the early stage of the rainfall event. As the rainfall system propagates southeastward, more latent heat is released due to the freezing of larger number of smaller cloud drops that are lofted above the freezing level, which is responsible for the increased updraft strength and convective invigoration over the downstream (southeast) area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

Simulation of Urban Climate with High-Resolution WRF Model: A Case Study in Nanjing, China

In this study, urban climate in Nanjing of eastern China is simulated using 1-km resolution Weather Research and Forecasting (WRF) model coupled with a single-layer Urban Canopy Model. Based on the 10-summer simulation results from 2000 to 2009 we find that the WRF model is capable of capturing the high-resolution features of urban climate over Nanjing area. Although WRF underestimates the total precipitation amount, the model performs well in simulating the surface air temperature, relative humidity, and precipitation frequency, diurnal cycle and inter-annual variability. We find that extremely hot events occur most frequently in urban area, with daily maximum (minimum) temperature exceeding 36ºC (28ºC) in around 40% (32%) of days. Urban Heat Island (UHI) effect at surface is more evident during nighttime than daytime, with 20% of cases the UHI intensity above 2.5ºC at night. However, The UHI affects the vertical structure of Planet Boundary Layer (PBL) more deeply during daytime than nighttime. Net gain for latent heat and net radiation is larger over urban than rural surface during daytime. Correspondingly, net loss of sensible heat and ground heat are larger over urban surface resulting from warmer urban skin. Because of different diurnal characteristics of urban-rural differences in the latent heat, ground heat and other energy fluxes, the near surface UHI intensity exhibits a very complex diurnal feature. UHI effect is stronger in days with less cloud or lower wind speed. Model results reveal a larger precipitation frequency over urban area, mainly contributed by the light rain events (<10 mm day-1). Consistent with satellite dataset, around 10-20% more precipitation occurs in urban than rural area at afternoon induced by more unstable urban PBL, which induces a strong vertical atmospheric mixing and upward moisture transport. A significant enhancement of precipitation is found in the downwind region of urban in our simulations in the afternoon.

Precipitation↗