Elevated ozone (O3) pollution in the warm season is an emerging environmental concern affecting global highly urbanized megacities. In southwestern China, full characterization of causes for O3 pollution has been stymied by limited observations and the dominant factors that influence O3 variability on a long-term basis still lack understanding. Herein, we identified O3 variations and inferred trends in precursor emissions in Chengdu over 2013–2020 based on extensive ambient measurements, emission inventory, and satellite data. Numerical models were used to investigate the changes in meteorological variability and biogenic emissions. Trends of O3 in urban areas show deterioration (+14.0% yr−1) between 2013 and 2016 followed by a slight decrease over 2017–2020, while O3 levels in rural areas generally show a downward trend (−2.9% yr−1) during 2014–2020. Both emission inventory (−3.7% yr−1) and OMI satellite columns (−4.5% yr−1) depict strong decline trends in NOx emissions, while satellite HCHO columns exhibit a flattened downward trend of VOC emissions (−1.8% yr−1), which caused rural areas shifted from VOCs-limited to transitional or NOx-limited regime since 2016. Considering metropolitan Chengdu remains VOCs-limited regime over time, the existing regulatory framework involving simultaneous NOx and VOCs control would result in evident O3 improvements in the near future. Despite benefits from anthropogenic emission reductions, we demonstrate that meteorological conditions and enhanced biogenic emissions over the warm season could partially or even fully offset effects attributed to emission changes, making the net effects obscure. This finding provides robust evidence of reductions in NOx and VOCs emission and informs effective O3 mitigation policies for megacities which undergo similar emission pathways in Chengdu.

Extreme solar particle events (ESPEs) are rare and the most potent known processes of solar eruptive activity. During ESPEs, a vast amount of cosmogenic isotopes (CIs) 10Be, 36Cl and 14C can be produced in the Earth’s atmosphere. Accordingly, CI measurements in natural archives allow us to evaluate particle fluxes during ESPEs. In this work, we present a new method of ESPE fluence (integral flux) reconstruction based on state-of-the-art modeling advances, allowing to fit together different CI data within one model. We represent the ESPE fluence as an ensemble of scaled fluence reconstructions for ground-level enhancement (GLE) events registered by the neutron monitor network since 1956 coupled with satellite and ionospheric measurements data. Reconstructed ESPE fluences appear softer in its spectral shape than earlier estimates, leading to significantly higher estimates of the low-energy (E<100 MeV) fluence. This makes ESPEs even more dangerous for modern technological systems than previously believed. Reconstructed ESPE fluences are fitted with a modified Band function, which eases the use of obtained results in different applications.

This study investigates the representation of stratocumulus (Sc) clouds, cloud variability, and precipitation statistics over the Southern Ocean (SO) to understand the dominant ice processes within the Icosahedral Nonhydrostatic (ICON) model at the kilometer scale using real case simulations. The simulations are evaluated using the shipborne observations as open-cell stratocumuli were continuously observed during two days (26th -27th of March 2016), south of Tasmania. The radar retrievals are used to effectively analyze the forward- simulated radar signatures from Passive and Active Microwave TRAnsfer (PAMTRA). We contrast cloud-precipitation statistics, and microphysical process rates between simulations performed with one-moment (1M) and two-moment (2M) microphysics schemes. We further analyze their sensitivity to primary and secondary ice-phase processes (Hallett–Mossop and collisional breakup). Both processes have previously been shown to improve the ice properties of simulated shallow mixed-phase clouds over the SO in other models. We find that only simulations with continuous formation, growth, and subsequent melting of graupel, and the effective riming of in-cloud rain by graupel, capture the observed cloud-precipitation vertical structure. In particular, the 2M microphysics scheme requires additional tuning for graupel processes in SO stratocumuli. Lowering the assumed graupel density and terminal velocity, in combination with secondary ice processes, enhances graupel formation in 2M microphysics ICON simulations. Overall, all simulations capture the observed intermittency of precipitation irrespective of the microphysics scheme used, and most of them sparsely distribute intense precipitation (>1mm h-1 ) events. Furthermore, the simulated clouds are too reflective as they are optically thick and/or have high cloud cover.

