This paper describes a method for analyzing the kinematic properties of ions composing the solar wind. The core technology is a velocity analysis performed by dual rotating electric field (REF) units arranged coaxially in tandem, where the electric field in the downstream unit is set in the opposite direction to the upstream one. When the solar wind flies freely through the REF units, ions diverge outwards in the upstream unit and converge inwards in the downstream unit. Since the degree of diversion and conversion correspond to each ion’s velocity, ions separate into multiple groups flying through the REF units, terminate their flights on an image sensor placed on the tail end, and create sorted patterns that exhibit the velocity distribution of ions. As the REF units act on ions only by dynamic lateral electric force, the initial velocity and charge state of the ions remain invariant during the analysis process, which can be advantageous for solar wind analysis. This paper introduces a proposed instrument equipped with multiple functions, including an energy analysis, an overall velocity analysis, and a detailed ion velocity analysis. The kinematic properties of the measured ions correspond to energy levels ranging from 300 eV to 20 keV, velocities from 20 km/s to 1,900 km/s, and mass numbers from 1 to 200.

Understanding the interaction between tectonic plates from geodetic data is relevant to the assessment of seismic hazard. To shed light on that prevalently slow aseismic interaction, we developed a new static-slip inversion strategy, the ELADIN (ELastostatic ADjoint INversion) method, that uses the adjoint elastostatic equations to compute the gradient of the cost function. To handle plausible slip constraints, ELADIN is a 2-step inversion algorithm. First it finds the slip that best explains the data without any constraint, and then refines the solution by imposing the constraints through a Gradient Projection Method. To obtain a selfsimilar, physically-consistent slip distribution that accounts for sparsity and uncertainty in the data, ELADIN reduces the model space by using a von Karman regularization function that controls the wavenumber content of the solution, and weights the observations according to their covariance using the data precision matrix. Since crustal deformation is the result of different concomitant interactions at the plate interface, ELADIN simultaneously determines the regions of the interface subject to both stressing (i.e., coupling) and relaxing slip regimes. For estimating the resolution, we introduce a mobile checkerboard that allows to determine lower-bound fault resolution zones for an expected slip-patch size and a given stations array. We systematically test ELADIN with synthetic inversions along the whole Mexican subduction zone and use it to invert the 2006 Guerrero Slow Slip Event (SSE), which is one of the most studied SSEs in Mexico. Since only 12 GPS stations recorded the event, careful regularization is thus required to achieve reliable solutions. We compared our preferred slip solution with two previously published models and found that our solution retains their most reliable features. In addition, although all three SSE models predict an upward slip penetration invading the seismogenic zone of the Guerrero seismic gap, our resolution analysis indicates that this penetration might not be a reliable feature of the 2006 SSE.

North African dust is known to be deposited in the Gulf of Mexico, but its deposition rate and associated supply of lithogenic dissolved metals, such as the abiotic metal thorium or the micronutrient metal iron, have not been well-quantified. 232Th is an isotope with similar sources as iron and its input can be quantified using radiogenic 230Th. By comparing dissolved 232Th fluxes at three sites in the northern Gulf of Mexico with upwind sites in the North Atlantic, we place an upper bound on North African dust contributions to 232Th and Fe in the Gulf of Mexico, which is about 30% of the total input. Precision on this bound is hindered by uncertainty in the relative rates of dust deposition in the North Atlantic and the northern Gulf of Mexico. Based on available radium data, shelf sources, including rivers, submarine groundwater discharge and benthic sedimentary releases are likely as important if not more important than dust in the budget of lithogenic metals in the Gulf of Mexico. In other words, it is likely there is no one dominant source of Th and Fe in the Gulf of Mexico. Finally, our estimated Fe input in the northern Gulf of Mexico implies an Fe residence time of less than 6 months, similar to that in the North Atlantic despite significantly higher supply rates in the Gulf of Mexico.

