Search

We thank Dr. Yurkin for his considered and constructive critique of our paper. We respond to his insightful comments below and hope that this discussion motivates further communication and cooperation and benefits to the scattering by small particles community. Key Points: Fenni et al. (2021, https://doi.org/10.1029/2020jd034172) main outcome is the enhancement in accuracy and computational cost of MIDAS with use of higher‐order quadrature schemesMajor conclusions of Fenni et al. (2021, https://doi.org/10.1029/2020jd034172) are validated on a large set of realistic arbitrarily‐shaped ice aggregates [ABSTRACT FROM AUTHOR]

This paper introduces the special collection in Geophysical Research Letters and Journal of Geophysical Research: Atmospheres on the exceptional stratospheric polar vortex in 2019/2020. Papers in this collection show that the 2019/2020 stratospheric polar vortex was the strongest, most persistent, and coldest on record in the Arctic. The unprecedented Arctic chemical processing and ozone loss in spring 2020 have been studied using numerous satellite and ground‐based data sets and chemistry‐transport models. Quantitative estimates of chemical loss are broadly consistent among the studies and show profile loss of about the same magnitude as in the Arctic in 2011, but with most loss at lower altitudes; column loss was comparable to or larger than that in 2011. Several papers show evidence of dynamical coupling from the mesosphere down to the surface. Studies of tropospheric influence and impacts link the exceptionally strong vortex to reflection of upward propagating waves and show coupling to tropospheric anomalies, including extreme heat, precipitation, windstorms, and marine cold air outbreaks. Predictability of the exceptional stratospheric polar vortex in 2019/2020 and related predictability of surface conditions are explored. The exceptionally strong stratospheric polar vortex in 2019/2020 highlights the extreme interannual variability in the Arctic winter/spring stratosphere and the far‐reaching consequences of such extremes. Plain Language Summary: The Arctic stratospheric polar vortex—a band of strong winds roughly encircling the pole at about 65°N latitude from about 15 to 50 km above the Earth's surface that forms every winter—was exceptionally strong during the 2019/2020 winter. The strong vortex in the stratosphere was linked to unusual conditions at both higher and lower altitudes. This collection of papers explores the far‐reaching consequences of the exceptionally strong stratospheric polar vortex in 2019/2020, including impacts on Arctic chemical ozone loss and on surface weather conditions. Chemical ozone loss in spring 2020 matched or exceeded the most previously on record (for 2011) and showed some features similar to the larger loss that occurs over the Antarctic every spring. The exceptionally strong stratospheric polar vortex was linked to weather extremes, including record heat, unusual patterns of precipitation, marine cold air outbreaks, and windstorms. Key Points: The stratospheric polar vortex in 2019/2020 was the strongest and longest‐lasting on record as described in this special collectionThis exceptionally strong and cold polar vortex led to unprecedented Arctic ozone loss, approaching that in some Antarctic wintersCirculation anomalies linked to the vortex spanned the mesosphere to the surface with implications for extreme weather and predictability [ABSTRACT FROM AUTHOR]

The Polar Mesospheric Cloud (PMC) Turbulence experiment performed optical imaging and Rayleigh lidar PMC profiling during a 6‐day flight in July 2018. A mosaic of seven imagers provided sensitivity to spatial scales from ∼20 m to 100 km at a ∼2‐s cadence. Lidar backscatter measurements provided PMC brightness profiles and enabled definition of vertical displacements of larger‐scale gravity waves (GWs) and smaller‐scale instabilities of various types. These measurements captured an interval of strong, widespread Kelvin‐Helmholtz instabilities (KHI) occurring over northeastern Canada on July 12, 2018 during a period of significant GW activity. This paper addresses the evolution of the KHI field and the characteristics and roles of secondary instabilities within the KHI. Results include the imaging of secondary KHI in the middle atmosphere and multiple examples of KHI "tube and knot" (T&K) dynamics where two or more KH billows interact. Such dynamics have been identified clearly only once in the atmosphere previously. Results reveal that KHI T&K arise earlier and evolve more quickly than secondary instabilities of uniform KH billows. A companion paper by Fritts et al. (2022), https://doi.org/10.1029/2021JD035834 reveals that they also induce significantly larger energy dissipation rates than secondary instabilities of individual KH billows. The expected widespread occurrence of KHI T&K events may have important implications for enhanced turbulence and mixing influencing atmospheric structure and variability. Key Points: First observation of unambiguous secondary Kelvin‐Helmholtz instabilities in high‐resolution images of the polar mesospheric cloud layerIdentification and quantification of Kelvin‐Helmholtz billow interactions leading to tubes and knots and accelerated billow breakdownEstimation of turbulence Reynolds number Returb ∼ 5,000 and νturb ∼ 3 times larger than the kinematic viscosity [ABSTRACT FROM AUTHOR]

The expression "lightning never strikes twice" is questioned in this paper because, among the randomness of lightning impacts, some spots are hit even more than twice year after year. This article introduces the recurrent lightning spots (RLS) concept, which are locations periodically impacted by cloud‐to‐ground lightning every consecutive year over a certain period. RLS are investigated in two regimes, with markedly different lightning climatology but similar orography, for 10 consecutive years: Catalonia (North East of Spain, Europe) and Barrancabermeja (North Central Colombia, South America). Results revealed 148 and 916 RLS in Catalonia and Barrancabermeja, respectively. RLS in both regions are typically found to be related to tall structures, mountain peaks, and steep terrain. The method allowed us to identify those tall towers and orographic relief frequently affected by lightning that are not detected with the mere computation of the ground flash density. In the case of Catalonia, some RLS are found offshore. Besides the scientific interest in understanding lightning, the new concept of RLS provides additional and valuable information applicable to lightning protection engineering. Plain Language Summary: The expression "lightning never strikes twice" is questioned in this paper because it shows that some spots are hit even more than twice: year after year. The RLSs (RLS) concept is presented, which are locations where lightning strike them periodically. The two regions of study with markedly different lightning climatology but similar orography are Catalonia (North East of Spain, Europe) and Barrancabermeja (in North Central Colombia, South America). RLS in both regions are typically found to be related to tall structures, mountain peaks and steep terrain. In the case of Catalonia, some RLS are found offshore. RLS new concept provides valuable information to lightning protection engineering. Key Points: Recurrent lightning spots (RLS) are targets with at least one cloud‐to‐ground lightning per year during a certain number of consecutive yearsIn Catalonia (Spain), 13% of the RLS are tall towers, 72% are mountainous peaks and 5% are offshore sitesThe method effectively finds tall structures affected by lightning in tropical climates [ABSTRACT FROM AUTHOR]

In a recent paper Fenni et al. (2021, https://doi.org/10.1029/2020jd034172) compared the code MIDAS, based on the direct solution of the volume‐integral equation combined with advanced cubatures for orientation averaging, to the code DDSCAT, a state‐of‐the‐art implementation of the discrete dipole approximation. This comment highlights methodological issues in this comparison and shows that the quantitative claims of Fenni et al. (2021, https://doi.org/10.1029/2020jd034172), related to superiority of MIDAS over DDSCAT, are based on very specific test cases with respect to particle symmetries or initial orientation, as well as to the selected scattering quantity of interest. Thus, these claims are not expected to hold for other similar particles. Moreover, the detailed discussion of these issues is relevant for all light‐scattering simulation methods, except those allowing analytical orientation averaging. Thus, the comment constructs general guidelines for fair evaluation of orientation‐averaging techniques in a wide range of light‐scattering methods and computer codes. Plain Language Summary: The paper discusses several issues that appear when one is comparing different orientation‐averaging techniques (cubatures) in combination with the same or different light‐scattering simulation methods. Fair evaluation of cubature performance in realistic general scenarios is important both for practitioners (to choose the most efficient combination of the existing codes and cubatures) and for code developers (to set their priorities on the new features with the largest expected benefits). Unfortunately, the performance of the cubatures is complexly interwoven with the internals of the simulation methods and depends on specific test particles and computed scattering quantities. This questions the generality of conclusions in some previous publications. Based on this discussion, the paper ends with general guidelines for fair evaluation of cubatures, allowing future studies to arrive at general conclusions, so that they can be directly used by other researchers. Key Points: Quantitative conclusions of Fenni et al. (2021) are based on very specific test casesOrientation‐averaging techniques should be compared on non‐symmetric particles, and not with a special initial orientationAny comparison of simulated results should consider their uncertainties accounting for all sources of errors [ABSTRACT FROM AUTHOR]

Wildfire‐related aircraft field campaigns frequently offer opportunities to validate remote‐sensing retrievals of aerosol properties and other quantities derived from satellite‐borne‐instrument observations. Satellite instruments often provide regional context‐imagery for more sparsely sampled aircraft and surface‐based measurements. However, aerosol amount, particle type, aerosol plume height and the associated wind vector products retrieved from the NASA Earth Observing System's Multi‐angle Imaging SpectroRadiometer (MISR) instrument have matured sufficiently that these quantities can also contribute substantially to a campaign data set, in regional context. This is especially useful when such measurements are not acquired at all from the suborbital platforms. During NOAA's California Fire Dynamics Experiment (CalFiDE), aircraft operations were coordinated with MISR overpasses on two occasions: for the Rum Creek fire on 30 August 2022, and for the Mosquito fire on 08 September. MISR‐retrieved aerosol properties show distinctly different patterns of black and brown smoke particle distributions and inferred plume evolution in the two cases. This paper presents the satellite‐retrieved results that complement the field observations, demonstrating what such measurements can offer, and contributing material for detailed fire dynamics and chemistry studies when combined with the CalFiDE suborbital observations and models in continuing studies. Plain Language Summary: A common use of aircraft field campaigns is to validate the radiances measured by space‐based instruments and the geophysical quantities derived from the satellite observations. However, satellite aerosol amount and properties derived from the NASA Earth Observing System's Multi‐angle Imaging SpectroRadiometer (MISR) instrument are sufficiently mature that they can also contribute directly to field‐campaign data sets. During NOAA's CalFiDE campaign is summer 2022, on two occasions the aircraft observed wildfire smoke plumes coordinated with MISR overpasses: for the Rum Creek fire on 30 August 2022, and for the Mosquito fire on 08 September. In addition to providing broad spatial context to the much more spatially limited aircraft measurements, the MISR results offer geometrically‐derived smoke‐plume height and plume‐level motion vectors from which smoke age can be estimated downwind along the plume. From MISR‐retrieved constraints on particle size, shape, and light‐absorption properties, the distribution of black and brown smoke can be inferred, along with the underlying processes responsible for plume‐particle evolution. This paper presents the satellite‐retrieved results that complement the suborbital data acquired for the CalFiDE campaign and associated modeling, for use in continuing studies of fire dynamics. Key Points: NASA's MISR multi‐angle imagery allows aerosol plume‐height, associated motion vector, and particle property retrievals from spaceTypically, field data are acquired to validate satellite data, but MISR data is mature enough to contribute directly to the CalFiDE campaignAircraft and MISR overflights were coordinated twice, yielding MISR aerosol context and detail for joint smoke‐plume dynamics and chemistry study [ABSTRACT FROM AUTHOR]

Volcanic ash clouds are carefully monitored as they present a significant hazard to humans and aircraft. The primary tool for forecasting the transport of ash from a volcano is dispersion modeling. These models make a number of assumptions about the size, sphericity and density of the ash particles. Few studies have measured the density of ash particles or explored the impact that the assumption of ash density might have on the settling dynamics of ash particles. In this paper, the raw apparent density of 23 samples taken from 15 volcanoes are measured with gas pycnometry, and a negative linear relationship is found between the density and the silica content. For the basaltic ash samples, densities were measured for different particle sizes, showing that the density is approximately constant for particles smaller than 100 μm, beyond which it decreases with size. While this supports the current dispersion model used by the London Volcanic Ash Advisory Centre (VAAC), where the density is held at a constant (2.3 g cm−3), inputting the measured densities into a numerical simulation of settling velocity reveals a primary effect from the silica content changing this constant. The VAAC density overestimates ash removal times by up to 18%. These density variations, including those varying with size beyond 100 μm, also impact short‐range particle‐size distribution measurements and satellite retrievals of ash. Plain Language Summary: Volcanic ash clouds are carefully monitored as they present a significant hazard to humans and aircraft. Dispersion modeling is a primary tool used to forecast ash flows from volcanoes. These models make a number of assumptions about the size, sphericity (roundness) and density of the ash particles. Few studies have measured the density of ash particles or explored the impact that the assumption of ash density might have on the dispersion forecasts. In this paper, the density of 23 samples taken from 15 volcanoes are measured, and a negative linear relationship is found between the density and the silica content. For the basaltic ash samples (the most common type of ash), densities were measured for different particle sizes, showing that the density is approximately constant for particles smaller than 100 μm, beyond which it decreases with size. This supports the London Volcanic Ash Advisory Centre keeping density constant in their current model, but in fact this constant changes with silica content, leading to an overestimation of ash removal times by up to 18%. These density deviations also impact short‐range particle‐size distribution measurements and satellite retrievals of ash. Key Points: The density of volcanic ash is measured as a function of particle size for a range of eruptionsSilica content and particle size negatively correlate with densityThe density of particles smaller than 100 μm is approximately constant but is dependent on silica content [ABSTRACT FROM AUTHOR]

Spatial climate data sets that extend back in time over many decades are an important resource for climate monitoring. The long‐term consistency of such data sets is, however, compromised by changes in the measurement systems over time. In this paper, we introduce a data set of monthly precipitation on a 5‐km grid over the European Alps that extends back to the late 19th century. In deriving the "long‐term Alpine precipitation reconstruction" (LAPrec), special care is taken of variations in the station network, in order for the data set to satisfy high standards in long‐term consistency. LAPrec builds on a reconstruction method that integrates the available information in two portions: The first is a set of high‐quality homogenized station series, taken from the HISTALP data archive, covering the entire period almost continuously. The second is a high‐resolution gridded precipitation analysis, taken from the "Alpine Precipitation Grid Data Set," constructed from thousands of rain‐gauges but covering a few decades only. We demonstrate how the reconstruction approach successfully introduces mesoscale structures that are not resolved by the available long‐term station series, more plausibly so than a predecessor data set using conventional interpolation. We also illustrate that LAPrec reveals long‐term precipitation trends that are spatially more consistent and more detailed than the trends in popular climate monitoring data sets. Over the period 1871–2017 a statistically significant increase is found in winter over the northern parts of the Alps (1%–2% per 10 years). LAPrec is available in two versions (back until 1871 and 1901 respectively) from the Copernicus climate data store. Plain Language Summary: Data sets extending back in time over many decades are an important resource for climate monitoring. In this paper, we introduce a data set on a regular grid of monthly precipitation over the European Alps that extends back to 1871. In deriving the "long‐term Alpine precipitation reconstruction" (LAPrec), special care is taken to keep the station network constant and make use of high‐quality homogenized time series, in order for the data set to be as consistent as possible over time. LAPrec builds on two sources of information: The first is a set of high‐quality station series from the HISTALP data archive, covering the entire period almost continuously. The second is a high‐resolution gridded precipitation analysis, the "Alpine Precipitation Grid Data Set," constructed from thousands of rain‐gauges but covering a few decades only. We demonstrate how this approach successfully introduces details that are not resolved by the available long‐term station series. We also illustrate that LAPrec reveals long‐term precipitation trends that are spatially more consistent and detailed than the trends in popular climate monitoring data sets. Over the period 1871–2017 a statistically significant increase in winter precipitation is found over the northern parts of the Alps. LAPrec is available from the Copernicus climate data store. Key Points: Long‐term spatial analyses of precipitation in the European Alps since 1871For applications such as climate monitoring and trend analysis, the data set meets high standards in long‐term consistency and homogeneityLong‐term Alpine precipitation reconstruction detects a significant trend in winter precipitation over the northern part of the Alps in the period starting from 1871 [ABSTRACT FROM AUTHOR]

This paper makes use of a combination of FLEXPART‐WRF simulations forced with ERA5 and the CESM2 model—incorporated in the CMIP6 project—to infer a series of changes over the present century in the behavior of the landfalling atmospheric rivers (ARs) arriving to the Iberian Peninsula. In addition, future changes in the intensity and position of their main moisture sources are studied. In overall terms, there is a noticeable increase in the amount of moisture transported by ARs in the study region, particularly accentuated by the end of the century. However, no significant changes in the number of events are observed. A northward shift of both the mean position of the ARs as well as their main sources of moisture is also detected, particularly for the end of the century, and in the summer and fall months. In relation to the latter, an increase in the contribution of moisture contribution is also observed, quantitatively compatible with Clausius‐Clapeyron amplification. Plain Language Summary: This paper makes use of a combination of simulations forced with reanalysis data and a climate model to infer a series of changes over the present century in the behavior of the landfalling atmospheric river—ARs, regions of intense moisture transport located in the lower layers of the atmosphere—arriving at the Iberian Peninsula. In addition, future changes in the intensity and position of their main moisture sources are studied. In overall terms, there is a noticeable increase in the amount of moisture transported by ARs in the study region, particularly accentuated by the end of the century. However, no significant changes in the number of events are observed. A northward shift of both the mean position of the ARs as well as their main sources of moisture is also detected, particularly for the end of the century, and in the summer and fall months. In relation to the latter, an increase in the contribution of moisture contribution is also observed, in a ratio similar to that expected. Key Points: FLEXPART‐WRF forced with CESM2 model has been able to reproduce the historical conditions of Atmospheric River over the Iberian PeninsulaA northward shift of the main source regions is projected, notable in summer and fall and particularly by the end of the centuryGradual strengthening in the intensity of Atmospheric Rivers is expected, observable from an increase in the amount of moisture transported [ABSTRACT FROM AUTHOR]

In this paper we analyze electric‐field and current measurements of upward leaders induced by a downward negative lightning flash that struck a residential building. The attachment process was recorded by two high‐speed cameras running at 37,800 and 70,000 images per second and the current measured in two lightning rods. Differently from previous works, here we show, for the first time, current measurements of multiple upward leaders that after initiation propagate to connect the negative downward moving leader. At the beginning of the propagation of the leaders that initiate on the instrumented lightning rods, current pulses appear superimposed to a steadily increasing DC current. The upward leader current pulses increase with the approach of the downward leader and are not synchronized but present an alternating pattern. All 2D leader speeds are approximately constant. The upward leaders are slower than the downward leader speed. The average time interval between current pulses in upward leaders is close to the interstep time interval found by optical or electric field sensors for negative cloud‐to‐ground stepped leaders. The upward leaders respond to different downward propagating branches and, as the branches alternate in propagation and intensity, so do the leaders accordingly. Right before the attachment process the alternating pattern of the leaders ceases, all downward leader branches intensify, and consequently upward leaders synchronize and pulse together. The average linear densities for upward leaders (49 and 82 μC/m) were obtained for the first time for natural lightning. Plain Language Summary: The effectiveness of a lightning protection system depends on its efficiency to intercept the down coming leader of a cloud‐to‐ground lightning flash. The interception is usually done by an upward connecting leader that initiates from grounded structures, humans, or living beings that protrude from nearby ground. The understanding of the upward connecting leader and of the attachment process with the downward leader plays an important role in the determination of the zone of protection and therefore in the improvement of a lightning protection system. Unconnected upward leaders, that is, upward leaders that fail to connect the downward leader, are also of great importance in lightning protection. They can be large enough to cause damage to equipment vulnerable to sparks or induced currents, and enough to injure people from who it initiates. In this paper we analyze electric‐field, speed, and current measurements of upward leaders induced by a downward negative lightning flash that struck a residential building. The attachment process was simultaneously recorded by two high‐speed cameras, an electric‐field sensor, and current sensors installed on two lightning rods. Differently from previous works we show, for the first time, current measurements of multiple upward leaders induced by the negative downward moving leader. Key Points: Current and charge density measurements of two upward leaders induced by the same downward leaderUpward leaders alternate their progression during initial propagationCurrent pulses of upward leaders increase intensity and synchronize right before attachment [ABSTRACT FROM AUTHOR]

Volcanic aerosols change the atmospheric composition and thereby affect weather and climate. Aerosol dynamic processes such as nucleation, condensation, and coagulation modify the shape, size, and mass of aerosol particles, which influence their atmospheric lifetime and radiative properties. Nevertheless, most models omit these processes for ash particles. In this work, we explore the ash aerosol aging and sulfate production during the first 4 days following the 2021 La Soufrière (St. Vincent) eruption with the ICON‐ART model (ICOsahedral Nonhydrostatic model with Aerosol and Reactive Trace gases). Online coupling of ICON‐ART with a one‐dimensional volcanic plume model calculates volcanic emission, which makes it possible to resolve the different eruption phases of the noncontinuous La Soufrière eruption. We compared our simulated aerosol distribution and composition with observations from the Cloud‐Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument, the Multiangle Imaging SpectroRadiometer (MISR) Research Aerosol (RA) Algorithm, and the Barbados Cloud Observatory (BCO). We show that online coupling is essential to adequately model the emissions and plume development close to the volcano. The modeled aerosol aging is in very good agreement with observations from MISR near the emission source and with CALIOP at larger distances. Furthermore, particle aging occurs faster in the troposphere than in the stratosphere due to the availability of water vapor and OH, but a layer of coated ash appears at the plume top due to faster oxidation of SO2 and lofting by aerosol‐radiation interaction. This paper gives the first direct comparison of aerosol aging in volcanic eruption plumes between simulations and observations. Plain Language Summary: Large volcanic eruptions can influence weather and climate, and endanger aviation and public health. To constrain these effects and risks, it is critical to reliably predict the volcanic plume dispersion. However, the atmospheric lifetime of the ash released during an eruption is influenced by many factors, such as emission height, meteorology, and aerosol dynamical processes. Aerosol dynamical processes lead to growth and aging of volcanic plume particles. They include the formation of new particles from precursor gases, condensation on existing particles, and coagulation of particles. This paper investigates the formation of aged ash particles in model simulations and observations following the 2021 La Soufrière eruption. We consider ash aging both close to the volcano and during further transport. We found that ash aging takes place already close to the volcano and the fraction of aged particles increases with distance from the source. During further transport, a layer of aged ash particles forms at the plume top due to interaction of these particles with radiation and subsequent warming of the plume. Key Points: Resolving individual eruption phases is essential for modeling the 2021 La Soufrière eruptionCombination of modeling and satellite data confirm, for the first time, the aging of volcanic ashAsh aging and sulfate production rates depend on distance from the source and altitude [ABSTRACT FROM AUTHOR]

This paper examines the empirical relationship between the Madden–Julian oscillation (MJO), the quasi‐biennial oscillation (QBO), and atmospheric river (AR) activity and precipitation in California on subseasonal time scales. We introduce an experimental forecast tool that uses observed anomaly patterns during a 38 yr period to predict the probability of above‐ and below‐normal AR activity and precipitation at lead times of 1–6 weeks based on the phase and amplitude of the MJO and QBO. The hindcast prediction skill of probabilistic AR activity and precipitation forecasts is evaluated for Northern, Central, and Southern California, as well as two sets of smaller geographical domains. These smaller domains are more relevant for water resource management and allow us to investigate the sensitivity of prediction skill to domain size. Consistent with previous studies, our results demonstrate that subseasonal AR activity and precipitation in California are strongly modulated by the MJO and QBO. The anomaly patterns of AR activity and precipitation vary considerably throughout the cool season, with a tendency toward below‐normal AR activity and precipitation during easterly QBO and above‐normal AR activity and precipitation during westerly QBO in JFM. The opposite patterns are generally observed in OND, but the anomaly signals are weaker and less coherent for AR activity. Certain combinations of MJO phase, QBO phase, lag time, and season yield notably higher skill scores, reinforcing the notion of "windows of opportunity" for skillful subseasonal‐to‐seasonal predictions. In California, these forecasts of opportunity are predominantly associated with easterly QBO in JFM and FMA. Plain Language Summary: This paper examines the relationship between the Madden–Julian oscillation (MJO), the quasi‐biennial oscillation (QBO), and atmospheric river (AR) activity and precipitation in California. We introduce an experimental forecast tool that shows the likelihood of above‐normal and below‐normal AR activity and precipitation based on the phases of the MJO and QBO. Consistent with previous studies, our results demonstrate that AR activity and precipitation in California are strongly influenced by the MJO and the QBO. There is a tendency for below‐normal AR activity and precipitation during easterly QBO conditions in January–March. The opposite pattern (above‐normal AR activity and precipitation) generally occurs during westerly QBO conditions. The results also suggest that our forecast tool has some potential to improve the prediction of AR activity and precipitation. The reliability and usefulness of this forecast tool depend on multiple factors, including the MJO phase, the QBO phase, and the time of year. In California, this forecast tool is likely to provide the most beneficial value during easterly QBO conditions in mid‐winter through early spring. Key Points: The modulation of atmospheric river (AR) activity and precipitation in California by the Madden‐Julian oscillation and quasi‐biennial oscillation is quantifiedA hindcast skill assessment of probabilistic AR activity and precipitation forecasts is presentedAn experimental subseasonal AR activity and precipitation forecast tool is introduced [ABSTRACT FROM AUTHOR]

This paper utilizes Aerosol Index (AI) data from the Total Ozone Mapping Spectrometer (TOMS) instrument, along with ERA5 reanalysis data, to identify atmospheric processes contributing to the uplift of dust from the Etosha Pan through the annual cycle. Etosha is one of the most prominent source areas in the Southern Hemisphere, although very little is known about its meteorology outside of the peak dust season (August–October). Emissions in December (AI = 1.6) are shown to be comparable to those in September (AI = 1.7), the dustiest month in the TOMS record. Unlike in September however, when a nocturnal low‐level jet is the primary emission mechanism, uplift in December is associated with an anabatic‐sea breeze that develops along the Namib coast, and propagates inland to reach Etosha during the evening. The system is a response to the thermal contrast between the elevated interior plateau and the adjacent waters of the cool Benguela Upwelling System, and so is at its strongest during austral summer, when the area of maximum diabatic heating shifts south over southern Africa. Topographic channeling of the flow through the east‐west orientated Hoanib River valley is shown to facilitate the inland propagation of the anabatic‐sea breeze, and explains the persistence of the system at Etosha's latitude. Evening surface winds at Etosha, associated with the anabatic‐sea breeze, are significantly stronger in the dustier December months, when diabatic heating over the subcontinent and hence the zonal thermal gradient are enhanced. Plain Language Summary: This paper uses satellite and meteorological data to identify the features of southern Africa's weather and climate that contribute to dust emission from the Etosha Pan throughout the year. Etosha is an important source of dust in the Southern Hemisphere, however very little work has been conducted there outside of the winter season, when emissions are at their highest. From the satellite data, it is shown that December is just as dusty as some of the winter months, however there is a difference in the low‐level winds between the two seasons; in winter, emissions are driven by a morning peak in surface wind speeds, whereas emissions in December are driven by maximum surface winds during the evening. This evening peak in surface winds is shown to coincide with the arrival of a sea breeze at Etosha, with the system a response to the strong heating of the southern African plateau at this time of year. Evening surface winds associated with the sea breeze are stronger during the dustiest December months at Etosha, and are driven by enhanced heating over the subcontinent. Key Points: An anabatic‐sea breeze helps to drive austral summer dust emissions from one of the Southern Hemisphere's most prominent source areasThe anabatic‐sea breeze is present throughout the year, however is at its strongest, and propagates furthest inland during austral summerVariability in the strength of the system is driven by changes in the pattern of diabatic heating over the interior of southern Africa [ABSTRACT FROM AUTHOR]

Following on our earlier single‐station, 2‐dimensional (2D) broadband interferometric mapping and polarization (BIMAP) observations of lightning discharges, we recently deployed two BIMAP stations at Los Alamos National Laboratory to map the lightning sources and their polarization in full 3‐dimensional (3D) space (BIMAP‐3D). The two stations are separated by 11.5‐km and each station consists of four antenna sets (instead of three for the original BIMAP) that form a Y‐shaped array for improved interferometric performance. In this paper, we report the BIMAP‐3D system design, a generalized and analytical 2D interferometry technique for noncoplanar antenna array, a two‐stage 3D mapping technique based on geometric triangulation and baseline‐based differential time of arrival, and a technique to reconstruct the polarization orientation in 3D space by combining the 2D polarization results from the two‐station observations. Along with description of the techniques, we demonstrate and discuss the initial lightning results, including 3D maps for a hybrid intracloud and cloud‐to‐ground flash and for a normal intracloud flash, development of abnormal K‐change leaders, and polarization signatures for a K‐change leader. We find that with the two‐stage 3D mapping techniques, the sources can be located to meters accuracy for a favorable event that occurs between the two stations. We also find the polarization vectors for the example K leader are mostly orthogonal to the leader channel after the full 3D polarization analysis. The main purpose of this paper is to report the BIMAP‐3D techniques and capabilities. Detailed analysis of more specific discharge processes will be reported in later studies. Plain Language Summary: A new 3‐dimensional broadband radio frequency interferometric mapping and polarization system (BIMAP‐3D) is developed and deployed at Los Alamos National Laboratory for lightning research. BIMAP‐3D provides an unprecedented capability in high‐resolution, time‐evolving 3D lightning source mapping and 3D source polarization detection for detailed study of lightning discharge physics. In this research, we described the BIMAP‐3D system, introduced a suite of advanced data processing techniques, and demonstrated BIMAP‐3D's capabilities with actual lightning observations. This new capability is expected to lead to a range of new understandings and discoveries for a variety of lightning discharge processes. Key Points: A new 3D broadband interferometric mapping and polarization system for lightning study is introducedA suite of new data process algorithms is reported, and 3D lightning results for overall flashes and K‐change leaders are demonstratedPolarization orientations for a K‐change leader are found mostly orthogonal to the leader channel [ABSTRACT FROM AUTHOR]

Based on directly measured triggered lightning currents, the characteristics of 70 continuing current (CC) waveforms and 106 M‐component waveforms were analyzed. The durations of CC without M‐component are all less than 10 ms, which are significantly less than those of CCs with M‐components. The first M‐component always appears at no more than 4 ms after the return stroke. The characteristics of the superimposed M‐components will be different for CC of different durations. The M‐components superimposed on CCs with duration greater than 10 ms are 3–4 times larger in terms of risetime, half‐peak width, duration, and transfer charge than that superimposed on CCs with duration less than 10 ms. The time difference (TD) between the peak electric field and the peak current of the M‐component is affected by both the continuing current level (ICC) and the magnitude of the M‐component (IM). TD and ICC (or IM) will not be large at the same time. The duration of a CC following a return stroke with a peak current greater than 21.2 kA will not exceed 46.0 ms. Key Points: The characteristics of continuing current and M‐component current waveforms from triggered lightning were analyzedThe characteristics of the M‐components superimposed on continuing currents with different durations are differentThe time difference between the peak electric field and the peak current of the M‐component is related to two current parameters [ABSTRACT FROM AUTHOR]

Lidar and radar observations of persistent atmospheric wave activity in the Antarctic atmosphere motivate investigation of generation of acoustic‐gravity waves (AGWs) by vibrations of ice shelves and exploiting their possible ionospheric manifestations as a source of information about the ice shelves' conditions and stability. A mathematical model of the waves radiated by vibrations of a finite area of the lower boundary of the atmosphere is developed in this paper by extending to AGWs an efficient, numerically exact approach that was originally developed in seismology and underwater acoustics. The model represents three‐dimensional wave fields as Fourier integrals of numerical or analytical solutions of a one‐dimensional wave equation and accounts for the source directionality, AGW refraction and diffraction, and the wind‐induced anisotropy of wave dissipation. Application of the model to the generation of atmospheric waves in Antarctica by free vibrations of the Ross Ice Shelf reveals a complex three‐dimensional structure of the AGW field and elucidates the impact of various environmental factors on the wave field. The intricate variation of the wave amplitude with altitude and in the horizontal plane is shaped by the spatial spectrum of the ice surface vibrations and the temperature and wind velocity stratification from the troposphere to the mesosphere. It is found that the waves due to the low‐order modes of the free oscillations of the Ross Ice Shelf, which have periods of the order of several hours, can transport energy to the middle and upper atmosphere in a wide range of directions from near‐horizontal to near‐vertical. Plain Language Summary: This research paper paves the way to infer the conditions and stability of ice shelves in Antarctica by looking at unusual wave activity in the atmosphere. The researchers have developed a mathematical model to understand how these waves, called acoustic‐gravity waves (AGWs), are created by the vibrations of ice shelves. The model is based on a method previously used in seismology and underwater acoustics and accounts for various factors that affect the waves, such as the properties of the source, the way the waves bend and spread due to wind, and how they dissipate. The researchers applied this model to study the atmospheric waves generated by vibrations of the Ross Ice Shelf in Antarctica. The results show a complex 3D structure of the AGW field, highlighting the impact of different environmental factors on the wave activity. The variation in wave amplitude depends on the ice surface vibrations and the temperature and wind conditions at different heights in the atmosphere. The study found that waves with periods of several hours can transfer energy from the ice shelf to the middle and upper atmosphere in various directions. This new approach could help scientists better understand the conditions and stability of ice shelves in the future. Key Points: Vibrations of large ice shelves radiate atmospheric waves that can travel in a wide range of directions from near‐horizontal to near‐zenithSimple, numerically efficient model is developed of atmospheric wave generation by finite sources at the ground levelWith the surface vibrations' spectrum as the input, the model quantifies atmospheric perturbations from the troposphere to the thermosphere [ABSTRACT FROM AUTHOR]

Multi‐resolution analysis methods can reveal the underlying physical dynamics of nonstationary signals, such as those from lightning. In this paper we demonstrate the application of two multi‐resolution analysis methods: Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) in a comparative way in the analysis of electric field change waveforms from lightning. EEMD and VMD decompose signals into a set of Intrinsic Mode Functions (IMFs). The IMFs can be combined using distance and divergence metrics to obtain noise reduction or to obtain new waveforms that isolate the physical processes of interest while removing irrelevant components of the original signal. We apply the EEMD and VMD methods to the observations of three close Narrow Bipolar Events (NBEs) that were reported by Rison et al. (2016, https://doi.org/10.1038/ncomms10721). The ΔE observations reveal the occurrence of complex oscillatory processes after the main NBE sferic. We show that both EEMD and VMD are able to isolate the oscillations from the main NBE, with VMD being more effective of the two methods since it requires the least user supervision. The oscillations are found to begin at the end of the NBEs' downward fast positive breakdown, and appear to be produced by a half‐wavelength standing wave within a weakly‐conducting resonant ionization cavity left behind in the wake of the streamer‐based NBE event. Additional analysis shows that one of the NBEs was likely initiated by an energetic cosmic ray shower, and also corrects a misinterpretation in the literature that fast breakdown is an artifact of NBE‐like events in interferometer observations. Plain Language Summary: This paper investigates the application of mode decomposition techniques to the analysis of Narrow Bipolar Events (NBEs). NBEs are high‐power discharges that often occur as the initiating event of lightning flashes, and are produced by streamer‐based activity called fast positive breakdown. We apply Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) methods to three NBEs that were observed at Langmuir Laboratory in New Mexico to extract and interpret oscillatory behavior that occurred following the NBEs. We show that both EEMD and VMD separate the oscillations from the strong electric field change of the parent NBE, with VMD being the preferred choice. The resulting waveforms are indicative of a shock‐excited residual process that lasts for tens of microseconds, beginning at the end of the downward fast breakdown activity. Although questions remain about the physical mechanism of the oscillations, they appear to be caused by the NBE's streamers creating a weakly‐conducting resonant cavity in its wake that supports half‐wavelength standing wave oscillations, analogous to the vibrations of a plucked guitar string. Key Points: Ensemble Empirical Mode Decomposition and Variational Mode Decomposition separate out Narrow Bipolar Event (NBE) oscillations for further studies of physical mechanisms responsible for oscillationsThe results show the oscillations initiate not during the NBE but at the end of the fast positive breakdown processThe oscillatory behavior indicates that residual ionic conductivity creates a resonant cavity which emits half‐wavelength radiation [ABSTRACT FROM AUTHOR]

A method to identify the turbulent mixing sources within the marine atmospheric boundary layer (MABL) based on the shipborne coherent Doppler lidar measurements is introduced in this paper. Combining with the coherent Doppler lidar signal‐to‐noise ratio, vertical velocity skewness, turbulence kinetic energy dissipation rate, and wind shear intensity, the categories of turbulent mixing sources and the specific turbulent mixing sources could be determined. The method is applied into two voyages of MABL observation during May 2021 in the South China Sea and during April 2022 in the Bohai Sea and Yellow Sea. The turbulent mixing processes are captured and the classification of the turbulence driven sources within the MABL are realized. The temporal‐spatial evolution characteristics of the turbulence mixing process in the MABL are investigated under different weather conditions containing clear‐sky day, cloudy‐sky day, and sea‐fog day. The convective mixing process is recognized in the daytime of the clear‐sky day and the intermittent cloud‐driven turbulence exists below the cloud layer. Additionally, the turbulent mixing is weak which could not act as the main driven source during the sea‐fog day. Furthermore, the dominant turbulence scale analyses of different turbulence sources are conducted based on the cospectra of the vertical velocity and the horizontal speed measurements. The turbulence parameters of different turbulence sources are statistical analyzed and compared in different sea areas. The classification method has the broad application prospects on the study of the air‐sea interaction. Plain Language Summary: The redistribution of substance and energy within the atmospheric boundary layer is achieved through the turbulent mixing in most cases. When the underlying surface is ocean, the complex turbulent mixing process within the marine atmospheric boundary layer (MABL) is the key topic in the air‐sea interaction research because it would affect the atmosphere circulation through changing the momentum, heat and water vapor distributions. Hence the quantitative measurements of the vertically resolved turbulence parameters and understanding of the main sources of the turbulent mixing are crucial. In this paper, the turbulent mixing driven sources are identified based on the high‐accuracy turbulence parameters measured by coherent Doppler lidar. Through the two voyages of MABL observation in the South China Sea, Bohai Sea, and Yellow Sea, the temporal‐spatial evolution characteristics of the turbulence mixing process are investigated under different weather conditions including clear‐sky day, cloudy‐sky day, and sea‐fog day. Additionally, the scale analysis focuses on the dominant turbulence and the statistical analysis of turbulence parameters of different turbulent driven sources are conducted in different sea areas. This classification method has the great potential and broad application prospects on the study of the turbulent mixing characteristics and air‐sea interaction. Key Points: A method to identify the turbulent mixing sources in marine atmospheric boundary layer based on coherent Doppler lidar is introducedThe classification method is first demonstrated with the lidar measurements over the South China Sea, Bohai Sea, and Yellow SeaThe temporal‐spatial evolution characteristics of the turbulence mixing process are investigated under different weather conditions [ABSTRACT FROM AUTHOR]

The Guangdong‐Hong Kong‐Macao Greater Bay Area (GBA), a cluster of world‐class cities, is undergoing rapid urbanization. However, the heterogeneity of the urban thermal environment resulting from the diversity of urban forms is not yet fully understood. This paper assesses the heterogeneity of the urban heat island (UHI) effect in the GBA using the coupled Weather Research and Forecasting (WRF) model/multi‐layer urban canopy and building energy model (BEP/BEM), with high‐resolution local climate zone (LCZ) map as urban land use/land cover data. The average UHI intensity is found to peak at 1.8 ± 0.4°C in the evening, when the average UHI intensity of LCZ 2 can reach a maximum of 2.4 ± 0.58°C. Properly setting air‐conditioning temperatures can effectively prevent the enhancement of the UHI phenomenon at night by the anthropogenic heat (AH) released from air‐conditioning. The UHI‐induced local circulations and enhanced surface roughness inhibit the penetration of sea breezes inland, and surface wind speed decreases in all LCZs, with a maximum change of more than 0.8 m s−1. However, the increased thermal difference between land and sea leads to enhanced sea breezes offshore, especially in the Pearl River estuary. In addition, a series of sensitivity experiments have been conducted in this paper on initial and boundary conditions, building drag coefficients and urban fractions, which paves the way for further analyzing urban climate in GBA using WRF model and LCZs. Plain Language Summary: With the rapid urbanization of the world, the demand for functional buildings has increased. Along with the diversification of urban forms, the differences in the thermal environment within cities are becoming more and more significant. This study therefore provides an in‐depth study of the urban thermal environment in the Guangdong‐Hong Kong‐Macao Greater Bay Area (GBA) based on numerical simulation and local climate zones (LCZs). It was found that the urban heat island (UHI) intensity in different urban forms has obvious differences, and may vary by 1°C. However, the daily variation trends are similar, all showing a stronger UHI intensity at night than during the day, and reasonable setting of air‐conditioning temperature can effectively mitigate the UHI intensity at night. The UHI‐induced local circulations and enhanced surface roughness weaken the surface wind speed and inhibit the penetration of sea breeze inland, but enhance the sea breezes offshore, especially in the Pearl River estuary. This study provides references for urban planning and future sustainable development, especially for areas located along the coast that are undergoing rapid development. In addition, a series of sensitivity experiments on initial and boundary conditions, building drag coefficients and urban fractions provide useful suggestions for numerical model configuration in the GBA. Key Points: Tests of initial and boundary conditions, building drag coefficient and urban fractions provide recommendations for Weather Research and Forecasting configurationUrban heat island (UHI) varies between different local climate zones, but all peak in evening, and proper setting of air‐conditioning temperatures can mitigate UHI at nightUrbanization weakens surface wind speeds and inhibits the penetration of sea breezes inland, but strengthens the sea breezes offshore [ABSTRACT FROM AUTHOR]

E. Williams et al. (2022, commented paper) questioned electron energy spectra derived from thunderstorm ground enhancements (TGEs) measured on Aragats; they concluded that "A more likely origin for any detected electrons at 3.2 km above sea level is Compton scattering and pair production activated by longer‐range bremsstrahlung gamma rays, themselves produced by runaway electron encounters with nuclei in the breakeven field at higher altitude." In this comment, we show that the selection criteria of "electron" TGEs unambiguously reject the assumption of the origination of TGE electrons measured on Aragats from the Compton and pair production processes. Thus, the strong accelerating electric field above the earth's surface can be significantly lower (25–150 m) than derived in the commented paper 500 m altitude. Plain Language Summary: Electron accelerators operate in the thunderous atmosphere, sending copious particles to the Earth's surface. To get inside the models of electron acceleration and multiplication by strong atmospheric fields, the critical problem is the measurement of electrons and their energies as they arrive at the earth's surface. It is rather tricky because electrons are fast attenuated in the air, and the flux of accompanied gamma rays is attenuated much less and reaches the ground in overwhelming amounts. We developed special hardware and software methods to prove electrons' existence in the vast particle fluxes reaching the ground and to measure their energies. Simulations and careful examination of the registered particle fluxes check these methods. Key Points: The contribution of the Compton scattered and pair‐production electrons to TGE flux is negligible and cannot "mimic" the TGE electron fluxThe criteria used in the energy spectrum recovery from Aragats Solar Neutron Telescope (ASNT) reliably select "electron" TGE events and reject TGE events with small electron contentIf the strong accelerating electric field terminates low above the earth's surface (25–100 m), electrons from the large RREAs reach ASNT, and their energy spectrum can be reliably recovered [ABSTRACT FROM AUTHOR]