Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research: Atmospheres”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources
USGS2023

VogCast: A framework for modeling volcanic air pollution and its application to the 2022 eruption of Mauna Loa Volcano, Hawai'i

Volcanic activity and the associated gas emissions into the atmosphere often result in adverse air quality conditions and present a hazard to human health and the environment. Building on a decade-long effort to provide operational surface sulfur dioxide and sulfate aerosol forecasts for the State of Hawai'i, we present an air quality modeling framework called VogCast. VogCast is designed to simplify ensemble air quality prediction on a regional scale by linking together multiple state-of-the-art models of meteorology, emissions, and dispersion. The framework is open-source and introduces a new dynamic plume-rise algorithm for distributing pollutants vertically. Using radar and satellite data, we demonstrate that VogCast reasonably captured the mean injection height, the location, and the general envelope of the vog plume during Mauna Loa's 2022 eruption. The results suggest that during the 12-day eruption period model performance varied between days with trade and non-trade wind conditions. Our findings also highlight the importance of sulfur dioxide emission rate and vent parameter inputs for improving forecast accuracy. The broad goal of this work is to better our understanding of vog dispersion and improve air quality prediction for impacted communities.

Hawaii↗
USGS2018

Investigating the origin of continual radio frequency impulses during explosive volcanic eruptions

Volcanic lightning studies have revealed that there is a relatively long‐lasting source of very high frequency radiation associated with the onset of explosive volcanic eruptions that is distinct from radiation produced by lightning. This very high frequency signal is referred to as “continual radio frequency (CRF)” due to its long‐lasting nature. The discharge mechanism producing this signal was previously hypothesized to be caused by numerous, small (10–100 m) leader‐forming discharges near the vent of the volcano. To test this hypothesis, a multiparametric data set of electrical and volcanic activity occurring during explosive eruptions of Sakurajima Volcano in Japan was collected from May to June 2015. Our observations show that a single CRF impulse has a duration on the order of 160 ns (giving an upper limit on discharge length of 10 m) and is distinct from near‐vent lightning discharges that were on the order of 30 m in length. CRF impulses did not produce discernible electric field changes and occurred in the absence of a net static electric field. Lightning mapping data and infrared video observations of the eruption column showed that CRF impulses originated from the gas thrust region of the column. These observations indicate that CRF impulses are not produced by small, leader‐forming discharges but rather are more similar to a streamer discharge, likely on the order of a few meters in length.

Sakurajima Volcano↗
USGS2011

The influence of the Atlantic Warm Pool on the Florida panhandle sea breeze

In this paper we examine the variations of the boreal summer season sea breeze circulation along the Florida panhandle coast from relatively high resolution (10 km) regional climate model integrations. The 23 year climatology (1979–2001) of the multidecadal dynamically downscaled simulations forced by the National Centers for Environmental Prediction–Department of Energy (NCEP‐DOE) Reanalysis II at the lateral boundaries verify quite well with the observed climatology. The variations at diurnal and interannual time scales are also well simulated with respect to the observations. We show from composite analyses made from these downscaled simulations that sea breezes in northwestern Florida are associated with changes in the size of the Atlantic Warm Pool (AWP) on interannual time scales. In large AWP years when the North Atlantic Subtropical High becomes weaker and moves further eastward relative to the small AWP years, a large part of the southeast U.S. including Florida comes under the influence of relatively strong anomalous low‐level northerly flow and large‐scale subsidence consistent with the theory of the Sverdrup balance. This tends to suppress the diurnal convection over the Florida panhandle coast in large AWP years. This study is also an illustration of the benefit of dynamic downscaling in understanding the low‐frequency variations of the sea breeze.

Florida↗
USGS2019

Laboratory experiments of volcanic ash resuspension by wind

