USGS Science⌕ Search

SEARCH · USGS Science

Results for “Episodes”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32Linked to original sources

Controls on methane released through ebullition in peatlands affected by permafrost degradation

Permafrost thaw in peat plateaus leads to the flooding of surface soils and the formation of collapse scar bogs, which have the potential to be large emitters of methane (CH 4 ) from surface peat as well as deeper, previously frozen, permafrost carbon (C). We used a network of bubble traps, permanently installed 20 cm and 60 cm beneath the moss surface, to examine controls on ebullition from three collapse bogs in interior Alaska. Overall, ebullition was dominated by episodic events that were associated with changes in atmospheric pressure, and ebullition was mainly a surface process regulated by both seasonal ice dynamics and plant phenology. The majority (>90%) of ebullition occurred in surface peat layers, with little bubble production in deeper peat. During periods of peak plant biomass, bubbles contained acetate-derived CH 4 dominated (>90%) by modern C fixed from the atmosphere following permafrost thaw. Post-senescence, the contribution of CH 4 derived from thawing permafrost C was more variable and accounted for up to 22% (on average 7%), in the most recently thawed site. Thus, the formation of thermokarst features resulting from permafrost thaw in peatlands stimulates ebullition and CH 4 release both by creating flooded surface conditions conducive to CH 4 production and bubbling as well as by exposing thawing permafrost C to mineralization.

Alaska↗

Characteristic length scales and time-averaged transport velocities of suspended sediment in the mid-Atlantic Region, USA

Watershed Best Management Practices (BMPs) are often designed to reduce loading from particle-borne contaminants, but the temporal lag between BMP implementation and improvement in receiving water quality is difficult to assess because particles are only moved downstream episodically, resting for long periods in storage between transport events. A theory is developed that describes the downstream movement of suspended sediment particles accounting for the time particles spend in storage given sediment budget data (by grain size fraction) and information on particle transit times through storage reservoirs. The theory is used to define a suspended sediment transport length scale that describes how far particles are carried during transport events, and to estimate a downstream particle velocity that includes time spent in storage. At 5 upland watersheds of the mid-Atlantic region, transport length scales for silt-clay range from 4 to 60 km, while those for sand range from 0.4 to 113 km. Mean sediment velocities for silt-clay range from 0.0072 km/yr to 0.12 km/yr, while those for sand range from 0.0008 km/yr to 0.20 km/yr, 4–6 orders of magnitude slower than the velocity of water in the channel. These results suggest lag times of 100–1000 years between BMP implementation and effectiveness in receiving waters such as the Chesapeake Bay (where BMPs are located upstream of the characteristic transport length scale). Many particles likely travel much faster than these average values, so further research is needed to determine the complete distribution of suspended sediment velocities in real watersheds.

Maryland, Pennsylvania, Virginia↗

The emergence of volcanic oceanic islands on a slow-moving plate: The example of Madeira Island, NE Atlantic

The transition from seamount to oceanic island typically involves surtseyan volcanism. However, the geological record at many islands in the NE Atlantic—all located within the slow-moving Nubian plate—does not exhibit evidence for an emergent surtseyan phase but rather an erosive unconformity between the submarine basement and the overlying subaerial shield sequences. This suggests that the transition between seamount and island may frequently occur by a relative fall of sea level through uplift, eustatic changes, or a combination of both, and may not involve summit volcanism. In this study, we explore the consequences for island evolutionary models using Madeira Island (Portugal) as a case study. We have examined the geologic record at Madeira using a combination of detailed fieldwork, biostratigraphy, and 40 Ar/ 39 Ar geochronology in order to document the mode, timing, and duration of edifice emergence above sea level. Our study confirms that Madeira's subaerial shield volcano was built upon the eroded remains of an uplifted seamount, with shallow marine sediments found between the two eruptive sequences and presently located at 320–430 m above sea level. This study reveals that Madeira emerged around 7.0–5.6 Ma essentially through an uplift process and before volcanic activity resumed to form the subaerial shield volcano. Basal intrusions are a likely uplift mechanism, and their emplacement is possibly enhanced by the slow motion of the Nubian plate relative to the source of partial melting. Alternating uplift and subsidence episodes suggest that island edifice growth may be governed by competing dominantly volcanic and dominantly intrusive processes.

