ALBERT

Description: Abstract Dust devils are convective vortices with a vertical axis of rotation made visible by lifted soil particles. Currently, there is great uncertainty about the extent to which dust devils contribute to the atmospheric aerosol input and thereby influence Earth's radiation budget. Past efforts to quantify the aerosol transport and study their formation, maintenance, and statistics using large‐eddy simulation (LES) have been of limited success. Therefore, some important features of dust devil‐like vortices simulated with LES still do not compare well with those of observed ones. One major difference is the simulated value of the core pressure drop, which is almost 1 order of magnitude smaller compared to the observed range of 250 to 450 Pa. However, most of the existing numerical simulations are based on highly idealized setups and coarse grid spacings. In this study, we investigate the effects of various factors on the simulated vortex strength with high‐resolution LES. For the fist time, we are able to reproduce observed core pressures by using a high spatial resolution of 2 m, a model setup with moderate background wind and a spatially heterogeneous surface heat flux. It is found that vortices mainly appear at the lines of horizontal flow convergence above the centers of the strongly heated patches, which is in contrast to some older observations in which vortices seemed to be created along the patch edges.

Description: Abstract Location information from super‐pressure balloons flown by Project Loon provide an unprecedented opportunity to analyze wind fields in the mid‐latitude stratosphere. Horizontal velocity spectra from the balloons' quasi‐intrinsic frame of reference show clear evidence of a persistent peak in the intrinsic wind spectrum around the inertial frequency. In the Southern Hemisphere mid‐latitudes, peak‐to‐peak amplitudes of horizontal velocity perturbations (on the order of 20 ms‐1) are larger than those seen in previous super‐pressure balloon campaigns in polar regions and similar to those observed in vertical soundings in the mid‐latitudes. A rotary spectral analysis shows that near‐circular anti‐cyclonic rotation of horizontal wind perturbations around the inertial frequency dominate at most times and locations. The strongest anti‐cyclonic rotation is more common in balloon flight segments with weak zonal winds and during the austral summer. Flight segments with strong eastward zonal velocities during austral winter and spring are more likely to have mixed cyclonic and anti‐cyclonic power around the inertial frequency. These results confirm previous model and radiosonde observations of the peak in horizontal kinetic energy at the inertial frequency, and demonstrate they are associated with increased anti‐cyclonic wave power indicative of near‐inertial oscillations or inertia‐gravity waves. Flight segments with mixed cyclonic and anti‐cyclonic power around the inertial frequency display a continuum of wave power from planetary to gravity‐wave scales. These results help explain the divergence of actual and modelled balloon trajectories in previous studies and provide a baseline against which reanalysis or meteorological model realizations of the intrinsic velocity field can be assessed.

Description: Abstract From 17–22 August 2017 simultaneous enhancements of ammonia (NH3), carbon monoxide (CO), hydrogen cyanide (HCN), and ethane (C2H6) were detected from ground‐based solar absorption Fourier transform infrared (FTIR) spectroscopic measurements at two high‐Arctic sites: Eureka (80.05°N, 86.42°W) Nunavut, Canada, and Thule (76.53°N, 68.74°W), Greenland. These enhancements were attributed to wildfires in British Columbia and the Northwest Territories of Canada using FLEXPART back‐trajectories and fire locations from Moderate Resolution Imaging Spectroradiometer (MODIS) and found to be the greatest observed enhancements in more than a decade of measurements at Eureka (2006–2017) and Thule (1999–2017). Observations of gas‐phase NH3 from these wildfires illustrate that boreal wildfires may be a considerable episodic source of NH3 in the summertime high Arctic. Comparisons of GEOS‐Chem model simulations using the Global Fire Assimilation System (GFASv1.2) biomass burning emissions to FTIR measurements and Infrared Atmospheric Sounding Interferometer (IASI) measurements showed that the transport of wildfire emissions to the Arctic was underestimated in GEOS‐Chem. However, GEOS‐Chem simulations showed that these wildfires contributed to surface layer NH3 and NH enhancements of 0.01–0.11 ppbv and 0.05–1.07 ppbv, respectively, over the Canadian Archipelago from 15–23 August 2017.

Description: Abstract Accurate estimates of NOx and SO2 emissions are important for air quality modeling and management. To incorporate chemical interactions of the two species in emission estimates, we develop a joint hybrid inversion framework to estimate their emissions in China and India (2005–2012). Pseudo observation tests and posterior evaluation with surface measurements demonstrate that joint assimilation of SO2 and NO2 can provide more accurate constraints on emissions than single‐species inversions. This occurs through synergistic change of O3 and OH concentrations, particularly in conditions where satellite retrievals of the species being optimized have large uncertainties. The percentage changes of joint posterior emissions from the single‐species posterior emissions go up to 242% at grid scales, although the national average of monthly emissions, seasonality, and interannual variations are similar. In China and India, the annual budget of joint posterior SO2 emissions is lower, but joint NOx posterior emissions are higher, because NOx emissions increase to increase SO2 concentration and better match Ozone Monitoring Instrument SO2 observations in high‐NOx regions. Joint SO2 posterior emissions decrease by 16.5% from 2008 to 2012, while NOx posterior emissions increase by 24.9% from 2005 to 2011 in China—trends which are consistent with the MEIC inventory. Joint NOx and SO2 posterior emissions in India increase by 15.9% and 19.2% from 2005 to 2012, smaller than the 59.9% and 76.2% growth rate using anthropogenic emissions from EDGARv4.3.2. This work shows the benefit and limitation of joint assimilation in emission estimates and provides an efficient framework to perform the inversion.

Description: Abstract By using 90 radiosonde stations with high vertical resolution data during the period 1998–2011, the latitudinal variation of the tropopause inversion layer (TIL) in different seasons and the interactions with the inertial gravity wave (IGW) activities in the region covering the Northern Hemispheric latitudes from 5° to 75° are studied. For the midlatitudes, the TIL features show obviously seasonal variations. In the Arctic region, TIL is strong and thick. The averaged Arctic TIL intensity peaks in summer. The intense interaction between the TIL and IGW is found in the region of 5°N to 75°N. The TIL could inhibit the upward propagation of IGWs from ~2 km below the tropopause in a larger region (40–75°N). It is found that for the middle‐latitude regions, the enhanced wind shear layer just above the tropopause could lead to instability and finally result in IGW breaking and intensive turbulence, which then leads to strong wave energy dissipation and a downward heat flux. The IGW‐induced cooling around the tropopause, which resulted from the downward heat flux, then makes a colder and sharper tropopause and finally form the TIL. The IGW‐associated strong downward heat flux is also found around the Arctic tropopause. However, there is no corresponding wind shear enhancement above the tropopause. This indicates that this strong heat flux may result from some other processes and then form the strong TIL in the Arctic.

Description: Abstract Proxy system models (PSMs) are an important bridge between climate simulations and climate records prior to the period where instrumental observations are available. PSMs help to interpret what proxies show and how they record climate. Although previous studies have evaluated PSMs for individual sites, their systematic evaluation on a global scale has not yet been conducted. This study evaluated the performance of PSMs for stable water isotopes in ice cores, corals, and tree‐ring cellulose for the period 1950–2007. Spatial distributions of the mean state were well simulated for all proxy types, albeit with a bias for tree‐ring cellulose. Interannual variability was well simulated for corals and tree‐ring cellulose. These results indicate that the models represent key mechanisms for the proxies. In contrast, the reproducibility of interannual variability in ice cores was markedly lower than that for the other proxies. Although the reproducibility was limited by the atmospheric forcing used to drive the model, the results suggest that the PSM may be missing post‐depositional processes, such as sublimation for ice cores on the interannual timescale.

Description: Abstract Three years of nighttime Cloud‐Aerosol Lidar and Infrared Pathfinder Satellite Observation data was used in synergy with CloudSat measurements to quantify how strongly aerosol type and aerosol load affect the cloud phase in low‐level clouds over the Arctic. Supercooled liquid layers were present in the majority of observed low‐level clouds (0.75 ≤ z ≤ 3.5 km) between −10 and −25 °C. Furthermore, based on the subset (6%) of data with high quality assurance for aerosol typing, ice formation is more common in the presence of dust or continental aerosols as opposed to marine or elevated smoke aerosols. With the first aerosol group, glaciated clouds were found at cloud top temperatures of 2 to 4 °C warmer than with the latter aerosol types. Further association of the aerosol concentration with the cloud phase showed that the aerosol concentration outweighs the aerosol type effect. Depending on the aerosol load, the temperature at which a cloud completely glaciates can vary by up to 6–10 °C. However, this behavior was most pronounced in stable atmospheric conditions and absent over open ocean with lower tropospheric stability values and probably less stratified clouds. Finally, more mixed‐phase clouds were associated with high aerosol load, suggesting that mixed‐phase clouds have an extended lifetime in the Arctic under high cloud condensation nuclei concentrations. This implies that in a pristine environment, with few or no local aerosol sources, and under the investigated conditions the amount of aerosol particles affects the cloud phase more than the aerosol type does.

Description: Abstract Aerosol pH is a useful diagnostic of aerosol chemistry for formation of secondary aerosol and has been hypothesized to be a key factor in specific chemical reaction routes producing sulfate and nitrate (Yue et al., 2009; Zhang et al., 2012; Hu et al., 2014). In this study, we measured hourly concentrations of water soluble ions (WS‐ions) in PM2.5, along with gaseous pollutants in Tianjin, China, from 4th January to 31st January 2015. The following source contributions to WS ions were estimated by PMF (Positive Matrix Factorization): secondary sulfate (13%), secondary nitrate (44%), coal (14%), vehicle (16%), and dust (13%). ISORROPIA‐II was used to investigate the complex relationships among aerosol pH, ammonia, and secondary aerosol formation. The estimated hourly aerosol pH varied from ‐0.3 to 7.7, with an average of 4.9 (±0.78); the median value was 4.89, and the interquartile range (IQR) was 0.72. During less polluted conditions, aerosol pH ranged from less than 0 to about 7; during heavily polluted conditions, pH was close to 5 (3.9‐7.9) despite large amounts of sulfate. Sufficient ammonia/ammonium were present to balance high sulfate and nitrate formation. NH4+/NH3 (g) helped stabilize pH while nonvolatile cations contributed less to decreasing aerosol acidity. High acidy (pH〈3), light pollution (Total water soluble ions (TWI)〈30 μg m‐3), and low water content (less than 5μg m‐3) were more correlated with higher rates of sulfate formation than nitrate formation in the winter.

Description: Abstract Understanding the mechanisms by which earthquake cycles produce folding and accommodate shortening is essential to quantify the seismic potential of active faults and integrate aseismic slip within our understanding of the physical mechanisms of the long‐term deformation. However, measuring such small deformation signals in mountainous areas is challenging with current space‐geodesy techniques, due to the low rates of motion relative to the amplitude of the noise. Here we successfully carry out a multitemporal Interferometric Synthetic Aperture Radar analysis over the North Qaidam fold‐thrust system in NE Tibet, where eight Mw〉 5.2 earthquakes occurred between 2003 and 2009. We report various cases of aseismic slip uplifting the thickened crust at short wavelengths. We provide a rare example of a steep, shallow, 13‐km‐long and 6‐km‐wide afterslip signal that coincides spatially with an anticline and that continues into 2011 in response to a Mw 6.3 event in 2003. We suggest that a buried seismic slip during the 2003 earthquake has triggered both plastic an‐elastic folding and aseismic slip on the shallow thrusts. We produce a first‐order two‐dimensional model of the postseismic surface displacements due to the 2003 earthquake and highlight a segmented slip on three fault patches that steepen approaching the surface. This study emphasizes the fundamental role of shallow aseismic slip in the long‐term and permanent deformation of thrusts and folds and the potential of Interferometric Synthetic Aperture Radar for detecting and characterizing the spatiotemporal behavior of aseismic slip over large mountainous regions.

Description: Abstract Structural details of the crust play an important role in controlling the distribution of volcanic activity in arc systems. In southwest Washington, several different regional structures associated with accretion and magmatism have been invoked to explain the broad distribution of Cascade volcanism in this region. In order to image these regional structures in the upper crust, Pg and Sg travel times from the imaging Magma Under St. Helens (iMUSH) active‐source seismic experiment are inverted for Vp, Vs, and Vp/Vs models in the region surrounding Mount St. Helens. Several features of these models provide new insights into the regional structure of the upper crust. A large section of the Southern Washington Cascades Conductor is imaged as a low Vp/Vs anomaly that is inferred to represent a broad sedimentary/metasedimentary sequence that composes the upper crust in this region. The accreted terrane Siletzia is imaged west of Mount St. Helens as north/south trending high Vp and Vp/Vs bodies. The Vp/Vs model shows relatively high Vp/Vs regions near Mount St. Helens and the Indian Heaven Volcanic Field, which could be related to the presence of magmatic fluids. Separating these two volcanic regions below 6‐km depth is a northeast trending series of high Vp and Vs bodies. These bodies have the same orientation as several volcanic/magmatic features at the surface, including Mount St. Helens and Mount Rainier, and it is argued that these high‐velocity features are a regional‐scale group of intrusive bodies associated with a crustal weak zone that focuses magma ascent.

Description: Abstract SO2 column densities from OMI provide important information on emission trends and missing sources, but there are discrepancies between different retrieval products. We employ three OMI SO2 retrieval products (NASA standard (SP), NASA prototype, and BIRA) to study the magnitude and trend of SO2 emissions. SO2 column densities from these retrievals are most consistent when viewing angles and solar zenith angles are small, suggesting more robust emission estimates in summer and at low latitudes. We then apply a hybrid 4D‐Var/mass balance emission inversion to derive monthly SO2 emissions from the NASA SP and BIRA products. Compared to HTAPv2 emissions in 2010, both posterior emission estimates are lower in US, India and Southeast China, but show different changes of emissions in North China Plain. The discrepancies between monthly NASA and BIRA posterior emissions in 2010 are less than or equal to 17% in China and 34% in India. SO2 emissions increase from 2005 to 2016 by 35% (NASA) ‐ 48% (BIRA) in India, but decrease in China by 23% (NASA) ‐ 33% (BIRA) since 2008. Compared to in‐situ measurements, the posterior GEOS‐Chem surface SO2 concentrations have reduced NMB in China, the US, and India but not in South Korea in 2010. BIRA posteriors have better consistency with the annual growth rate of surface SO2 measurement in China and spatial variability of SO2 concentration in China, South Korea and India, whereas NASA SP posteriors have better seasonality. These evaluations demonstrate the capability to recover SO2 emissions using OMI observations.

Description: Abstract Seismicity of several intraplate seismic zones in the North American midcontinent is believed to be related to reactivation of ancient faults in Precambrian continental rifts by the contemporary stress field. Existence of such a rift system beneath the Wabash Valley Seismic Zone (WVSZ) is not clear. Here we obtained a crustal structural image along a 300‐km‐long profile across WVSZ using a dense linear seismic array. We first calculated teleseismic receiver functions of stations and applied the Common‐Conversion‐Point stacking method to image crustal interfaces and the Moho. We then used ambient noise cross correlation to obtain phase and group velocities of Rayleigh and Love waves. Finally, we jointly inverted the receiver function and surface wave dispersion data to determine shear wave velocity structure along the profile. The results show a thick (50‐ to 60‐km) crust with a typical Proterozoic crustal layering: a 1‐ to 2‐km thick Phanerozoic sedimentary layer, an upper crust ∼15 km thick, and a 30‐ to 40‐km‐thick lower crust. The unprecedented high‐resolution image also reveals a 50‐km‐wide high‐velocity body above an uplifted Moho and several velocity anomalies in the upper and middle crust beneath the La Salle Deformation Belt. We interpreted them as features produced by magmatic intrusions in a failed, immature continental rift during the end of Precambrian. Current seismicity in WVSZ is likely due to reactivation of ancient faults of the rift system by a combination of stress fields from the far‐field plate motion and prominent crustal and upper mantle heterogeneities in the region.