Fresh volcanic eruption deposits tend to be loose, bare, and readily resuspended by wind. Major resuspension events in Patagonia, Iceland, and Alaska have lofted ash clouds with potential to impact aircraft, infrastructure, and downwind communities. However, poor constraints on this resuspension process limit our ability to model this phenomenon. Here, we present laboratory experiments measuring threshold shear velocities and emission rates of resuspended ash under different environmental conditions, including relative humidity of 25–75% and simulated rainfall with subsequent drying. Eruption deposits were replicated using ash collected from two major eruptions: the 18 May 1980 eruption of Mount St. Helens and the 1912 eruption of Novarupta, in Alaska's Valley of Ten Thousand Smokes. Samples were conditioned in a laboratory chamber and prepared with bulk deposit densities of 1,300–1,500 kg/m 3 . A control sample of dune sand was included for comparison. The deposits were subjected to different wind speeds using a modified PI‐SWERL® instrument. Under a constant relative humidity of 50% and shear velocities 0.4–0.8 m/s, PM 10 emission by resuspension ranged from 10 to >100 mg·m −2 ·s −1 . Addition of liquid water equivalent to 5 mm of rainfall had little lasting effect on Mount St. Helens wind erosion potential, while the Valley of Ten Thousand Smokes deposits exhibited lower emissions for at least 12 days. The results indicate that particle resuspension due to wind erosion from ash deposits potentially exceeds that of most desert surfaces and approaches some of the highest emissions ever measured.

Journal of Geophysical Research - Atmospheres↗
USGS2023

Flow development and entrainment in turbulent particle-laden jets

Explosive eruptions expel volcanic gases and particles at high pressures and velocities. Within this multiphase fluid, small ash particles affect the flow dynamics, impacting mixing, entrainment, turbulence, and aggregation. To examine the role of turbulent particle behavior, we conducted an analogue experiment using a particle-laden jet. We used compressed air as the carrier fluid, considering turbulent conditions at Reynolds numbers from approximately 5,000 to 20,000. Two different particles were examined: 14-μm diameter solid nickel spheres and 13-μm diameter hollow glass spheres. These resulted in Stokes numbers between 1 and 35 based on the convective scale. The particle mass percentage in the mixture is varied from 0.3% to more than 20%. Based on a 1-D volcanic plume model, these Stokes numbers and mass loadings corresponded to millimeter-scale particle diameters at heights of 4–8 km above the vent during large, sustained eruptions. Through particle image velocimetry, we measured the mean flow behavior and the turbulence statistics in the near-exit region, primarily focusing on the dispersed phase. We show that the flow behavior is dominated by the particle inertia, with high Stokes numbers reducing the entrainment by more than 40%. When applied to volcanic plumes, these results suggest that high-density particles can greatly increase the probability of column collapse.

Journal of Geophysical Research: Atmospheres↗
USGS1999

Carbon cycling in boreal wetlands: A comparison of three approaches

Three independent methods were used to measure net ecosystem production (NEP) in four wetlands near Thompson, Manitoba, Canada. The first method calculated NEP by subtracting heterotrophic respiration from net primary productivity, using both measurements and estimates derived from the literature. The second method used radiocarbon data from cores to derive long-term NEP averaged over the past several decades. The third method used direct measurement of NEP combined with a model to fill in for days with no data. The three methods, with their independently derived uncertainties, all show the same magnitude and pattern of NEP variation across four different wetland types. However, direct measurement yielded distinctly lower estimates of NEP in the most productive sites. Highest NEP (31–180 gC m −2 yr −1 ) was observed in the two wetlands with the highest proportion of sedge vegetation. A bog collapse scar and a nutrient-rich fen had NEP values not significantly different from zero. The maximum NEP at sites with intermediate nutrient status is due to slower overall decomposition and is likely associated with greater allocation of production below ground by sedges. The three methods for estimating NEP differ in the effort required, the sources of error, and in the timescale over which they apply. Used in combination, they allow estimation of parameters such as below- ground production and the contribution of heterotrophic decomposition to total soil respiration. Using the radiocarbon method, we also derived estimates of the rate of N accumulation in the four wetland types.

Manitoba↗
USGS2008

Reconstructed historical land cover and biophysical parameters for studies of land-atmosphere interactions within the eastern United States

