Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “One Earth”

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 1,315 records · Page 73Linked to original sources
USGS2012

Seismic constraints on the nature of lower crustal reflectors beneath the extending Southern Transition Zone of the Colorado Plateau, Arizona

We determine the reflection polarity and exploit variations in P and S wave reflectivity and P wave amplitude versus offset (AVO) to constrain the origin of lower crustal reflectivity observed on new three-component seismic data recorded across the structural transition of the Colorado Plateau. The near vertical incidence reflection data were collected by Stanford University in 1989 as part of the U.S. Geological Survey Pacific to Arizona Crustal Experiment that traversed the Arizona Transition Zone of the Colorado Plateau. The results of independent waveform modeling methods are consistent with much of the lower crustal reflectivity resulting from thin, high-impedance layers. The reflection polarity of the cleanest lower crustal events is positive, which implies that these reflections result from high-velocity contrasts, and the waveform character indicates that the reflectors are probably layers less than or approximately equal to 200 m thick. The lower crustal events are generally less reflective to incident S waves than to P waves, which agrees with the predicted behavior of high-velocity mafic layering. Analysis of the P wave AVO character of lower crustal reflections demonstrates that the events maintain a constant amplitude with offset, which is most consistent with a mafic-layering model. One exception is a high-amplitude (10 dB above background) event near the base of lower crustal reflectivity which abruptly decreases in amplitude at increasing offsets. The event has a pronounced S wave response, which along with its negative AVO trend is a possible indication of the presence of fluids in the lower crust. The Arizona Transition Zone is an active but weakly extended province, which causes us to discard models of lower crustal layering resulting from shearing because of the high degree of strain required to create such layers. Instead, we favor horizontal basaltic intrusions as the primary origin of high-impedance reflectors based on (1) The fact that most xenoliths in eruptive basalts of the Transition Zone are of mafic igneous composition, (2) indications that a pulse of magmatic activity crossed the Transition Zone in the late Tertiary period, and (3) the high regional heat flow observed in the Transition Zone. The apparent presence of fluids near the base of the reflective zone may indicate a partially molten intrusion. We present a mechanism by which magma can be trapped and be induced to intrude horizontally at rheologic contrasts in extending crust.

Arizona↗
USGS1986

Radarclinometry

A mathematical theory and a corresponding algorithm have been developed to derive topographic maps from radar images as photometric arrays. Thus, as radargrammetry is to photogrammetry, so radarclinometry is to photoclinometry. Photoclinometry is endowed with a fundamental indeterminacy principle even for terrain homogeneous in normal albedo. This arises from the fact that the geometric locus of orientations of the local surface normal that is consistent with a given reflected specific-intensity of radiation is more complicated than a fixed line in space. For a radar image, the locus is a cone whose half-angle is the incidence angle and whose axis contains the radar. The indeterminacy is removed throughout a region if one possesses a control profile as a boundary-condition. In the absence of such ground-truth, a point-boundary-condition will suffice only in conjunction with a heuristic assumption, such as that the strike-line runs perpendicularly to the line-of-sight. In the present study I have implemented a more reasonable assumption which I call 'the hypothesis of local cylindricity'. Firstly, a general theory is derived, based solely on the implicit mathematical determinacy. This theory would be directly indicative of procedure if images were completely devoid of systematic error and noise. The theory produces topography by an area integration of radar brightness, starting from a control profile, without need of additional idealistic assumptions. But we have also theorized separately a method of forming this control profile, which method does require an additional assumption about the terrain. That assumption is that the curvature properties of the terrain are locally those of a cylinder of inferable orientation, within a second-order mathematical neighborhood of every point of the terrain. While local strike-and-dip completely determine the radar brightness itself, the terrain curvature determines the brightness-gradient in the radar image. Therefore, the control profile is formed as a line integration of brightness and its local gradient starting from a single point of the terrain where the local orientation of the strike-line is estimated by eye. Secondly, and independently, the calibration curve for pixel brightness versus incidence-angle is produced. I assume that an applicable curve can be found from the literature or elsewhere so that our problem is condensed to that of properly scaling the brightness-axis of the calibration curve. A first estimate is found by equating the average image brightness to the point on the brightness axis corresponding to the complement of the effective radar depression-angle, an angle assumed given. A statistical analysis is then used to correct, on the one hand, for the fact that the average brightness is not the brightness that corresponds to the average incidence angle, as a result of the non-linearity of the calibration curve; and on the other hand, we correct for the fact that the average incidence angle is not the same for a rough surface as it is for a flat surface (and therefore not the complement of the depression angle). Lastly, the practical modifications that were interactively evolved to produce an operational algorithm for treating real data are developed. They are by no means considered optimized at present. Such a possibility is thus far precluded by excessive computer-time. Most noteworthy in this respect is the abandonment of area integration away from a control profile. Instead, the topography is produced as a set of independent line integrations down each of the parallel range lines of the image, using the theory for control-profile formation. An adaptive technique, which now appears excessive, was also employed so that SEASAT images of sand dunes could be processed. In this, the radiometric calibration was iterated to force the endpoints of each profile to zero elevation. A secondary algorithm then employed line-averages of appropriate quantities to adjust the mean t

