skip to main content


Title: Atmospheric and Ionospheric Responses to Hunga‐Tonga Volcano Eruption Simulated by WACCM‐X
Abstract

High‐resolution Whole Atmosphere Community Climate Model with thermosphere/ionosphere extension is used to simulate the responses to the Hunga‐Tonga volcano eruption on 15 January 2022. Global propagation of the Lamb wave L’0and L’1pseudomodes are reproduced in the simulation, with the exponential growth of wave amplitudes with altitudes. The wavefront is vertical up to the lower thermosphere, and tilts outward above. These features are consistent with theoretical results. With simulated surface pressure perturbation agreeing with observations (∼100–250 Pa), thermospheric wind perturbations over 100 ms−1are comparable with reported satellite and ground‐based observations. Traveling ionospheric disturbances in the total electron contents from the simulation show good agreement with observations, including magnitude and propagating speed and evidence of conjugacy in the first 1–2 hr after eruption. Conjugacy inE × Bdrift, on the other hand, is more persistent.

 
more » « less
NSF-PAR ID:
10415582
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Volume:
50
Issue:
10
ISSN:
0094-8276
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    We report the first lidar observations of vertical fluxes of sensible heat and meteoric Na from 78 to 110 km in late May 2020 at McMurdo, Antarctica. The measurements include contributions from the complete temporal spectrum of gravity waves and demonstrate that wave‐induced vertical transport associated with atmospheric mixing by non‐breaking gravity waves, Stokes drift imparted by the wave spectrum, and perturbed chemistry of reactive species, can make significant contributions to constituent and heat transport in the mesosphere and lower thermosphere (MLT). The measured sensible heat and Na fluxes exhibit downward peaks at 84 km (−3.0 Kms−1and −5.5 × 104 cm−2s−1) that are ∼4 km lower than the flux peak altitudes observed at midlatitudes. This is likely caused by the strong downwelling over McMurdo in late May. The Na flux magnitude is double the maximum at midlatitudes, which we believe is related to strong persistent gravity waves in the MLT at McMurdo. To achieve good agreement between the measured Na flux and theory, it was necessary to infer that a large fraction of gravity wave energy was propagating downward, especially between 80 and 95 km where the Na flux and wave dissipation were largest. These downward propagating waves are likely secondary waves generated in‐situ by the dissipation of primary waves that originate from lower altitudes. The sensible heat flux transitions from downward below 90 km to upward from 97 to 106 km. The observations are explained with the fully compressible solutions for polarization relations of primary and secondary gravity waves withλz > 10 km.

     
    more » « less
  2. Abstract

    We perform ensemble simulations of radiation belt electron acceleration using the quasi‐linear approach during the storm on 9 October 2012, where chorus waves dominated electron acceleration atL = 5.2. Based on a superposed epoch analysis of 11 similar storms when both multi‐MeV electron flux enhancements and chorus wave activities were observed by Van Allen Probes, we use percentiles to sample the normalized input distributions for the four key inputs to estimate their relative perturbations. Using 11 points in each input parameter including chorus wave amplitudeBw, chorus wave peak frequencyfm, background magnetic fieldB0, and electron densityNe, we ran 114simulations to quantify the impact of uncertainties in the input parameters on the resulting simulated electron acceleration by chorus. By comparing the simulations to observations, our ensemble simulations reveal that inaccuracies in all four input parameters significantly affect the simulated electron acceleration, with the largest simulation errors attributed to the uncertainties inBw,Ne, andfm. The simulation can deviate from the observations by four orders of magnitude, while members with largest probability density (smallest perturbations in the input) provide reasonable estimations of output fluxes with log accuracy errors concentrated between ∼−2.0 and 0.5. Quantifying the uncertainties in our study is a prerequisite for the validation of our radiation belt electron model and improvements of accurate electron flux predictions.

     
    more » « less
  3. Abstract

    Quasi‐random vertical displacement fluctuations, caused by the spectrum of non‐breaking gravity waves, mix the atmosphere, similar to turbulence, which induces significant vertical transport of heat and constituents in the upper atmosphere. Multi‐decade observations of temperature, made between 85 and 100 km with a Na lidar at Colorado State University (CSU, 40.6°N, 105.1°W), are used to derive the seasonal variations of the wave‐induced thermal (KH) and constituent (KWave) diffusivities. Both show strong annual oscillations with maxima in winter, which increase with increasing altitude.KHandKWaveexhibit summer minima of ∼40 and ∼70 m2s−1, respectively, that are approximately constant with altitude. In winter,KHvaries from ∼50 at 85 to ∼180 m2s−1at 100 km, whileKWavevaries from ∼110 at 85 to ∼340 m2s−1at 100 km. These values are much larger than the eddy diffusivity (Kzz∼ 35 m2s−1) predicted for this site by the Whole Atmosphere Community Climate Model. The CSU diffusivities are comparable to similar measurements made at other mid‐latitude mountain sites in both hemispheres, and derived from global observations of atomic O. However, the seasonal variations differ from the O observations, which may reflect differences in wave sources at these sites and the different approaches employed to derive the wave diffusivities. Even so, the CSU results demonstrate that heat and constituent transport by unresolved, non‐breaking gravity waves are important processes that need to be incorporated in global chemistry models to properly characterize the thermal and constituent structure of the upper atmosphere.

     
    more » « less
  4. Abstract

    We present the analysis of ∼100 pc scale compact radio continuum sources detected in 63 local (ultra)luminous infrared galaxies (U/LIRGs;LIR≥ 1011L), using FWHM ≲ 0.″1–0.″2 resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8–170 pc, which are classified into four main categories—“AGN” (active galactic nuclei), “AGN/SBnuc” (AGN-starburst composite nucleus), “SBnuc” (starburst nucleus), and “SF” (star-forming clumps)—based on ancillary data sets and the literature. We find that “AGN” and “AGN/SBnuc” more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with “SBnuc” and “SF,” which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities (ΣSFR) are measured for “SF” and “SBnuc” using both the total 33 GHz continuum emission (SFR ∼ 0.14–13Myr−1, ΣSFR∼ 13–1600Myr−1kpc−2) and the thermal free–free emission from Hiiregions (median SFRth∼ 0.4Myr−1,ΣSFRth44Myr−1kpc−2). These values are 1–2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have a much flatter median 15–33 GHz spectral index (∼−0.08) compared with “SBnuc” and “SF” (∼−0.46), which may reflect higher nonthermal contribution from supernovae and/or interstellar medium densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100 pc scales.

     
    more » « less
  5. Abstract

    We report observations of the ground state transitions of12CO,13CO, C18O, HCN, and HCO+at 88–115 GHz in the inner region of the nearby galaxy IC 342. These data were obtained with the 16 pixel spectroscopic focal plane array Argus on the 100 m Robert C. Byrd Green Bank Telescope (GBT) at 6″–9″ resolution. In the nuclear bar region, the intensity distributions of12CO(1–0) and13CO(1–0) emission trace moderate densities, and differ from the dense gas distributions sampled in C18O(1–0), HCN(1–0), and HCO+(1–0). We observe a constant HCN(1–0)-to-HCO+(1–0) ratio of 1.2 ± 0.1 across the whole ∼1 kpc bar. This indicates that the HCN(1–0) and HCO+(1–0) lines have intermediate optical depth, and that the correspondingnH2of the gas producing the emission is of order 104.5−6cm−3. We show that HCO+(1–0) is thermalized and HCN(1–0) is close to thermalization. The very tight correlation between the HCN(1–0) and HCO+(1–0) intensities across the 1 kpc bar suggests that this ratio is more sensitive to the relative abundance of the two species than to the gas density. We confirm an angular offset (∼10″) between the spatial distribution of molecular gas and the star formation sites. Finally, we find a breakdown of theLIRLHCNcorrelation at high spatial resolution due to the effect of incomplete sampling of star-forming regions by HCN emission in IC 342. The scatter of theLIRLHCNrelation decreases as the spatial scale increases from 10″ to 30″ (170–510 pc), and is comparable to the scatter of the global relation at a scale of 340 pc.

     
    more » « less