Corridors of size-selected crescent-shaped dunes, known as barchans, are commonly found in water, air, and other planetary environments. The growth of barchans results from the interplay between a fluid flow and a granular bed, but their size regulation involves intricate exchanges between different barchans within a field. One size-regulating mechanism is the binary interaction between nearby dunes, when two dunes exchange mass via the near flow field or by direct contact (collision). In a recent Letter (Assis & Franklin, GRL, 2020), we identified five different patterns arising from binary interactions of subaqueous barchans, and proposed classification maps. In this paper, we further inquire into binary exchanges by investigating the motion of individual grains while barchans interact with each other. The experiments were conducted in a water channel where the evolution of pairs of barchans in both aligned and off-centered configurations was recorded by conventional and high-speed cameras. Based on image processing, we obtained the morphology of dunes and motion of grains for all interaction patterns. We present the trajectories of individual grains, from which we show the origin and destination of moving grains, and their typical lengths and velocities. We also show that grains from the impacting dune spread with a diffusion-like component over the target barchan, and we propose a diffusion length. Our results provide new insights into the size-regulating mechanisms of barchans and barchanoid forms found on Earth and other planets.

Hydroperiod, or the amount of time a lentic waterbody contains water, shapes communities of aquatic organisms. Precise measurement of hydroperiod features such as inundation timing and duration can help predict community dynamics and ecosystem stability. In areas defined by high spatial and temporal variability, fine-scale temporal variation in inundation timing and duration may drive community structure, but that variation may not be captured using common approaches including remote sensing technology. Here, we provide methods to accurately capture inundation timing by fitting hidden Markov models to measurements of daily temperature standard deviation collected from temperature loggers. We describe a rugged housing design to protect loggers from physical damage and apply our methods to a group of intermittent ponds in southeastern Arizona, showing that initial pond inundation timing is highly variable across a small geographic scale (~50km2). We also compare a 1-logger (pond only) and 2-logger (pond + control) design and show that, although a single logger may be sufficient to capture inundation timing in most cases, a 2-logger design can increase confidence in results. These methods are cost-effective and show promise in capturing variation in intraregional inundation timing that may have profound effects on aquatic communities, with implications for how these communities may respond to hydroperiod alteration from a changing climate.

Ionospheric perturbations induced by tsunamis and earthquakes can be used for tsunami early warning and remote sensing of earthquakes, provided the perturbations are characterized properly to distinguish them from the ones caused by other sources. The ionospheric perturbations are increasingly being obtained from Global Positioning System (GPS) based Total Electron Content (TEC) measurements sampled at uniform time intervals. However, the sampling is not uniform in space. The nonuniform spatial sampling along the GPS satellite tracks introduces aliasing if it is not accounted while computing the ionospheric perturbations. All the methods hitherto used to detect the co-seismic and tsunamigenic ionospheric perturbations did not account the nonuniform spatial sampling while computing these perturbations. In addition, the residual approach used to obtain the perturbations by detrending the TEC time series using high-order polynomial fit introduces artifacts. These aliasing and artifacts corrupt amplitude, Signal-to-Noise Ratio (SNR), phase, and frequency of ionospheric perturbations which are vital to distinguish the perturbations induced by tsunamis and earthquakes from the rest. We show that Spatio-Periodic Leveling Algorithm (SPLA) successfully removes such aliasing and artifacts. The efficiency of SPLA in removing the aliases and artifacts is validated under two simulated scenarios, and using GPS observations carried out during two natural disasters – the 2004 Indian Ocean tsunami and the 2015 Nepal-Gorkha earthquake. We, further, studied the severity of aliasing and artifacts on co-seismic and tsunamigenic perturbations by analyzing its characteristics employing SNR, spatiotemporal, and wavelet analyses. The results reveal that removal of aliasing and artifacts using SPLA i) increases the SNR up to ~149% compared to the residual method and ~39% compared to the differential method, ii) distinctly resolves signals from sharp static variations, and iii) detects 50% more co-seismic ionospheric perturbations and 25% more tsunami-induced ionospheric perturbations in the two events studied. Cross-correlation of the perturbation time series obtained using the residual method and SPLA reveals that aliasing and artifacts shift the time of occurrence by -7.64 minutes to +4.21 minutes. Further, the results show that the SPLA efficiently detects the ionospheric perturbations at low elevation angles, thereby removes the need of applying elevation cut-off and increases the area of ionospheric exploration of a GPS receiver.