Earth, Moon and Planets↗
USGS1998

Suppression of large earthquakes by stress shadows: A comparison of Coulomb and rate-and-state failure

Stress shadows generated by California's two most recent great earthquakes (1857 Fort Tejon and 1906 San Francisco) substantially modified 19th and 20th century earthquake history in the Los Angeles basin and in the San Francisco Bay area. Simple Coulomb failure calculations, which assume that earthquakes can be modeled as static dislocations in an elastic half-space, have done quite well at approximating how long the stress shadows, or relaxing effects, should last and at predicting where subsequent large earthquakes will not occur. There has, however, been at least one apparent exception to the predictions of such simple models. The 1911 M >6.0 earthquake near Morgan Hill, California, occurred at a relaxed site on the Calaveras fault. We examine how the more complex rate-and-state friction formalism based on laboratory experiments might have allowed the 1911 earthquake. Rate-and-state time-to-failure calculations are consistent with the occurrence of the 1911 event just 5 years after 1906 if the Calaveras fault was already close to failure before the effects of 1906. We also examine the likelihood that the entire 78 years of relative quiet (only four M ≥6 earthquakes) in the bay area after 1906 is consistent with rate-and-state assumptions, given that the previous 7 decades produced 18 M ≥6 earthquakes. Combinations of rate-and-state variables can be found that are consistent with this pattern of large bay area earthquakes, assuming that the rate of earthquakes in the 7 decades before 1906 would have continued had 1906 not occurred. These results demonstrate that rate-and-state offers a consistent explanation for the 78-year quiescence and the 1911 anomaly, although they do not rule out several alternate explanations.

California↗
USGS2007

Paleomagnetic and mineral magnetic constraints on Zn-Pb ore genesis in the Pend Oreille Mine, Metaline district, Washington, USA

Zinc-lead mineralization in the Metaline mining district of northeastern Washington, USA, is hosted by the Cambrian Metaline Formation and is classified into Yellowhead-type (YO) and Josephine-type (JO) ore based on texture and mineralogy. Paleomagnetic results are reported for four Cambrian Metaline Formation sites, one Ordovician Ledbetter slate site, 12 YO and 13 JO (including two breccia sites) mineralization sites in the Pend Oreille Mine, and eight sites from the nearby Cretaceous Kaniksu granite batholith. Thermal and alternating field step demagnetization, saturation isothermal remanence analysis, and synthetic specimen tests show that the remanence in the host carbonates and Zn-Pb mineralization is carried mostly by pseudosingle (PSD) to single domain (SD) pyrrhotite and mostly by PSD to SD magnetite in the Kaniksu granite. Based on thermomagnetic measurements, sphalerite and galena concentrates and tailings from the mine's mill contain hexagonal and monoclimc pyrrhotite. The postfolding characteristic remanent magnetization (ChRM), known thermal data, and paleoarc method of dating suggest that the Zn-Pb mineralization carries a primary chemical remanent magnetization (CRM), and Metaline Formation carbonates a secondary CRM that were acquired during the Middle Jurassic (166 ??6 Ma) during the waning stages of the Nevadan orogeny. A paleomagnetic breccia test favours a solution-collapse origin for the Josephine breccia. Finally, the Kaniksu paleopole is concordant with the North American Cretaceous reference paleopole, suggesting the Kootenay terrane has not been rotated since emplacement of the batholith at ???94 Ma. ?? 2007 NRC Canada.

