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 685 records · Page 38

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

Abstract. 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.

54 ENVIRONMENTAL SCIENCES↗

Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport

In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Manipulating electromagnetic responses of metal wires at the deep subwavelength scale via both near- and far-field couplings

A hybrid coupling model containing both near- and far-field couplings is developed for radiating two-resonator structures. We demonstrate that the near- and far-field couplings make distinguished contributions to electromagnetic responses. Compared to the classical electromagnetically induced transparency configurations, the presence of far-field coupling provides more flexibility in tuning lineshapes. Planar metamaterials composed of metal wires are designed based on this model, and various electromagnetic responses are experimentally observed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗

A coupled equilibrium boundary layer model with stable water isotopes and its application to local water recycling

The contribution of evapotranspiration (ET) to regional precipitation, known as “local water recycling”, is a key process in the water cycle that can affect water management. However, the amount of planetary boundary layer (PBL) moisture that arises from ET is highly uncertain due to complex atmosphere and land surface conditions. In this study, an idealized two-layer equilibrium boundary layer model was coupled with a stable water isotope module including HDO and H 2 18 O to constrain PBL growing season water transport processes. The model was validated using turbulent heat fluxes and isotope ratios of water vapor (δ m ) and precipitation (δ P ) measured at a cropland site and a nearby tall tower in the Upper Midwest, United States. The results show that the PBL equilibrium features of δ m and δ P are well-constrained by thermal and moisture equilibrium in the PBL. For this study region, the summer values of rain evaporation fraction (f) and local water recycling ratio (LRR) are estimated to be 0.09 and 0.29 ± 0.12, respectively. The summer LRR values for the years 2006–2010 were 0.35, 0.36, 0.17, 0.29, and 0.29, respectively. The small value of LRR in 2008 corresponded to a drought condition with the lowest precipitation and second lowest ET among the five years. The summer magnitude of the amount effect is –2.8‰ (mm day –1 ) –1 and –0.8‰ (mm day –1 ) –1 for HDO and H 2 18 O, respectively. The local water recycling is identified as a significant factor influencing the continental isotope effect. Cropland has likely changed the regional LRR by –7.6 to 19.5% under different pre-agriculture land use scenarios. In conclusion, the feedback processes revealed here indicate that local water recycling is expected to be weakened under drought conditions, but it will be enhanced if irrigation is applied more intensely with more frequent drought events as the climate continues to warm.

54 ENVIRONMENTAL SCIENCES↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