Precipitation amount (A), frequency (F), intensity (I), and duration (D) are important properties of precipitation, but their estimates are sensitive to data resolution. This study investigates this resolution dependence, and the influences of different model physics, by analyzing simulations by the Community Atmospheric Model (CAM) version 4 (CAM4) and version 5 (CAM5) with varying grid sizes from ~0.25 to 2.0°. Results show that both CAM4 and CAM5 greatly overestimate F and D but underestimate I at all resolutions, despite realistic A. These biases partly result from too much parameterized (convective) precipitation with high F and D but low I. Different cloud microphysics schemes contribute to the precipitation differences between CAM4 and CAM5. The A, F, I, and D of convective and nonconvective precipitation react differently to grid‐size decreases, leading to the large decreases in F and D but increases in the I for total precipitation as model resolution increases. This resolution dependence results from the increased probability of precipitation over a larger area (area aggregation effect, which is smaller than in observations) and the varying performance of model physics under changing resolution (model adjustment effect), which roughly enhances the aggregation‐induced dependence. Finer grid sizes not only increase resolved precipitation, which has higher intensity and thus improves overall precipitation intensity in CAM, but also reduce the area aggregation effect. Thus, the long‐standing
Advances in high‐performance computing make it possible to run atmospheric general circulation models (AGCMs) over an increasingly wider range of grid resolutions, using either globally uniform or variable‐resolution grids. In principle, this is an exciting opportunity to resolve atmospheric process and scales in a global model and in unprecedented detail, but in practice this grid flexibility is incompatible with the non‐ or weakly converging solutions with increasing horizontal resolution that have long characterized AGCMs. In the the Community Atmosphere Model (CAM), there are robust sensitivities to horizontal resolution that have persisted since the model was first introduced over thirty years ago; the atmosphere progressively dries and becomes less cloudy with resolution, and parametrized deep convective precipitation decreases at the expense of stratiform precipitation. This study documents a convergence experiment using CAM, version 6, and argues that a unifying cause, the sensitivity of resolved dynamical modes to native grid resolution, feeds back into other model components and explains these robust sensitivities to resolution. The increasing magnitudes of resolved vertical velocities with resolution are shown to fit an analytic scaling derived for the equations of motion at hydrostatic scales. This trend in vertical velocities results in an increase in resolved upward moisture fluxes at cloud base, balanced by an increase in stratiform precipitation rates with resolution. Compensating, greater magnitude subsiding motion with resolution has previously been shown to dry out the atmosphere and reduce cloud cover. Here, it is shown that both the increase in condensational heating from stratiform cloud formation and greater subsidence drying contribute to an increase in atmospheric stability with resolution, reducing the activity of parametrized convection. The impact of changing the vertical velocity field with native grid resolution cannot be ignored in any effort to recover convergent solutions in AGCMs, and, in particular, the development of scale‐aware physical parametrizations.
more » « less- NSF-PAR ID:
- 10371242
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Quarterly Journal of the Royal Meteorological Society
- Volume:
- 146
- Issue:
- 733
- ISSN:
- 0035-9009
- Page Range / eLocation ID:
- p. 3789-3807
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract drizzling problem in climate models may be mitigated by increasing model resolution and modifying model physics to suppress parameterized convective precipitation and enhance resolved nonconvective precipitation. -
more » « less
Abstract A new configuration of the Community Earth System Model (CESM)/Community Atmosphere Model with full chemistry (CAM‐chem) supporting the capability of horizontal mesh refinement through the use of the spectral element (SE) dynamical core is developed and called CESM/CAM‐chem‐SE. Horizontal mesh refinement in CESM/CAM‐chem‐SE is unique and novel in that pollutants such as ozone are accurately represented at human exposure relevant scales while also directly including global feedbacks. CESM/CAM‐chem‐SE with mesh refinement down to ∼14 km over the conterminous US (CONUS) is the beginning of the Multi‐Scale Infrastructure for Chemistry and Aerosols (MUSICAv0). Here, MUSICAv0 is evaluated and used to better understand how horizontal resolution and chemical complexity impact ozone and ozone precursors over CONUS as compared to measurements from five aircraft campaigns, which occurred in 2013. This field campaign analysis demonstrates the importance of using finer horizontal resolution to accurately simulate ozone precursors such as nitrogen oxides and carbon monoxide. In general, the impact of using more complex chemistry on ozone and other oxidation products is more pronounced when using finer horizontal resolution where a larger number of chemical regimes are resolved. Large model biases for ozone near the surface remain in the Southeast US as compared to the aircraft observations even with updated chemistry and finer horizontal resolution. This suggests a need for adding the capability of replacing sections of global emission inventories with regional inventories, increasing the vertical resolution in the planetary boundary layer, and reducing model biases in meteorological variables such as temperature and clouds.
-
more » « less
Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement.
Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.
-
more » « less
Within the pycnocline, where diapycnal mixing is suppressed, both the vertical movement (uplift) of isopycnal surfaces and upward motion along sloping isopycnals supply nutrients to the euphotic layer, but the relative importance of each of these mechanisms is unknown. We present a method for decomposing vertical velocity w into two components in a Lagrangian frame: vertical velocity along sloping isopycnal surfaces [Formula: see text] and the adiabatic vertical velocity of isopycnal surfaces [Formula: see text]. We show that [Formula: see text], where [Formula: see text] is the isopycnal slope and [Formula: see text] is the geometric aspect ratio of the flow, and that [Formula: see text] accounts for 10%–25% of the total vertical velocity w for isopycnal slopes representative of the midlatitude pycnocline. We perform the decomposition of w in a process study model of a midlatitude eddying flow field generated with a range of isopycnal slopes. A spectral decomposition of the velocity components shows that while [Formula: see text] is the largest contributor to vertical velocity, [Formula: see text] is of comparable magnitude at horizontal scales less than about 10 km, that is, at submesoscales. Increasing the horizontal grid resolution of models is known to increase vertical velocity; this increase is disproportionately due to better resolution of [Formula: see text], as is shown here by comparing 1- and 4-km resolution model runs. Along-isopycnal vertical transport can be an important contributor to the vertical flux of tracers, including oxygen, nutrients, and chlorophyll, although we find weak covariance between vertical velocity and nutrient anomaly in our model.
-
more » « less
Abstract Characteristics of tropical cyclones (TCs) in global climate models (GCMs) are known to be influenced by details of the model configurations, including horizontal resolution and parameterization schemes. Understanding model-to-model differences in TC characteristics is a prerequisite for reducing uncertainty in future TC activity projections by GCMs. This study performs a process-level examination of TC structures in eight GCM simulations that span a range of horizontal resolutions from 1° to 0.25°. A recently developed set of process-oriented diagnostics is used to examine the azimuthally averaged wind and thermodynamic structures of the GCM-simulated TCs. Results indicate that the inner-core wind structures of simulated TCs are more strongly constrained by the horizontal resolutions of the models than are the thermodynamic structures of those TCs. As expected, the structures of TC circulations become more realistic with smaller horizontal grid spacing, such that the radii of maximum wind (RMW) become smaller, and the maximum vertical velocities occur off the center. However, the RMWs are still too large, especially at higher intensities, and there are rising motions occurring at the storm centers, inconsistently with observations. The distributions of precipitation, moisture, and radiative and surface turbulent heat fluxes around TCs are diverse, even across models with similar horizontal resolutions. At the same horizontal resolution, models that produce greater rainfall in the inner-core regions tend to simulate stronger TCs. When TCs are weak, the radial gradient of net column radiative flux convergence is comparable to that of surface turbulent heat fluxes, emphasizing the importance of cloud–radiative feedbacks during the early developmental phases of TCs.