As the largest natural source of sulfur-containing gases into the atmosphere, ocean organism-derived dimethyl sulfide (DMS) has been considered to play a critical role in the Earth’s climate system. Yet there are great uncertainties in modeling the spatiotemporal variations of DMS and incomplete knowledge of influencing factors in different oceanic regions. Moreover, little is known about the future change of global DMS, which limits our understanding of the feedback of marine ecosystem to climate change. Here we develop an artificial neural network model and combine data mining approaches to address these issues. Phytoplankton biomass and salinity are currently predominant factors associated with DMS variability in the coastal and Arctic regions, respectively. In the mid- and low-latitude open oceans, nutrients and temperature are also crucial factors in addition to radiation and mixed layer depth, and their relationships with DMS show reversals when passing certain thresholds. Although the global average DMS concentration and emission slightly decline from 2005 to 2100, they may change considerably in specific regions. In contrast to the DMS decreases in the low-latitudes mainly related with phosphate reduction and temperature rise and in the North Atlantic subpolar gyre attributed to salinity decline, warming will cause DMS increase in the Southern Ocean and sea ice loss will dramatically enhance DMS emission in the Arctic. Although the global negative feedback loop between oceanic DMS and climate may not operate, the future spatial redistribution of DMS may lead to the change in cloud cover pattern and significantly affect regional climate.

We simulated the Nov 4, 2021 geomagnetic storm event penetrating electric field using the Multiscale Atmosphere-Geospace Environment (MAGE) and compared with the NASA ICON observation. The ICON observation showed enhancement of the vertical ion drift when the penetrating electric field arrived at the equatorial region. The simulated vertical ion drifts are consistent with ICON observation. Hence, we are able to verify the MAGE simulation with ICON observation. On the dusk side, the MAGE simulation showed strong pre-reversal enhancement (PRE), whereas the ICON observation did not display any sign of the PRE. The MAGE simulation did show that PRE amplitude decreases as altitude increase. Because the ICON orbital height is above the model upper boundary, it could be a factor for the discrepancy. Instrumental issue cannot be ruled out at this moment. GOLD UV image at the same time exhibits multiple plasma bubbles, which seem to suggest the existence of the PRE.

Searches for phosphine in Venus’ atmosphere have sparked a debate. Cordiner et al. 2022 analyse spectra from the Stratospheric Observatory For Infrared Astronomy (SOFIA) and infer <0.8 ppb of PH3. We noticed that spectral artefacts arose mainly from inessential calibration-load signals. By-passing these signals allows simpler post-processing, and 6.5σ detection of 1 ppb of PH3 at ~75 km altitude (just above the clouds). Compiling six phosphine results would suggest the abundance inverts: decreasing above the clouds but rising again in the mesosphere from some unexplained source. However, no such extra source is needed if phosphine is undergoing destruction by sunlight (photolysis), as it does on Earth. Low values/limits were found where the viewed part of the super-rotating Venusian atmosphere had passed through sunlight, while the high values are from views moving into sunlight. We suggest Venusian phosphine is indeed present, and so merits further work on models of its origins.

Dynamic influences on summertime seasonal United States rainfall variability are not well understood. A major cause of moisture transport is the Great Plains low-level jet (LLJ). Using observations and a dry atmospheric general circulation model, this study explored the distinct and combined impacts of two prominent atmospheric teleconnections - the East Asian monsoon (EAM) and North Atlantic subtropical high (NASH) - on the Great Plains LLJ in the summer. Separately, a strong EAM and strong western NASH are linked to a strengthened LLJ and positive rainfall anomalies in the Plains/ Midwest. Overall, NASH variability is more important for considering the LLJ impacts, but strong EAM events amplify western NASH-related Great Plains LLJ strengthening and associated rainfall signals. This occurs when the EAM-forced Rossby wave pattern over North America constructively interferes with low-level wind field, providing upper-level support for the LLJ and increasing mid- to upper-level divergence.