Geochemistry, Geophysics, Geosystems↗

Ambient tremors in a collisional orogenic belt

Deep-seated tectonic tremors have been regarded as an observation tied to interconnected fluids at depth, which have been well documented in worldwide subduction zones and transform faults but not in a collisional mountain belt. In this study we explore the general features of collisional tremors in Taiwan and discuss the possible generation mechanism. In the 4 year data, we find 231 ambient tremor episodes with durations ranging from 5 to 30 min. In addition to a coseismic slip-induced stress change from nearby major earthquake, increased tremor rate is also highly correlated with the active, normal faulting earthquake swarms at the shallower depth. Both the tremor and earthquake swarm activities are confined in a small, area where the high attenuation, high thermal anomaly, the boundary between high and low resistivity, and localized veins on the surfaces distributed, suggesting the involvement of fluids from metamorphic dehydration within the orogen.

Geophysical Research Letters↗

Explosive dome eruptions modulated by periodic gas-driven inflation

Volcan Santiaguito (Guatemala) “breathes” with extraordinary regularity as the edifice's conduit system accumulates free gas, which periodically vents to the atmosphere. Periodic pressurization controls explosion timing, which nearly always occurs at peak inflation, as detected with tiltmeters. Tilt cycles in January 2012 reveal regular 26 ± 6 min inflation/deflation cycles corresponding to at least ~10 1 kg/s of gas fluxing the system. Very long period (VLP) earthquakes presage explosions and occur during cycles when inflation rates are most rapid. VLPs locate ~300 m below the vent and indicate mobilization of volatiles, which ascend at ~50 m/s. Rapid gas ascent feeds pyroclast-laden eruptions lasting several minutes and rising to ~1 km. VLPs are not observed during less rapid inflation episodes; instead, gas vents passively through the conduit producing no infrasound and no explosion. These observations intimate that steady gas exsolution and accumulation in shallow reservoirs may drive inflation cycles at open-vent silicic volcanoes.

Geophysical Research Letters↗

Laboratory constraints on models of earthquake recurrence

In this study, rock friction ‘stick-slip’ experiments are used to develop constraints on models of earthquake recurrence. Constant-rate loading of bare rock surfaces in high quality experiments produces stick-slip recurrence that is periodic at least to second order. When the loading rate is varied, recurrence is approximately inversely proportional to loading rate. These laboratory events initiate due to a slip rate-dependent process that also determines the size of the stress drop [Dieterich, 1979; Ruina, 1983] and as a consequence, stress drop varies weakly but systematically with loading rate [e.g., Gu and Wong, 1991; Karner and Marone, 2000; McLaskey et al., 2012]. This is especially evident in experiments where the loading rate is changed by orders of magnitude, as is thought to be the loading condition of naturally occurring, small repeating earthquakes driven by afterslip, or low-frequency earthquakes loaded by episodic slip. As follows from the previous studies referred to above, experimentally observed stress drops are well described by a logarithmic dependence on recurrence interval that can be cast as a non-linear slip-predictable model. The fault’s rate dependence of strength is the key physical parameter. Additionally, even at constant loading rate the most reproducible laboratory recurrence is not exactly periodic, unlike existing friction recurrence models. We present example laboratory catalogs that document the variance and show that in large catalogs, even at constant loading rate, stress drop and recurrence co-vary systematically. The origin of this covariance is largely consistent with variability of the dependence of fault strength on slip rate. Laboratory catalogs show aspects of both slip and time predictability and successive stress drops are strongly correlated indicating a ‘memory’ of prior slip history that extends over at least one recurrence cycle.

Journal of Geophysical Research↗

Seismic source dynamics of gas-piston activity at Kı̄lauea Volcano, Hawai‘i

