skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: A Method for Estimating Global Subgrid‐Scale Orographic Gravity‐Wave Temperature Perturbations in Chemistry‐Climate Models

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 amplitudesconsistent 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 ofbetween 0.6 and 1 based on comparisons to altitude‐dependentdistributions 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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Advances in Modeling Earth Systems
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    TIMED/GUVI limb measurements and first‐principles simulations from the Thermosphere Ionosphere Electrodynamics Global Circulation Model (TIEGCM) are used to investigate thermospheric atomic oxygen (O) and molecular nitrogen (N2) responses in the middle thermosphere on a constant pressure surface (∼160 km) to the November 20 and 21, 2003 superstorm. The consistency between GUVI observations and TIEGCM simulated composition changes allows us to utilize TIEGCM outputs to investigate the storm‐time behaviors of O and N2systematically. Diagnostic analysis shows that horizontal and vertical advection are the two main processes that determine the storm‐induced perturbations in the middle thermosphere. Molecular diffusion has a relatively smaller magnitude than the two advection processes, acting to compensate for the changes caused by the transport partly. Contributions from chemistry and eddy diffusion are negligible. During the storm initial and main phases, composition variations at high latitudes are determined by both horizontal and vertical advection. At middle‐low latitudes, horizontal advection is the main driver for the composition changes where O mass mixing ratiodecreases (N2mass mixing ratioincreases); whereas horizontal and vertical advection combined to dominate the changes in the regions whereincreases (decreases). Over the entire storm period, horizontal advection plays a significant role in transporting high‐latitude composition perturbations globally. Our results also demonstrate that storm‐time temperature changes are not the direct cause of the composition perturbations on constant pressure surfaces.

    more » « less
  2. Abstract

    We explore the response of northeastern Pacific sea surface temperature (SST) to deglacial (16–7 ka) climate variability as recorded in‐based SST reconstructions spanning 65°N to 10°S. Included in the analysis is a new 23 kyr SST record from core NH8P from the northwest Mexican Margin. We isolate spatiotemporal patterns in regional SSTs with trend empirical orthogonal function (TEOF) analysis. The dominant TEOF mode reflects deglacial warming associated with rising. Tropical and subtropical SSTs correlated most strongly with this mode, suggesting that the thermodynamic response of the tropical eastern Pacific to greenhouse gas forcing was the dominant driver of regional SST change during deglaciation. The second TEOF mode reflects millennial‐scale variability and is most strongly expressed in subpolar SSTs. The synchronous timing between North Pacific and North Atlantic SST oscillations is evidence for the rapid transmission of millennial‐scale climate perturbations between the basins, likely through an atmospheric teleconnection. SSTs at NH8P have no correlation with either leading TEOF mode as there is minimal change in SST at this site after20 ka. A model simulation of the LGM indicates that glacial cooling was muted in much of the Eastern Pacific Warm Pool (EPWP), in which NH8P lies, due to reductions in latent heat flux. This suggests that the wind‐evaporation‐SST feedback was responsible for the attenuation of EPWP cooling. Overall, this study highlights the distinct latitudinal trends in the Pacific's response to deglaciation.

    more » « less
  3. Abstract

    We provide evidence that midlatitude postsunrise traveling ionospheric disturbances (TIDs) are comprised of electrified waves with an eastward propagation component. The post‐sunrise gravity wave (GW) wind‐induced dynamo action effectively generated periodic meridional polarization electric fields (PEFs), facilitating TID zonal propagation in a similar fashion as GW‐driven neutral perturbations. A combination of near‐simultaneous eastward and upward observations using the Millstone Hill incoherent scatter radar along with 2‐dimensional total electron content maps allowed resolution of TID vertical and horizontal propagation as well as zonal ion drifts(meridional PEFs). In multiple observations,oscillated in the early morning during periods when TIDs exhibited downward phase progression, 30–60 min period,140 m/s eastward speed, and 70 km vertical wavelength. Inside these TIDs, multiple flow vortexes occurred in a vertical‐zonal plane spanning the ionospheric topside and bottomside. Subsequently, PEFs weakened after a few hours as TID horizontal wavefronts rotated clockwise.

    more » « less
  4. Abstract

    Silicon stable isotope ratios (30Si) of over 150 stream water samples were measured during seven storm events in six small critical zone observatory (CZO) catchments spanning a wide range in climate (sub‐humid to wet, tropical) and lithology (granite, volcanic, and mixed sedimentary). Here we report a cross‐site analysis of this dataset to gain insight into stream30Si variability across low‐order catchments and to identify potential climate (i.e., runoff), hydrologic, lithologic, and biogeochemical controls on observed stream Si chemical and isotopic signatures. Event‐based30Si exhibit variability both within and across sites (−0.22‰ to +2.27‰) on the scale of what is observed globally in both small catchments and large rivers. Notably, each site shows distinct30Si signatures that are preserved even after normalization for bedrock composition. Successful characterization of observed cross‐site behavior requires the merging of two distinct frameworks in a novel combined model describing both non‐uniform fluid transit time distributions and multiple fractionating pathways in application to low‐order catchments. The combined model reveals that site‐specific architecture (i.e., biogeochemical reaction pathways and hydrologic routing) regulates stream silicon export signatures even when subject to extreme precipitation events.

    more » « less
  5. Abstract

    In this study, we report on turbulent mixing observed during the annual stratification cycle in the hypolimnetic waters of Lake Michigan (USA), highlighting stratified, convective, and transitional mixing periods. Measurements were collected using a combination of moored instruments and microstructure profiles. Observations during the stratified summer showed a shallow, wind‐driven surface mixed layer (SML) with locally elevated dissipation rates in the thermocline () potentially associated with internal wave shear. Below the thermocline, turbulence was weak () and buoyancy‐suppressed (< 8.5), with low hypolimnetic mixing rates () limiting benthic particle delivery. During the convective winter period, a diurnal cycle of radiative convection was observed over each day of measurement, where temperature overturns were directly correlated with elevated turbulence levels throughout the water column (;). A transitional mixing period was observed for spring conditions when surface temperatures were near the temperature of maximum density (TMD3.98) and the water column began to stably stratify. While small temperature gradients allowed strong mixing over the transitional period (), hypolimnetic velocity shear was overwhelmed by weakly stable stratification (;), limiting the development of the SML. These results highlight the importance of radiative convection for breaking down weak hypolimnetic stratification and driving energetic, full water column mixing during a substantial portion of the year (>100 days at our sample site). Ongoing surface water warming in the Laurentian Great Lakes is significantly reducing the annual impact of convective mixing, with important consequences for nutrient cycling, primary production, and benthic‐pelagic coupling.

    more » « less