Within the Charlotte, North Carolina, to Atlanta, Georgia, megaregion (Charlanta), the Atlanta metropolitan area has been shown to augment proximal cloud-to-ground (CG) lightning occurrence. Although numerous studies have documented this “urban lightning effect” (ULE) with regard to CG lightning, relatively few have investigated urban effects on distributions of total lightning (TL). Moreover, there has yet to be a study of the ULE using TL observations from the Geostationary Lightning Mapper (GLM). In an effort to fill this gap, we investigated spatial distributions of TL around the cities of Atlanta, GA, Greenville, SC, and Charlotte, NC, using GLM data collected during the warm seasons of 2018–2021. Analyses reveal augmentation of total lightning intensity and frequency over the major cities of Atlanta and Charlotte, with a diminished urban signal over the smaller city of Greenville. This work also demonstrated the potential efficacy of the emerging satellite-based TL climatology in ULE studies.

The Hunga Tonga-Hunga Ha’apai (HTHH) volcanic eruption in January 2022 injected extreme amounts of water vapor (H2O) and a moderate amount of the aerosol precursor (SO2) into the Southern Hemisphere (SH) stratosphere. The H2O and aerosol perturbations have persisted and resulted in large-scale SH stratospheric cooling, equatorward shift of the Antarctic polar vortex, and slowing of the Brewer-Dobson circulation associated with a substantial ozone reduction in the SH winter midlatitudes. Chemistry-climate model simulations forced by realistic HTHH inputs of H2O and SO2 reproduce the observed stratospheric cooling and circulation effects, demonstrating the observed behavior is due to the volcanic influences. Furthermore, the combination of aerosol transport to polar latitudes and a cold polar vortex enhances springtime Antarctic ozone loss, consistent with observed polar ozone behavior in 2022.

Natural aerosols and their interactions with clouds remain an important uncertainty within climate models, especially at the poles. Here, we study the behavior of sea salt aerosols (SSaer) in the Arctic and Antarctic within 12 climate models from CMIP6. We investigate the driving factors that control SSaer abundances and show large differences based on the choice of the source function, and the representation of aerosol processes in the atmosphere. Close to the poles, the CMIP6 models do not match observed seasonal cycles of surface concentrations, likely due to the absence of wintertime SSaer sources such as blowing snow. Further away from the poles, simulated concentrations have the correct seasonality, but have a positive mean bias of up to one order of magnitude. SSaer optical depth is derived from the MODIS data and compared to modeled values, revealing good agreement, except for winter months. Better agreement for AOD than surface concentration may indicate a need for improving the vertical distribution, the size distribution and/or hygroscopicity of modeled polar SSaer. Source functions used in CMIP6 emit very different numbers of small SSaer, potentially exacerbating cloud-aerosol interaction uncertainties in these remote regions. For future climate scenarios SSP126 and SSP585, we show that SSaer concentrations increase at both poles at the end of the 21st century, with more than two times mid-20th century values in the Arctic. The pre-industrial climate CMIP6 experiments suggest there is a large uncertainty in the polar radiative budget due to SSaer.

In a context of climate change, the stakes surrounding water availability are getting higher. Decomposing and quantifying the effects of climate on discharge allows to better understand their impact on water resources. We propose a methodology to separate the effect of change in annual mean of climate variables from the effect of intra-annual distribution of precipitations. It combines the Budyko framework with outputs from a Land Surface Model (LSM). The LSM is used to reproduces the behavior of 2134 reconstructed watersheds over Europe between 1902 and 2010, with climate inputs as the only source of change. We fit to the LSM outputs a one parameter approximation to the Budyko framework. It accounts for the evolution of annual mean in precipitation (P) and potential evapotranspiration (PET). We introduce a time-varying parameter in the equation which represents the effect of long-term variations in the intra-annual distribution of P and PET. To better assess the effects of changes in annual means or in intra-annual distribution of P, we construct synthetic forcings fixing one or the other. The results over Europe show that the changes in discharge due to climate are dominated by the trends in the annual averages of P. The second main climate driver is PET, except over the Mediterranean area where changes in intra-annual variations of P have a higher impact on discharge than trends in PET. Therefore the effects of changes in intra-annual distribution of climate variables are not to be neglected when looking at changes in annual discharge.