Since 2008, eruptive activity at the summit of Kı̄lauea Volcano, Hawai‘i has been confined to the new Overlook pit crater within the Halema‘uma‘u Crater. Among the broad range of magmatic processes observed in the new pit are recurring episodes of gas pistoning. The gas-piston activity is accompanied by seismic signals that are recorded by a broadband network deployed in the summit caldera. We use raw data recorded with this network to model the source mechanism of representative gas-piston events in a sequence that occurred on 20–25 August 2011 during a gentle inflation of the Kı̄lauea summit. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–10,000 s. Most of the seismic wavefield produced by gas-pistoning originates in a source region ∼1 km below the eastern perimeter of the Halema‘uma‘u pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks featuring an east striking crack (dike) dipping 80°to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each gas-piston event is marked by a similar rapid inflation lasting a few minutes, trailed by a slower deflation ramp extending up to 15 min, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes accompanying the growth and collapse of a layer of foam at the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes : foam thickness (10–50 m), foam cell diameter (0.04–0.10 m), and gas-injection velocity (0.01–0.06 m s −1 ). Gas-piston activity occurs in a narrow pipe with diameter of 6 m connecting the Halema‘uma‘u pit crater to the subjacent dike system. The height of the magma column is estimated at ∼104 m at the start of the sequence based on the period of very long period (VLP) oscillations accompanying the onset of the gas-piston signal. Based on the change in the period of VLP oscillations and tilt evidence, the height of the magma column is inferred to have risen by up to ∼23 m by the end of the 5 day long sequence. A penny-shaped crack model of the dike geometry yields effective diameters of ∼1.2–2.9 km for the east dike and 0.7 km for the north dike. The shallower north dike segment is embedded in a relatively weak medium, compatible with expected mechanical properties in the hydrothermal environment of this dike.

Hawaii↗

Abrupt termination of Marine Isotope Stage 16 (Termination VII) at 631.5 ka in Santa Barbara Basin, California

The Marine Isotope Stage 16–15 boundary (Termination VII) is the first deglacial warming step of the late Quaternary following the mid-Pleistocene transition (MPT), when 41 kyr climatic cycles shifted to strong 100 kyr cycles. The detailed structure of this important climatic event has remained unknown until now. Core MV0508-19JPC from Santa Barbara Basin, California, contains a decadal-scale climatic and geochemical sediment record of 4000 years duration that includes the early part of this deglacial episode. This record reveals that the climatic shift during the early deglacial occurred rapidly (<700 years), in a progression of three abrupt warming steps. The onset of Marine Isotope Stage (MIS) 15 was remarkably abrupt with 4–5°C sea surface warming in ~50 years. The deglacial sequence contains the well-dated Lava Creek tephra (631.3 ± 4 ka) from Yellowstone Caldera used to date the onset of Termination VII at 631.5 ka. The late MIS 16 and early MIS 15 interval exhibits multiple decadal-scale negative excursions in δ 13 C of planktic foraminifera, likely the result of repeated discharges of methane from methane hydrates associated with both ocean warming and low sea level. A warm interstadial that interrupts late MIS 16 is marked by elevated concentrations of redox-sensitive elements indicating sulfidic, oxygen-deficient bottom and pore-waters, and elevated concentrations of total organic carbon and Cd, reflecting increased surface productivity. Unlike younger sediments on the California margin, these indicators of increased productivity and low dissolved oxygen do not consistently correspond with each other or with preserved laminations, possibly reflecting instability of a still evolving ocean-atmosphere system following the MPT.

California↗

Dynamics of the Yellowstone hydrothermal system

The Yellowstone Plateau Volcanic Field is characterized by extensive seismicity, episodes of uplift and subsidence, and a hydrothermal system that comprises more than 10,000 thermal features, including geysers, fumaroles, mud pots, thermal springs, and hydrothermal explosion craters. The diverse chemical and isotopic compositions of waters and gases derive from mantle, crustal, and meteoric sources and extensive water-gas-rock interaction at variable pressures and temperatures. The thermal features are host to all domains of life that utilize diverse inorganic sources of energy for metabolism. The unique and exceptional features of the hydrothermal system have attracted numerous researchers to Yellowstone beginning with the Washburn and Hayden expeditions in the 1870s. Since a seminal review published a quarter of a century ago, research in many fields has greatly advanced our understanding of the many coupled processes operating in and on the hydrothermal system. Specific advances include more refined geophysical images of the magmatic system, better constraints on the time scale of magmatic processes, characterization of fluid sources and water-rock interactions, quantitative estimates of heat and magmatic volatile fluxes, discovering and quantifying the role of thermophile microorganisms in the geochemical cycle, defining the chronology of hydrothermal explosions and their relation to glacial cycles, defining possible links between hydrothermal activity, deformation, and seismicity; quantifying geyser dynamics; and the discovery of extensive hydrothermal activity in Yellowstone Lake. Discussion of these many advances forms the basis of this review.