Canadian Journal of Earth Sciences↗
USGS2025

An automated compositing method for producing annual clear images from Landsat Collection 2 for annual NLCD production

Quality image input is fundamental to the quality of derived land cover products. Substantial time and effort are usually required to prepare images. Here, we present a novel and streamlined compositing algorithm that ingests Landsat Collection 2 Analysis Ready Data (ARD) and outputs cloud-free and gap-free composite imagery, which can be directly used for classification. This method leverages and improves the previous National Land Cover Database (NLCD) Virtual Median Value Point (VMVP) compositing method, the first part of the image preparation for NLCD 2019 operational production. The NLCD 2019 image preparation approach includes a second part, a residual cloud and cloud shadow detection and gap-filling method, to produce final cloud-free and gap-free composite imagery. The second part requires one clear reference image for each target year. Additional reference images are needed for producing reasonable observations for perennial ice/snow areas because Pixel QA (Quality Assessment) from ARD has difficulties differentiating ice/snow areas from clouds. Unlike the NLCD 2019 image preparation approach, our new compositing method, which is referred to as Automated VMVP (AVMVP), uses Landsat ARD as the only input and does not require reference images and extra steps. In this method, we developed new spectral filter criteria coupled with counts of clear observations using Pixel QA to identify potential cloud and cloud shadow observations on initially selected observations from the NLCD VMVP compositing algorithm. We also automate “gap-filling” using clear observations retrieved from a maximum of ±2 years around the target year when needed. Finally, a percentile-filtered compositing method was developed for the perennial ice/snow areas. All these steps are streamlined, pixel-based, and directly run on Landsat Collection 2 ARD. We have run successful tests on the conterminous United States (CONUS). Composite images derived from our innovative method were used to produce the CONUS Annual NLCD Collection 1 product suite that covers the period from 1985 to 2023.

conterminous United States↗
USGS2018

An updated method for estimating landslide‐event magnitude

Summary statistics derived from the frequency–area distribution (FAD) of inventories of triggered landslides allows for direct comparison of landslides triggered by one event (e.g. earthquake, rainstorm) with another. Such comparisons are vital to understand links between the landslide‐event and the environmental characteristics of the area affected. This could lead to methods for rapid estimation of landslide‐event magnitude, which in turn could lead to estimates of the total triggered landslide area. Previous studies proposed that the FAD of landslides follows an inverse power‐law, which provides the basis to model the size distribution of landslides and to estimate landslide‐event magnitude ( mLS ), which quantifies the severity of the event. In this study, we use a much larger collection of earthquake‐induced landslide (EQIL) inventories ( n =45) than previous studies to show that size distributions are much more variable than previously assumed. We present an updated model and propose a method for estimating mLS and its uncertainty that better fits the observations and is more reproducible, robust, and consistent than existing methods. We validate our model by computing mLS for all of the inventories in our dataset and comparing that with the total landslide areas of the inventories. We show that our method is able to estimate the total landslide area of the events in this larger inventory dataset more successfully than the existing methods.

Earth Surface Processes and Landforms↗
USGS2018

The influence of neap-spring tidal variation and wave energy on sediment flux in salt marsh tidal creeks