We present an analysis of the magnetospheric response to the two Coronal Mass Ejection (CME) impacts which led to the destruction of 38 out of 49 Starlink satellites in early February 2022. We employ the Expanding-Contracting Polar Cap model to analyse the variation in the size of the ionospheric polar caps and thereby quantify the Universal Time (UT) effect of the diurnal motions of the geomagnetic poles in a geocentric frame of reference. The results show that use of quasi-steady convection model predicts a very similar global power deposition into the thermosphere as that inferred here, but does not give the same division of that power between the northern and southern hemispheres. We demonstrate that, through the combined effects of the Russell-McPherron dipole-tilt mechanism on solar-wind magnetosphere coupling and of the diurnal polar cap motions in a geocentric frame, the power deposited varies significantly with the arrival UT of the CMEs at Earth. We show that in the events of early February 2022, both CMEs arrived at almost the optimum UT to cause maximum thermospheric heating.

Previous studies suggest that boreal summer intraseasonal variations along the subtropical westerly jet (SJ), featuring quasi-biweekly periodicity, frequently modulate downstream subseasonal variations over East Asia (EA). Based on subseasonal hindcasts from six dynamical models, this study discovered that the leading two-three-week prediction skills for surface air temperature (SAT) are improved significantly in summer when the SJ has strengthened intraseasonal signals, which are best demonstrated over the eastern Tibetan Plateau, Southwest Basin, and North China. The reasons are that the enhanced quasi-biweekly wave and the associated energy dispersion along the SJ cause more regular quasi-biweekly periodic variations of downstream SAT, which potentially increase regional predictability. This study suggests not only that intraseasonal variations along the SJ could provide a window of opportunity for achieving better subseasonal prediction over EA, but also that intraseasonal waves along the SJ are crucial for improving EA subseasonal prediction.

Predicting secondary organic aerosol (SOA) formation relies either on extremely detailed, numerically expensive models accounting for the condensation of individual species or on extremely simplified, numerically affordable models parameterizing SOA formation for large-scale simulations. In this work, we explore the possibility of creating a random forest to reproduce the behavior of a detailed atmospheric organic chemistry model at a fraction of the numerical cost. A comprehensive dataset was created based on thousands of individual detailed simulations, randomly initialized to account for the variety of atmospheric chemical environments. Recurrent random forests were trained to predict organic matter formation from dodecane and toluene precursors, and the partitioning between gas and particle phases. Validation tests show that the random forests perform well without any divergence over 10 days of simulations. The distribution of errors shows that the sampling of initial conditions for the training simulations needs to focus on chemical regimes where SOA production is the most sensitive. Sensitivity tests show that specializing multiple random forests for a specific chemical regime is not more efficient than training a single general random forest for the entire dataset. The most important predictors are those providing information about the chemical regime, oxidants levels and existing organic mass. The choice of predictors is crucial as using too many unimportant predictors reduces the performances of the random forests.