Changes in the concentration of greenhouse gases within the atmosphere lead to changes in radiative fluxes throughout the atmosphere. The value of this change, called the instantaneous radiative forcing, varies across climate models, due partly to differences in the distribution of clouds, humidity, and temperature across models, and partly due to errors introduced by approximate treatments of radiative transfer. This paper describes an experiment within the Radiative Forcing Model Intercomparision Project that uses benchmark calculations made with line-by-line models to identify parameterization error in the representation of absorption and emission by greenhouse gases. The clear-sky instantaneous forcing by greenhouse gases to which the world has been subject is computed using a set of 100 profiles, selected from a re-analysis of present-day conditions, that represent the global annual mean forcing with sampling errors of less than 0.01 \si{\watt\per\square\meter}. Six contributing line-by-line models agree in their estimate of this forcing to within 0.025 \si{\watt\per\square\meter} while even recently-developed parameterizations have typical errors four or more times larger, suggesting both that the samples reveal true differences among line-by-line models and that parameterization error will be readily resolved. Agreement among line-by-line models is better in the longwave than in the shortwave where differing treatments of the water vapor vapor continuum affect estimates of forcing by carbon dioxide and methane. The impacts of clouds on instantaneous radiative forcing are roughly estimated, as are adjustments due to stratospheric temperature change. Adjustments are large only for ozone and for carbon dioxide, for which stratospheric cooling introduces modest non-linearity.

Previous studies have interpreted Last Interglacial (LIG; ~129-116 ka) sea-level estimates in multiple different ways to calibrate projections of future Antarctic ice-sheet (AIS) mass loss and associated sea-level rise. This study systematically explores the extent to which LIG constraints could inform future Antarctic contributions to sea-level rise. We develop a Gaussian process emulator of an ice-sheet model to produce continuous probabilistic projections of Antarctic sea-level contributions over the LIG and a future high-emissions scenario. We use a Bayesian approach conditioning emulator projections on a set of LIG constraints to find associated likelihoods of model parameterizations. LIG estimates inform both the probability of past and future ice-sheet instabilities and projections of future sea-level rise through 2150. Although best-available LIG estimates do not meaningfully constrain Antarctic mass loss projections or physical processes until 2060, they become increasingly informative over the next 130 years. Uncertainties of up to 50 cm remain in future projections even if LIG Antarctic mass loss is precisely known (+/-5 cm), indicating there is a limit to how informative the LIG could be for ice-sheet model future projections. The efficacy of LIG constraints on Antarctic mass loss also depends on assumptions about the Greenland ice sheet and LIG sea-level chronology. However, improved field measurements and understanding of LIG sea levels still have potential to improve future sea-level projections, highlighting the importance of continued observational efforts.

We demonstrate scintillation analysis from a network of geodetic Global Positioning System (GPS) receivers which provide data at 1-second resolution. We introduce proxy phase ($\sigma_{TEC}$) and amplitude ($SNR_4$) scintillation indices and validate them against the rate of change of TEC index (ROTI), and $S_4$. Additionally, we validate scintillation observations against a CASES scintillation receiver. We develop receiver dependent scintillation event thresholding using hardware-dependent noise variance. We analyze six-days adjacent to the 7-8 September 2017 geomagnetic storm, using 169 receivers covering magnetic latitudes between 15$^\circ$ and 65$^\circ$ in the American longitude sector. We leverage the available spatial sampling coverage to construct 2-D maps of scintillation and present episodic evolution of scintillation intensifications during the storm. We show that low-latitude and high-latitude scintillation morphology match well-established scintillation climatology patterns. At mid-latitudes, spatiotemporal evolution of scintillation partially agrees with known scintillation patterns. Additionally, the results reveal previously undocumented mid-latitude scintillation producing structures. The results provide an unprecedented view into the spatiotemporal development of scintillation-producing plasma irregularities and provide a resource to further exploit scintillation evolution at large spatial scales.