Description: Abstract The Charlevoix Seismic Zone (CSZ) is located along the early Paleozoic St. Lawrence rift zone in southeastern Quebec at the location of a major Devonian impact structure. The impact structure superimposed major, steeply dipping basement faults trending approximately N35°E. Approximately 250 earthquakes are recorded each year and are concentrated within and beneath the impact structure. Most M4+ earthquakes associated with the rift faults occurred outside the impact structure. Apart from the unique distribution of earthquakes, stress inversion of focal mechanisms shows stress rotations within the CSZ, and in the CSZ relative to the stress orientation determined from borehole breakouts. The primary goal of this research is to investigate the combined effects of the preexisting structures and regional stresses on earthquake activity and stress rotations in the CSZ. We approach this using PyLith, a finite‐element code for simulations of crustal deformation. Adopting the results from recent hypocenter relocation and 3‐D tomography studies, we modify the locations and dips of the rift faults and assess the effect of the new fault geometries on stress distributions. We also discuss the effects of resolved velocity anomalies. We find that the observed stress rotation is due to the combined effect of the rift faults and the impact structure. One‐dimensional velocity models of the CSZ with an embedded impact structure and a combination of 65°‐40°‐40° and constant 70° fault dip models with a very low friction coefficient of 0.3 and cohesion of 0 MPa can explain the observed seismicity and more than 50% of the stress rotations.

Description: Abstract Seismic anisotropy provides important information on the structure and geodynamics of the Earth. The forearc mantle wedge in subduction zones mainly exhibits trench‐parallel azimuthal anisotropy globally, which is inconsistent with the model of olivine a axis aligning with the slab‐driven corner flow. Its formation mechanism is currently unclear. Here we present high‐resolution 3‐D P wave anisotropic tomography of the Tohoku subduction zone. We suggest that ductile deformation of the forearc lithospheric mantle of the overriding plate induces the trench‐parallel azimuthal anisotropy and positive radial anisotropy (i.e., horizontal velocity 〉 vertical velocity) in Tohoku. Our results provide the first seismic anisotropic evidence for the slab‐mantle decoupling at a common depth of ~70 km. On the basis of the high‐resolution seismic images, we propose a geodynamic model suggesting that the forearc mantle wedge anisotropy is produced via ductile deformation of dry olivine or hydrous antigorite lithospheric mantle, which accords well with the trench‐parallel shear wave splitting measurements dominant in subduction zones globally.