Changes in atmospheric iron (Fe) deposition to the open ocean affect net primary productivity, nitrogen fixation, and carbon uptake rates. We investigate the changes in soluble Fe (SFe) deposition from the pre-industrial period to the late 21st century using the EC-Earth3-Iron Earth System model, which stands out for its comprehensive representation of the atmospheric oxalate, sulfate, and Fe cycles. We show how anthropogenic activity has modified the magnitude and spatial distribution of SFe deposition by increasing combustion Fe emissions along with atmospheric acidity and oxalate levels. We find that SFe deposition has doubled since the early Industrial Era using the Coupled Model Intercomparison Project Phase 6 (CMIP6) emission inventory, with acidity being the main solubilization pathway for dust Fe, and ligand-promoted (oxalate) processing dominating the solubilization of combustion Fe. We project a global SFe deposition increase of 40% by the late 21st century relative to present day under Shared Socioeconomic Pathway (SSP) 3-7.0, which assumes weak climate change mitigation policies. In contrast, sustainable and business-as-usual SSPs (1-2.6 and 2-4.5) result in 35% and 10% global decreases, respectively. Despite these differences, SFe deposition consistently increases and decreases across SSPs over the (high nutrient low chlorophyl) equatorial Pacific and Southern Ocean (SO), respectively. Future changes in dust and wildfires with climate remains a key challenge for constraining SFe projections. We show that the equatorial Pacific and the SO would be sensitive not only to changes in Australian or South American dust emissions, but also to those in North Africa.

Convective cold pools (CPs) are known to mediate the interaction between convective rain cells and thereby help organize thunderstorm clusters, in particular mesoscale convective systems and extreme rainfall events. Unfortunately, the observational detection of CPs on a large scale has so far been hampered by the lack of relevant large-scale nearsurface data. Unlike numerical studies, where high-resolution near-surface fields of relevant quantities such as virtual temperature and winds are available and frequently used to detect cold pools, observational studies mainly identify CPs based on surface time series. Since research vessels or weather stations measure these time series locally, the characterization of cold pools from observations is limited to regional or station-based studies. To eventually enable studies on a global scale, we here develop and evaluate a methodology for the detection of CPs that relies only on data that (i) is globally available and (ii) has high spatio-temporal resolution. We trained convolutional neural networks to segment CPs in cloud and rainfall fields from high-resolution cloud resolving simulation output. Such data is not only available from simulations, but also from geostationary satellites that fulfill both (i) and (ii). The networks make use of a U-Net architecture, a common choice for image segmentation due to its strength in learning spatial correlations at different scales. Based on cloud and rainfall fields only, the trained networks systematically identify CP pixels in the simulation output. Our methodology may thus open for reliable global CP detection from space-borne sensors. As it also provides information on the spatial extent and the relative positioning of CPs over time, our method may offer new insight into the role of CPs in convective organization.

The imprint of marine atmospheric boundary layer (MABL) dynamical structures on sea surface roughness, as seen from Sentinel-1 Synthetic Aperture Radar (SAR) acquisitions, is investigated. We focus on February 13th, 2020, a case study of the EUREC4A (Elucidating the role of clouds-circulation coupling in climate) field campaign. For suppressed conditions, convective rolls imprint on sea surface roughness is confirmed through the intercomparison with MABL turbulent organization deduced from airborne measurements. A discretization of the SAR wide swath into 25 x 25 km$^2$ tiles then allows us to capture the spatial variability of the turbulence organization varying from rolls to cells. Secondly, we objectively detect cold pools within the SAR image and combine them with geostationary brightness temperature. The geometrical or physically-based metrics of cold pools are correlated to cloud properties. This provides a promising methodology to analyze the dynamics of convective systems as seen from below and above.

In this paper, we applied a variety of statistical methods to study gravity waves in the troposphere and lower stratosphere in the Brazilian sector, using an unprecedented large database from Instituto de Controle do Espaço Aéreo (ICEA) of radiosonde measurements carried out in 2014 on 32 locations in the Brazilian territory totaling 49,652 wind profiles. The average wind profiles were computed and classified by means of a hierarchical cluster analysis. The kinetic and potential energy densities of the gravity waves were determined using a detrending technique based on the least squares method and the Fast Fourier Transform. The time series of the energy densities were analyzed in detail and some persistent and seasonal behavior was found in some cases. A systematic search for quasi monochromatic waves was carried out and the main characteristics of such waves propagating in the troposphere and the lower stratosphere were found. The correlation analysis between the troposphere and the lower ionosphere based on parameters observed on both layers was used to investigate the wave coupling between the two layers. The results we found have implications in the so-called seeding problem of the equatorial ionospheric irregularities.