The Kelvin-Helmholtz instability (KHI) at Earth’s magnetopause and associated turbulence are suggested to play a role in the transport of mass and momentum from the solar wind into Earth’s magnetosphere. We investigate electromagnetic turbulence observed in KH vortices encountered at the dusk flank magnetopause by the Magnetospheric Multiscale (MMS) spacecraft under northward interplanetary magnetic field (IMF) conditions in order to reveal its generation process, mode properties, and role. A comparison with another MMS event at the dayside magnetopause with reconnection but no KHI signatures under a similar IMF condition indicates that while high-latitude magnetopause reconnection excites a modest level of turbulence in the dayside low-latitude boundary layer, the KHI further amplifies the turbulence, leading to magnetic energy spectra with a power-law index –5/3 at magnetohydrodynamic scales even in its early nonlinear phase. The mode of the electromagnetic turbulence is analyzed with a single-spacecraft method based on Ampère’s law, developed by Bellan (2016), for estimating wave vectors as a function of spacecraft-frame frequency. The results suggest that the turbulence does not consist of propagating normal-mode waves, but is due to interlaced magnetic flux tubes advected by plasma flows in the vortices. The turbulence at sub-ion scales in the early nonlinear phase of the KHI may not be the cause of the plasma transport across the magnetopause, but rather a consequence of three-dimensional vortex induced reconnection, the process that can cause an efficient transport by producing tangled reconnected field lines.

The southern Granite-Rhyolite province contains a comprehensive record of lithospheric evolution in North America. During the last decade, increased seismicity along with improved seismic monitoring installations in Oklahoma provided a rich catalog of local earthquakes. The source-receiver geometry of this dataset is well posed to illuminate the middle and lower crust through long offset recordings of the Pg phase. We present a 3-D P-wave velocity model for central and north Oklahoma developed through a non-standard processing scheme applied to local earthquake waveforms recorded from 2010-2017, focusing on the deeper crust. We employed common-mid-point sorting, stacking, and inversion of Pg-phases which resulted in a set of localized velocity-depth functions up to depths of 40 km. Using this methodology, we significantly increased the S/N ratio for far offset (~250 km) local earthquake waveforms which led to the increase in depth of investigation for our final 3-D velocity model. We find high velocity (> 7 km/s) lower crust throughout the investigated area which suggests a mafic lower crust. The high velocities support previously established models which state that the lower crust of the Granite-Rhyolite province was derived from melting of older crust. We further relate shallow and middle crustal velocity anomalies to other data sets such as gravimetric and magnetic anomalies, and the spatial distribution of earthquakes. We interpret the Nemaha Fault system as a deep-rooted discontinuity which separates two crustal domains. On the contrary, we do not find clear evidence for the existence of the Midcontinent rift (MCR) in northern Oklahoma.

The Lower St. Lawrence Seaway (LSLS), in eastern Canada, is an important habitat for several species of endangered baleen whale. As we seek to reduce the hazards that these endangered species face from human activity, there is increasing demand for detailed knowledge of their habitat use. Only a sparse network of hydrophones exists in the LSLS to remotely observe whales. However, there is also a network of onshore seismometers, designed to monitor earthquakes, that have sufficiently high sample rates to record fin and blue whale calls. We present a simple method for detecting band-limited, regularly repeating calls, such as the 20 Hz calls of fin and blue whales, and apply the method to build a catalog of fin and blue whale detections at 14 onshore seismometers across the LSLS, over approximately a four-year period. The resulting catalog contains >600000 fin whale calls and >60000 blue whale calls. Individual calls are rarely detected at more than one seismometer. Fin whale calls recorded onshore appear to travel mainly through solid earth, rather than only entering the earth at the shoreline, and they often have a complex ~2 s sequence of P-like and S-like phases. Onshore seismometers provide a valuable, previously unused source of data for monitoring baleen whales. However, in the LSLS, the current seismometer network cannot provide high-precision whale tracking alone, so a denser deployment of onshore and/or offshore seismometers is required.