Idaho, Montana, Wyoming↗

Coupling centennial-scale shoreline change to sea-level rise and coastal morphology in the Gulf of Mexico using a Bayesian network

Predictions of coastal evolution driven by episodic and persistent processes associated with storms and relative sea-level rise (SLR) are required to test our understanding, evaluate our predictive capability, and to provide guidance for coastal management decisions. Previous work demonstrated that the spatial variability of long-term shoreline change can be predicted using observed SLR rates, tide range, wave height, coastal slope, and a characterization of the geomorphic setting. The shoreline is not suf- ficient to indicate which processes are important in causing shoreline change, such as overwash that depends on coastal dune elevations. Predicting dune height is intrinsically important to assess future storm vulnerability. Here, we enhance shoreline-change predictions by including dune height as a vari- able in a statistical modeling approach. Dune height can also be used as an input variable, but it does not improve the shoreline-change prediction skill. Dune-height input does help to reduce prediction uncer- tainty. That is, by including dune height, the prediction is more precise but not more accurate. Comparing hindcast evaluations, better predictive skill was found when predicting dune height (0.8) compared with shoreline change (0.6). The skill depends on the level of detail of the model and we identify an optimized model that has high skill and minimal overfitting. The predictive model can be implemented with a range of forecast scenarios, and we illustrate the impacts of a higher future sea-level. This scenario shows that the shoreline change becomes increasingly erosional and more uncertain. Predicted dune heights are lower and the dune height uncertainty decreases.

Earth's Future↗

Age of the Lava Creek supereruption and magma chamber assembly at Yellowstone based on 40 Ar/ 39 Ar and U-Pb dating of sanidine and zircon crystals

The last supereruption from the Yellowstone Plateau formed Yellowstone caldera and ejected the >1000 km 3 of rhyolite that composes the Lava Creek Tuff. Tephra from the Lava Creek eruption is a key Quaternary chronostratigraphic marker, in particular for dating the deposition of mid Pleistocene glacial and pluvial deposits in western North America. To resolve the timing of eruption and crystallization history for the Lava Creek magma, we performed (1) 40 Ar/ 39 Ar dating of single sanidine crystals to delimit eruption age and (2) ion microprobe U-Pb and trace-element analyses of the crystal faces and interiors of single zircons to date the interval of zircon crystallization and characterize magmatic evolution. Sanidines from the two informal members composing Lava Creek Tuff yield a preferred 40 Ar/ 39 Ar isochron date of 631.3&thinsp;&plusmn;&thinsp;4.3 ka. Crystal faces on zircons from both members yield a weighted mean 206 Pb/ 238 U date of 626.5&thinsp;&plusmn;&thinsp;5.8 ka, and have trace element concentrations that vary with the eruptive stratigraphy. Zircon interiors yield a mean 206 Pb/ 238 U date of 659.8&thinsp;&plusmn;&thinsp;5.5 ka, and reveal reverse and/or oscillatory zoning of trace element concentrations, with many crystals containing high U concentration cores that likely grew from highly evolved melt. The occurrence of distal Lava Creek tephra in stratigraphic sequences marking the Marine Isotope Stage 16&ndash;15 transition supports the apparent eruption age of &sim;631 ka. The combined results reveal that Lava Creek zircons record episodic heating, renewed crystallization, and an overall up-temperature evolution for Yellowstone's subvolcanic reservoir in the 10 3 &minus;10 4 year interval before eruption.

Wyoming↗

Surface melt dominates Alaska glacier mass balance

Mountain glaciers comprise a small and widely distributed fraction of the world's terrestrial ice, yet their rapid losses presently drive a large percentage of the cryosphere's contribution to sea level rise. Regional mass balance assessments are challenging over large glacier populations due to remote and rugged geography, variable response of individual glaciers to climate change, and episodic calving losses from tidewater glaciers. In Alaska, we use airborne altimetry from 116 glaciers to estimate a regional mass balance of &minus;75&thinsp;&plusmn;&thinsp;11&thinsp;Gt&thinsp;yr &minus;1 (1994&ndash;2013). Our glacier sample is spatially well distributed, yet pervasive variability in mass balances obscures geospatial and climatic relationships. However, for the first time, these data allow the partitioning of regional mass balance by glacier type. We find that tidewater glaciers are losing mass at substantially slower rates than other glaciers in Alaska and collectively contribute to only 6% of the regional mass loss.