Over the past 350 years, the eastern half of the United States experienced extensive land cover changes. These began with land clearing in the 1600s, continued with widespread deforestation, wetland drainage, and intensive land use by 1920, and then evolved to the present-day landscape of forest regrowth, intensive agriculture, urban expansion, and landscape fragmentation. Such changes alter biophysical properties that are key determinants of land-atmosphere interactions (water, energy, and carbon exchanges). To understand the potential implications of these land use transformations, we developed and analyzed 20-km land cover and biophysical parameter data sets for the eastern United States at 1650, 1850, 1920, and 1992 time slices. Our approach combined potential vegetation, county-level census data, soils data, resource statistics, a Landsat-derived land cover classification, and published historical information on land cover and land use. We reconstructed land use intensity maps for each time slice and characterized the land cover condition. We combined these land use data with a mutually consistent set of biophysical parameter classes, to characterize the historical diversity and distribution of land surface properties. Time series maps of land surface albedo, leaf area index, a deciduousness index, canopy height, surface roughness, and potential saturated soils in 1650, 1850, 1920, and 1992 illustrate the profound effects of land use change on biophysical properties of the land surface. Although much of the eastern forest has returned, the average biophysical parameters for recent landscapes remain markedly different from those of earlier periods. Understanding the consequences of these historical changes will require land-atmosphere interactions modeling experiments.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m-2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

Journal of Geophysical Research D: Atmospheres↗
USGS2008

Northerly surface winds over the eastern North Pacific Ocean in spring and summer

Persistent spring and summer northerly surface winds are the defining climatological feature of the western coast of North America, especially south of the Oregon coast. Northerly surface winds are important for upwelling and a vast array of other biological, oceanic, and atmospheric processes. Intermittence in northerly coastal surface wind is characterized and wind events are quantitatively defined using coastal buoy data south of Cape Mendocino on the northern California coast. The defined wind events are then used as a basis for composites in order to explain the spatial evolution of various atmospheric and oceanic processes. Wind events involve large-scale changes in the three-dimensional atmospheric circulation including the eastern North Pacific subtropical anticyclone and southeast trade winds. Composites of QSCAT satellite scatterometer wind estimates from 1999 to 2005 based on a single coastal buoy indicate that wind events typically last 72-96 h and result in anomalies in surface wind and Ekman pumping that extend over 1000 kin from the west coast of North America. It may be useful to consider ocean circulation and dependent ecosystem dynamics and the distribution of temperature, moisture, and aerosols in the atmospheric boundary layer in the context of wind events defined herein. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2008

Diurnal and vertical variability of the sensible heat and carbon dioxide budgets in the atmospheric surface layer

The diurnal and vertical variability of heat and carbon dioxide (CO2) in the atmospheric surface layer are studied by analyzing measurements from a 213 in tower in Cabauw (Netherlands). Observations of thermodynamic variables and CO2 mixing ratio as well as vertical profiles of the turbulent fluxes are used to retrieve the contribution of the budget terms in the scalar conservation equation. On the basis of the daytime evolution of turbulent fluxes, we calculate the budget terms by assuming that turbulent fluxes follow a linear profile with height. This assumption is carefully tested and the deviation ftom linearity is quantified. The budget calculation allows us to assess the importance of advection of heat and CO2 during day hours for three selected days. It is found that, under nonadvective conditions, the diurnal variability of temperature and CO2 is well reproduced from the flux divergence measurements. Consequently, the vertical transport due to the turbulent flux plays a major role in the daytime evolution of both scalars and the advection is a relatively small contribution. During the analyzed days with a strong contribution of advection of either heat or carbon dioxide, the flux divergence is still an important contribution to the budget. For heat, the quantification of the advection contribution is in close agreement with results from a numerical model. For carbon dioxide, we qualitatively corroborate the results with a Lagrangian transport model. Our estimation of advection is compared with, traditional estimations based on the Net Ecosystem-atmosphere Exchange (NEE). Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Oxygen-18 concentrations in recent precipitation and ice cores on the Tibetan Plateau