To help understand and determine the driver of jovian auroral X-rays, we present the first statistical study to focus on the morphology and dynamics of the jovian northern hot spot (NHS) using Chandra data. The catalogue we explore dates from 18 December 2000 up to and including 8 September 2019. Using a numerical criterion, we characterize the typical and extreme behaviour of the concentrated NHS emissions across the catalogue. The mean power of the NHS is found to be 1.91 GW with a maximum brightness of 2.02 Rayleighs (R), representing by far the brightest parts of the jovian X-ray spectrum. We report a statistically significant region of emissions at the NHS center which is always present, the averaged hot spot nucleus (AHSNuc), with mean power of 0.57 GW and inferred average brightness of ∼ 1.2 R. We use a flux equivalence mapping model to link this distinct region of X-ray output to a likely source location and find that the majority of mappable NHS photons emanate from the pre-dusk to pre-midnight sector, coincident with the dusk flank boundary. A smaller cluster maps to the noon magnetopause boundary, dominated by the AHSNuc, suggesting that there may be multiple drivers of X-ray emissions. On application of timing analysis techniques (Rayleigh, Monte Carlo, Jackknife), we identify several instances of statistically significant quasi-periodic oscillations (QPOs) in the NHS photons ranging from ∼ 2.3-min to 36.4-min, suggesting possible links with ultra-low frequency activity on the magnetopause boundary (e.g. dayside reconnection, Kelvin-Helmholtz instabilities).

A thermodynamic model to calculate the “Sulfide Content at Sulfide Saturation” or SCSS of basaltic and intermediate composition silicate melts has been built from four independently measurable thermodynamic entities, namely the standard state Gibbs free energy of the saturation reaction, the “sulfide capacity”, and the activities of FeO in the silicate melt and of FeS in the coexisting sulfide: ln [S2-]SCSS = ∆G(FeO-FeS)/RT + ln C(S2- )- ln a(FeO)(sil melt) + ln a(FeS)sulf The model was calibrated for silicate melts of basic and intermediate composition from published experimental results as a function of temperature, silicate melt composition, and sulfide matte composition in the system Fe-Ni-Cu-S-O at 1 bar. The likely effects of pressure and H2O content on SCSS were included in an exploratory way. The model was tuned against the large dataset of S contents in OFB glasses of Jenner and O’Neill (2012), giving it a precision comparable to that of the S analyses themselves, which is ~ 5%. All but 3% the OFB glasses were found to be sulfide saturated within uncertainty; these 3% have lost S by devolatization, revealed by their low S/Se. Applying the model to other OFB datasets suggests sulfide saturation is ubiquitous, including olivine-hosted melt inclusions proposed previously to be sulfide undersaturated. The sulfur fugacity (fS2) of undegassed Ocean Floor Basalts varies proportionally to fO2, with log10fS2 typically within the range -0.6 to +0.4.

Advances in machine learning (ML) techniques and computational capacity have yielded state-of-the-art methodologies for processing, sorting, and analyzing large seismic data sets. In this work, we consider an application of ML for automatically identifying dominant types of impulsive seismicity contained in observations from a 34-station broadband seismic array deployed on the Ross Ice Shelf (RIS), Antarctica from 2014 to 2017. The RIS seismic data contain signals and noise generated by many glaciological processes that are useful for monitoring the integrity and dynamics of ice shelves. Deep clustering was employed to efficiently investigate these signals. Deep clustering automatically groups signals into hypothetical classes without the need for manual labeling, allowing for comparison of their signal characteristics and spatial and temporal distribution with potential source mechanisms. The method uses spectrograms as input and encodes their salient features into a lower-dimensional latent representation using an autoencoder, a type of deep neural network. For comparison, two clustering methods are applied to the latent data: a Gaussian mixture model (GMM) and deep embedded clustering (DEC). Eight classes of dominant seismic signals were identified and compared with environmental data such as temperature, wind speed, tides, and sea ice concentration. The greatest seismicity levels occurred at the RIS front during the 2016 El Niño summer, and near grounding zones near the front throughout the deployment. We demonstrate the spatial and temporal association of certain classes of seismicity with seasonal changes at the RIS front, and with tidally driven seismicity at Roosevelt Island.