Sediment flux in marsh tidal creeks is commonly used to gage sediment supply to marshes. We conducted a field investigation of temporal variability in sediment flux in tidal creeks in the accreting tidal marsh at China Camp State Park adjacent to northern San Francisco Bay. Suspended-sediment concentration (SSC), velocity, and depth were measured near the mouths of two tidal creeks during three six-to-ten-week deployments: two in winter and one in summer. Currents, wave properties and SSC were measured in the adjacent shallows. All deployments spanned the largest spring tides of the season. Results show that tidally-averaged suspended-sediment flux (SSF) in the tidal creeks decreased with increasing tidal energy, and SSF was negative (bayward) for tidal cycles with maximum water surface elevation above the marsh plain. Export during the largest spring tides dominated the cumulative SSF measured during the deployments. During ebb tides following the highest tides, velocities exceeded 1 m/s in the narrow tidal creeks, resulting in negative tidally-averaged water flux, and mobilizing sediment from the creek banks or bed. Storm surge also produced negative SSF. Tidally-averaged SSF was positive in wavey conditions with moderate tides. Spring-tide sediment export was about 50% less at a station 130 m further up the tidal creek than at the creek mouth. The negative tidally-averaged water flux near the creek mouth during spring tides indicates that in the lower marsh, some of the water flooding directly across the bay--marsh interface drains through the tidal creeks, and suggests that this interface may be a pathway for sediment supply to the lower marsh as well.

California↗
USGS2012

Re–Os geochronology of the lacustrine Green River Formation: Insights into direct depositional dating of lacustrine successions, Re–Os systematics and paleocontinental weathering

Lacustrine sedimentary successions provide exceptionally high-resolution records of continental geological processes, responding to tectonic, climatic and magmatic influences. These successions are therefore essential for correlating geological and climatic phenomena across continents and furthermore the globe. Producing accurate geochronological frameworks within lacustrine strata is challenging because the stratigraphy is often bereft of biostratigraphy and directly dateable tuff horizons. The rhenium–osmium (Re–Os) geochronometer is a well-established tool for determining precise and accurate depositional ages of marine organic-rich rocks. Lake systems with stratified water columns are predisposed to the preservation of organic-rich rocks and thus should permit direct Re–Os geochronology of lacustrine strata. We present Re–Os systematics from one of the world's best documented lacustrine systems, the Eocene Green River Formation, providing accurate Re–Os depositional dates that are supported by Ar–Ar and U–Pb ages of intercalated tuff horizons. Precision of the Green River Formation Re–Os dates is controlled by the variation in initial 187 Os/ 188 Os and the range of 187 Re/ 188 Os ratios, as also documented in marine systems. Controls on uptake and fractionation of Re and Os are considered to relate mainly to depositional setting and the type of organic matter deposited, with the need to further understand the chelating precursors of Re and Os in organic matter highlighted. In addition to geochronology, the Re–Os data records the 187Os/188Os composition of lake water (1.41–1.54) at the time of deposition, giving an insight into continental runoff derived from weathering of the geological hinterland of the Green River Formation. Such insights enable us to evaluate fluctuations in continental climatic, tectonic and magmatic processes and provide the ability for chemostratigraphic correlation combined with direct depositional dates. Furthermore, initial 187Os/188Os values can be used as a diagnostic tool to distinguish between lacustrine and marine depositional settings when compared to known oceanic 187 Os/ 188 Os values.

Colorado;Utah;Wyoming↗
USGS1989

Style of extensional tectonism during rifting, Red Sea and Gulf of Aden

