Subkilometer processes are critical to the physics of aerosol‐cloud interaction (ACI) but have been dependent on parameterizations in global model simulations. We thus report the strength of ACI in the Ultra‐Parameterized Community Atmosphere Model (UPCAM), a multiscale climate model that uses coarse exterior resolution to embed explicit cloud‐resolving models with enough resolution (250 m horizontal, 20 m vertical) to quasi‐resolve subkilometer eddies. To investigate the impact on ACIs, UPCAM's simulations are compared to a coarser multiscale model with 4 km horizontal resolution. UPCAM produces cloud droplet number concentrations (
Model intercomparison studies often report a large spread in simulation results, but quantifying the causes of these differences is hindered by the fact that several processes contribute to the model spread simultaneously. Here we use the Multi‐Scale Infrastructure for Chemistry and Aerosols (MUSICA) version 0 to investigate the model resolution dependencies of simulated chemical species, with a focus on the differences between global uniform grid and regional refinement grid simulations with the same modeling framework. We construct two global (ne30 [∼112 km] and ne60 [∼56 km]) and two regional refinement grids over Korea (ne30x8 [∼14 km] and ne30x16 [∼7 km]). The grid resolution can change chemical concentrations by an order of magnitude in the boundary layer, and the importance increases as the species' reactivity increases (e.g., up to 50% and 1,000% changes for ethane and xylenes, respectively). The diurnal cycle of oxidants (OH, O3, and NO3) also varies with the grid resolution, which leads to different oxidation pathways of volatile organic compounds (e.g., the fraction of monoterpenes reacting with NO3in Seoul around midnight is 90% for ne30, but 65% for ne30x16). The models with high‐resolution grids usually do a better job at reproducing aircraft observations during the KORUS‐AQ campaign, but not always, implying compensating errors in the coarse grid simulations. For example, ozone is better reproduced by the coarse grid due to the artificial mixing of NOx. When developing new chemical mechanisms and evaluating models over urban areas, the uncertainties associated with model resolution should be considered.
more » « less- NSF-PAR ID:
- 10435590
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 15
- Issue:
- 7
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract N dN d -
more » « less
SUMMARY Precisely constraining the source parameters of large earthquakes is one of the primary objectives of seismology. However, the quality of the results relies on the quality of synthetic earth response. Although earth structure is laterally heterogeneous, particularly at shallow depth, most earthquake source studies at the global scale rely on the Green's functions calculated with radially symmetric (1-D) earth structure. To avoid the impact of inaccurate Green's functions, these conventional source studies use a limited set of seismic phases, such as long-period seismic waves, broad-band P and S waves in teleseismic distances (30° < ∆ < 90°), and strong ground motion records at close-fault stations. The enriched information embedded in the broad-band seismograms recorded by global and regional networks is largely ignored, limiting the spatiotemporal resolution. Here we calculate 3-D strain Green's functions at 30 GSN stations for source regions of 9 selected global earthquakes and one earthquake-prone area (California), with frequency up to 67 mHz (15 s), using SPECFEM3D_GLOBE and the reciprocity theorem. The 3-D SEM mesh model is composed of mantle model S40RTS, crustal model CRUST2.0 and surface topography ETOPO2. We surround each target event with grids in horizontal spacing of 5 km and vertical spacing of 2.0–3.0 km, allowing us to investigate not only the main shock but also the background seismicity. In total, the response at over 210 000 source points is calculated in simulation. The number of earthquakes, including different focal mechanisms, centroid depth range and tectonic background, could further increase without additional computational cost if they were properly selected to avoid overloading individual CPUs. The storage requirement can be reduced by two orders of magnitude if the output strain Green's functions are stored for periods over 15 s. We quantitatively evaluate the quality of these 3-D synthetic seismograms, which are frequency and phase dependent, for each source region using nearby aftershocks, before using them to constrain the focal mechanisms and slip distribution. Case studies show that using a 3-D earth model significantly improves the waveform similarity, agreement in amplitude and arrival time of seismic phases with the observations. The limitations of current 3-D models are still notable, dependent on seismic phases and frequency range. The 3-D synthetic seismograms cannot well match the high frequency (>40 mHz) S wave and (>20 mHz) Rayleigh wave yet. Though the mean time-shifts are close to zero, the standard deviations are notable. Careful calibration using the records of nearby better located earthquakes is still recommended to take full advantage of better waveform similarity due to the use of 3-D models. Our results indicate that it is now feasible to systematically study global large earthquakes using full 3-D earth response in a global scale.
-
more » « less
Abstract Many chemical processes depend non‐linearly on temperature. Gravity‐wave‐induced temperature perturbations have been shown to affect atmospheric chemistry, but accounting for this process in chemistry‐climate models has been a challenge because many gravity waves have scales smaller than the typical model resolution. Here, we present a method to account for subgrid‐scale orographic gravity‐wave‐induced temperature perturbations on the global scale for the Whole Atmosphere Community Climate Model. Temperature perturbation amplitudes
consistent with the model's subgrid‐scale gravity wave parameterization are derived and then used as a sinusoidal temperature perturbation in the model's chemistry solver. Because of limitations in the parameterization, we explore scaling of between 0.6 and 1 based on comparisons to altitude‐dependent distributions of satellite and reanalysis data, where we discuss uncertainties. We probe the impact on the chemistry from the grid‐point to global scales, and show that the parameterization is able to represent mountain wave events as reported by previous literature. The gravity waves for example, lead to increased surface area densities of stratospheric aerosols. This increases chlorine activation, with impacts on the associated chemical composition. We obtain large local changes in some chemical species (e.g., active chlorine, NOx, N2O5) which are likely to be important for comparisons to airborne or satellite observations, but the changes to ozone loss are more modest. This approach enables the chemistry‐climate modeling community to account for subgrid‐scale gravity wave temperature perturbations interactively, consistent with the internal parameterizations and are expected to yield more realistic interactions and better representation of the chemistry. -
more » « less
Abstract A severe derecho impacted the Midwestern United States on 10 August 2020, causing over $12 billion (U.S. dollars) in damage, and producing peak winds estimated at 63 m s−1, with the worst impacts in Iowa. The event was not forecast well by operational forecasters, nor even by operational and quasi-operational convection-allowing models. In the present study, nine simulations are performed using the Limited Area Model version of the Finite-Volume-Cubed-Sphere model (FV3-LAM) with three horizontal grid spacings and two physics suites. In addition, when a prototype of the Rapid Refresh Forecast System (RRFS) physics is used, sensitivity tests are performed to examine the impact of using the Grell–Freitas (GF) convective scheme. Several unusual results are obtained. With both the RRFS (not using GF) and Global Forecast System (GFS) physics suites, simulations using relatively coarse 13- and 25-km horizontal grid spacing do a much better job of showing an organized convective system in Iowa during the daylight hours of 10 August than the 3-km grid spacing runs. In addition, the RRFS run with 25-km grid spacing becomes much worse when the GF convective scheme is used. The 3-km RRFS run that does not use the GF scheme develops spurious nocturnal convection the night before the derecho, removing instability and preventing the derecho from being simulated at all. When GF is used, the spurious storms are removed and an excellent forecast is obtained with an intense bowing echo, exceptionally strong cold pool, and roughly 50 m s−1surface wind gusts.
-
more » « less
Abstract Prior research indicates that land use and land cover change (LULCC) in the central United States has led to significant changes in surface climate. The spatial resolution of simulations is particularly relevant in this region due to its influence on model skill in capturing mesoscale convective systems (MCSs) and on representing the spatial heterogeneity. Recent advances in Earth system models (ESMs) make it feasible to use variable resolution (VR) meshes to study regional impacts of LULCC while avoiding inconsistencies introduced by lateral boundary conditions typically seen in limited area models. Here, we present numerical experiments using the Community Earth System Model version 2–VR to evaluate (1) the influence of resolution and land use on model skill and (2) impacts of LULCC over the central United States at different resolutions. These simulations are configured either on the 1° grid or a VR grid with grid refinement to 1/8° over the contiguous United States for the period of 1984–2010 with two alternative land use data sets corresponding to the preindustrial and present day states. Our results show that skill in simulating precipitation over the central United States is primarily dependent on resolution, whereas skill in simulating 2‐m temperature is more dependent on accurate land use. The VR experiments show stronger LULCC‐induced precipitation increases over the Midwest in May and June, corresponding to an increase in the number of MCS‐like features and a more conductive thermodynamic environment for convection. Our study demonstrates the potential of using VR ESMs for hydroclimatic simulations in regions with significant LULCC.