Alaska, British Columbia, Yukon Territory↗

Hydrothermal response to a volcano-tectonic earthquake swarm, Lassen, California

The increasing capability of seismic, geodetic, and hydrothermal observation networks allows recognition of volcanic unrest that could previously have gone undetected, creating an imperative to diagnose and interpret unrest episodes. A November 2014 earthquake swarm near Lassen Volcanic National Park, California, which included the largest earthquake in the area in more than 60&thinsp;years, was accompanied by a rarely observed outburst of hydrothermal fluids. Although the earthquake swarm likely reflects upward migration of endogenous H 2 O-CO 2 fluids in the source region, there is no evidence that such fluids emerged at the surface. Instead, shaking from the modest sized (moment magnitude 3.85) but proximal earthquake caused near-vent permeability increases that triggered increased outflow of hydrothermal fluids already present and equilibrated in a local hydrothermal aquifer. Long-term, multiparametric monitoring at Lassen and other well-instrumented volcanoes enhances interpretation of unrest and can provide a basis for detailed physical modeling.

California↗

Cascadia subducting plate fluids channelled to fore-arc mantle corner: ETS and silica deposition

In this study we first summarize the constraints that on the Cascadia subduction thrust, there is a 70 km gap downdip between the megathrust seismogenic zone and the Episodic Tremor and Slip (ETS) that lies further landward; there is not a continuous transition from unstable to conditionally stable sliding. Seismic rupture occurs mainly offshore for this hot subduction zone. ETS lies onshore. We then suggest what does control the downdip position of ETS. We conclude that fluids from dehydration of the downgoing plate, focused to rise above the fore-arc mantle corner, are responsible for ETS. There is a remarkable correspondence between the position of ETS and this corner along the whole margin. Hydrated mineral assemblages in the subducting oceanic crust and uppermost mantle are dehydrated with downdip increasing temperature, and seismic tomography data indicate that these fluids have strongly serpentinized the overlying fore-arc mantle. Laboratory data show that such fore-arc mantle serpentinite has low permeability and likely blocks vertical expulsion and restricts flow updip within the underlying permeable oceanic crust and subduction shear zone. At the fore-arc mantle corner these fluids are released upward into the more permeable overlying fore-arc crust. An indication of this fluid flux comes from low Poisson's Ratios (and V p / V s ) found above the corner that may be explained by a concentration of quartz which has exceptionally low Poisson's Ratio. The rising fluids should be silica saturated and precipitate quartz with decreasing temperature and pressure as they rise above the corner.

Journal of Geophysical Research↗

Differential uplift and incision of the Yakima River terraces, central Washington State

The fault-related Yakima folds deform Miocene basalts and younger deposits of the Columbia Plateau in central Washington State. Geodesy implies ~2&thinsp;mm/yr of NNE directed shortening across the folds, but until now the distribution and rates of Quaternary deformation among individual structures has been unclear. South of Ellensburg, Washington, the Yakima River cuts a ~600&thinsp;m deep canyon across several Yakima folds, preserving gravel-mantled strath terraces that record progressive bedrock incision and related rock uplift. Here we integrate cosmogenic isochron burial dating of the strath terrace gravels with lidar analysis and field mapping to quantify rates of Quaternary differential incision and rock uplift across two folds transected by the Yakima River: Manastash and Umtanum Ridge. Isochron burial ages from in situ produced 26 Al and 10 Be at seven sites across the folds date episodes of strath terrace formation over the past ~2.9&thinsp;Ma. Average bedrock incision rates across the Manastash (~88&thinsp;m/Myr) and Umtanum Ridge (~46&thinsp;m/Myr) anticlines are roughly 4 to 8 times higher than rates in the intervening syncline (~14&thinsp;m/Myr) and outside the canyon (~10&thinsp;m/Myr). These contrasting rates demonstrate differential bedrock incision driven by ongoing Quaternary rock uplift across the folds at rates corresponding to ~0.13 and ~0.06&thinsp;mm/yr shortening across postulated master faults dipping 30&thinsp;&plusmn;&thinsp;10&deg;S beneath the Manastash and Umtanum Ridge anticlines, respectively. The reported Quaternary shortening across the anticlines accounts for ~10% of the ~2&thinsp;mm/yr geodetic budget, suggesting that other Yakima structures actively accommodate the remaining contemporary deformation.