Models describing the development of the Red Sea and the Gulf of Aden, prior to the present periods of sea-floor spreading, include those that use block faulting on steep normal faults, uniform diffuse shear in continental crust, simple shear on large detachment faults that cut the entire lithosphere, combinations involving detachment faults/ductile deformation/plutonic inflation, and ones that minimize the role of mechanical extension in favor of an earlier stage of sea-floor spreading. Geologic and geophysical studies from the Arabian continental margin in the southern Red Sea and LANDSAT analysis of the northern Somalia margin in the Gulf of Aden suggest that the early continental rifts were long narrow features that formed by extension on closely spaced normal faults above moderate- to shallow-dipping detachments with break-away zones defining one rift flank and root zones under the opposing rift flank. The rift flanks presently form the opposing continental margins across each ocean basin. The detachment on the Arabian margin dips gently to the west, with a breakaway zone now eroded above the deeply dissected terrain of the Arabian escarpment. The Arabian detachment projects westward to middle crustal levels beneath the sediment of the southern Red Sea coastal plain. Strata in the upper plate dip as steeply as 60° to the west, and the beds are repeated by numerous planar and listric normal faults that dip to the east. Most of the faults truncate downward at the detachment. Thus, the upper plate is highly extended and the rocks in its eastern part have been translated about 20 km westward and 21/2- to 5-km downward relative to the rest of Arabia. A prominent detachment surface, with a north dip, is evident in northernmost Somalia where it breaks away north of the Somalian escarpment in an otherwise undeformed section of cratonic strata of Jurassic to Eocene age. The upper plate of the Somalian detachment consists of a highly faulted collage of the cratonic strata. This fault projects to middle crustal levels in the opposing Arabian margin to the northeast.

Journal of African Earth Sciences↗
USGS2009

Diverse lavas from closely spaced volcanoes drawing from a common parent: Emmons Lake Volcanic Center, Eastern Aleutian Arc

Emmons Lake Volcanic Center (ELVC) on the lower Alaskan Peninsula is one of the largest and most diverse volcanic centers in the Aleutian Arc. Since the Middle Pleistocene, eruption of ~ 350 km 3 of basalt through rhyolite has produced a 30 km, arc front chain of nested calderas and overlapping stratovolcanoes. ELVC has experienced as many as five major caldera-forming eruptions, the most recent, at ~ 27 ka, produced ~ 50 km 3 of rhyolitic ignimbrite and ash fall. These violent silicic events were interspersed with less energetic, but prodigious, outpourings of basalt through dacite. Holocene eruptions are mostly basaltic andesite to andesite and historically recorded activity includes over 40 eruptions within the last 200 yr, all from Pavlof volcano, the most active site in the Aleutian Arc. Geochemical and geophysical observations suggest that although all ELVC eruptions derive from a common clinopyroxene + spinel + plagioclase fractionating high-aluminum basalt parent in the lower crust, magma follows one of two closely spaced, but distinct paths to the surface. Under the eastern end of the chain, magma moves rapidly and cleanly through a relatively young (~ 28 ka), hydraulically connected dike plexus. Steady supply, short magma residence times, and limited interaction with crustal rocks preserve the geochemistry of deep crustal processes. Below the western part of the chain, magma moves haltingly through a long-lived (~ 500 ka) and complex intrusive column in which many generations of basaltic to andesitic melts have mingled and fractionated. Buoyant, silicic melts periodically separate from the lower parts of the column to feed voluminous eruptions of dacite and rhyolite. Mafic lavas record a complicated passage through cumulate zones and hydrous silicic residues as manifested by disequilibrium phenocryst textures, incompatible element enrichments, and decoupling of REEs and HFSEs ratios. Such features are absent in mafic lavas from the younger part of the chain, highlighting the importance of plumbing architecture and longevity in creating petrologic diversity. Supplemental Data include 156 major element (XRF) and 128 trace element (ICP-MS) whole-rock analyses, 23 new 40 Ar/ 39 Ar ages, a generalized geologic map with associated unit descriptions and field photographs, and photomicrographs of key petrographic features.

Alaska↗
USGS2003

Geochronology and eruptive history of the Katmai volcanic cluster, Alaska Peninsula