Description: Abstract We investigate 3‐D seismic structures (Vp, Vs, and Poisson's ratio) and Vp azimuthal anisotropy in the source area of the 2018 Eastern Iburi earthquake (M 6.7) in Hokkaido, Japan. Its mainshock occurred at the edge of a high‐Vp (2–4%) seismogenic zone. Significant low‐Vs (−1% to −3%) and high Poisson's ratio (2–7%) anomalies are imaged in and below the source zone and extend to the upper surface of the subducting Pacific slab, most likely reflecting ascending fluids released by the slab dehydration. A high consistency between the fault plane and the low‐Vs and high Poisson's ratio anomalies indicates that the fluids may have entered the fault and affected the rupture nucleation. A high‐V (1–3%) anomaly is revealed in the fore‐arc mantle wedge and connects with the high‐V seismogenic zone, probably reflecting a lithospheric fragment and contributing to cool down the mantle wedge. Complex seismic anisotropy is revealed in the crust in and around the source area, which may reflect complicated stress regime and strong structural heterogeneities there.

Description: Abstract Deep convective clouds similar to those arising in the TC eyewall are simulated using a parcel model and 2D slab symmetric cloud model with spectral bin microphysics (the Hebrew University Cloud Model, HUCM). The size distribution of sea spray particles (SSP) at cloud base is calculated using the Lagrangian‐Eulerian bin‐microphysics model (LEM). The model describes the SSP production, advection and formation of the size distribution of SSP in the hurricane atmospheric boundary layer at different strong wind speeds. The SSP distributions calculated by the LEM are used in the parcel model and the HUCM to investigate the microphysical and dynamical effects of SSP on clouds. The SSP ascending in cloud updrafts dramatically increase the number concentration of cloud drops within a wide range of drop sizes. As a result, sea spray creates clouds with unique property combinations of both maritime and continental types. These clouds have droplet size distributions characterized by a high drop concentration and a low effective radius, as in continental clouds. At the same time, the presence of SSP of a few hundred microns in radii triggers intense rain just above the cloud base, which is typical of extreme maritime clouds. In the presence of large sea spray drops, the smallest cloud condensational nuclei, including the smallest SSP, are activated, giving rise to the permanent in‐cloud nucleation of small droplets which produce a high concentration of small ice crystals above the level of homogeneous freezing. We showed that the SSP substantially increased the maximum vertical velocity, cloud water content and mass contents of ice particles. The results are compared with available observed data.

Description: Abstract Terrestrial gamma ray flashes (TGFs) are very short bursts of gamma radiation associated to thunderstorm activity and are the manifestation of the highest‐energy natural particle acceleration phenomena occurring on Earth. Photon energies up to several tens of megaelectronvolts are expected, but the actual upper limit and high‐energy spectral shape are still open questions. Results published in 2011 by the AGILE team proposed a high‐energy component in TGF spectra extended up to ≈100 MeV, which is difficult to reconcile with the predictions from the Relativistic Runaway Electron Avalanche (RREA) mechanism at the basis of many TGF production models. Here we present a new set of TGFs detected by the AGILE satellite and associated to lightning measurements capable to solve this controversy. Detailed end‐to‐end Monte Carlo simulations and an improved understanding of the instrument performance under high‐flux conditions show that it is possible to explain the observed high‐energy counts by a standard RREA spectrum at the source, provided that the TGF is sufficiently bright and short. We investigate the possibility that single high‐energy counts may be the signature of a fine‐pulsed time structure of TGFs on time scales ≈4 μs, but we find no clear evidence for this. The presented data set and modeling results allow also for explaining the observed TGF distribution in the (Fluence × duration) parameter space and suggest that the AGILE TGF detection rate can almost be doubled.

Description: Abstract The processes that control the isotopic composition of precipitation in the mid‐latitudes are complicated, but can provide valuable insights into precipitation‐generating processes and are critical for interpreting stable isotope‐based paleoclimate records. In this study, we investigated the controls on changes in the isotopic composition of rainwater in central Texas using a combination of existing monthly stable isotope data from the global network of isotopes in precipitation (GNIP) and 20 months of event‐based rainwater collection from Austin, TX. We find that the strongest control on the isotopic composition of precipitation is the varying proportion of convective and stratiform rainfall, with other factors such as precipitation amount, temperature, storm track playing a secondary role. Isotopic values are generally lower in the cold season than the warm season precipitation because cold season precipitation is predominantly stratiform often associated with a northerly storm track. However, the majority of the precipitation in the south‐central United States (US) occurs during the warm season in association with mesoscale convective systems (MCS) that are fed with moisture by the southerly winds. MCS are characterized by a combination of a leading edge of organized deep convection and trailing stratiform precipitation. Stronger MCS tend to contain higher proportions of stratiform rainfall and as a result, have more isotopically depleted values. Therefore, changes in the stable isotopic composition of rainfall may be interpreted as reflecting changes in the intensity of MCS.

Description: Abstract Atmospheric nitrate (NO3− = particulate NO3− + gas‐phase nitric acid [HNO3]) and sulfate (SO42−) are key molecules that play important roles in numerous atmospheric processes. Here, the seasonal cycles of NO3− and total suspended particulate sulfate (SO42−(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 NO3− was largely affected by snowpack emissions in which [NO3−] and δ15N(NO3−) were highest (49.3 ± 21.4 ng/m3, n = 8) and lowest (−47.0 ± 11.7‰, n = 5), respectively, during periods of sunlight in the interior of Antarctica. The seasonal cycle of Δ17O(NO3−) reflected tropospheric chemistry year‐round with lower values observed during sunlight periods and higher values observed during dark periods, reflecting shifts from HOx‐ to O3‐dominated oxidation chemistry. SO42−(TSP) concentrations were highest during austral summer and fall (86.7 ± 73.7 ng/m3, n = 18) and are indicated to be derived from dimethyl sulfide (DMS) emissions, as δ34S(SO42−)(TSP) values (18.5 ± 1.0‰, n = 10) were similar to literature δ34S(DMS) values. The seasonal cycle of Δ17O(SO42−)(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 HOx‐ to O3‐dominated chemistry in the atmospheric derived SO42− component. Overall, the budgets of NO3− and SO42−(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.

Description: Abstract Understanding global soil moisture‐air temperature (θ‐Ta) coupling is needed to improve the representation of land‐atmosphere interactions in Earth system models. Most studies on θ‐Ta coupling have focused on hot extremes, where precipitation‐related indices and model‐derived soil moisture products are commonly used. In this study, global θ‐Ta coupling is examined based on monthly air temperature anomalies and the Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage (TWS). A discrete wavelet decomposition is used to partition the TWS into different components. The results show that TWS is useful in revealing the spatial patterns of θ‐Ta coupling. Decomposed GRACE TWS shows improved skill compared to raw TWS in explaining temporal variability of monthly air temperature, which likely reflects different roles of soil moisture at different depths in the θ‐Ta coupling. The explanatory power improves further by using a combination of decomposed GRACE TWS and precipitation. Such improvement is observed particularly in places where vegetation tends to have a deeper rooting system, such as eastern region of South America, the southern tip of Africa, and north of the Tropic of Capricorn in Australia. The occurrence of θ‐Ta coupling is mainly constrained by the coupling of root zone moisture and land surface temperature. In addition to deeper rooting systems, clear wet and dry season alternation is another favorable factor for developing significant monthly θ‐Ta coupling.

Description: Abstract This study investigates heat wave variability over Korea during 1979‐2017. It is found that most of heat waves in Korea can be classified into two distinct types based on the spatial patterns of atmospheric circulation anomalies: the zonal wave (Z‐wave) type and the meridional wave (M‐wave) type. The Z‐wave type is accompanied by large‐scale atmospheric waves across the Eurasian continent, while the M‐wave type is associated with convective activities over the subtropical western North Pacific. The Z‐wave type occurs when the high‐pressure node of eastward propagating wave located around Korea and it seems that the associated wave energy could originate from North Atlantic Ocean. The M‐wave type, on the other hand, is driven by northward propagating wave train from subtropical western North Pacific to East Asia, which is triggered by anomalous convective activity over the subtropical western North Pacific. By analyzing thermodynamical as well as dynamical variables, detailed descriptions on the physical characteristics of two types of heat wave are presented in this study along with the possible implications for summer climate variability over Korea.

Description: Abstract Recent laboratory evidence shows that compaction creep in porous rocks may develop through stages of acceleration, especially if the material is susceptible to strain localization. This paper provides a mechanical interpretation of compaction creep based on viscoplasticity and nonlinear dynamics. For this purpose, a constitutive operator describing the evolution of compaction creep is defined to evaluate the spontaneous accumulation of pore collapse within an active compaction band. This strategy enables the determination of eigenvalues associated with the stability of the response, i.e. able to differentiate decelerating from accelerating strain. This mathematical formalism was linked to a constitutive law able to simulate compaction localization. Material point simulations were then used to identify the region of the stress space where unstable compaction creep is expected, showing that accelerating strains correspond to pulses of inelastic strain rate. Such pulses were also found in full‐field numerical analyses of delayed compaction, revealing that they correspond to stages of inception and propagation of new bands across the volume of the simulated sample. These results illustrate the intimate relation between the spatial patterns of compaction and their temporal dynamics, showing that while homogeneous compaction develops with decaying rates of accumulation, localized compaction occurs through stages of accelerating deformation caused by the loss of strength taking place during the formation of a band. In addition, they provide a predictive modeling framework to simulate and explain the spatiotemporal dynamics of compaction in porous sedimentary formations.

Description: ABSTRACT Detailed P wave velocity and anisotropy structure of the uppermost mantle below the central United States is presented based on a tomographic inversion of Pn traveltimes for earthquakes in the range 2 to 14°. Dense raypath coverage throughout the northern Mississippi Embayment is obtained using the Northern Embayment Lithosphere Experiment and U.S. Transportable Array data sets. A detailed analysis of the trade‐off between velocity and anisotropy variations demonstrates that both are well resolved over most of the study area. Anomalously fast Pn velocities are identified below the northern Mississippi Embayment, centered on the New Madrid seismic zone. A prominent region of low velocity coincides with the southwestern margin of the Illinois basin. Pn anisotropy displays complex patterns and differs from absolute plate motion directions and SKS splitting directions. A circular pattern of fast anisotropy directions is centered on the New Madrid seismic zone and may be related to the presence of the mafic “rift pillow.”

Description: Abstract Sedimentary relative paleointensity (RPI) records are often carried by complex magnetic mineral mixtures, including detrital and biogenic magnetic minerals. Recent studies have demonstrated that magnetic inclusions within larger detrital silicate particles can make significant contributions to sedimentary paleomagnetic records. However, little is known about the role such inclusions play in sedimentary paleomagnetic signal recording. We analyzed paleomagnetic and mineral magnetic data for marine sediment core MD01‐2421 from the North Pacific Ocean, offshore of central Japan, to assess how magnetic inclusions and other detrital magnetic minerals record sedimentary paleomagnetic signals. Stratigraphic intervals in which abundant magnetic inclusions dominate the magnetic signal are compared with other intervals to assess quantitatively their contribution to sedimentary RPI signals. The normalized remanence record from core MD01‐2421 does not correlate clearly with global RPI stacks, which we attribute to a demonstrated lower paleomagnetic recording efficiency of magnetic inclusions compared to other detrital magnetic minerals. We also carried out the first laboratory redeposition experiments under controlled Earth‐like magnetic fields for particles with magnetic inclusions using material from core MD01‐2421. Our results confirm that such particles can be aligned by ambient magnetic fields but with a lower magnetic recording efficiency compared to other detrital magnetic minerals, which is consistent with normalized remanence data from core MD01‐2421. Our demonstration of the role of sedimentary magnetic inclusions should have wide applicability for understanding sedimentary paleomagnetic recording.

Description: Abstract In this study, we estimate atmospheric turbulence in the free atmosphere in terms of the Thorpe scale (LT) and eddy dissipation rate (ε) using U.S. high vertical‐resolution radiosonde data over 4 years (September 2012 to August 2016) at 68 operational stations. In addition, same calculations are conducted for 12 years (October 2005 to September 2017) at four stations among the 68 stations. These high vertical‐resolution radiosonde data have a vertical resolution of approximately 5 m and extend to an altitude of approximately 33 km, and thus, turbulence can be retrieved in the entire troposphere and lower stratosphere. There are thicker and stronger turbulent layers in the troposphere than in the stratosphere, with mean ε values of 1.84 × 10−4 and 1.37 × 10−4 m2/s3 in the troposphere and stratosphere, respectively. The vertical structure of ε exhibits strong seasonal variations, especially in the upper troposphere and lower stratosphere, with the largest ε values in summer and the smallest in winter. In the horizontal distribution of ε, large ε is seen mainly above the mountainous region in the troposphere, but this pattern is not seen in the stratosphere. Although ε is estimated by the square of LT multiplied by the cube of the Brunt‐Väisälä frequency (N), the regions of large ε are matched with large LT regions where N is relatively small. For the time series of ε near the tropopause for 12 years at four stations, an annual variation is prominent at all stations without significant interannual variations. There is, however, a slightly increasing trend of ε at two stations.

Description: Abstract The mechanical dynamics of volcanic systems can be better understood with detailed knowledge on strength of a volcanic edifice and subsurface. Previous work highlighting this on Mt. Etna has suggested that its carbonate basement could be a significant zone of widespread planar weakness. Here, we report new deformation experiments to better quantify such effects. We measure and compare key deformation parameters using Etna basalt (EB), which is representative of upper edifice lava flows, and Comiso limestone (CL), which is representative of the carbonate basement, under upper crustal conditions. These data are then used to derive empirical constitutive equations describing changes in rocks strength with pressure, temperature and strain rate. At a constant strain rate of 10‐5 s‐1 and an applied confining pressure of 50 MPa the brittle to ductile transitions were observed at 975 °C (EB) and 350 °C (CL). For the basaltic edifice of Mt. Etna, the strength is described with a Mohr‐coulomb failure criterion with μ ~0.704, C = 20 MPa. For the carbonate basement, strength is best described by a power law‐type flow in two regimes: a low‐T regime with stress exponent n ~5.4 and an activation energy Q ~ 170.6 kJ/mol and a high‐T regime with n~ 2.4 and Q ~ 293.4 kJ/mol. We show that extrapolation of these data to Etna's basement predicts a brittle to ductile transition that corresponds well with the generally observed trends of the seismogenic zone underneath Mt. Etna. This in turn may be useful for future numerical simulations of volcano‐tectonic deformation of Mt. Etna, and other volcanoes with limestone basements.

Description: Abstract On 20 April 2013, an Mw 6.6 Lushan earthquake occurred on the southwestern segment of the Longmen Shan fault belt, which is the tectonic block boundary between the eastern Tibetan plateau and the Sichuan basin. Seismic reflection profiles and aftershock relocation indicate that there exists a back thrust fault in the source region but whether it is ruptured during the Lushan earthquake remains controversial. Here the precise leveling data are firstly used together with Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), and strong motion data to invert for the fault geometry and slip distribution associated with the earthquake. The joint inversion result shows that the Lushan earthquake occurred on a blind thrust fault with strike N208.5 °E and dip 42.1° to the NW and did not rupture the back reverse fault. The coseismic slip model reveals the Lushan earthquake involves the rupture of one major asperity. The coseismic slip is mainly concentrated on a steeply dipping fault plane. The coseismic rupture terminates on the southwestern side of the seismic gap between the Wenchuan and Lushan earthquakes. Topographic stress may be the dominant mechanism of coseismic rupture termination.

Description: Abstract Experimental data show that inelastic straining occurs even at very low pressure before and during “brittle” fracturing. This process is therefore investigated within the framework of elastoplasticity using 2D, 3‐layer FD modeling. The constitutive model includes both tensile and shear failure mechanisms coupled at the level of the strain softening law. The modeling results show that sets of parallel joints initiate as pure dilation bands, the narrow σ3‐normal bands of localized dilatant damage (inelastic deformation). The band thickness, length, and the initial strain softening degree within it are proportional to the ductility of the material, which increases with the effective stress level (σ1) or pressure. The strength reduction within the bands is accelerated at a certain stage, and the strength locally reaches zero resulting in fracture initiation. The initial fracture then propagates in mode I following the propagating band. The fracture (joint) appears thus as a band of damaged material with the increased porosity, which is maximum along the axial zone of the band where the material is completely broken. The damage is due to both tensile and shear mechanisms. The role of shear failure increases with the ductility (pressure) increase, which also leads to the band thickness increase. These processes can result in small (band thickness)‐scale shear fractures within the band, causing the increase in the roughness of fracture walls organized in plumose patterns typical of both natural and experimentally generated joints.

Description: Abstract The second Radiative Heating in Underexplored Bands Campaign (RHUBC‐II) was conducted in 2009 by the U.S. Department of Energy Atmospheric Radiation Measurement program to improve water vapor spectroscopy in the far‐infrared spectral region. RHUBC‐II was located in an extremely dry region of Chile to ensure very low opacities in this spectral region. Spectrally resolved measurements by a far‐infrared spectrometer and a submillimeter interferometer from RHUBC‐II are compared with line‐by‐line radiative transfer model calculations. Water vapor amounts and temperatures used in the calculations come from collocated radiosondes, with extensive adjustments to correct for issues due to the campaign's dry conditions and mountainous terrain. A reanalysis is also performed of far‐infrared measurements taken at the Atmospheric Radiation Measurement North Slope of Alaska site before and during the first RHUBC campaign. These analyses determine that differences between the measurements and model calculations using existing spectroscopic parameters are significant in the far‐infrared and submillimeter regions, leading to the derivation of improved water vapor continuum absorption coefficients and air‐broadened widths of 74 water vapor lines. The foreign continuum is increased by more than 50% in part of the far‐infrared and the widths of more than 20 lines are changed by more than 10%. The uncertainty in the foreign continuum coefficients is estimated as greater than 20% in some spectral regions, primarily a consequence of the uncertainty in the specification of water vapor. The improved far‐infrared spectroscopic parameters have a notable impact on calculated spectral radiances and a modest impact on broadband radiative fluxes and heating rates.

Description: Abstract Deciphering the relationship between lateral growth of faults and along‐strike deformation (i.e., shortening and uplift) in the Earth's upper crust remains a challenge. Here we gain insight into the relation between these processes by studying the Kashi anticline, an asymmetric, doubly plunging thrust‐fault‐related fold located in the southwest Tian Shan, China. We use seismic interpretation and field observations, together with 2‐D trishear and excess area methods, to quantify the distribution of shortening along this structure. The shortening distribution along strike of the Kashi anticline is nonlinear and has a peaked, asymmetric, bell shape, with a maximum value of 5.9 ± 0.2 km. After comparing the 3‐D structural model of the Kashi anticline and our trishear models, we propose that lateral propagation‐to‐maximum shortening ratio, initiation fault length, and lateral propagation rate control the lateral fault propagation process and the fold terminations. Moreover, the 3‐D fault morphology and the ages of the growth strata suggest that the Kashi anticline experienced two stages of lateral growth with propagation rates of 60 km/Ma between 1.4 ± 0.2 Myr and 0.9 ± 0.3 Ma, and ~67 km/Myr from 0.9 ± 0.3 Ma to present. These observations highlight the relation between the evolution of lateral fault growth and the along‐strike shortening distribution, allowing us to use the latter (which we can measure) to infer the former (which we cannot). These novel insights from the Kashi anticline can be used to understand lateral growth of thrust and normal faults worldwide.

Description: Abstract Seismically detected ultralow velocity zones (ULVZs) at the the core‐mantle boundary (CMB) reflect the dynamical state and geological evolution of the silicate‐metal frontier of Earth's deep interior. However, modeling the dynamical context of ULVZs is hampered by challenges, such as the necessity of fine scale resolution and the accurate treatment of large viscosity contrasts. Here we extend the treatment of ULVZs using a lubrication theory approach and apply it to numerical and analytical models relevant for mantle convection in the CMB region. A generic model of a thin and dense low viscosity ULVZ layer embedded between an overlying convecting viscous mantle and an underlying inviscid core can explain several features that are consistent with seismic inferences, such as the absence of ULVZs in some regions and a tabular shape where they are concentrated. The model explains how the topography of a ULVZ layer tends to saturate and flatten as it becomes thicker, due to a non‐linear feedback between viscous aggregation beneath upwelling mantle currents and gravitational spreading/relaxation. Implementation of the ULVZ equation in thermal convection models indicates that ULVZs are preferentially gathered beneath long‐lived plumes, and may not exist beneath newly formed plume roots where there is no source of layer material. The presence/absence of ULVZs and their detailed shapes may provide important insights into the dynamical state and convective instability of the lowermost mantle thermal boundary layer.

Description: Abstract We implement a Coulomb rate‐and‐state approach to explore the nonlinear relation between stressing rate and seismicity rate in the Groningen gas field. Coulomb stress rates are calculated, taking into account the 3‐D structural complexity of the field and including the poroelastic effect of the differential compaction due to fault offsets. The spatiotemporal evolution of the Groningen seismicity must be attributed to a combination of both (i) spatial variability in the induced stressing rate history and (ii) spatial heterogeneities in the rate‐and‐state model parameters. Focusing on two subareas of the Groningen field where the observed event rates are very contrasted even though the modeled seismicity rates are of similar magnitudes, we show that the rate‐and‐state model parameters are spatially heterogeneous. For these two subareas, the very low background seismicity rate of the Groningen gas field can explain the long delay in the seismicity response relative to the onset of reservoir depletion. The characteristic periods of stress perturbations, due to gas production fluctuations, are much shorter than the inferred intrinsic time delay of the earthquake nucleation process. In this regime the modeled seismicity rate is in phase with the stress changes. However, since the start of production and for two subareas of our analysis, the Groningen fault system is unsteady and it is gradually becoming more sensitive to the stressing rate.

Description: Abstract From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T‐3). The resulting measurements, which are unevenly distributed on Alpha‐Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T‐3 heat flow is ~54.7 ± 11.3 mW m‐2, and Nautilus Basin is the only well‐surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m‐2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near‐linear, implying that conductive heat transport dominates and that near‐seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A numerical model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T‐3 dataset.

Description: Abstract The Early Cretaceous Ontong Java Plateau (OJP) in the southwestern Pacific Ocean is the largest oceanic plateau by volume on Earth, and a broad range of observations has been conducted to reveal its formation and evolution. However, because seafloor seismic observations of the OJP and surrounding areas have been insufficient so far, such experiments are capable of generating additional information regarding the crustal and mantle structure of the OJP. To image seismic velocity discontinuities from the crust to the uppermost mantle, we applied receiver function (RF) analysis to seismic records acquired by 17 broadband ocean bottom seismometers deployed across the region in and around the OJP and 3 broadband stations located on ocean islands in Micronesia (one: permanent, two: temporary). The results revealed mid‐crustal discontinuities and the Moho at depths of 10–20 km and 30–40 km (from the top of the basement), respectively, in the central OJP. Moreover, a mantle discontinuity was also imaged at the depth of 55–60 km (from the top of the basement) in the central OJP. These boundaries were not imaged outside the OJP, implying they are characteristic features of the OJP. In addition, RF images showed Moho signals at the depth of 20 km in the eastern OJP, where few previous seismic exploration surveys have been conducted. This depth is comparable with that found in the Manihiki and Hikurangi plateaus that were potentially separated from the OJP.

Description: Abstract Previous compilation of crustal structure in South America had large unsampled areas including the thin crust in the Sub‐Andean lowlands, largely estimated by gravity data, and the sparsely sampled Amazon Craton. A deployment of 35 seismic stations in Brazil, Bolivia, Paraguay, Argentina and Uruguay improved the coverage of the Pantanal Basin in Western Brazil, the intracratonic Paraná and the Chaco basins. Crustal thicknesses and Vp/Vs ratios were estimated with a modified H‐k method by producing three stacked traces to enhance the three Moho conversions (the direct Ps and the two multiples Ppps and Ppss). This modified method gives lower uncertainties than previous studies and shows more regional consistency between nearby stations. The temporary stations and the Brazilian network (RSBR) have characterized the crustal structure as follows. The Paraná Basin has a thick crust 40‐45 km, and average Vp/Vs ratio (1.71‐1.77), while the Chaco Basin has a slightly thinner crust (35‐40 km) and higher Vp/Vs ratio (1.75‐1.79). This confirms the lack of widespread magmatic underplating in the Paraná Basin that could be related to the origin of the flood basalts during the South Atlantic opening. A belt of thin crust (30‐35 km) with low Vp/Vs (〈1.74) is confined to the eastern edge of the Pantanal Basin. Normal crust (38‐43 km) is observed along the western edge of the Pantanal, from the southern part of the Amazon craton to the Rio Apa cratonic block. This study, combined with other published data, provides an updated crustal thickness map of South America.

Description: Abstract We consider fluid‐induced seismicity and present closed‐form expressions for the elastic displacements, strains and stresses resulting from injection into or production from a reservoir with displaced faults. We apply classic inclusion theory to two‐dimensional finite‐width and infinite‐width reservoir models. First we simplify the fault model to the bare minimum while still maintaining its essential features: a vertical fault in a homogeneous reservoir of infinite width in an infinite domain. We confirm and sharpen findings from earlier numerical studies and furthermore conclude that the development of infinitely large elastic shear stresses in a displaced fault, at the internal and external reservoir/fault corners, implies that even small amounts of injection or production will result in some amount of slip or other non‐elastic deformation. Another finding is that there is a marked difference between the shear stress patterns resulting from injection and production in a reservoir with a displaced fault. In both situations two slip patches emerge but at the start of injection some amount of slip occurs immediately in the overburden and underburden, whereas during production the slip may remain inside the reservoir region. Next we derive similar, but more complicated expressions for displaced inclined (normal or reverse) faults and conclude that our findings for vertical faults also apply to inclined faults.

Description: Abstract Loess is an important dust component of airborne particles in Northwestern China. Knowledge of the chemical composition, mixing state, and processing of loess particles in urban plumes is still limited. Urban loess particles were characterized using a single‐particle aerosol mass spectrometer. To understand sources and processing of loess particles, source samples from the road, urban background, soil, construction, and biomass burning ash were collected in the urban areas and characterized. Loess particles were determined as a kind of calcium‐silicate‐rich ones, which were internally mixed with calcium, silicates, potassium, elemental carbon, organics, ammonium, sulfate, and nitrate. Road and soil were major sources of loess particles. Among the aged loess particles, the average peak areas of taken‐up nitrate and sulfate were comparable to that of (Fe+Ca+Al). Diurnal uptake profiles of chloride, sulfate, oxalate, and nitrate on loess particles were analyzed. The nocturnal elevation of chloride occurred significantly due to the uptake of HCl (g). Nighttime nitrate formation occurred prevalently under high relative humidity conditions via the heterogeneous hydrolysis of N2O5. The nighttime enrichment of oxalate, which is a marker for aqueous‐formatted secondary organic aerosol, was also found. Besides the nighttime chemistry, the daytime photochemical activities were also a drive for the elevations of sulfate, nitrate, and ammonium. Conclusively, the processing of loess particles in polluted urban plumes significantly altered their chemical composition and mixing state.

Description: Abstract As one of the typical midlatitude synoptic‐scale disturbances, extra‐tropical cyclones (ECs) can exert significant impacts on the atmospheric general circulation through its interaction with the time‐mean flow. Under the background of global warming, Eurasian continent exhibits evident non‐uniform warming, which has the potential to alter the atmospheric baroclinicity by changing the meridional temperature gradient and further affect the ECs activity. In this study, we investigated the possible connection between the land surface thermal anomaly over Eurasia and the summer ECs activity over East Asia together with relevant mechanisms. We found that the land surface warming (cooling) near 50°N of East Asia is associated with anomalous weak (strong) summer ECs activity over East Asia. Warm (cool) land surface usually reduces (increases) the meridional temperature gradient and further the atmospheric baroclinicity in the key area of cyclone activity, resulting in low (high) frequency of the extratropical cyclogenesis and weakened (intensified) ECs activity. The land surface warming (cooling) can also depress (benefit) the associated baroclinic conversion between the time‐mean effective potential energy and eddy effective potential energy, resulting in decrease (increase) in the eddy kinetic energy. As a result, the energy obtained by the synoptic‐scale eddy from the time‐mean flow has been reduced (increased), which favors (hampers) the extra‐tropical cyclogenesis, causing weak (strong) cyclone activity in the middle latitude of East Asia.

Description: Abstract This paper reports a study to understand the radio spectrum of thunderstorm narrow bipolar events (NBEs) or compact intracloud discharges, which are powerful sources of high frequency (HF) and very high frequency (VHF) electromagnetic radiation. The radio spectra from 10 kHz to about 100 MHz are obtained for three NBEs, including one caused by fast positive breakdown and two by fast negative breakdown. The results indicate that the two polarities of fast breakdown have similar spectra, with a relatively flat spectrum in the HF and VHF band. The ratio of energy spectral densities in the very low frequency and high frequency bands is (0.9‐5)× 105. We develop a statistical modeling approach to investigate if a system of streamers can explain the main features of fast breakdown. Assuming that the current moment peak and charge moment of individual streamers vary in the ranges of 5‐10 A‐m and 5‐20 μC‐m, respectively, the modeling results indicate that a system of 107‐108 streamers can reproduce the current moment, charge transfer and radio spectrum of fast breakdown. The rapid current variation on a timescale of nanoseconds required for fast breakdown to produce strong HF/VHF emissions is provided by exponentially accelerating and expanding streamers. Our study therefore supports the hypothesis that fast breakdown is a system of streamers. Finally, suggestions are given regarding future streamer simulations and NBE measurements in order to further develop our understanding of NBEs and lightning initiation.

Description: Abstract The long‐term trend in dust loading over East Asia remains under debate and is dependent on the study period chosen. In this study, the long‐term trends in springtime dust over East Asia and the North Pacific Ocean (EA&NPO) during 1980–2017 were examined based on the Modern‐Era Retrospective Analysis for Research and Applications version 2 (MERRA‐2) reanalysis. Results showed that there was a spatial gradient in dust aerosol loadings, with decreases from western China eastward towards the NPO. This pattern was corroborated by Cloud‐Aerosol Lidar with Orthogonal Polarization (CALIOP) observations. Furthermore, the empirical orthogonal function (EOF) method was used to reveal the leading modes of springtime dust aerosol optical depth (AOD) over EA&NPO. An abrupt shift occurred in the dust AOD trend in 2010 for the EOF1 mode. The dust AOD increased at a rate of approximately 2×10‐4 per year during 1999–2009, and then decreased more sharply (around 5×10‐4 per year) afterwards. This trend reversal of dust AOD was closely associated with a decrease in 10‐m wind velocity, which induces reduced dust emission. Compared with 10‐m wind, the soil moisture is less correlated with the trend reversal in dust AOD. Additionally, the trends of dry (wet) deposition were closely associated with the trends of the dust AOD, especially for the period 2010–2016. Overall, our findings add new insights to the long‐term nonlinear variability of dust.

Description: Abstract The collision of the Indian plate with Eurasia has played a major role in controlling the dynamics of central Asia leading to the world's largest continental deformation zone. In order to study the deformation within the Indian plate as well as the India‐Eurasia collision zone, we model the lithospheric stress field by calculating the two primary sources of stress, one arising due to topography and shallow lithospheric structure estimated by gravitational potential energy (GPE) differences and the other arising from basal tractions derived from density driven mantle convection. We use several tomography models to calculate horizontal tractions using the convection code HC for two radially varying viscosity structures. We also take into account lateral viscosity variations in the lithosphere model arising from stiff cratons, weak plate boundaries and strength variations due to old and young oceanic lithosphere. We do a quantitative comparison of our predicted deviatoric stresses, strain rates and plate velocities with surface observables and find that the regional tomography model of (A. Singh, Mercier, Ravi Kumar, Srinagesh, & Chadha, 2014) embedded in the global S‐wave model S40RTS does a remarkable job of fitting the observations of GPS velocities and strain rates as well as intraplate stress field from the World Stress Map.

Description: Abstract The eastern and northeastern Tibetan plateau is a key region to study the growth and expansion of the plateau and associated extrusion tectonics. We studied the seismic anisotropic structure in this region by shear‐wave splitting analysis of teleseismic records from a dense linear seismic array, to constrain the lithospheric deformation and processes. We detected small‐scale variations in anisotropy, including changes of splitting parameters around major faults and different anisotropy patterns among individual tectonic blocks and units but with consistent interior features. Our results combined with previous observations suggest that, in addition to the dominant effects of lateral extrusion induced by the India‐Eurasia collision, major faults and tectonic heterogeneity may have also exerted significant impacts on the deformation and thus anisotropic structure of the lithosphere. In particular, we constructed two‐layer anisotropy models for both the Longmenshan sub‐block in the easternmost Songpan‐Ganzi terrane and the Western Qinling orogen, indicating crust‐mantle decoupling in these areas. The lower anisotropic layer of both areas shows a general NW‐SE fast polarization direction (FPD). We attribute this feature to the large‐scale mantle deformation, due to the lateral extrusion of Tibet associated with the India‐Eurasia collision. The upper‐layer anisotropy in both areas features an optimal NEE‐SWW FPD. While in the Longmenshan sub‐block it may stem from crustal deformation under the combined effects of mid‐lower crustal flow, faulting and tectonic heterogeneity, that in the Western Qinling Orogen is probably resulted from shearing caused by upper‐crustal displacement along a mid‐crustal detachment.

Description: Abstract Cross‐correlation of fully diffuse wavefields averaged over time should converge to the Green's function; however, the ambient seismic field in the real Earth is not fully diffuse, which interferes with that convergence. We apply blind signal separation to reduce the effect of spurious non‐diffuse components on the cross‐correlation tensor of the ambient seismic field. We describe the diffuse component as having uncorrelated neighboring frequencies and equal intensity at all azimuths, and an independent (i.e., statistically uncorrelated) non‐diffuse component arising from a spatially isolated point source for which neighboring frequencies are correlated. Under the assumption of linear independence of the spurious non‐diffuse wave outside the stationary phase zone and the constructive interference of noise waves within that zone, we can suppress the spurious non‐diffuse component from the noise interferometry. Our numerical simulations show good separation of one spurious non‐diffuse noise source component for either non‐diffuse Rayleigh or Love waves. We apply this separation to the Rayleigh‐wave component of the Green's function for 136 cross‐correlation pairs from 17 stations in Southern California. We perform beamforming over different frequency bands for the cross‐correlations before and after the separation, and find that the reconstructed Rayleigh waves are more coherent. We also estimate the bias in Rayleigh wave phase velocity for each receiver pair due to the spurious non‐diffuse contribution.

Description: No abstract is available for this article.

Description: Abstract The Monin‐Obukhov similarity theory (MOST) is widely used for the surface turbulence flux‐gradient relations in modeling and data analysis. Here we quantify multi‐scale turbulence processes by applying our newly developed analysis technique to large‐eddy simulation data, and find that in the unstable surface layer, large convective eddies (with the scaling of boundary layer depth) and local free convection exist in addition to small eddies. An empirical MOST function (considering the last two processes only) is found to underestimate the surface friction velocity and heat flux both by about 30%. Much better results can be obtained using a function that explicitly considers all three processes. Generally, the non‐dimensional wind shear exhibits larger scatter and deviates more from the MOST than the temperature gradient. Based on these results, we propose the revised Sorbjan (1986) function (with coefficients determined from this study) for wind shear and MOST function for temperature gradient, for estimating surface fluxes in the unstable surface layer. The three‐dimensional multi‐scale analysis method we develop in this study is of general nature and can be of interest for problems of three‐dimensional multi‐scale process description in other disciplines.

Description: Abstract This paper investigates the relationship between long‐term trends (1980‐2017) in intensity and wind evolution for tropical cyclones (TCs) within the Western Tropical Atlantic (WTA) and Central/Eastern Tropical Atlantic (CETA) sub‐basins. Long‐term TC trends in intensity, intensification time, and wind variability for the CETA were generally more significant than, and in some cases opposite to those for the WTA. Both the TC intensity levels, as measured by the power dissipation index (PDI) normalized by storm hours and proportion of rapid intensification (RI) intervals (defined as a 12‐hour wind speed increase of 20 knots or more), exhibits no long‐term trends in either sub‐basin. A TC wind variability index (WVI) calculated over 72‐h intervals of the TC lifecycle decreases for the WTA over the decades, while the CETA has the 72‐h intervals with the greatest wind speed fluctuations. The average period of intensification before the peak in TC intensity increases ~0.97 h per year for the CETA. TC maximum intensity exhibits no trend, suggesting that TCs in the tropical North Atlantic have a trend favoring a longer intensification period to reach their lifetime maximum intensity. A correlation analysis suggests that warmer sea surface temperatures (SST) and greater moisture favors longer intensification and greater WVI. In contrast, greater 850‐200 hPa vertical wind shear is associated with shorter intensification periods and less WVI.

Description: ABSTRACT Decreases in stratospheric NOx associated with enhanced aerosol have been observed after large volcanic eruptions, e.g., after the eruption of Mt. Pinatubo in 1991. While the 1991 Mt. Pinatubo eruption was the last large explosive eruption, recent studies have shed light on the impacts of moderate‐sized eruptions since the year 2000 on the global stratospheric aerosol budget. We use an ensemble of simulations from a coupled climate‐chemistry model to quantify and analyze changes in NO and NO2 (NOx), N2O5, HNO3, ClO, and ClONO2 during periods of increased stratospheric volcanic aerosol concentrations since 2000. By using an ensemble approach, we are able to distinguish forced responses from internal variability. We also compare the model ensemble results to satellite measurements of these changes in atmospheric composition, including measurements from the Optical Spectrograph and Infrared Imaging Spectrometer on the Odin satellite and the Aura Microwave Limb Sounder. We find decreases in stratospheric NOx concentrations up to 20 hPa, consistent with increases in stratospheric HNO3 concentrations. The HNO3 perturbations also extend higher, up to 5 hPa associated with periods of increased volcanic aerosol concentrations in both model simulations and observations, though correlations with volcanic aerosol are considerably higher in the model simulations. The model simulates increases in ClO at altitudes and magnitudes similar to the NOx reductions, but this response is below the detectable limit in the available observations (100 pptv). We also demonstrate the value of accounting for transport‐related anomalies of atmospheric trace gases by regression onto N2O anomalies.

Description: Abstract Information regarding dust concentrations and size distributions is very important for determining air quality and aerosol–radiation–cloud interactions. Only by using a correct erosion database can the sectional dust emission scheme resolve detailed size distributions and determine where and how dust will be emitted. In this paper, the bias and reasons for dust emission in China Meteorological Administration Unified Atmospheric Chemistry Environment – Dust (CUACE/Dust), an operational dust forecasting model, are analyzed using a heavy sand and dust storm (SDS) episode. We used 18 years of routine SDS phenomena recorded at meteorological stations to retrieve and update the desertification details in the MBA sectional dust emission scheme adopted in CUACE/Dust. New desertification details include decreased erodibility in the area adjacent to Uzbekistan, Turkmenistan, and southern Kazakhstan, where Kyzylkum, Karakum, and Aralkum are located in Central Asia, and in the Chinese deserts of Onqin Daga, Mu Us, and Gurbantungut. New desertification also results in increased erodibility in northern Mongolia. Comparisons show that the new desertification database decreases overestimation of dust emission in Central Asia, including western Mongolia. It improves the underestimation of dust emission in northern Mongolia and the Gobi Desert in southeast Mongolia, and the Taklimakan Desert in China. Consequently, it corrects the overestimated dust cloud in the source area and in areas impacted by dust transportation. The timing, quantitative mass concentrations, and dust size distributions determined here are all more reasonable and rational than those of the original case.

Description: Abstract Groundwater (GW) and recharge as the main drivers of the water budget are challenging to quantify, due to the complexity of hydrological processes and limited observations. Understanding these processes in relation to climate is crucial for evaluating future water availability of Tibetan Plateau (TP). By computing storage changes in Gravity Recovery and Climate Experiment (GRACE) Terrestrial Water Storage (TWS) and Global Land Surface Data Assimilation System (GLDAS) land surface state variables and water balance approach, we calculated GW storage changes and recharges. In the Qaidam Basin (northern TP), TWS from the GRACE revealed a significant increasing trend of 25.5 mm/year during 2002‐2012. However, an obviously turning point was found around 2012 and TWS revealed a significant decreasing rate of 37.9 mm/year during 2013‐2016. Similarly, GW (recharge) had a significant increasing trend of 21.2 (4.5) mm/year before 2012 and a decreasing rate of 32.1 (10.9) mm/year after 2012. Domain‐averaged difference (P‐ET) between precipitation (P) and evapotranspiration (ET) also exhibited an increasing trend of 4.4 mm/year during 2002‐2012 and a decreasing rate of 9.0 mm/year during 2013‐2016. Precipitation followed dissimilar pattern with an increasing rate of 5.3 mm/year during 2002‐2012 while no significant trend during 2013‐2016. However, ET had a consistent increasing trend over the basin during the past 15 years (0.9 mm/year before 2012 and 9.0 mm/year thereafter). This study concluded that GW amount and distribution is mainly controlled by precipitation and ET. Decrease in precipitation at high elevations and increase in ET may impact future groundwater availability in this region.

Description: Abstract In this work we investigate interannual variations in lower stratospheric ozone from 1984 to 2016 based on a satellite‐derived data set and simulations from a chemical transport model. An empirical orthogonal function (EOF) analysis of ozone variations between 2000 and 2016 indicates that the first, second, and third EOF modes are related to the quasi‐biennial oscillation (QBO), canonical El Niño–Southern Oscillation (ENSO), and ENSO Modoki events, respectively; these three leading EOFs capture nearly 80% of the variance. However, for the period 1984–2000, the first, second, and third modes are related to the QBO, ENSO Modoki, and canonical ENSO events, respectively. The explained variance of the second mode in relation to ENSO Modoki is nearly twice that of the third mode for canonical ENSO. Since the frequency of ENSO Modoki events was higher from 1984 to 2000 than after 2000, the Brewer‐Dobson circulation anomalies related to ENSO Modoki were stronger during 1984–2000, which caused ENSO Modoki events to have a greater effect on lower stratospheric ozone before 2000 than after. Ozone anomalies associated with QBO, ENSO Modoki, and canonical ENSO events are largely caused by dynamic processes, and the effect of chemical processes on ozone anomalies is opposite to that of dynamic processes. Ozone anomalies related to dynamic processes are 3–4 times greater than those related to chemical processes.

Description: Abstract Oxides of nitrogen are critical trace gases in the troposphere and are precursors for nitrate aerosol and ozone, which is an important pollutant and greenhouse gas. Lightning is the major source of NOx (NO + NO2) in the mid‐ to upper troposphere. We estimate the production efficiency (PE) of lightning NOx (LNOx) using satellite data from the Ozone Monitoring Instrument (OMI) and the ground‐based World Wide Lightning Location Network (WWLLN) in three northern midlatitude, primarily continental regions that include much of North America, Europe and East Asia. Data were obtained over 5 boreal summers, 2007 – 2011 and comprise the largest number of midlatitude convective events to date for estimating the LNOx PE with satellite NO2 and ground‐based lightning measurements. In contrast to some previous studies, the algorithm assumes no minimum flash‐rate threshold and estimates freshly produced LNOx by subtracting a background of aged NOx estimated from the OMI dataset itself. We infer an average value of 180 ± 100 moles LNOx produced per lightning flash. We also show evidence of a dependence of PE on lightning flash rate and find an approximate empirical power function relating moles LNOx to flashes. PE decreases by an order of magnitude for a 2‐order of magnitude increase in flash rate. This phenomenon has not been reported in previous satellite LNOx studies but is consistent with ground‐based observations suggesting an inverse relationship between flash rate and size.

Description: Abstract We study the statistical properties of tidal weather (variability period 〈30 days) of DW1 amplitude using the extended Canadian Middle Atmospheric Model (eCMAM) and Sounding of the Atmosphere using Broadband Emission Radiometry (SABER). A hierarchy of statistical models, for example, the autoregressive (AR), vector AR, and parsimonious vector AR models, are built to predict tidal weather. The quasi 23‐day oscillation found in the tidal weather is a key parameter in the statistical models. Comparing to the more complex vector AR and parsimonious vector AR models, which consider the spatial correlations of tidal weather, the simplest AR model can predict one‐day tidal weather with an accuracy of 89% (R2: correlation coefficient squared). In the AR model, 23 coefficients at each latitude and height are obtained from seven years of eCMAM data. Tidal weather is predicted via a linear combination of 23 days of tidal weather data prior to the prediction day. Different sensitivity tests are performed to prove the robustness of these coefficients. These coefficients obtained from eCMAM are in very good agreement with those from SABER. SABER tidal weather is predicted with an accuracy of 86% and 87% at one day by the AR models with coefficients from eCMAM and SABER, respectively. The five‐day forecast accuracy is between 60 and 65%.

Description: Abstract Top‐side reverberations off mantle discontinuities are commonly observed at long periods, but their interpretation is complicated because they include both near‐source and near‐receiver reflections. We have developed a method to isolate the station‐side reflectors in large data sets with many sources and receivers. Analysis of USArray transverse‐component data from 3200 earthquakes, using direct S as a reference phase, shows clear reflections off the 410‐ and 660‐km discontinuities, which can be used to map the depth and brightness of these features. Because our results are sensitive to the impedance contrast (velocity and density), they provide a useful complement to receiver‐function studies, which are primarily sensitive to the S velocity jump alone. In addition, reflectors in our images are more spread out in time than in receiver functions, providing good depth resolution. Our images show strong discontinuities near 410 and 660 km across the entire USArray footprint, with intriguing reflectors at shallower depths in many regions. Overall, the discontinuities in the east appear simpler and more monotonous with a uniform transition zone thickness of ~250 km compared to the western United States. In the west, we observe more complex discontinuity topography and small‐scale changes below the Great Basin and the Rocky Mountains, and a decrease in transition‐zone thickness along the western coast. We also observe a dipping reflector in the west that aligns with the top of the high‐velocity Farallon slab anomaly seen in some tomography models, but which also may be an artifact caused by near‐surface scattering of incoming S waves.

Description: Abstract Inelastic rheological behavior, such as viscoelasticity, is increasingly utilized in the modeling of volcanic ground deformation, as elevated thermal regimes induced by magmatic systems may necessitate the use of a mechanical model containing a component of time‐dependent viscous behavior. For the modeling of a given amplitude and footprint of ground deformation, incorporating a viscoelastic regime has been shown to reduce the magma reservoir overpressure requirements suggested by elastic models. This phenomenon, however, is restricted to pressure‐based analyses and the associated creep behavior. Viscoelastic materials exhibit additional constitutive time‐dependent behaviors, determined by the stress and strain states, that are yet to be analyzed in the context of volcanic ground deformation. By utilizing a mechanically homogeneous model space and distinct reservoir evolutions, we provide a comparison of three viscoelastic rheological models, including the commonly implemented Maxwell and Standard Linear Solid configurations, and their time‐dependent behaviors from a fundamental perspective. We also investigate the differences between deformation time series resulting from a pressurization or volume change, two contrasting approaches that are assumed to be equivalent through elastic modeling. Our results illustrate that the perceived influence of viscoelasticity is dependent on the mode of deformation, with stress‐based pressurization models imparting enhanced deformation relative to the elastic models, thus reducing pressure requirements. Strain‐based volumetric models, however, exhibit reduced levels of deformation and may produce episodes of apparent ground subsidence induced by source inflation or vice versa, due to the relaxation of crustal stresses, dependent on whether the reservoir is modeled to be expanding or contracting, respectively.

Description: Abstract Layer 2A, the porous and permeable uppermost igneous oceanic crust, permits the circulation of fluid within the crust, the exchange of dissolved mineral species between the ocean and crust, and the convective dissipation of heat from the crust. We examine the presence, temporal extent, thickness, and evolution of layer 2A using multichannel seismic data collected at 30°S in the South Atlantic across crustal age ranges of 0–70 Ma and half spreading rates of 12–31 mm/year. We observe the layer 2A/2B boundary in 0–48 Myr old crust but not in crust older than ~48 Ma. The thickness of layer 2A in the South Atlantic has substantial variability, with a mean of 760 m and a standard deviation of 290 m. Layer 2A has no systematic change in thickness with age in the South Atlantic, and thickness does not correlate with spreading rate. The crust in the South Atlantic is never fully sealed by sediment cover, which implies that the fluid circulation system in the upper crust never becomes fully closed and the thickness of layer 2A can work as a proxy for the depth at which significant circulation can occur. The disappearance of the layer 2A/2B boundary in older crust implies that fluid circulation within the upper crust continues to occur for at least ~48 Myr after crustal formation in the South Atlantic, after which layer 2A becomes indistinguishable from layer 2B in reflection images.

Description: Abstract There is growing evidence that outgassing through transient fracture networks exerts an important control on conduit processes and explosive‐effusive activity during silicic eruptions. Indeed, the first modern observations of rhyolitic eruptions have revealed that degassed lava effusion may depend upon outgassing during simultaneous pyroclastic venting. The outgassing is thought to occur as gas and pyroclastic debris are discharged through shallow fracture networks within otherwise low‐permeability, conduit‐plugging lava domes. However, this discharge is only transient, as these fractures become clogged and eventually blocked by the accumulation and sintering of hot, melt‐rich pyroclastic debris, drastically reducing their permeability and creating particle‐filled tuffisites. In this study we present the first published permeability measurements for rhyolitic tuffisites, using samples from the recent rhyolitic eruptions at Chaitén (2008‐2009) and Cordón Caulle (2011‐2012) in Chile. To place constraints on tuffisite permeability evolution, we combine (1) laboratory measurements of the porosity and permeability of tuffisites that preserve different degrees of sintering, (2) theoretical estimates on grainsize‐ and temperature‐dependent sintering timescales, and (3) H2O diffusion constraints on pressure‐time paths. The inferred timescales of sintering‐driven tuffisite compaction and permeability loss, spanning seconds (in the case of compaction‐driven sintering) to hours (surface tension‐driven sintering), coincide with timescales of diffusive degassing into tuffisites, observed vent pulsations during hybrid rhyolitic activity (extrusive behaviour coincident with intermittent explosions) and, more broadly, timescales of pressurisation accompanying silicic lava dome extrusion. We discuss herein the complex feedbacks between fracture opening, closing, and sintering, and their role in outgassing rhyolite lavas and mediating hybrid explosive‐effusive activity.

Description: Abstract Summertime surface‐level ozone (O3) is known to vary with temperature, but the relative roles of different processes responsible for causing the O3‐temperature relationship are not well quantified. In this study we use simulations of NASA's Global Modeling Initiative (GMI) chemical transport model (CTM) to isolate and assess the relative impact of atmospheric transport, chemistry, and emissions on large‐scale O3 variability, events, and its covariance with temperature. Using observations from CASTNet in the contiguous United States, we show that the GMI CTM reproduces the spatiotemporal variability of O3 and its relationship with temperature during the summer. We use the change in O3 given a change in temperature (dO3/dT) along with other metrics to understand differences between our simulations. In regions with moderate to strong positive correlations between temperature and O3 such as the Northeast, Great Lakes, and Great Plains, temperature's association with transport yields a majority of the total O3‐temperature relationship (∼60%) while temperature‐dependent chemistry and anthropogenic NO emissions play smaller roles (∼30% and ∼10%, respectively). There are regions, however, with insignificant correlations between temperature and O3, and our findings suggest that transport is still an important driver of O3 variability in these regions, albeit not correlated with temperature. Transport is not directly dependent on temperature but rather is linked through an indirect association, and it is therefore important to understand the exact mechanisms that link transport to O3 and how these mechanisms will change in a warming world.

Description: Abstract Snow cover in mountainous terrain plays an important role in regional and global water and energy balances, climate change, and ecosystems. Blowing snow is a frequent and important weather phenomenon over the Tibetan Plateau (TP); however, this process is neglected in most current land surface models, despite the consequential role it plays in the land surface and atmospheric water and energy budgets. In this paper, we present a blowing snow model PIEKTUK coupled with the Community Land Model (CLM4.5) to provide a better estimate of the snow dynamics for the consideration of snow redistribution induced by wind. Two simulations with a 0.065° spatial resolution were performed in 2010 over the TP, namely, a sensitivity experiment with the inclusion of blowing snow effects (CLM_BS) and a control run with the original model (CLM). A specific objective of this study was to evaluate the improvements in the simulations of snow dynamics and other key variables in surface energy partitioning provided by the coupled model, such as the surface albedo and land surface temperature (LST). Compared with a variety of remote‐sensing observations, the results show that the surface snow cover, snow depth, and surface albedo can be better reproduced in most of the TP region by CLM_BS than by the original CLM, particularly in the Kunlun Mountains, Hoh Xil area, and the southwestern TP. In areas with reduced bias, variations in the monthly mean snow cover fraction can be reflected particularly well by CLM_BS. For LST, however, a significant decrease in the nighttime LST bias was detected in CLM_BS, while the bias in the daytime LST increases. The results show considerable potential for the inclusion of the blowing snow process to promote the modeling of snow dynamics and land‐atmosphere interactions on the TP.

Description: Abstract The Lightning Cluster Filter Algorithm (LCFA) in the Geostationary Lightning Mapper (GLM) ground system identifies lightning flashes from the stream of event detections. It excels at clustering simple flashes, but experiences anomalies with complex flashes that last longer than 3 s or contain more than 100 groups, leading to flashes being artificially split. We develop a technique that corrects these anomalies and apply it to the 2018 GLM data to document all lightning across the Americas. We produce statistics describing the characteristics and frequencies of “reclustered” GLM flashes as well as thunderstorm “area” features. The average GLM Americas flash rate in 2018 was 11.7 flashes s‐1 with the greatest flash rate densities occurring over Lake Maracaibo (157 flashes km‐2 year‐1). Lloró, Chocó, Colombia had the most thunderstorm activity with 256 thunder days. The longest GLM flash spanned 673 km, the largest flash covered 114,997 km2, and the longest‐lasting flash had a 13.496 s duration. The first case occurred over Rio Grande do Sul in Brazil, while the other two cases occurred in the central United States. All three extreme flashes are located in the stratiform regions of Mesoscale Convective Systems (MCSs). The highest flash rate for a thunderstorm area feature was 17.6 flashes s‐1 while the largest thunderstorm was 216,865 km2 in size. Both storms occurred in South America. These initial results demonstrate the value that the development of a reprocessed GLM “science” product would offer and how such a product might be created at a reduced computational cost.

Description: Abstract We present the results of tomographic studies using seismic velocity and attenuation in the area of the Colima Volcanic complex (CVC). Our dataset comprises body waves from local earthquakes recorded by the temporary seismic stations of the CODEX network in the Colima area and a few stations of the regional Mapping the Rivera Subduction Zone (MARS) networks, both deployed in 2006–2008. We obtain three‐dimensional distributions of seismic velocities and attenuation in the crust beneath the CVC area. At shallow depths, we observe a large negative anomaly to the south of CVC, coinciding with the location of the Central Colima Graben. This anomaly may represent debris avalanche deposits, as well as shallow magma reservoirs feeding the eruptions of the presently active Volcán de Colima. In contrast, the volcano edifice of Nevado de Colima, which is built of rigid igneous rocks, is associated with high‐velocity and low‐attenuation anomalies at shallow depths. In the deeper section, a major anomaly with high Vp/Vs, low Vs, and high S wave attenuation corresponds to the location of the regional Tamazula fault. As this represents a mechanically weakened zone of the crust, it may form the pathway that feeds CVC. Both velocity and attenuation models show that the fault‐associated conduit brought magma from the mantle through the lower crust to a depth of 15 km. Then, a light fraction of magma may continue to ascend, forming shallow reservoirs beneath the southern flank of CVC.

Description: Abstract At extensional volcanic arcs, faulting often acts to localize magmatism. Santorini is located on the extended continental crust of the Aegean microplate and is the most active volcano of the Hellenic arc, but the relationship between tectonism and magmatism remains poorly constrained. As part of the PROTEUS experiment, seismic data were acquired across the Santorini caldera and the surrounding region using a dense amphibious array of 〉14,300 marine sound sources and 156 short period seismometers, covering an area 120 km by 45 km. Here, a P‐wave velocity model of the shallow, upper‐crustal structure (〈3 km depth), obtained using travel‐time tomography, is used to delineate fault zones, sedimentary basins, and tectono‐magmatic lineaments. Our interpretation of tectonic boundaries and regional faults are consistent with prior geophysical studies, including the location of basin margins and E‐W oriented basement faults within the Christiana basin west of Santorini. Reduced seismic velocities within the basement east of Santorini, near the Anydros and Anafi basins, are coincident with a region of extensive NE‐SW faulting and active seismicity. The structural differences between the eastern and western sides of Santorini are in agreement with previously proposed models of regional tectonic evolution. Additionally, we find regional magmatism has been localized in NE‐SW trending basin‐like structures that connect the Christiana, Santorini, and Kolumbo volcanic centers. At Santorini itself, we find that magmatism has been localized along NE‐SW trending lineaments that are subparallel to dikes, active faults, and regional volcanic chains. These results show strong interaction between magmatism and active deformation.

Description: Abstract The central and western Sahara (CWS) is the largest source of mineral aerosols during boreal summer, but observed ground‐based data are extremely scarce and typically distant from key source regions. Knowledge of dust emission mechanisms has therefore been mostly limited to short‐term observations from a point or model approximations. To address this deficiency, dust plumes from the CWS are classified according to emission mechanism for June, July and August of 2004‐2017 using an automated inference method which accurately tracks the timing, convective association and geometry of plumes observed with the Spinning Enhanced Visible and Infrared Imager (SEVIRI) aboard Meteosat Second Generation satellites. From these characteristics, plumes are classified as either low‐level jet or cold pool outflow events. The extensive data set is used to generate the largest available climatology of dust emission sources and Saharan emission mechanisms. Automated inference compares well with ground‐based measurements from the Fennec Campaign (76% accuracy) as well as with an entirely manual approach (88% accuracy). Cold pool activity accounts for 82% of total observed dust and 88% at the point of emission. Dust from cold pools evolves seasonally from hotspots around the Mali‐Niger‐Algeria border triple point towards the central Sahara to the northwest, while dust from low‐level jets is organised along the axis of the northeasterly Harmattan, and dominates emission within the Tidihelt Depression of central Algeria. The widespread importance of cold pool outflows in this research supports the findings of the Fennec Campaign, but low‐level jets remain highly significant in certain isolated hotspots.

Description: Abstract Intense forest fires in western North America during August 2017 caused smoke plumes that reached the stratosphere. While this phenomenon has often been observed, this particular event caused increases in stratospheric aerosol extinction at higher altitudes with greater magnitude than previously observed in the satellite record. Here we use multiple satellite limb sounding observations, which provide high sensitivity to thin aerosol layers and good vertical resolution, to show that enhancements in aerosol extinction from the fires reached as high as 23 km in altitude and persisted for more than 5 months. Within 1 month, the aerosol is observed to cover latitudes from 20°N to 60°N, which is essentially the northernmost limit of the observations. At midlatitudes between 15‐ and 20‐km altitudes, the sustained level of median aerosol extinction measured at 750 nm increased by almost an order of magnitude, from approximately 10−4 km−1 to nearly 10−3 km−1. Agreement between limb scatter and occultation measurements is generally within 20% despite potential bias due to modified aerosol shape and composition.

Description: Abstract Entrainment rate (λ) in convective parameterizations remains a sensitive parameter with much uncertainty in model simulations. This study estimates λ using carbon monoxide (CO) measurements jointly from the Microwave Limb Sounder (MLS) and Tropospheric Emission Spectrometer (TES) onboard the Aura satellite, associated with deep convective cases identified by CloudSat and CALIPSO observations. CO is treated as a conserved quantity over convective transport time scales and a simple entraining‐plume model is used to derive entrainment rates. The relationships of the observational estimates of λ as a function of convective height, environmental relative humidity and convective available potential energy (CAPE) derived from Atmospheric Infrared Sounder (AIRS) data are compared with those from Goddard Earth Observing System Model (GEOS‐5) simulations. Bulk statistics of λ show that the values of λ are predominately below 20 % km‐1 for deep convection and the occurrence frequency of any λ decreases with increasing λ. Composite λ values are generally lower in the tropics compared to northern mid‐latitudes in both observations and the GEOS‐5 model. A decrease of λ with increasing convective height is found in both observations and model simulations. We also find that λ tends to decrease with increasing CAPE in the observation‐based λ’s and plume‐based GEOS‐5 λ’s, although the model given λ’s have a non‐monotonic relation with CAPE. The observed λ’s have a weak relation with lower‐to‐mid tropospheric RH, while both GEOS‐5 plume‐based and given λ increases with increasing RH.

Description: Abstract 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/m3. 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, PM10 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.

Description: Abstract Clouds and aerosols play essential roles in regulating surface incident solar radiation (Rs). It has been suggested that the increased aerosol loading over China is a key factor for the decadal variability in Rs and can explain the bias in its trend from reanalyses because the reanalyses do not include the interannual variability of aerosols. In this study, we compare the Rs derived from sunshine duration at 2,400 weather stations in China and that from five reanalyses from 1980 to 2014. The determining factors for the biases in the mean values and trends of Rs from the reanalyses are examined, with the help of Rs and the cloud fraction (CF), from satellite and 2,400 weather stations. Our results show that all reanalyses overestimate the multiyear Rs by 24.10–40.00 W/m2 due to their underestimations of CF, which is more obvious in southern China. The biases in the simulated CF in the reanalyses can explain the biases in Rs by 55–41%, and the bias in clear‐sky surface solar radiation (Rc), which is primarily due to biases in aerosol loading, can explain 32–9% of the bias in Rs. The errors in the trend of the simulated CF can explain the errors in the Rs trends in the reanalyses by 73–12%, and the trend errors in the Rc can explain 43–30% of the trend error in Rs. Our study suggests that more work is needed to improve the simulation of aerosols, clouds, and aerosol‐cloud‐radiation interactions in the reanalyses.

Description: Abstract Tropopause‐penetrating convection is a frequent seasonal feature of the Central United States climate. This convection presents the potential for consistent transport of water vapor into the upper troposphere and lower stratosphere (UTLS) through the lofting of ice, which then sublimates. Water vapor enhancements associated with convective ice lofting have been observed in both in situ and satellite measurements. These water vapor enhancements can increase the probability of sulfate aerosol‐catalyzed heterogeneous reactions that convert reservoir chlorine (HCl and ClONO2) to free radical chlorine (Cl and ClO) that leads to catalytic ozone loss. In addition to water vapor transport, lofted ice may also scavenge nitric acid and further impact the chlorine activation chemistry of the UTLS. We present a photochemical model that resolves the vertical chemical structure of the UTLS to explore the effect of water vapor enhancements and potential additional nitric acid removal. The model is used to define the response of stratospheric column ozone to the range of convective water vapor transported and the temperature variability of the lower stratosphere currently observed over the Central United States in conjunction with potential nitric acid removal and to scenarios of elevated sulfate aerosol surface area density representative of possible future volcanic eruptions or solar radiation management. We find that the effect of HNO3 removal is dependent on the magnitude of nitric acid removal and has the greatest potential to increase chlorine activation and ozone loss under UTLS conditions that weakly favor the chlorine activation heterogeneous reactions by reducing NOx sources.

Description: Abstract With the highlight of environmental problems over the Tibetan Plateau (TP), aerosol pollution and the influence of this pollution on cloud properties are becoming a new area of research. Based on the aerosol index and cloud property parameters derived from satellite observations, in this study, the inconsistent effects of aerosols on ice cloud properties between daytime and nighttime over the TP are investigated. The results indicate that ice clouds are mainly distributed over the TP margin area, especially over the north slope, during both daytime and nighttime. The occurrence frequency of ice cloud is higher during the daytime than during the nighttime over the margin areas of the TP. Similarly, aerosols are mainly concentrated over the northern margin of the TP. A potential relationship may exist between the aerosol index and ice cloud properties. When the aerosol index increases from 0.05 to 0.17, the ice cloud droplet radius (ICDR) during the daytime decreases from 32.1 to 27.9 μm, while the ICDR during the nighttime remains almost constant (approximately 25 μm); furthermore, the ice water path (IWP) during the daytime decreases slightly due to the saturation effect, while the nocturnal IWP increases significantly. The changes in ice cloud optical depth (ICOD) during daytime and nighttime show significant and completely opposite trends. The removal of the influence of meteorological factors showed that aerosols have a more dominant influence than meteorological conditions on ice cloud properties (except for the nocturnal ICDR and IWP during the daytime).

Description: Abstract The Royal Netherlands Meteorological Institute (KNMI) operates a three‐dimensional microbarometer array at the Cabauw Experimental Site for Atmospheric Research observatory. The array consists of five microbarometers on a meteorological tower up to an altitude of 200 m. Ten ground‐based microbarometers surround the tower with an array aperture of 800 m. This unique setup allows for the study of infrasound propagation in three dimensions. The added value of the vertical dimension is the sensitivity to wind and temperature in the atmospheric boundary layer over multiple altitudes. In this study, we analyze infrasound generated by an accidental chemical explosion at the Moerdijk petrochemical plant on 3 June 2014. The recordings of the tower microbarometers show two sequential arrivals, whereas the recordings on the ground show one wavefront. This arrival structure is interpreted to be the upgoing and downgoing wavefronts. The observations are compared with propagation modeling results using global‐scale and mesoscale atmospheric models. Independent temperature and wind measurements, which are available at the Cabauw Experimental Site for Atmospheric Research, are used for comparison with model output. The modeling results explain the signal arrival times; however, the tower wavefront arrivals are not explained. This study is important for understanding the influence of the atmospheric boundary layer on infrasound detections and propagation.

Description: Abstract The aftershock productivity is known to strongly vary for different mainshocks of the same magnitude, which cannot be simply explained by random fluctuations. In addition to variable source mechanisms, different rheological properties might be responsible for the observed variations. Here we show, for the subduction zone of northern Chile, that the aftershock productivity is linearly related to the degree of mechanical coupling along the subduction interface. Using the earthquake catalog of Sippl et al. (2018, https://doi.org/10.1002/2017JB015384), which consists of more than 100,000 events between 2007 and 2014, and three different coupling maps inferred from interseismic geodetic deformation data, we show that the observed aftershock numbers are significantly lower than expected from the Båth's law. Furthermore, the productivity decays systematically with depth in the uppermost 80 km, while the b value increases. We show that this lack of aftershocks and the observed depth dependence can be simply explained by a linear relationship between the productivity and the coupling coefficient, leading to Båth law only in the case of full coupling. Our results indicate that coupling maps might be useful to forecast aftershock productivity and vice versa.

Description: Abstract Spectral induced polarization spectra were carried out on three graphitic schists and two graphitic sandstones. The microstructural arrangement of graphite of two graphitic schists was studied with thin sections using transmitted and reflected light optical and electron microscopic methods. Chemical maps of selected areas confirm the presence of carbon. The complex conductivity spectra were measured in the frequency range 10 mHz to 45 kHz and in the temperature range +20 °C down to −15 °C. The measured spectra are fitted with a double Cole‐Cole complex conductivity model with one component associated with the polarization of graphite and the second component associated with the Maxwell‐Wagner polarization. The Cole‐Cole exponent and the chargeability are observed to be almost independent of temperature including in freezing conditions. The conductivity and relaxation time are dependent on the temperature in a predictable way. As long as the temperature decreases, the electrical conductivity decreases and the relaxation time increases. A finite element model is able to reproduce the observed results. In this model, we consider an intragrain polarization mechanism for the graphite and a change of the conductivity of the background material modeled with an exponential freezing curve. One of the core sample (a black schist), very rich in graphite, appears to be characterized by a very high conductivity (approximately 30 S/m). Two induced polarization profiles are discussed in the area of Thorens. The model is applied to the chargeability data to map the volumetric content of graphite.

Description: Abstract Receiver function analysis is widely used to image sharp structures in the Earth, such as the Moho or transition zone discontinuities. Standard procedures either rely on the assumption that underlying discontinuities are horizontal (common conversion point stacking) or are computationally expensive and usually limited to 2‐D geometries (reverse time migration and generalized Radon transform). Here, we develop a teleseismic imaging method that uses fast 3‐D traveltime calculations with minimal assumption about the underlying structure. This allows us to achieve high computational efficiency without limiting ourselves to 1‐D or 2‐D geometries. In our method, we apply acoustic Kirchhoff migration to transmitted and reflected teleseismic waves (i.e., receiver functions). The approach expands on the work of Cheng et al. (2016, https://doi.org/10.1093/gji/ggw062) to account for free surface multiples. We use an Eikonal solver based on the fast marching method to compute traveltimes for all scattered phases. Three‐dimensional scattering patterns are computed to correct the amplitudes and polarities of the three component input signals. We consider three different stacking methods (linear, phase weighted, and second root) to enhance the structures that are most coherent across scattering modes and find that second‐root stack is the most effective. Results from synthetic tests show that our imaging principle can recover scattering structures accurately with minimal artifacts. Application to real data from the Multidisciplinary Experiments for Dynamic Understanding of Subduction under the Aegean Sea experiment in the Hellenic subduction zone yields images that are similar to those obtained by 2‐D generalized Radon transform migration at no additional computational cost, further supporting the robustness of our approach.

Description: Abstract Forecasting the onset of a volcanic eruption from a closed system requires understanding its stress state and failure potential, which can be investigated through numerical modeling. However, the lack of constraints on model parameters, especially rheology, may substantially impair the accuracy of failure forecasts. Therefore, it is essential to know whether large variations and uncertainties in rock properties will preclude the ability of models to predict reservoir failure. A series of two‐dimensional, axisymmetric models are used to investigate sensitivities of brittle failure initiation to assumed rock properties. The numerical experiments indicate that the deformation and overpressure at failure onset simulated by elastic models will be much lower than the viscoelastic models, when the timescale of pressurization exceeds the viscoelastic relaxation time of the host rock. Poisson's ratio and internal friction angle have much less effect on failure forecasts than Young's modulus. Variations in Young's modulus significantly affect the prediction of surface deformation before failure onset when Young's modulus is 〈 40 GPa. Longer precursory volcano‐tectonic events may occur in weak host rock (E 〈 40 GPa) due to well‐developed Coulomb failure prior to dike propagation. Thus, combining surface deformation with seismicity may enhance the accuracy of eruption forecast in these situations. Compared to large and oblate magma systems, small and prolate systems create far less surface uplift prior to failure initiation, suggesting that more frequent measurements are necessary.

Description: Abstract Core‐mantle boundary (CMB) topography may provide useful hints on the deep mantle thermochemical structure, as clusters of thermal plumes and piles of chemically differentiated material, which are usually proposed as end‐member explanations for the large low shear‐wave velocity regions observed in the deep mantle, have different actions on this topography. CMB topography is further sensitive to several parameters, including mantle viscosity and its variations with thermal and compositional changes. Here we assess the influence of the postperovskite (pPv) phase viscosity on deep mantle dynamics and on CMB topography. We perform numerical simulations of thermal and thermochemical convection in spherical geometry, varying the ratio between pPv and bridgmanite viscosities, ΔηpPv, between 1 (regular pPv) and 10−3 (weak pPv). Thermochemical structures are dominated by smaller‐scale wavelengths (spherical harmonic degrees 3 to 6) and are more stable in weak than in regular pPv models. The amplitude of CMB topography is reduced by about a factor of 2 as ΔηpPv changes from 1 to 10−3, mostly due to a sharp drop in the depressions induced by downwellings reaching the CMB. By contrast, the topographies induced by plumes clusters and thermochemical piles are mostly unaffected. For all the values of ΔηpPv we tested, long‐wavelength CMB topography and reconstructed shear‐wave tomography are anticorrelated in purely thermal models, and correlated in thermochemical models with strong chemical density contrast (ΔρC = 140 kg/m3). In models with smaller density contrast (ΔρC = 90 kg/m3), topography and tomography are anticorrelated at ΔηpPv = 1, but correlated at ΔηpPv = 10−3.

Description: Abstract In this study, the micromechanical interparticle contact behavior of “De NoArtri” (DNA‐1A) grains is investigated, which is a lunar regolith simulant, using a custom‐built micromechanical loading apparatus, and the results on the DNA‐1A are compared with Ottawa sand which is a standard quartz soil. Material characterization is performed through several techniques. Based on microhardness intender and surface profiler analyses, it was found that the DNA‐1A grains had lower values of hardness and higher values of surface roughness compared to Ottawa sand grains. In normal contact micromechanical tests, the results showed that the DNA‐1A had softer behavior compared with Ottawa sand grains and that cumulative plastic displacements were observed for the DNA‐1A simulant during cyclic compression, whereas for Ottawa sand grains elastic displacements were dominant in the cyclic sequences. In tangential contact micromechanical tests, it was shown that the interparticle friction values of DNA‐1A were much greater than that of Ottawa sand grains, which was attributed to the softer contact response and greater roughness of the DNA‐1A grains. Widely used theoretical models both in normal and tangential directions were fitted to the experimental data to obtain representative parameters, which can be useful as input in numerical analyses which use the discrete element method.

Description: Abstract Volcanic plumes from small and moderate eruptions represent a challenge in the study of plume morphology due to eruption source parameter uncertainties and atmospheric influence. Sakurajima volcano, Japan, features such activity and due to its continuous eruptions in the recent years provides an ideal natural laboratory. A data set of 896 eruptions between 2009 and 2016 with well‐constrained plume heights, estimated erupted mass, and associated atmospheric conditions has been compiled. Plume heights ranged between 1,500 and 5,000 m and mainly developed under stable atmospheric stratification and low background wind speeds. The eruptions presented in the database were used to drive FPLUME, a 1‐D integral volcanic plume model, to study the simulated plume morphology. FPLUME was seen to provide consistent results under stable atmospheric stratification. A method for the real‐time monitoring of erupted mass used in the Sakurajima observatory was seen to provide appropriate first guess estimates for the eruptions, showing agreement with analytical and simulated mass flow rate calculations. Volcanic plumes from Sakurajima show significant influence by the atmospheric environment. The plume scaling parameter (Π) was used to characterize the expected degree of plume bending with results correlating well against modeled plume angles. The vertical wind profile was seen to have a significant impact on the resolved plume. Wind shear characteristics were seen to have a mechanical effect on the plume, aiding or inhibiting bending. Finally, potential issues were identified in simulations under unstable atmospheric conditions as the model either failed to provide a solution or overestimated the plume height.

Description: Abstract The sudden stratospheric warming (SSW) of 12 February 2018 was not forecast by any extended‐range model beyond 12 days. From early February, all forecast models that comprise the subseasonal‐to‐seasonal (S2S) database abruptly transitioned from indicating a strong stratospheric polar vortex (SPV) to a high likelihood of a major SSW. We demonstrate that this forecast evolution was associated with the track and intensity of a cyclone in the northeast Atlantic, with an associated anticyclonic Rossby wave break, which was not well forecast. The wave break played a pivotal role in building the Ural high, which existing literature has shown was a precursor of the 2018 SSW. The track of the cyclone built an anomalously strong sea level pressure dipole between Scandinavia and Greenland (termed the S‐G dipole), which we use as a diagnostic of the wave break. Forecasts that did not capture the magnitude of this event had the largest errors in the SPV strength and did not show enhanced vertical wave activity. A composite of 49 similarly strong wintertime (November–March) S‐G dipoles in reanalysis shows associated anticyclonic wave breaking leading to significantly enhanced vertical wave activity and a weakened SPV in the following days, which occurred in 35% of the 15‐day periods preceding observed major SSWs. Our results indicate a particular transient trigger for weakening the SPV, complementing existing results on the importance of tropospheric blocking for disruptions to the Northern Hemisphere extratropical stratospheric circulation.

Description: Abstract The physical mechanism of intermediate‐depth earthquakes is still uncertain. Dehydration embrittlement and thermal shear heating mechanisms are the leading hypotheses, and each has been supported both by observations and experiments. Slab character is likely to affect either mechanism. We apply uniform processing to data sets from the two main subduction zones in Japan: the older, colder, and faster‐subducting Pacific plate and the younger, warmer, and slower‐subducting Philippine Sea plate. We compare the stress drops and radiated efficiencies of intermediate‐depth earthquakes in these settings and find no significant differences between the scaling of source properties. In particular, we find both an increase of stress drop and apparent stress with increasing moment for the Pacific Plate subduction in Hokkaido and for the Philippine Sea Plate subduction in Kyushu. We suggest that this, along with apparent invariance of radiated efficiency, suggests that an embrittlement process is more important in these regions than a thermal shear mechanism.

Description: Abstract Temperature distribution at depth is of key importance for characterizing the crust, defining its mechanical behavior and deformation. Temperature can be retrieved by heat flow measurements in boreholes that are sparse, shallow, and have limited reliability, especially in active and recently active areas. Laboratory data and thermodynamic modeling demonstrate that temperature exerts a strong control on the seismic properties of rocks, supporting the hypothesis that seismic data can be used to constrain the crustal thermal structure. We use Rayleigh wave dispersion curves and receiver functions, jointly inverted with a transdimensional Monte Carlo Markov Chain algorithm, to retrieve the VS and VP/VS within the crust in the Italian peninsula. The high values (〉1.9) of VP/VS suggest the presence of filled‐fluid cracks in the middle and lower crust. Intracrustal discontinuities associated with large values of VP/VS are interpreted as the α−β quartz transition and used to estimate geothermal gradients. These are in agreement with the temperatures inferred from shear wave velocities and exhibit a behavior consistent with the known tectonic and geodynamic setting of the Italian peninsula. We argue that such methods, based on seismological observables, provide a viable alternative to heat flow measurements for inferring crustal thermal structure.

Description: Abstract Low‐δ26Mg basalts are commonly interpreted to represent melts derived from carbonated mantle sources. The mantle domain feeding low‐δ26Mg Cenozoic basalts in eastern China overlaps the so‐called Big Mantle Wedge (BMW) above the stagnant Pacific slab in the mantle transition zone, which indicates that the BMW is an important carbon reservoir generated by the slab. However, Mg isotopic composition in the nearby mantle beyond the BMW and, thus, the spatial extent of carbonated components in the mantle beneath eastern Asia have not yet been extensively characterized. Therefore, it remains largely unconstrained if additional or alternative carbon reservoirs exist. Here we carried out a geochemical study on Cenozoic Huihe nephelinites, which crop out ~500 km west of the present‐day BMW. These rocks are characterized by negative K, Zr, Hf, and Ti anomalies, high Zr/Hf, Ca/Al ratios, and low δ26Mg values, which suggest that they are derived from a carbonated mantle source. The composition of the nephelinites demonstrates that low δ26Mg mantle components exist at significant distances from the present‐day BMW, which highlights that in addition to the stagnant Pacific slab, other oceanic slab(s) also contribute(s) carbonate‐bearing crustal materials to the mantle sources of Cenozoic volcanism in eastern Asia.

Description: Abstract The 2016–2017 Central Apennines earthquake sequence is a recent example of how damages from subsequent aftershocks can exceed those caused by the initial mainshock. Recent studies reveal that physics‐based aftershock forecasts present comparable skills to their statistical counterparts, but their performance remains a controversial subject. Here we employ physics‐based models that combine the elasto‐static stress transfer with rate‐and‐state friction laws, and short‐term statistical Epidemic Type Aftershock Sequence (ETAS) models to describe the spatiotemporal evolution of the earthquake cascade. We then track the absolute and relative model performance using log‐likelihood statistics for a 1‐year horizon after the 24 August 2016 Mw = 6.0 Amatrice earthquake. We perform a series of pseudoprospective experiments by producing seven classes of Coulomb rate‐state (CRS) forecasts with gradual increase in data input quality and model complexity. Our goal is to investigate the influence of data quality on the predictive power of physics‐based models and to assess the comparative performance of the forecasts in critical time windows, such as the period following the 26 October Visso earthquakes leading to the 30 October Mw = 6.5 Norcia mainshock. We find that (1) the spatiotemporal performance of the basic CRS models is poor and progressively improves as more refined data are used, (2) CRS forecasts are about as informative as ETAS when secondary triggering effects from M3+ earthquakes are included together with spatially variable slip models, spatially heterogeneous receiver faults, and optimized rate‐and‐state parameters. After the Visso earthquakes, the more elaborate CRS model outperforms ETAS highlighting the importance of the static stress transfer for operational earthquake forecasting.

Description: Abstract We present a theoretical study focusing on exploring the possibility of controlling anthropogenic and natural seismicity. We actively control the pressure of injected fluids using a negative‐feedback control system. Our analysis is based on the spring‐slider model for modeling the earthquake instability. We use a general Coulomb‐type rheology for describing the frictional behavior of a fault system. This model leads to a nonautonomous system, whose steady state and stability are studied using a double‐scale asymptotic analysis. This approach renders the dominant order of the system time invariant. Established tools from the classical mathematical theory of control are used for designing a proper stabilizing controller. We show that the system is stabilizable by controlling fluid pressure. This is a central result for industrial operations. A stabilizing controller is then designed and tested. The controller regulates in real time the applied pressure in order to assure stability, avoid unwanted seismicity, and drive the system from unstable states of high potential energy, to stable ones of low energy. The controller performs well even in the absence of complete knowledge of the frictional properties of the system. Finally, we present two numerical examples (scenarios) and illustrate how anthropogenic and natural earthquakes could be, in theory, prevented.

Description: Abstract We integrate paleoseismic datasets along the Mt. Vettore‐Mt. Bove normal fault‐system (VBFS) rupturing at surface in the 30 October 2016 Norcia earthquake. Through the analysis of new trenches from this work and a review of the pre‐existing data, we correlate events among trench sites along antithetic and synthetic fault splays. We recognize seven M6.5, 2016 Norcia‐type (or larger) surface‐faulting events in the last ~22 kyr, including 2016. Before 2016, one event occurred in the past two millennia (260‐575 CE), and possibly corresponds to the event damaging Rome in 443 CE or 484/508 CE. Three previous events occurred between 10590 BCE and 415 BCE, whereas the two oldest ones date between 19820 BCE and 16540 BCE. The average recurrence time is 3360–3640 yrs for the last ~22 kyr, and 1220‐1970 yrs for the last ~4 kyr. We infer a minimum dip‐slip rate of 0.26‐0.38 mm/yr on the master fault in the central portion of the VBFS, and a dip‐slip rate of at least 0.10 mm/yr on the southernmost portion. We infer a Middle‐Late Pleistocene inception of the long‐term scarp of the investigated splays. The along‐strike variation of slip rates well reproduces the trend of the 2016 surface slip, thus the time window exposed in the trenches is representative for the present fault activity. Based on trenching data, different earthquake rupture scenarios should be also considered for local hazard assessment.

Description: Abstract Soot particles are generally considered to be poor ice nucleating particles. Involvement of soot in clouds and their release back into the atmosphere can form residual particles with altered cloud forming potential. The impact and extent of such different cloud processing scenarios on ice nucleation is however not well understood. In this work, we present the impact of cloud processing of soot aerosols on subsequent ice nucleation cycles at T ≤ 233 K. Coupling of two continuous flow diffusion chambers allows the simulation of different cloud processing scenarios and investigation of subsequent ice nucleation activity of the processed particles. The processing scenarios presented here encompass contrail, cirrus and mixed‐phase cloud processing, mimicking typical pathways that soot particles can be exposed to in the atmosphere. For all scenarios tested, the processed particles showed an enhanced ice active fraction for T 〈 233 K. The relative humidity with respect to water for the ice nucleation onset was observed to be on average approximately 10% (relative humidity with respect to ice, RHi ≈ 16 %) lower for the cloud processed particles compared to the unprocessed soot, for which ice nucleation was observed close to or at homogeneous freezing conditions of solution droplets. We attribute the enhanced ice nucleation abilities of the cloud processed soot to a pore condensation and freezing mechanism and have identified key parameters governing these changes. Enhanced ice nucleation abilities of soot in cirrus clouds can have significant impacts, given the importance of the atmospheric ice phase for precipitation formation and global climate.

Description: Abstract The mobile south flank of Kīlauea Volcano hosts two normal fault systems, the Koa'e fault system (KFS) and the Hilina fault system (HFS). In historical time, at least three M〉6.5 earthquakes have occurred on the basal detachment of the Kīlauea Volcano's south flank, with the most recent being the May 4, 2018 M6.9 earthquake. Here we analyze kinematic GPS data collected from 2001 to 2017, and InSAR data before, during and after the 2018 M6.9 earthquake to determine the crustal motion across the HFS and KFS faults. Our results indicate that the HFS faults did not significantly slip during the interseismic period from 2007 to 2011. Despite its substantial magnitude, InSAR shows that the 2018 M6.9 earthquake triggered sub‐cm level slip along sections of the previously mapped HFS branches. Up to 20 cm of offset occurred on what appears to be a newly formed (or previously unknown) fault near the eastern end of the HFS. During the 3 months following the M6.9 earthquake, up to more than 30 cm of slip occurred along the KFS, which helps accommodate rapid large‐scale subsidence of Kīlauea's summit region as large volumes of summit reservoir magma fed the lower East Rift Zone eruption. The HFS appears to activate only in concert with large earthquakes on the basal detachment. The KFS, on the other hand, moves both seismically during small local earthquakes, and aseismically in response to nearby earthquakes and caldera subsidence.

Description: Abstract Correlations within and between Precambrian basins are heavily reliant on precise dating of volcanic units (i.e., tuff beds and lava flows) in the absence of biostratigraphy. However, felsic tuffs and lavas are rare or absent in many basins and direct age determinations of Precambrian basaltic lavas have proven to be challenging. In this paper, we report the first successful application of 40Ar/39Ar dating to pyroxene from a Neoproterozoic basalt unit, the Keene Basalt in the Officer Basin of central Australia. 40Ar/39Ar analyses of igneous pyroxene crystals yielded an age of 752 ± 4 Ma (MSWD = 0.69, probability = 72%), which is underpinned by 40Ar/39Ar plagioclase age (753.04 ± 0.84 Ma) from the basalt. This age is significant because the Keene Basalt is one of the very few extrusive igneous rocks identified within the Neoproterozoic successions of central Australia, and is potentially an important time marker for correlating the Neoproterozoic stratigraphy within, and beyond, the central Australian basins. Our geochronological and geochemical data show that the Keene Basalt, which is characterized by enriched elemental and Nd‐Pb isotopic signatures, is strikingly similar to, and coeval with, the 755 ± 3 Ma Mundine Well Dolerite in northwestern Australia. Here, we suggest that both are part of the same large igneous province (~6.5 × 105 km2) related to breakup of the supercontinent Rodinia. This study demonstrates the potential of pyroxene 40Ar/39Ar geochronology to date ancient flood basalts, and to provide pivotal time‐constraints for stratigraphic correlations of Precambrian basins.

Description: Abstract Brucite, Mg (OH)2, is an important analog for studying the thermodynamics of hydrous silicate minerals in the deep Earth, as well as H/D isotope fractionation between minerals and water. In this study, we measured in situ Raman and Fourier transform infrared spectra for the natural and deuterated brucite samples, at high temperatures to 650 K, just before the dehydration of brucite at ambient pressure. All of the optical modes systematically shift to lower frequencies at elevated temperature, while deuterium substitution reduces the magnitudes of the temperature dependence. The isobaric mode Grüneisen parameters (γiP), as well as the intrinsic anharmonic parameters (ai), have been evaluated for the vibrational modes between Mg (OH)2 and Mg (OD)2. The anharmonic contribution to the thermodynamic properties (such as internal energy, isochoric and isothermal heat capacities, and entropy) is negative and severe at high temperature. The difference in the heat capacity is up to ~7% at 700 K due to the anharmonic effect. The deuterium isotopic effect on the thermodynamics is positive, and the magnitude of the isotopic effect is comparable to that from the anharmonic effect. On the other hand, the anharmonicity significantly increases the magnitude of the positive pressure dependence of the D/H fractionation β factor for brucite, and this correction could be more important at elevated temperature. At the temperature of 800 K, 103·(∂lnβ/∂P)T increases from +0.23 GPa−1 (for quasi‐harmonic approximation) to +0.44 GPa−1, due to the anharmonic correction.

Description: No abstract is available for this article.

Description: Abstract To evaluate the effect of melt viscosity on bubble nucleation, we formulated the homogeneous nucleation rate of water bubbles to explicitly include melt viscosity. The viscosity coefficient appears in the preexponential factor of the nucleation rate in terms of the Péclet number: the ratio of the bubble growth timescale by molecular diffusion and the viscous relaxation timescale. The preexponential factor is almost constant when viscosity is low (or a high Péclet number), whereas it linearly decreases with increasing viscosity (or a decreasing Péclet number) exceeding the crossover value of viscosity, under a given supersaturation. The crossover point depends on whether homogeneous or heterogeneous nucleation takes place. We numerically solved the evolution of bubble nucleation and growth processes in ascending magmas by using the new nucleation rate formula and a precise approximation of moment equations of the bubble size distribution function. The resultant bubble number density has two regimes, similar to the previous study, but the transition point between the diffusion‐controlled regime and the viscosity‐controlled regime moves to higher viscosity or higher decompression rates by 0.6 log units at the maximum. In the viscosity‐controlled regime, the effect of the better approximation of bubble size distribution moment equations reduces bubble number density by a few orders of magnitude compared with the previous study. As a result of compiling the past laboratory experimental data, it turned out that all the experiments are conducted under the conditions equivalent to the diffusion‐controlled regime. We propose an experimental condition to confirm the presence of the viscosity‐controlled regime.

Description: Abstract Seismic observations suggest (1) significant accumulation of subducted slabs above the 670‐km discontinuity in many subduction zones, (2) possible structure change at ~1,000‐km depth, and (3) the large low shear wave velocity provinces above the core‐mantle boundary in the African and Pacific lower mantle be associated with chemical heterogeneity. Global mantle convection models with realistic plate motion history reproduce most of these structures. However, it remains unclear how the convection models compare with seismic models at different spatial wavelengths and depths. By conducting quantitative analysis between mantle convection and seismic models, we found that mantle convective structures show significant correlations with seismic structures in the upper mantle and mantle transition zone for wavelengths up to spherical harmonic degree 20. However, the global correlation is weak at intermediate to short wavelengths (for degrees 4 and higher) in the lower mantle below ~1,000‐km depth. A weak layer beneath the spinel‐to‐postspinel phase change help consistently reproduce stagnant slabs in the western Pacific, while having insignificant effects elsewhere, that is, the large low shear wave velocity province structures. The cold slab structures and their correlations with the seismically fast anomalies are nearly identical for our convection models with and without the plumes, indicating that seismically fast anomalies in the mantle mainly result from the subducted slabs. Models with viscosity increase at 1,000‐km depth and the 670‐km depth phase change may reproduce seismic slab structures including the stagnant slabs in the mantle transition zone equally well as models with a thin weak layer below the 670‐km phase boundary.

Description: Abstract Based on a successful cloud‐resolving simulation with the Weather Research and Forecasting Model, this study examines the evolution and the role of midtropospheric mesoscale cyclonic vortex in the formation of Super Typhoon Nepartak (2016). The midtropospheric vortex is correlated with the convective activity in pre‐Nepartak. Once the deep convection outbreaks, the midtropospheric vortex intensifies first via the vertical advection associated with the severe updrafts and then through the midlevel convergence associated with stratiform precipitation. As the stratiform precipitation dissipates, the midlevel vortex weakens slightly in the following shallow convection phase. The above‐described processes recur sequentially during the pregenesis of Nepartak, and the midtropospheric vortex demonstrates diurnal variations. Its intensification usually corresponds to the weakening of low‐level cyclonic circulation except for the deep convection phase, indicating that the development of midtropospheric vortex can inhibit the development of self‐sustained low‐level cyclonic circulation. Although the midtropospheric vortex is not always a quasi‐balanced perturbation, a cold core can be found in the lower troposphere below it during the most of the pregenesis stage. The appearance of the cold core enhances the low‐level temperature gradient around it, which favors convection burst. In addition, the closed cyclonic circulation associated with the midlevel vortex can serve as a pouch protecting the vorticity, moisture, and convection inside from the vertical wind shear and dry air intrusion when the low‐level and midlevel vortices are overlapped in the late pregenesis stage, which facilitates the sustained deep convection and the formation of Nepartak.

Description: Abstract The dynamical behavior of the mesosphere and lower thermosphere (MLT) region during strongly disturbed wintertime conditions commonly known as polar‐night jet oscillations (PJOs) is described in detail and compared to other wintertime conditions. For this purpose, wind measurements provided by two specular meteor radars located at Andenes (69°N, 16°E) and Juliusruh (54°N, 13°E) are used to estimate horizontal mean winds and tides as an observational basis. Winds and tidal main features are analyzed and compared for three different cases: major sudden stratospheric warming (SSW) with (a) strong PJO event, (b) non‐PJO event, and (c) no major SSWs. We show that the distinction into strong PJOs, non‐PJOs, and winters with no major SSWs is better suited to identify differences in the behavior of the mean winds and tides during the boreal winter. To assess the impact of the stratospheric disturbed conditions on the MLT region, we investigate the 30‐year nudged simulation by the Extended Canadian Middle Atmosphere Model. Analysis of geopotential height disturbances suggests that changes in the location of the polar vortex at mesospheric heights are responsible for the jets observed in the MLT mean winds during strong PJOs, which in turn influence the evolution of semidiurnal tides by increasing or decreasing their amplitudes depending on the tidal component.

Description: Abstract Forced by Pacific Decadal Oscillation‐related sea surface temperature (SST) anomalies with the same pattern but opposite signs in the western‐central North Pacific, nonlinear wintertime atmospheric responses are produced by a state‐of‐the‐art atmospheric general circulation model (GFDL AM2.1); that is, an obvious equivalent barotropic geopotential low appears over the cold SST forcing (“CSST”), whereas a weak baroclinic structure shows up corresponding to the warm SST forcing (“WSST”), and both of them have similar characteristics in the lower troposphere. Specifically, because of the relatively dry environment in the central North Pacific, nonlinear responses of moisture process including latent heat flux and low‐level atmosphere moisture advection induce asymmetric diabatic heating (Qd): in WSST, Qd tends to increase in the middle‐lower troposphere but decrease in the middle‐upper level, whereas it always increases in the whole troposphere in CSST. Thus, Qd has the same low‐level positive vertical gradient in both CSST and WSST, which produces similar atmospheric circulation anomalies in the lower troposphere. In turn, the asymmetric responses of low‐level temperature advection further modify air temperature meridional gradient as well as atmospheric baroclinicity in the lower troposphere, significantly shifting the transient eddy activities southward in CSST and greatly weakening their intensity in WSST, respectively. Accordingly, the transient eddy vorticity forcing primarily determines the upper‐level atmospheric responses in CSST, but it has unsystematic effects in WSST that are overtaken by Qd. Therefore, the dominance of diabatic heating in WSST and transient eddy forcing in CSST over the central North Pacific lead to the asymmetric atmospheric responses among which the asymmetry of moisture plays an essential role.

Description: Abstract Thirty southerly low‐level jet (LLJ) events were observed during the Plains Elevated Convection at Night (PECAN) field campaign in the Great Plains region of the United States during summer 2015. Here we present Doppler lidar wind data from three PECAN instrumentation sites to explore characteristics of LLJs and the boundary layer as well as some of the heterogeneities possible within the wind field of a LLJ. Southerly LLJs were observed on 66% of nights at the southwestern site (Greensburg, KS) but only 56% and 53% of nights at the eastern and northern sites, respectively (Hesston and Ellis, KS). The northernmost site had a relative abundance of weaker jets or nonjet conditions due to fronts or convective systems that only affected part of the observation domain. Plotting mean wind fields of each LLJ type reveals that the strongest LLJs tend to develop under very similar conditions but begin to show variability in wind profile evolution after several hours. A robust mixed layer height retrieval algorithm is used to investigate the interplay between the jets and the turbulent convective boundary layer, showing that stronger LLJs are preceded by deeper afternoon mixed layers and often have a later decoupling of mixing between the upper convective mixed layer and the near‐surface layer. Only the strongest LLJs generated a shallow mixing layer overnight. Comparing jet strength and direction to pristine nocturnal convection initiation shows that the strongest southerly LLJs yielded the most pristine nocturnal convection initiation events per night, and the pristine nocturnal convection initiation occurred farther north.

Description: Abstract In this study, the interannual variations in the tropical cyclone (TC) over the western North Pacific (WNP) and the influences of regional sea surface temperature (SST) anomalies are documented by separating the WNP into four quadrants considering nonuniform SST‐induced environmental changes. Our analysis shows that the TC variations in the northwest and southeast quadrants are related to both equatorial central‐eastern Pacific Ocean (EPO) and tropical Indian Ocean (TIO) SST anomalies. The TC variation in the northeast quadrant is mainly related to tropical North Atlantic Ocean SST anomalies. The main environmental variables differ for the TC variations in the four quadrants. Lower‐level (850‐hPa) vorticity is important for the TC variations in the northwest, southwest, and southeast quadrants. Midlevel (700‐hPa) humidity contributes to the TC variations in the northwest, northeast, and southeast quadrants. The vertical shear has a supplementary contribution to the TC variation in the southeast quadrant. The potential intensity (PI) negatively affects the TC variations in the southwest and southeast quadrants. The remote SST anomalies modulate different environmental variables over the WNP. The TIO SST influence is manifested in the lower‐level vorticity and vertical motion. The tropical North Atlantic SST impact occurs through the lower‐level vorticity change. The EPO SST effect occurs via changing the lower‐level vorticity and vertical motion as well as the midlevel moisture and vertical shear. The environmental variables experience more prominent changes when SST anomalies coexist in two remote regions. Numerical experiments confirm the EPO and TIO SST anomaly impacts on the environmental conditions affecting the WNP TC variations.

Description: Abstract An anomalous “north‐south” dipole mode of the snow water equivalent (SWE) persisting from winter to spring is detected over the Eurasian mid‐to‐high latitudes in this study. Using observational data sets and numerical experiments of the Community Atmospheric Model (5.0), we show that this mode contributes to prolonged winter‐springtime coldness in midlatitude Eurasia and is closely linked to the declining November Arctic sea ice concentration. The decline in the sea ice concentration over the Barents‐Laptev Seas can induce a teleconnection pattern over the mid‐to‐high latitudes in the following winter, accompanied by an anomalous ridge over the Ural Mountains and an anomalous trough over Europe and East Asia. Such changes in the large‐scale circulation lead to more cold surges and heavy snowfall in the midlatitudes and light snowfall in the high latitudes, forming an anomalous north‐south dipole mode of the SWE, which further reduces the temperature through thermodynamic feedback. Due to seasonal memory, this SWE pattern can persist into the following spring and can lead to springtime midlatitude coldness via thermodynamic and dynamic processes. For the thermodynamic process, the anomalous SWE condition can lead to anomalous wet soil, reduced incoming surface solar radiation, and cooling air in the midlatitudes. This phenomenon induces an enhanced Siberian High and a deepened East Asian trough via the snow‐Siberian high‐feedback mechanism, which favors a cold spring in northern East Asia. Further analysis suggests that an empirical seasonal prediction model based on the SWE can reasonably predict East Asian spring temperature anomalies.

Description: Abstract A complete and quantitative understanding of cumulus entrainment remains elusive, in part due to the difficulty of directly observing cloud entrainment rates. Multiple approaches to ground‐based observational retrieval of bulk fractional entrainment rates (ε) within cumuli have been developed, such as the parcel model by Jensen and Del Genio (JDG, 2006, https://doi.org/10.1175/JCLI3722.1) and Entrainment Rate In Cumulus Algorithm (ERICA) by Wagner et al. (2013, https://doi.org/10.1175/JTECH-D-12-00187.1). In this paper, a new cumulus entrainment retrieval based on a turbulent kinetic energy (TKE) similarity theory is presented. This method estimates ε based on only the environmental and subcloud conditions. By conducting large‐eddy simulations of a range of continental and maritime shallow cumulus convection cases as Observing System Simulations Experiments, the first numerical verification of the three retrieval methods is produced. These simulations consider a broad range of shallow cumulus environments along with variations of the numerical configuration. The diagnosed ε from these simulations is found to be robustly larger in cumuli over the ocean than in cumuli over land. For continental cumuli, the experiments also reveal a diurnal cycle with increasing ε in the late afternoon. These diagnosed ε serve as the “truth” against which the pseudo‐retrieved entrainment rates from several different implementations of each retrieval are verified. Overall, the simpler JDG and TKE retrievals outperform the more sophisticated ERICA method and better capture the sensitivity to continentality. Only the TKE method reproduces the diurnal variations in ε within continental cumuli. The mean error in the ε retrievals are between 20% and 30% for the TKE and JDG methods, but 50% for ERICA.

Description: Abstract The present study documents intraseasonal snow cover variations over western Eurasia and associated atmospheric processes using the latest Moderate Resolution Imaging Spectroradiometer/Terra daily snow cover product and National Centers for Environmental Prediction/National Center for Atmospheric Research atmospheric reanalysis. It is found that 9‐ to 30‐day variation dominates total intraseasonal snow cover variations over western Siberia. Composite analysis based on 69 positive snow events over western Siberia reveals that atmospheric circulation anomalies control the 9‐ to 30‐day snow variation over western Siberia. A zonal wave train associated with the North Atlantic Oscillation leads to the development of an anomalous cyclone over western Siberia. The associated anomalous ascending motion, anomalous water vapor convergence, and water vapor increase in the atmosphere provide a favorable condition for snowfall. The snowfall starts when anomalous ascending motion reaches the strongest. The maximum snow cover appears about 1 day after the peak of anomalous descending motion and water vapor flux divergence. The surface air temperature tends to vary out of phase with snow cover over western Siberia. Surface air temperature anomalies over western Siberia are contributed by horizontal advection and diabatic heating. The adiabatic heating has a damping effect in surface air temperature variation.

Description: Abstract With a six‐year (2009–2014) summer climate simulation using the Weather Research and Forecasting model at convection‐permitting resolution (4‐km grid spacing), the effects of microphysics parameterization (MP) schemes on precipitation characteristics are investigated in this study. The convection‐permitting simulations employ three popular MP schemes, namely, Lin (single‐moment bulk MP), Weather Research and Forecasting Single‐Moment 5‐class (one‐moment and mixed‐phased MP), and Thompson (two‐moment and mixed‐phase MP) scheme. By evaluating the simulations against the CMORPH, rain gauge (Station), and ERA‐Interim data, it is found that the convection‐permitting model reproduce well the summer precipitation amount and the associated large‐scale atmospheric circulations, which are insensitive to the choice of MP schemes. The simulations with three MP schemes overestimate the precipitation amount, especially over the Yangtze‐Huaihe River Valley. The overestimations may be due to the systematic biases, and cannot be significantly reduced by using different MP schemes. Moreover, all simulations capture well the major features of precipitation diurnal variations and their transition characteristics, but they significantly overestimate the precipitation frequency while underestimate the precipitation intensity. The analysis on the microphysical hydrometeors shows that the model‐simulated precipitation amount is considerably affected by the vertical profiles of solid hydrometeors, especially the snow and graupel particles. The Thompson scheme creates more snow particles and less graupel than the other schemes, while produces the least precipitation amount that best matches the CMORPH.

Description: Abstract Carbonaceous matter in the atmosphere has an important influence on climate change. Currently, the deposition of carbonaceous matter is one of the largest uncertainties in the climate system. This phenomenon is common in remote regions, such as the Himalayas and Tibetan Plateau. In this study, for the first time, we reported in situ measurements of wet and dry deposition rates of carbonaceous matter at three remote stations: Nam Co, Lulang, and Everest. The results showed that the annual wet deposition rates of water‐insoluble organic carbon (WIOC) and black carbon (BC) were 60.2 and 5.8 mg·m−2·year−1, 330 and 34.6 mg·m−2·year−1, and 47.0 and 2.6 mg·m−2·year−1 at the Nam Co, Lulang, and Everest stations, respectively. Seasonal variations in the wet deposition rates of WIOC and BC were controlled by precipitation amount and their atmospheric concentrations. In addition, the wet scavenging ratios of WIOC and BC at Nam Co Station were close to those observed in other remote areas. The total BC deposition at Nam Co Station (15.3 mg·m−2·year−1) was higher than that from chemical transport models, implying a dominant role of dry deposition of BC in the total deposition at this station and an urgent need to improve the aerosol deposition in models for the Himalayas and Tibetan Plateau. It was found that the deposition rates of carbonaceous matter in the Himalayas and Tibetan Plateau had large spatial variation; thus, high‐resolution models need to be applied in the future.