A detailed study of the climatic significance of δ 18 O in precipitation was completed on a 1500 km southwest-northeast transect of the Tibetan Plateau in central Asia. Precipitation samples were collected at four meteorological stations for up to 9 years. This study shows that the gradual impact of monsoon precipitation affects the spatial variation of δ 18 O-T relationship along the transect. Strong monsoon activity in the southern Tibetan Plateau results in high precipitation rates and more depleted heavy isotopes. This depletion mechanism is described as a precipitation “amount effect” and results in a poor δ 18 O-T relationship at both seasonal and annual scales. In the middle of the Tibetan Plateau, the effects of the monsoon are diminished but continue to cause a reduced correlation of δ 18 O and temperature at the annual scale. At the monthly scale, however, a significant δ 18 O-T relationship does exist. To the north of the Tibetan Plateau beyond the extent of the effects of monsoon precipitation, δ 18 O in precipitation shows a strong temperature dependence. δ 18 O records from two shallow ice cores and historic air temperature data were compared to verify the modern δ 18 O-T relationship. δ 18 O in Dunde ice core was positively correlated with air temperature from a nearby meteorological station in the north of the plateau. The δ 18 O variation in an ice core from the southern Plateau, however, was inversely correlated with precipitation amount at a nearby meteorological station and also the accumulation record in the ice core. The long-term variation of δ 18 O in the ice core record in the monsoon regions of the southern Tibetan Plateau suggest past monsoon seasons were probably more expansive. It is still unclear, however, how changes in large-scale atmosphere circulation might influence summer monsoon precipitation on the Tibetan Plateau.

Journal of Geophysical Research D: Atmospheres↗
USGS2003

Multiple oxygen and sulfur isotopic analyses on water-soluble sulfate in bulk atmospheric deposition from the southwestern United States

Sulfate is a major component of bulk atmospheric deposition (including dust, aerosol, fog, and rain). We analyzed sulfur and oxygen isotopic compositions of water-soluble sulfate from 40 sites where year-round dust traps collect bulk atmospheric deposition in the southwestern United States. Average sulfur and oxygen isotopic compositions (δ 34 S and δ 18 O) are 5.8 ± 1.4 (CDT) and 11.2 ± 1.9 (SMOW) (n = 47), respectively. Samples have an oxygen 17 anomaly (Δ 17 O), with an average value of 1.0 ± 0.6‰. Except for a weak positive correlation between δ 18 O and Δ 17 O values (r 2 ≈ 0.4), no correlation exists for δ 18 O versus δ 34 S, Δ 17 O versus δ 34 S, or any of the three isotopic compositions versus elevation of the sample site. Exceptional positive Δ 17 O values (up to 4.23‰) are found in samples from sites in the vicinity of large cities or major highways, and near-zero Δ 17 O values are found in samples close to dry lakes. Comparison of isotopic values of dust trap sulfate and desert varnish sulfate from the region reveals that varnish sulfate has average isotopic values that are ∼4.8‰ lower for δ 18 O, ∼2.1‰ higher for δ 34 S, and ∼0.3‰ lower for Δ 17 O than those of the present-day bulk deposition sulfate. Although other factors could cause the disparity, this observation suggests a possibility that varnish sulfate may have recorded a long-term atmospheric sulfate deposition during the Holocene or Pleistocene, as well as the differences between sulfur and oxygen isotopic compositions of the preindustrial bulk deposition sulfate and those of the industrial era.

Journal of Geophysical Research D: Atmospheres↗
USGS1997

Winter CO2 fluxes in a boreal forest