In the Katmai district of the Alaska Peninsula, K–Ar and 40 Ar/ 39 Ar ages have been determined for a dozen andesite–dacite stratocones on the arc front and for 11 rear-arc volcanoes, 10 of which are monogenetic. Tied to mapping and stratigraphic studies, our dating emphasized proximal basal lavas that rest on basement rocks, in order to estimate ages of inception of each polygenetic cone. Oldest among arc-front cones is Alagogshak Volcano (690–43 ka), succeeded in the Holocene by the active Mount Martin cone. Mount Mageik consists of four overlapping subedifices, basal lavas of which give ages of 93, 71, and 59 ka, and Holocene. The three small prehistoric cones of Trident Volcano yield ages of 143, 101–58, and 44 ka. Falling Mountain and Mount Cerberus, dacite domes near the 1912 Novarupta vent, are related compositionally to the Trident group and give ages of 70 ka and 114 ka. Mount Katmai, which underwent caldera collapse in 1912, consists of two subedifices that overlapped in space and time, and is the only arc-front center here to include basalt and rhyolite; one cone began by 90 ka, the other by 47 ka. Snowy Mountain also consists of two contiguous cones, which started around 200 and 171 ka, respectively, the younger remaining active into the Holocene. Devils Desk, the only mafic cone on the arc front, was short-lived at about 245 ka. In the rear-arc, (1) Mount Griggs produced mafic-to-silicic andesite in several episodes between 292 ka and the Holocene; (2) the Savonoski River cluster includes a Pliocene dacite dome and five small mafic cones (390–88 ka); (3) Gertrude Creek cone (49.8% SiO 2 ) yields an age of 500 ka; and (4) the Saddlehorn Creek cluster includes five Pliocene basalt-to-andesite remnants. Eruptive volumes were reconstructed, permitting estimates of average eruption rates for edifice lifetimes. Since the mid Pleistocene, total volume erupted along the arc front here is 210±47 km 3 and in the rear-arc 39±6 km 3 , of which Mount Griggs alone accounts for 35±5 km 3 . Most productive has been Mount Katmai at 70±18 km 3 , yielding a rate of ∼1 km 3 /kyr, followed by Mount Mageik (0.33 km 3 /kyr) and Mount Griggs (0.3 km 3 /kyr since 50 ka).

Alaska↗
USGS1988

The surface integral approach to Radarclinometry

Because radarclinometry is fundamentally describable in terms of a nonlinear, first-order, partial differential equation, one expects that it can, in principle, be carried out by direct deterministic integration beginning at a given threshold profile along the azimuthal coordinate. Such a boundary condition could be provided by the altimetry profile obtained on a preceding or succeeding orbital revolution of the radar-bearing spacecraft. Notwithstanding the mismatched resolutions of the radar altimeter and the radar imaging system as planned for the Megallan mission to Venus, there are fundamental considerations, not involving system noise, that influence the possibility of success of this approach. From the topographic map of the Lake Champlain West quadrangle in the Adirondack Mountains of the U.S., a radar image is synthesized. Radarclinometry, in surface integral form, recaptures the topographic map when the applicable radar reflectance function is weakly variable over the range of application, but it diverges beyond a certain point for nominally variable reflectance functions. The effect can be understood by using results from the "shape-from-shading" literature. (This literature is produced by a group within the artificial intelligence community who have been independently attacking, for all practical purposes, photoclinometry, except that they have not given primacy to images of terrain.) The ubiquity of the instability suggests that the value of the surface integral approach is much in doubt. ?? 1988 Kluwer Academic Publishers.

Earth, Moon and Planets↗
USGS2014

Source and progression of a submarine landslide and tsunami: The 1964 Great Alaska earthquake at Valdez