A month-long data assimilation experiment is carried out to assess the impact of CrIS and IASI Transformed Retrievals (TRs) on the accuracy of analyses and forecasts from a 3-h Weather Research and Forecasting (WRF) cycling system implemented over the central North Pacific Ocean. Conventional observations and satellite MicroWave (MW) radiance data are assimilated along with TRs in comparative experiments. Both the NCEP Global Forecasting System (GFS) and the European Centre for Medium-Range Weather Forecasts (ECMWF) analyses are used in the evaluation process. The results show that the assimilation of TRs, both alone, and in combination with MW radiance assimilation, have the greatest impact on the characterization of the moisture field in the middle atmospheric levels (800 to 300 hPa), and particularly in the lower portion (800 to 600 hPa). The latter improvement is likely due to a refinement in the vertical definition of the trade-wind inversion.

The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) is a scientific mission that collects high spatio-temporal resolution (~70 m, 1-5 days average revisit time) thermal images since its launch on 29 June 2018. As a predecessor of future missions, one of the main objectives of ECOSTRESS is to retrieve and understand the spatio-temporal variations in terrestrial evapotranspiration (ET) and its responses to soil water availability. In the European ECOSTRESS Hub (EEH), by taking advantage of land surface temperature retrievals, we generated ECOSTRESS ET products over Europe and Africa using three structurally contrasting models, namely Surface Energy Balance System (SEBS) and Two Source Energy Balance (TSEB) parametric models, as well as the non-parametric Surface Temperature Initiated Closure (STIC) model. A comprehensive evaluation of the EEH ET products was conducted with respect to flux measurements from 19 eddy covariance sites over 6 different biomes with diverse aridity levels. Results revealed comparable performances of STIC and SEBS (RMSE of ~70 W m-2). However, the relatively complex TSEB model produced a higher RMSE of ~90 W m-2. Comparison between STIC ET estimate and the operational ECOSTRESS ET product from NASA PT-JPL model showed a difference in RMSE between the two ET products around 50 W m-2. Substantial overestimation (>80 W m-2) was noted in PT-JPL ET estimates over shrublands and savannas presumably due to the weak constraint of LST in the model. Overall, the EEH is promising to serve as a support to the Land Surface Temperature Monitoring (LSTM) mission.

On 12 October 2020, the NASA’s Global-scale Observations of the Limb and Disk (GOLD) mission observed three differently shaped EPBs within a 12o longitude range, near the subsatellite point. One is straight aligned to the magnetic field line, whereas the poleward extensions of the others are tilted eastward and westward from the magnetic field line resembling a C-shape and reversed C-shape structures. These EPBs were inside the GOLD imager’s field-of-view for a period of ~3 hours. This allowed us to compute their zonal motion and determine their drift velocities. EPBs’ drift velocities were derived from measuring their longitudinal shifts at the magnetic equator and at both EIA crests. This unique observation, showing three morphologies in a narrow longitude sector, indicates the occurrence of strong longitudinal gradients in the typical parameters associated with the dynamics of EPBs, like neutral winds, electric fields, or other parameters within such a narrow longitude range.