We measured soil respiration during two winters in three different ecotypes of the BOREAS northern study area. The production of CO 2 was continuous throughout the winter and, when totaled for the winter of 1994–1995, was equivalent to the release of ∼40–55 g C/m 2 from the soil surface. As soils cooled in the early winter, the CO 2 production rate decreased in a manner that appeared to be exponentially related to shallow soil temperatures. This exponential relationship was not observed when soils began to warm, possibly indicating that there may be additional or different processes responsible for increased CO 2 production during winter warming events. We also measured CO 2 concentrations in soil gas and the Δ 14 C of the soil CO 2 . These measurements show that the CO 2 produced in winter is not simply the return to the atmosphere of the carbon fixed during the previous growing season. We suggest that the wintertime production of CO 2 originates, at least in part, from the decomposition of old organic carbon stored at depth in the soil.

Journal of Geophysical Research D: Atmospheres↗
USGS2012

Improving the accuracy of S0 2 column densities and emission rates obtained from upward-looking UV-spectroscopic measurements of volcanic plumes by taking realistic radiative transfer into account

Sulfur dioxide (SO 2 ) is monitored using ultraviolet (UV) absorption spectroscopy at numerous volcanoes around the world due to its importance as a measure of volcanic activity and a tracer for other gaseous species. Recent studies have shown that failure to take realistic radiative transfer into account during the spectral retrieval of the collected data often leads to large errors in the calculated emission rates. Here, the framework for a new evaluation method which couples a radiative transfer model to the spectral retrieval is described. In it, absorption spectra are simulated, and atmospheric parameters are iteratively updated in the model until a best match to the measurement data is achieved. The evaluation algorithm is applied to two example Differential Optical Absorption Spectroscopy (DOAS) measurements conducted at Kilauea volcano (Hawaii). The resulting emission rates were 20 and 90% higher than those obtained with a conventional DOAS retrieval performed between 305 and 315 nm, respectively, depending on the different SO 2 and aerosol loads present in the volcanic plume. The internal consistency of the method was validated by measuring and modeling SO2 absorption features in a separate wavelength region around 375 nm and comparing the results. Although additional information about the measurement geometry and atmospheric conditions is needed in addition to the acquired spectral data, this method for the first time provides a means of taking realistic three-dimensional radiative transfer into account when analyzing UV-spectral absorption measurements of volcanic SO 2 plumes.

Journal of Geophysical Research D: Atmospheres↗
USGS2016

Influence of land-atmosphere feedbacks on temperature and precipitation extremes in the GLACE-CMIP5 ensemble

We examine how soil moisture variability and trends affect the simulation of temperature and precipitation extremes in six global climate models using the experimental protocol of the Global Land-Atmosphere Coupling Experiment of the Coupled Model Intercomparison Project, Phase 5 (GLACE-CMIP5). This protocol enables separate examinations of the influences of soil moisture variability and trends on the intensity, frequency, and duration of climate extremes by the end of the 21st century under a business-as-usual (Representative Concentration Pathway 8.5) emission scenario. Removing soil moisture variability significantly reduces temperature extremes over most continental surfaces, while wet precipitation extremes are enhanced in the tropics. Projected drying trends in soil moisture lead to increases in intensity, frequency, and duration of temperature extremes by the end of the 21st century. Wet precipitation extremes are decreased in the tropics with soil moisture trends in the simulations, while dry extremes are enhanced in some regions, in particular the Mediterranean and Australia. However, the ensemble results mask considerable differences in the soil moisture trends simulated by the six climate models. We find that the large differences between the models in soil moisture trends, which are related to an unknown combination of differences in atmospheric forcing (precipitation, net radiation), flux partitioning at the land surface, and how soil moisture is parameterized, imply considerable uncertainty in future changes in climate extremes.

Journal of Geophysical Research D: Atmospheres↗
USGS2019

Assessing the seasonal dynamics of nitrate and sulfate aerosols at the South Pole utilizing stable isotopes