Like many subduction zone earthquakes, the deadliest aspects of the 1964 M = 9.2 Alaska earthquake were the tsunamis it caused. The worst of these were generated by local submarine landslides induced by the earthquake. These caused high runups, engulfing several coastal towns in Prince William Sound. In this paper, we study one of these cases in detail, the Port Valdez submarine landslide and tsunami. We combine eyewitness reports, preserved film, and careful posttsunami surveys with new geophysical data to inform numerical models for landslide tsunami generation. We review the series of events as recorded at Valdez old town and then determine the corresponding subsurface events that led to the tsunami. We build digital elevation models of part of the pretsunami and posttsunami fjord‐head delta. Comparing them reveals a ~1500 m long region that receded 150 m to the east, which we interpret as the primary delta landslide source. Multibeam imagery and high‐resolution seismic reflection data identify a ~400 m wide chute with hummocky deposits at its terminus, which may define the primary slide path. Using these elements we run hydrodynamic models of the landslide‐driven tsunamis that match observations of current direction, maximum inundation, and wave height at Valdez old town. We speculate that failure conditions at the delta front may have been influenced by manmade changes in drainage patterns as well as the fast retreat of Valdez and other glaciers during the past century.

Alaska↗
USGS2020

Seismic and geodetic progression of the 2018 summit caldera collapse of Kīlauea Volcano

The 2018 eruption of Kīlauea volcano, Hawaiʻi, resulted in a major collapse of the summit caldera along with an effusive eruption in the lower East Rift Zone. The caldera collapse comprised 62 highly similar collapse cycles of strong ground deformation and earthquake swarms that ended with a magnitude 5 collapse event and one partial cycle that did not end with a collapse event. We analyzed geodetic and seismic data to better understand how the caldera collapse progressed over 3 months of activity, focusing on the cyclical activity. We identified 3 main phases of collapse: initial ring-fault activation and small explosions (Phase 1), an eastward shift in activity and freeing of the central piston (Phase 2), and a recoupling of the piston to the reservoir followed by relatively steady behavior until the eruption’s end (Phase 3). Additionally, we observed geodetic evidence of tangential motion from the localization of the main ring fault (Phase 2) and the formation of a major peripheral ring fault on the eastern side of the collapse caldera during Phase 3. Both geodetic and seismic parameters suggest that the collapse may have had an eastward-component of motion after the ring fault system had formed. The cyclical seismic and geodetic parameters show no obvious signs that the collapse was coming to an end, with the only notable change being a significant increase in the ratio of cyclical displacement to co-collapse displacement observed during the last complete cycle on GNSS stations outside the caldera region.

Hawaii↗
USGS2001

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Journal of Geophysical Research B: Solid Earth↗
USGS2015

Validation of eDNA markers for New Zealand mudsnail surveillance and initial eDNA monitoring at Mississippi River Basin sites

The performance of newly developed New Zealand mudsnail ( Potamopyrgus antipodarum ; NZMS) genetic markers for environmental (eDNA) analysis of water were compared across two laboratories. The genetic markers were tested in four quantitative polymerase chain reaction assays targeting two regions of the NZMS mitochondrial genome, specifically the cytochrome c oxidase subunit 1 (coi) and cytochrome b (cytb) genes. In a blind study, analysts tested each sample eight times with each assay. There were 10 expected-negative samples from the Black River in La Crosse, Wisconsin, 10 expected-positive samples from the Black Earth Creek in Black Earth, Wisconsin, and 10 known-positive samples from the Black River spiked with NZMS DNA. Previously extracted samples, kept at the Upper Midwest Environmental Sciences Center, were pooled by sample location and then equal quantities were distributed between the Upper Midwest Environmental Sciences Center and the Molecular Conservation Genetics Laboratory at the University of Wisconsin-Stevens Point for analysis. The assays tested were (1) the assay targeting cytb with a minor groove binder probe described by Goldberg and others (2013), (2) the cytb assay with a modified double-quenched probe, (3) an assay targeting coi with a double-quenched probe, and (4) a duplex reaction combining the modified cytb assay and the coi assay. Samples were considered positive for the presence of NZMS DNA when quantitative polymerase chain reaction amplification and probe signal was higher than the normalized threshold value above baseline fluorescence. For the duplex assay, samples were considered positive only when both probe signals were higher than the normalized threshold value above baseline fluorescence. Positive results were then confirmed by sequencing the products. All four assays detected the DNA of NZMS in all expected-positive and known-positive samples in both labs. The modified cytb assay, the coi assay, and the duplex assay all failed to detect the DNA of NZMS in all expected-negative samples in both labs. The cytb assay, as described by Goldberg and others (2013), failed to detect the DNA of NZMS in all expected-negative samples for the Molecular Conservation Genetics Laboratory, but some reactions resulted in positive detection in late cycles for 9 of the 10 expected-negative samples at the Upper Midwest Environmental Sciences Center. Amplicons for expected-negative samples with positive reactions were sent for sequencing, and none were confirmed as NZMS. Six amplicons failed to give readable sequences, and three gave sequences without similarity to any known sequence in GenBank. Amplicons from each assay for one representative positive sample were sequenced and identified as NZMS with greater than 99 percent identity. The duplex assay was chosen as the most efficient assay and was used at the Upper Midwest Environmental Sciences Center to analyze triplicate samples from 29 streams in Wisconsin, 8 streams in Illinois, and 8 streams in Iowa. In order to verify results, additional triplicate samples were collected from two of the streams in Iowa and two of the streams in Wisconsin for analysis at the Molecular Conservation Genetics Laboratory. All samples at all sites were negative for NZMS DNA.

