Large impacts simultaneously reset both the surface age and the magnetization of the entire depth of crust over areas comparable to the final size of the resulting craters. These properties make ...large impact craters (>300 km in diameter) ideal “magnetic markers” for constraining the history of the Martian core dynamo. However, the relationship between crustal magnetization and magnetic field measured in orbit is nonunique, making the measured magnetic field signature of an impact crater only a proxy for the magnetization (or lack thereof) below. Using Monte Carlo Fourier domain modeling of subsurface magnetization, we calculate probability distributions of the magnetic field signatures of partially and completely demagnetized craters. We compare these distributions to measured magnetic field signatures of 41 old impact craters on Mars larger than 300 km in diameter and calculate probabilities of their magnetization state. We compare these probabilities to cratering densities and absolute model ages and in this manner arrive at a robust time history of Martian large‐crater magnetization and hence of the Martian dynamo. We conclude that the most likely scenario was a Mars dynamo active when the oldest detectable basins formed, ceasing before the Hellas and Utopia impacts, between 4.0 and 4.1 Ga (in model age) and not thereafter restarting. The Mars atmosphere was thereafter exposed directly to erosion by the solar wind, significantly altering the path of climate evolution. Further improvements to the history of the Martian dynamo will require better crater age estimates and lower altitude magnetic field data.
Key Points
Magnetization probability distributions are calculated with statistical methods.
The Martian dynamo ceased before the Prometheus impact at a model age of ~4 Ga.
The dynamo was still inactive after the Antoniadi impact (model age of 3.8 Ga).
On the 7 August 2019, a 195 km2 raft of andesitic pumice was produced at 200 m below sea level at an unnamed submarine volcano in the Tonga Islands (Southwest Pacific Ocean). Drifting chiefly ...westward, the raft reached the Fiji Islands on the 19 September. Yachts that crossed the raft as early as 2 days post‐eruption provided an outstanding data set of raft characteristics and pristine samples. Further, exceptional tracking of raft dispersal by satellite images allows us to contrast virtual particle tracking methods with ocean model currents to explore the relative influence of surface currents, wind, and wave action on pumice flotsam dispersal over up to 2 years. Attenuation of ocean waves by large and compact pumice rafts appears to reduce the effect of Stokes drift. The coupling of real‐time satellite observations with oceanographic Lagrangian simulations allows near‐real time forecasting for global maritime hazard mitigation.
Plain Language Summary
Although 70% of volcanism occurs underwater, submarine eruptions are seldom witnessed owing to their remoteness and absence of technologies enabling their detection. On the 7 August 2019, a submarine volcanic eruption 200 m below sea level in the Tonga Islands formed a 195 km2 pumice raft. The raft was first encountered and sampled by yacht crews, and is visible on satellite images. Pumice rafts are dispersed by ocean currents, wind, and waves, and can cross entire oceans or get stranded on coasts. After 7 weeks of dispersal, the pumice raft reached the Fiji Islands, and a fraction continued its westward route toward Vanuatu and eastern Australia. Excellent imaging of the raft by satellites permitted reconstruction of daily raft dispersal for the first 8 weeks. Here we show that drift calculations, including components of ocean current, wind, and wave action, can usefully forecast raft dispersal. We tested and tuned these drift calculations by comparing the simulated drift with daily satellite images of the raft. Further, a combination of satellite images and drift calculations based on oceanographic models were used for maritime hazard mitigation in near‐real time.
Key Points
A >30 million m3 pumice raft formed at an unnamed volcano in Tonga midday on the 7 August 2019 and was dispersed by ocean currents and wind
Forecasts of pumice raft dispersal, based on particle drift with model currents, winds, and waves, are evaluated with daily satellite images
Stokes drift does not improve short‐term predictions on large and compact pumice rafts because these rafts attenuate ocean waves
The effects of earthquakes on groundwater and aquifer properties can be quantified and monitored using water-level changes produced by tides and barometric pressure. Tidal and barometric responses ...are particularly useful in evaluating the impacts of unexpected events, such as earthquakes, because the signals are continuously generated and recorded over large areas of the Earth’s surface. The techniques for the extraction of tidal and barometric signals from the water-level time series are described in many excellent papers, here, we focus on reviewing the hydrogeologic interpretations of, and earthquake impacts on, these responses. We review how hydrogeology and earthquakes impact the groundwater response to Earth tides, and changes in barometric pressure and barometric tides. Next, we review the current understanding of the mechanisms responsible for earthquake-induced changes in aquifer confinement and permeability. We conclude with a summary of open questions and topics for future research, notably the value in long-term monitoring and analysis of the earthquake response at multiple tidal and barometric frequencies.
Hydrologic responses to earthquakes and their mechanisms have been widely studied. Some responses have been attributed to increases in the vertical permeability. However, basic questions remain: How ...do increases in the vertical permeability occur? How frequently do they occur? Is there a quantitative measure for detecting the occurrence of aquitard breaching? We try to answer these questions by examining data from a dense network of ∼50 monitoring stations of clustered wells in a sedimentary basin near the epicenter of the 1999 M7.6 Chi‐Chi earthquake in western Taiwan. While most stations show evidence that confined aquifers remained confined after the earthquake, about 10% of the stations show evidence of coseismic breaching of aquitards, creating vertical permeability as high as that of aquifers. The water levels in wells without evidence of coseismic breaching of aquitards show tidal responses similar to that of a confined aquifer before and after the earthquake. Those wells with evidence of coseismic breaching of aquitards, on the other hand, show distinctly different postseismic tidal response. Furthermore, the postseismic tidal response of different aquifers became strikingly similar, suggesting that the aquifers became hydraulically connected and the connection was maintained many months thereafter. Breaching of aquitards by large earthquakes has significant implications for a number of societal issues such as the safety of water resources, the security of underground waste repositories, and the production of oil and gas. The method demonstrated here may be used for detecting the occurrence of aquitard breaching by large earthquakes in other seismically active areas.
Key Points:
Earthquakes can enhance vertical permeability and by disrupt the confinement of aquifers
∼10% of wells in the near field show increased permeability that persists for months
Impacts on the safety of water resources and underground waste repositories
The cause of the Lusi mud eruption remains controversial. The review by Miller and Mazzini (2017) firmly dismisses a role of drilling operations at the adjacent Banjarpanji-1 well and argues that the ...eruption was triggered by the M6.3 Yogyakarta earthquake some 254 km away. We disagree with these conclusions. We review drilling data and the daily drilling reports, which clearly confirm that the wellbore was not intact and that there was a subsurface blowout. Downhole pressure data from Lusi directly witness the birth of Lusi at the surface on the 29th of May 2006, indicating a direct connection between the well and the eruption. Furthermore, the daily drilling reports specifically state that Lusi activity was visibly altered on three separate occasions by attempts to kill the eruption by pumping dense fluid down the BJP-1 well, providing further evidence of a connection between the wellbore and Lusi. By comparison with other examples of newly initiated mud eruptions elsewhere by other earthquakes, the Yogyakarta earthquake was far away given its magnitude. The seismic energy density of the Yogyakarta earthquake was only 0.0043 J/m3, which is less than a quarter of the minimum 0.019 J/m3 seismic energy density that has ever been inferred to trigger other mud eruptions. We show that the Lusi area had previously experienced other shallow earthquakes with similar frequencies and stronger ground shaking that did not trigger an eruption. Finally, the data from the BJP-1 well indicates that there was no prior hydrodynamic connection between deep overpressured hydrothermal fluids and the shallow Kalibeng clays, and that there was no evidence of any liquefaction or remobilization of the Kalibeng clays induced by the earthquake. We thus strongly favor initiation by drilling and not an earthquake.
•Cause of the Lusi eruption remains controversial, with proposals that it was triggered by either an earthquake or drilling.•Drilling reports and data confirm that the wellbore was not intact, there was a blowout, and a direct surface connection.•We show that the Yogyakarta earthquake was too far away given its magnitude to have initiated a new eruption.•We strongly favor initiation of the eruption by drilling and not an earthquake.
The development of discrete volcanic centers reflects a focusing of magma ascending from the source region to the surface. We suggest that this organization occurs via mechanical interactions between ...magma chambers, volcanic edifices, and dikes and that the stresses generated by these features may localize crustal magma transport before the first eruption occurs. We develop a model for the focusing or “lensing” of rising dikes by magma chambers beneath a free surface, and we show that chambers strongly modulate dike focusing by volcanic edifices. We find that the combined mechanical effects of chambers, edifice loading, and dike propagation are strongly coupled. Chambers deeper than ∼20 km below the surface with magmatic overpressure in the range of 20–100 MPa should dominate dike focusing, while more shallow systems are affected by both edifice and chamber focusing.
Volcanic clasts in many pyroclastic density current deposits are notably more round than their counterparts in corresponding fall deposits. This increase in roundness and sphericity reflects ...different degrees of comminution, abrasion and breakup during transport. We performed experimental measurements to determine an
empirical
relationship between particle shape and mass loss caused by particle–particle interactions. We consider, as examples, pumice from four volcanoes: Medicine Lake, California; Lassen, California; Taupo, New Zealand; Mount St Helens, Washington. We find that average sample roundness reaches a maximum value once particles lose between 15% and 60% of their mass. The most texturally homogeneous clasts (Taupo) become the most round. Crystal-rich pumice abrades more slowly than crystal-free pumice of similar density. Abrasion rates also decrease with time as particles become less angular. We compare our experimental measurements with the shapes of clasts in one of the May 18, 1980 pyroclastic density current units at Mount St Helens, deposited 4–8 km from the vent. The measured roundness of these clasts is close to the experimentally determined maximum value. For a much smaller deposit from the 1915 Lassen eruption, clast roundness is closer to the value for pumice in fall deposits and suggests that only a few volume percent of material was removed from large clasts. In neither field deposit do we see a significant change in roundness with increasing distance from the vent. We suggest that this trend is recorded because much of the rounding and ash production occur in proximal regions where the density currents are the most energetic. As a result, all clasts that are deposited have experienced similar amounts of comminution in the proximal region, and similar amounts of abrasion as they settle through the dense, near-bed region prior to final deposition.
Young flood channels emanate from Cerberus Fossae, Mars. In order to determine whether subsurface aquifers can discharge water sufficiently fast to explain these flood features, I develop a model ...that couples groundwater flow in a sub‐cryosphere aquifer with flow through a fissure that penetrates from the surface to the aquifer. The model is constrained by estimates of peak discharge, and the distance water flowed uphill from the fissure. Aquifer permeabilities similar to those of young basalt aquifers on the Earth, ∼10−9 m2, can produce inferred discharges of 106 m3/s Burr et al., 2002a but probably require multiple floods to create the channels.