Atmospheric nitrate (NO 3 − = particulate NO 3 − + gas‐phase nitric acid [HNO 3 ]) and sulfate (SO 4 2− ) are key molecules that play important roles in numerous atmospheric processes. Here, the seasonal cycles of NO 3 − and total suspended particulate sulfate (SO 4 2− (TSP) ) were evaluated at the South Pole from aerosol samples collected weekly for approximately 10 months (26 January to 25 October) in 2002 and analyzed for their concentration and isotopic compositions. Aerosol NO 3 − was largely affected by snowpack emissions in which [NO 3 − ] and δ 15 N(NO 3 − ) were highest (49.3 ± 21.4 ng/m 3 , n = 8) and lowest (−47.0 ± 11.7‰, n = 5), respectively, during periods of sunlight in the interior of Antarctica. The seasonal cycle of Δ 17 O(NO 3 − ) reflected tropospheric chemistry year‐round with lower values observed during sunlight periods and higher values observed during dark periods, reflecting shifts from HO x ‐ to O 3 ‐dominated oxidation chemistry. SO 4 2− (TSP) concentrations were highest during austral summer and fall (86.7 ± 73.7 ng/m 3 , n = 18) and are indicated to be derived from dimethyl sulfide (DMS) emissions, as δ 34 S(SO 4 2− ) (TSP) values (18.5 ± 1.0‰, n = 10) were similar to literature δ 34 S(DMS) values. The seasonal cycle of Δ 17 O(SO 4 2− ) (TSP) exhibited minima during austral summer (0.9 ± 0.1‰, n = 5) and maxima during austral fall (1.3 ± 0.3‰, n = 6) and austral spring (1.6 ± 0.1‰, n = 5), indicating a shift from HO x ‐ to O 3 ‐dominated chemistry in the atmospheric derived SO 4 2− component. Overall, the budgets of NO 3 − and SO 4 2− (TSP) at the South Pole were complex functions of transport, localized chemistry, biological activity, and meteorological conditions, and these results will be important for interpretations of oxyanions in ice core records in the interior of Antarctica.

Journal of Geophysical Research D: Atmospheres↗
USGS2002

Fire dynamics and implications for nitrogen cycling in boreal forests

We used a dynamic, long-term mass balance approach to track cumulative carbon (C) and nitrogen (N) losses to fire in boreal Manitoba over the 6500 years since deglaciation. Estimated C losses to decomposition and fire, combined with measurements of N pools in mature and burned forest floors, suggest that loss of N by combustion has likely resulted in a long-term loss that exceeds the amount of N stored in soil today by 2 to 3 times. These estimates imply that biological N fixation rates could be as high as 5 to 10 times atmospheric deposition rates in boreal regions. At the site scale, the amount of N lost is due to N content of fuels, which varies by stand type and fire severity, which in turn vary with climate and fire dynamics. The interplay of fire frequency, fire severity, and N partitioning during regrowth are important for understanding rates and sustainability of nutrient and carbon cycling over millenia and over broad regions.

Manitoba↗
USGS2003

GCIP water and energy budget synthesis (WEBS)

As part of the World Climate Research Program's (WCRPs) Global Energy and Water-Cycle Experiment (GEWEX) Continental-scale International Project (GCIP), a preliminary water and energy budget synthesis (WEBS) was developed for the period 1996-1999 fromthe "best available" observations and models. Besides this summary paper, a companion CD-ROM with more extensive discussion, figures, tables, and raw data is available to the interested researcher from the GEWEX project office, the GAPP project office, or the first author. An updated online version of the CD-ROM is also available at http://ecpc.ucsd.edu/gcip/webs.htm/. Observations cannot adequately characterize or "close" budgets since too many fundamental processes are missing. Models that properly represent the many complicated atmospheric and near-surface interactions are also required. This preliminary synthesis therefore included a representative global general circulation model, regional climate model, and a macroscale hydrologic model as well as a global reanalysis and a regional analysis. By the qualitative agreement among the models and available observations, it did appear that we now qualitatively understand water and energy budgets of the Mississippi River Basin. However, there is still much quantitative uncertainty. In that regard, there did appear to be a clear advantage to using a regional analysis over a global analysis or a regional simulation over a global simulation to describe the Mississippi River Basin water and energy budgets. There also appeared to be some advantage to using a macroscale hydrologic model for at least the surface water budgets. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres↗