Illinois, Iowa, Wisconsin↗
USGS2003

Crystallisation ages in coeval silicic magma bodies: 238U-230Th disequilibrium evidence from the Rotoiti and earthquake flat eruption deposits, Taupo volcanic zone, New Zealand

The timescales over which moderate to large bodies of silicic magma are generated and stored are addressed here by studies of two geographically adjacent, successive eruption deposits in the Taupo Volcanic Zone, New Zealand. The earlier, caldera-forming Rotoiti eruption (>100 km3 magma) at Okataina volcano was followed, within months at most, by the Earthquake Flat eruption (??? 10 km3 magma) from nearby Kapenga volcano; both generated nonwelded ignimbrite and coeval widespread fall deposits. The Rotoiti and Earthquake Flat deposits are both crystal-rich high-silica rhyolites, with sparse glass-bearing granitoid fragments also occurring in Rotoiti lag breccias generated during caldera collapse. Here we report 238U-230Th disequilibrium data on whole rocks and mineral separates from representative Rotoiti and Earthquake Flat pumices and the co-eruptive Rotoiti granitoid fragments using TIMS and in situ zircon analyses by SIMS. Multiple-grain zircon-controlled crystallisation ages measured by TIMS from the Rotoiti pumice range from 69??3 ka (<63 ?? size fraction) to 76??6 (125-250 ?? fraction), with a weighted mean of all size fractions of 71??2 ka. SIMS model ages from single zircons in pumice range from 50??24 ka to > 350 ka, with a pronounced peak at 70-90 ka. The weighted mean of isochrons is 83??14 ka, in accord with the TIMS data. One glass-bearing Rotoiti granitoid clast yielded an age of 57??8 ka by TIMS (controlled by Th-rich phases that, however, are not apparently present in the juvenile pumices). Another glass-bearing Rotoiti granitoid yielded SIMS zircon model ages peaking at 60-90 ka, having a similar age distribution to the pumice. Age data from pumices are consistent with a published 64??4 ka eruptive age (now modified to 62??2 ka), but chemical and/or mineralogical data imply that the granitoid lithics are not largely crystalline Rotoiti rhyolite, but instead represent contemporaneous partly molten intrusions reflecting different sources in their chemistries and mineralogies. Similarly, although the Earthquake Flat eruption immediately followed (and probably was triggered by) the Rotoiti event, age data from juvenile material are significantly different. A multiple-grain zircon-controlled crystallisation age measured by TIMS from a representative pumice is 173??5 ka, while SIMS model ages range from 70-26+34 ka to >350 ka, with a peak at 105 ka. These age data coupled with previously published geochemical and isotopic data show that the Rotoiti and Earthquake Flat deposits were erupted from independent, unconnected magma bodies. Published by Elsevier Science B.V.

Earth and Planetary Science Letters↗
USGS2019

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