Washington↗

The geologic history of Margaritifer basin, Mars

In this study, we investigate the fluvial, sedimentary, and volcanic history of Margaritifer basin and the Uzboi-Ladon-Morava (ULM) outflow channel system. This network of valleys and basins spans more than 8000 km in length, linking the fluvially dissected southern highlands and Argyre Basin with the northern lowlands via Ares Vallis. Compositionally, thermophysically, and morphologically distinct geologic units are identified and are used to place critical relative stratigraphic constraints on the timing of geologic processes in Margaritifer basin. Our analyses show that fluvial activity was separated in time by significant episodes of geologic activity, including the widespread volcanic resurfacing of Margaritifer basin and the formation of chaos terrain. The most recent fluvial activity within Margaritifer basin appears to terminate at a region of chaos terrain, suggesting possible communication between surface and subsurface water reservoirs. We conclude with a discussion of the implications of these observations on our current knowledge of Martian hydrologic evolution in this important region.

Journal of Geophysical Research E: Planets↗

Fitful and protracted magma assembly leading to a giant eruption, Youngest Toba Tuff, Indonesia

The paroxysmal eruption of the 74 ka Youngest Toba Tuff (YTT) of northern Sumatra produced an extraordinary 2800 km 3 of non-welded to densely welded ignimbrite and co-ignimbrite ash-fall. We report insights into the duration of YTT magma assembly obtained from ion microprobe U-Th and U-Pb dates, including continuous age spectra over >50% of final zircon growth, for pumices and a welded tuff spanning the compositional range of the YTT. A relatively large subpopulation of zircon crystals nucleated before the penultimate caldera-related eruption at 501 ka, but most zircons yielded interior dates 100-300 ka thereafter. Zircon nucleation and growth was likely episodic and from diverse conditions over protracted time intervals of >100 to >500 ka. Final zircon growth is evident as thin rim plateaus that are in Th/U chemical equilibrium with hosts, and that give crystallization ages within tens of ka of eruption. The longevity and chemical characteristics of the YTT zircons, as well as evidence for intermittent zircon isolation and remobilization associated with magma recharge, is especially favored at the cool and wet eutectoid conditions that characterize at least half of the YTT, wherein heat fluxes could dissolve major phases but have only a minor effect on larger zircon crystals. Repeated magma recharge may have contributed to the development of compositional zoning in the YTT but, considered together with limited allanite, quartz, and other mineral dating and geospeedometry, regular perturbations to the magma reservoir over >400 ka did not lead to eruption until 74 ka ago.

Lake Toba↗

Ephemerality of discrete methane vents in lake sediments

Methane is a potent greenhouse gas whose emission from sediments in inland waters and shallow oceans may both contribute to global warming and be exacerbated by it. The fraction of methane emitted by sediments that bypasses dissolution in the water column and reaches the atmosphere as bubbles depends on the mode and spatiotemporal characteristics of venting from the sediments. Earlier studies have concluded that hot spots&mdash;persistent, high-flux vents&mdash;dominate the regional ebullitive flux from submerged sediments. Here the spatial structure, persistence, and variability in the intensity of methane venting are analyzed using a high-resolution multibeam sonar record acquired at the bottom of a lake during multiple deployments over a 9 month period. We confirm that ebullition is strongly episodic, with distinct regimes of high flux and low flux largely controlled by changes in hydrostatic pressure. Our analysis shows that the spatial pattern of ebullition becomes homogeneous at the sonar's resolution over time scales of hours (for high-flux periods) or days (for low-flux periods), demonstrating that vents are ephemeral rather than persistent, and suggesting that long-term, lake-wide ebullition dynamics may be modeled without resolving the fine-scale spatial structure of venting.

Geophysical Research Letters↗