Prior investigations have attempted to characterize the longitudinal variability of the column number density ratio of atomic oxygen to molecular nitrogen (ΣO/N2) in the context of non-migrating tides. The retrieval of thermospheric ΣO/N2 from far ultra-violet (FUV) emissions assumes production is due to photoelectron impact excitation on O and N2. Consequently, efforts to characterize the tidal variability in O/N2 have been limited by ionospheric contamination from O+ radiative recombination at afternoon local times (LT) around the equatorial ionization anomaly. The retrieval of ΣO/N2 from FUV observations by the Ionospheric Connection Explorer (ICON) provides an opportunity to address this limitation. In this work, we derive modified ΣO/N2 datasets to delineate the response of thermospheric composition to non-migrating tides as a function of LT in the absence of ionospheric contamination. We assess estimates of the ionospheric contribution to 135.6 nm emission intensities based on either Global Ionospheric Specification (GIS) electron density, International Reference Ionosphere (IRI) model output, or observations from the Extreme Ultra-Violet imager (EUV) onboard ICON during March and September equinox conditions in 2020. Our approach accounts for any biases between the ionospheric and airglow datasets. We found that the ICON-FUV dataset, corrected for ionospheric contamination based on GIS, uncovered a previously obscured diurnal eastward wavenumber 2 tide in a longitudinal wavenumber 3 pattern at March equinox in 2020. This finding demonstrates not only the necessity of correcting for ionospheric contamination of the FUV signals but also the utility of using GIS for the correction.

Insight and other observations of the Martian surface at different locations have recorded the diurnal variation in surface pressure (Ps) with two rapid fluctuations that occur at dawn and dusk (around LT0800 and LT2000). These short-period surface pressure perturbations at specific local times are typically observed near Martian equinox. Similar phase-locked surface pressure fluctuations over most areas of the middle and low latitudes are simulated by the Martian General Circulation Model at the Dynamic Meteorology Laboratory (LMD). This phenomenon is thus likely to be global rather than local. By reconstructing the surface pressure variation from the horizontal mass flux, the pressure fluctuations in a sol can be attributed to the diurnal variation in the horizontal wind divergence and convergence in the Martain tropical troposphere in the GCM simulations. The background diurnal variation in Ps is related to the diurnal migrating tidal wind, while the enhanced convergence due to the overlap of the 4-hour and 6-hour tides before LT0800 and LT2000 is responsible for the Ps peaks occurring at dawn and twilgith. Although the amplitudes of the 4-hour and 6-hour tides are smaller than those of diurnal tides, the phases of these tides remain similar in the Martain troposphere, which suggests that the convergences and divergences due to 4 h/6 h tidal winds at different altitudes are in phase and together create a mass flux comparable to that induced by diurnal/semidiurnal components and lead to rapid pressure fluctuations.

We examine the distribution of aerosol optical depth (AOD) across 27,707 northern hemisphere (NH) midlatitude cyclones for 2005-2018 using retrievals from the Moderate Resolution Spectroradiometer (MODIS) sensor on the Aqua satellite. Cyclone-centered composites show AOD enhancements of 20-45% relative to background conditions in the warm conveyor belt (WCB) airstream. Fine mode AOD (fAOD) accounts for 68% of this enhancement annually. Relative to background conditions, coarse mode AOD (cAOD) is enhanced by more than a factor of two near the center of the composite cyclone, co-located with high surface wind speeds. Within the WCB, MODIS AOD maximizes in spring, with a secondary maximum in summer. Cyclone-centered composites of AOD from the Modern Era Retrospective analysis for Research and Applications, version 2 Global Modeling Initiative (M2GMI) simulation reproduce the magnitude and seasonality of the MODIS AOD composites and enhancements. M2GMI simulations show that the AOD enhancement in the WCB is dominated by sulfate (37%) and organic aerosol (25%), with dust and sea salt each accounting for 15%. MODIS and M2GMI AOD are 60% larger in North Pacific WCBs compared to North Atlantic WCBs and show a strong relationship with anthropogenic pollution. We infer that NH midlatitude cyclones account for 355 Tg yr-1 of sea salt aerosol emissions annually, or 60% of the 30-80oN total. We find that deposition within WCBs is responsible for up to 35% of the total aerosol deposition over the NH ocean basins. Furthermore, the cloudy environment of WCBs leads to efficient secondary sulfate production.