Meaningful analysis of uranium-series isotopic disequilibria in basaltic lavas relies on the use of complex forward numerical models like dynamic melting (McKenzie, 1985) and equilibrium porous flow (Spiegelman and Elliott, 1993). Historically, such models have either been solved analytically for simplified scenarios, such as constant melting rate or constant solid/melt trace element partitioning throughout the melting process, or have relied on incremental or numerical calculators with limited power to solve problems and/or restricted availability. The most public numerical solution to reactive porous flow, UserCalc (Spiegelman, 2000) was maintained on a private institutional server for nearly two decades, but that approach has been unsustainable in light of modern security concerns. Here we present a more long-lasting solution to the problems of availability, model sophistication and flexibility, and long-term access in the form of a cloud-hosted, publicly available Jupyter notebook. Similar to UserCalc, the new notebook calculates U-series disequilibria during time-dependent, equilibrium partial melting in a one-dimensional porous flow regime where mass is conserved. In addition, we also provide a new disequilibrium transport model which has the same melt transport model as UserCalc, but approximates rate-limited diffusive exchange of nuclides between solid and melt using linear kinetics. The degree of disequilibrium during transport is controlled by a Damköhler number, allowing the full spectrum of equilibration models from complete fractional melting (Da = 0) to equilibrium transport (Da = ∞).

The Weddell Gyre’s variability on seasonal and interannual timescales is investigated using an ocean-sea ice model at three different horizontal resolutions. The model is evaluated against available observations to demonstrate that the highest resolution configuration (0.1$^{\circ}$ in the horizontal) best reproduces observed features of the region. The simulations suggest that the gyre is subject to large variability in its circulation that is not captured by summer-biased or short-term observations. The Weddell Gyre’s seasonal cycle consists of a summer minimum and a winter maximum and accounts for changes that are between one third and a half of its mean transport . On interannual time scales we find that the gyre’s strength is correlated with the local Antarctic easterlies and that extreme events of gyre circulation are associated with changes in sea ice concentration and the characteristics of warm inflow at the eastern boundary.

Wind erosion is integral to the evolution of arid landscapes on Earth and Mars, but the nature of long-term wind erosion of bedrock is poorly understood. Here we describe the Salina del Fraile (SdF) depression in the Puna Plateau of the Central Andes, NW Argentina, as a landform excavated by wind over several million years. New structural cross-sections and a compilation of chronostratigraphic ages rule out the hypothesis that the depression was created by transtensional tectonics. Dated remnant lacustrine and alluvial deposits in the floor of the depression constrain the rate and timing of erosion. Late Oligocene–Miocene compressional folding uplifted friable strata that were preferentially eroded, resulting in the high-relief (900 m) depression. Up to 1.95 km and an average of 1.05 km of strata were eroded during the last 8.2 to 17 Ma, at rates of 0.06 to 0.23 mm/yr. These rates are similar to long-term average wind erosion rates reported in other regions. Coarse-grained eolian megaripples, yardangs, and elongated ridges indicate ongoing eolian abrasion and deflation, aided by salt weathering, of the floor of the depression. Megaripple migration across stony lag surfaces exposes fresh bedrock to continued erosion. The SdF also contains kilometerscale mesas and ridges that we interpret as erosional remnants. These landforms are similar to megayardangs and erosional topography identified on the lower flanks of Mount Sharp, Gale crater, Mars. In such hyperarid landscapes characterized by lithologic heterogeneities, high relief landforms can be generated and sustained by wind erosion, without significant fluvial or glacial incision.