SEARCH · Search USGS

Results for “Journal of Time Series Analysis”

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

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8Linked to original sources

USGS2007

Swarms of repeating long-period earthquakes at Shishaldin Volcano, Alaska, 2001-2004

During 2001–2004, a series of four periods of elevated long-period seismic activity, each lasting about 1–2 months, occurred at Shishaldin Volcano, Aleutian Islands, Alaska. The time periods are termed swarms of repeating events , reflecting an abundance of earthquakes with highly similar waveforms that indicate stable, non-destructive sources. These swarms are characterized by increased earthquake amplitudes, although the seismicity rate of one event every 0.5–5 min has remained more or less constant since Shishaldin last erupted in 1999. A method based on waveform cross-correlation is used to identify highly repetitive events, suggestive of spatially distinct source locations. The waveform analysis shows that several different families of similar events co-exist during a given swarm day, but generally only one large family dominates. A network of hydrothermal fractures may explain the events that do not belong to a dominant repeating event group, i.e. multiple sources at different locations exist next to a dominant source. The dominant waveforms exhibit systematic changes throughout each swarm, but some of these waveforms do reappear over the course of 4 years indicating repeatedly activated source locations. The choked flow model provides a plausible trigger mechanism for the repeating events observed at Shishaldin, explaining the gradual changes in waveforms over time by changes in pressure gradient across a constriction within the uppermost part of the conduit. The sustained generation of Shishaldin's long-period events may be attributed to complex dynamics of a multi-fractured hydrothermal system: the pressure gradient within the main conduit may be regulated by temporarily sealing and reopening of parallel flow pathways, by the amount of debris within the main conduit and/or by changing gas influx into the hydrothermal system. The observations suggest that Shishaldin's swarms of repeating events represent time periods during which a dominant source is activated.

Alaska↗

USGS2010

Subevents of long-period seismicity: implications for hydrothermal dynamics during the 2004-2008 eruption of Mount St. Helens

One of the most striking aspects of seismicity during the 2004–2008 eruption of Mount St. Helens (MSH) was the precise regularity in occurrence of repetitive long-period (LP) or “drumbeat” events over sustained time periods. However, this precise regularity was not always observed, and at times the temporal occurrence of LP events became more random. In addition, accompanying the dominant LP class of events during the 2004–2008 MSH eruption, there was a near-continuous, randomly occurring series of smaller seismic events. These subevents are not always simply small-amplitude versions of the dominant LP class of events but appear instead to result from a separate random process only loosely coupled to the main LP source mechanism. We present an analysis of the interevent time and amplitude distributions of the subevents, using waveform cross correlation to separate LP events from the subevents. We also discuss seismic tremor that accompanied the 8 March 2005 phreatic explosion event at MSH. This tremor consists of a rapid succession of LPs and subevents triggered during the explosion, in addition to broadband noise from the sustained degassing. Immediately afterward, seismicity returned to the pre-explosion occurrence pattern. This triggering in relation to the rapid ejection of steam from the system, and subsequent return to pre-explosion seismicity, suggests that both seismic event types originated in a region of the subsurface hydrothermal system that was (1) in contact with the reservoir feeding the 8 March 2005 phreatic explosion but (2) not destroyed or drained by the explosion event. Finally, we discuss possible thermodynamic conditions in a pressurized hydrothermal crack that could give rise to seismicity. Pressure drop estimates for typical LP events are not generally large enough to perturb pure water in a shallow hydrothermal crack into an unstable state. However, dissolved volatiles such as CO 2 may lead to a more unstable system, increasing the seismogenic potential of a hydrothermal crack subject to rapid heat flux. The interaction of hydrothermal and magmatic systems beneath MSH in 2004–2008 thus appears able to explain a wide range of observed phenomena, including subevents, LP events, larger (M d > 2) events, and phreatic explosions.

Washington↗

USGS2007

Mangrove forest distributions and dynamics (1975–2005) of the tsunami-affected region of Asia

Aim We aimed to estimate the present extent of tsunami-affected mangrove forests and determine the rates and causes of deforestation from 1975 to 2005. Location Our study region covers the tsunami-affected coastal areas of Indonesia, Malaysia, Thailand, Burma (Myanmar), Bangladesh, India and Sri Lanka in Asia. Methods We interpreted time-series Landsat data using a hybrid supervised and unsupervised classification approach. Landsat data were geometrically corrected to an accuracy of plus-or-minus half a pixel, an accuracy necessary for change analysis. Each image was normalized for solar irradiance by converting digital number values to the top-of-the atmosphere reflectance. Ground truth data and existing maps and data bases were used to select training samples and also for iterative labelling. We used a post-classification change detection approach. Results were validated with the help of local experts and/or high-resolution commercial satellite data. Results The region lost 12% of its mangrove forests from 1975 to 2005, to a present extent of c . 1,670,000 ha. Rates and causes of deforestation varied both spatially and temporally. Annual deforestation was highest in Burma ( c . 1%) and lowest in Sri Lanka (0.1%). In contrast, mangrove forests in India and Bangladesh remained unchanged or gained a small percentage. Net deforestation peaked at 137,000 ha during 1990–2000, increasing from 97,000 ha during 1975–90, and declining to 14,000 ha during 2000–05. The major causes of deforestation were agricultural expansion (81%), aquaculture (12%) and urban development (2%). Main conclusions We assessed and monitored mangrove forests in the tsunami-affected region of Asia using the historical archive of Landsat data. We also measured the rates of change and determined possible causes. The results of our study can be used to better understand the role of mangrove forests in saving lives and property from natural disasters such as the Indian Ocean tsunami, and to identify possible areas for conservation, restoration and rehabilitation.

Journal of Biogeography↗

USGS2018

Landsat time series analysis of fractional plant cover changes on abandoned energy development sites

Oil and natural gas development in the western United States has increased substantially in recent decades as technological advances like horizontal drilling and hydraulic fracturing have made extraction more commercially viable. Oil and gas pads are often developed for production, and then capped, reclaimed, and left to recover when no longer productive. Understanding the rates, controls, and degree of recovery of these reclaimed well sites to a state similar to pre-development conditions is critical for energy development and land management decision processes. Here we use a multi-decadal time series of satellite imagery (Landsat 5, 1984–2011) to assess vegetation regrowth on 365 abandoned well pads located across the Colorado Plateau in Utah, Colorado, and New Mexico. We developed high-frequency time series of the Soil-Adjusted Total Vegetation Index (SATVI) for each well pad using the Google Earth Engine cloud computing platform. BFAST time-series models were used to fit temporal trends, identifying when vegetation was cleared from the site and the magnitudes and rates of vegetation change after abandonment. The time series metrics are used to calculate the relative fractional vegetation cover (RFVC) of each pad, a measure of post-abandonment vegetation cover relative to pre-drilling condition. Mean and median RFVC were 36% (s.d. 33%) and 26%, respectively, five years after abandonment, with one third of well pads having RFVC greater than 50%. Statistical analyses suggest that much of the high vegetation cover is associated with weedy invasive annual species such as cheatgrass (Bromus tectorum) and Russian thistle (Salsola spp.). Climate conditions and the year of abandonment also play a role, with increased cover in later years associated with a wetter period. Non-linear change at many pads suggests longer recovery times than would be estimated by linear extrapolation. New techniques implemented here address a complex response of cover change to soils, management, and climate over time, and can be extended to the operational monitoring of energy development across large areas.

International Journal of Applied Earth Observation↗

USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗

USGS2023

Machine-learning model to delineate sub-surface agricultural drainage from satellite imagery

Knowing subsurface drainage (tile-drain) extent is integral to understanding how landscapes respond to precipitation events and subsequent days of drying, as well as how soil characteristics and land management influence stream response. Consequently, a time series of tile-drain extent would inform one aspect of land management that complicates our ability to explain streamflow and water-quality as a function of climate variability or conservation management. We trained a UNet machine-learning model, a convolutional neural network designed to highlight objects of interest within an image, to delineate tile-drain networks in panchromatic satellite imagery without additional data on soils, topography, or historical tile-drain extent. This was done by training the model to match the accuracy of human experts manually tracing the surface representation of tile drains in satellite imagery. Our approach began with a library of images that were used to train and quantify the accuracy of the model, with model performance tested on imagery from two areas that were not used to train the model. Satellite imagery included acquisition dates from 2008 to 2020. Training imagery was from agricultural areas within the US Great Lakes basin. Validation imagery was from the upper Maumee River, tributary to western Lake Erie, and an Indiana, Ohio-River headwater tributary. Our analysis of the satellite imagery paired with meteorological and soil data found that during spring, a combination of relatively high solar radiation, intermediate soil-water content and bare fields enabled the best model performance. Each area of interest was heavily tile-drained, where better understanding the movement of water, nutrients, and sediment from fields to downstream water bodies is key to managing harmful algal blooms and hypoxia. The trained UNet model successfully identified tile drains visible in the validation imagery with an accuracy of 93%–96% and balanced accuracy of 52%–54%, similar to performance for training data (95% and 63%, respectively). Model performance will benefit from ongoing contributions to the training library.

Indiana, Michigan, Ohio, Wisconsin↗

USGS1997

Southern California Permanent GPS Geodetic Array: Continuous measurements of regional crustal deformation between the 1992 Landers and 1994 Northridge earthquakes

The southern California Permanent GPS Geodetic Array (PGGA) was established in 1990 across the Pacific-North America plate boundary to continuously monitor crustal deformation. We describe the development of the array and the time series of daily positions estimated for its first 10 sites in the 19-month period between the June 28, 1992 ( M w =7.3), Landers and January 17, 1994 ( M w =6.7), Northridge earthquakes. We compare displacement rates at four site locations with those reported by Feigl et al. [1993], which were derived from an independent set of Global Positioning System (GPS) and very long baseline interferometry (VLBI) measurements collected over nearly a decade prior to the Landers earthquake. The velocity differences for three sites 65–100 km from the earthquake's epicenter are of order of 3–5 mm/yr and are systematically coupled with the corresponding directions of coseismic displacement. The fourth site, 300 km from the epicenter, shows no significant velocity difference. These observations suggest large-scale postseismic deformation with a relaxation time of at least 800 days. The statistical significance of our observations is complicated by our incomplete knowledge of the noise properties of the two data sets; two possible noise models fit the PGGA data equally well as described in the companion paper by Zhang et al . [this issue]; the pre-Landers data are too sparse and heterogeneous to derive a reliable noise model. Under a fractal white noise model for the PGGA data we find that the velocity differences for all three sites are statistically different at the 99% significance level. A white noise plus flicker noise model results in significance levels of only 94%, 43%, and 88%. Additional investigations of the pre-Landers data, and analysis of longer spans of PGGA data, could have an important effect on the significance of these results and will be addressed in future work.

Journal of Geophysical Research B: Solid Earth↗

USGS2017

Riverine discharges to Chesapeake Bay: Analysis of long-term (1927–2014) records and implications for future flows in the Chesapeake Bay basin

The Chesapeake Bay (CB) basin is under a total maximum daily load (TMDL) mandate to reduce nitrogen, phosphorus, and sediment loads to the bay. Identifying shifts in the hydro-climatic regime may help explain observed trends in water quality. To identify potential shifts, hydrologic data (1927–2014) for 27 watersheds in the CB basin were analyzed to determine the relationships among long-term precipitation and stream discharge trends. The amount, frequency, and intensity of precipitation increased from 1910 to 1996 in the eastern U.S., with the observed increases greater in the northeastern U.S. than the southeastern U.S. The CB watershed spans the north-to-south gradient in precipitation increases, and hydrologic differences have been observed in watersheds north relative to watersheds south of the Pennsylvania—Maryland (PA-MD) border. Time series of monthly mean precipitation data specific to each of 27 watersheds were derived from the Precipitation-elevation Regression on Independent Slopes Model (PRISM) dataset, and monthly mean stream-discharge data were obtained from U.S. Geological Survey streamgage records. All annual precipitation trend slopes in the 18 watersheds north of the PA-MD border were greater than or equal to those of the nine south of that border. The magnitude of the trend slopes for 1927–2014 in both precipitation and discharge decreased in a north-to-south pattern. Distributions of the monthly precipitation and discharge datasets were assembled into percentiles for each year for each watershed. Multivariate correlation of precipitation and discharge within percentiles among the groups of northern and southern watersheds indicated only weak associations. Regional-scale average behaviors of trends in the distribution of precipitation and discharge annual percentiles differed between the northern and southern watersheds. In general, the linkage between precipitation and discharge was weak, with the linkage weaker in the northern watersheds compared to those in the south. On the basis of simple linear regression, 26 of the 27 watersheds are projected to have higher annual mean discharge in 2025, the target date for implementation of the TMDL for the CB basin.

Chesapeake Bay Basin↗

USGS1967

The ternary system K2SO4 MgSO4 CaSO4

Melting and subsolidus relations in the system K 2 SO 4 MgSO 4 CaSO 4 were studied using heating-cooling curves, differential thermal analysis, optics, X-ray diffraction at room and high temperatures and by quenching techniques. Previous investigators were unable to study the binary MgSO 4 CaSO 4 system and the adjacent area in the ternary system because of the decomposition of MgSO 4 and CaSO 4 at high temperatures. This problem was partly overcome by a novel sealed-tube quenching method, by hydrothermal synthesis, and by long-time heating in the solidus. As a result of this study, we found: (1) a new compound, CaSO 4 ·3MgSO 4 (m.p. 1201°C) with a field extending into the ternary system; (2) a high temperature form of MgSO 4 with a sluggishly reversible inversion. An X-ray diffraction pattern for this polymorphic form is given; (3) the inversion of β-CaSO 4 (anhydrite) to α-CaSO 4 at 1195°C, in agreement with grahmann ; (1) (4) the melting point of MgSO 4 is 1136°C and that of CaSO 4 is 1462°C (using sealed tube methods to prevent decomposition of the sulphates); (5) calcium langbeinite (K 2 SO 4 ·2CaSO 4 ) is the only compound in the K 2 SO 4 CaSO 4 binary system. This resolved discrepancies in the results of previous investigators; (6) a continuous solid solution series between congruently melting K 2 SOP 4 ·2MgSO 4 (langbeinite) and incongruently melting K 2 SO 4 ·2CaSO 4 (calcium langbeinite); (7) the liquidus in the ternary system consists of primary phase fields of K 2 SO 4 , MgSO 4 , CaSO 4 , langbeinite-calcium langbeinite solid solution, and CaSO 4 ·3MgSO 4 . The CaSO 4 field extends over a large portion of the system. Previously reported fields for the compounds (K 2 SO 4 ·MgSO 4 ·nCaSO 4 ), K 2 SO 4 ·3CaSO 4 and K 2 SO 4 ·CaSO 4 were not found; (8) a minimum in the ternary system at: 740°C, 25% MgSO 4 , 6% CaSO 4 , 69% K 2 SO 4 ; and ternary eutectics at 882°C, 49% MgSO 4 , 19% CaSO 4 , 32% K 2 SO 4 ; and 880°, 67·5% MgSO 4 , 5% CaSO 4 , 27·5% K 2 SO 4 .

Journal of Inorganic and Nuclear Chemistry↗

USGS1995

Graphical and PC-software analysis of volcano eruption precursors according to the Materials Failure Forecast Method (FFM)

The Materials Failure Forecasting Method for volcanic eruptions (FFM) analyses the rate of precursory phenomena. Time of eruption onset is derived from the time of “failure” implied by accelerating rate of deformation. The approach attempts to fit data, Ω, to the differential relationship Ω Ω ¨ = A Ω ˙ "> ¨ = AΩ˙ , where the dot superscript represents the time derivative, and the data Ω may be any of several parameters describing the accelerating deformation or energy release of the volcanic system. Rate coefficients, A and α, may be derived from appropriate data sets to provide an estimate of time to “failure”. As the method is still an experimental technique, it should be used with appropriate judgment during times of volcanic crisis. Limitations of the approach are identified and discussed. Several kinds of eruption precursory phenomena, all simulating accelerating creep during the mechanical deformation of the system, can be used with FFM. Among these are tilt data, slope-distance measurements, crater fault movements and seismicity. The use of seismic coda, seismic amplitude-derived energy release and time-integrated amplitudes or coda lengths are examined. Usage of cumulative coda length directly has some practical advantages to more rigorously derived parameters, and RSAM and SSAM technologies appear to be well suited to real-time applications. One graphical and four numerical techniques of applying FFM are discussed. The graphical technique is based on an inverse representation of rate versus time. For α = 2, the inverse rate plot is linear; it is concave upward for α < 2 and concave downward for α > 2. The eruption time is found by simple extrapolation of the data set toward the time axis. Three numerical techniques are based on linear least-squares fits to linearized data sets. The “linearized least-squares technique” is most robust and is expected to be the most practical numerical technique. This technique is based on an iterative linearization of the given rate-time series. The hindsight technique is disadvantaged by a bias favouring a too early eruption time in foresight applications. The “log rate versus log acceleration technique”, utilizing a logarithmic representation of the fundamental differential equation, is disadvantaged by large data scatter after interpolation of accelerations. One further numerical technique, a nonlinear least-squares fit to rate data, requires special and more complex software. PC-oriented computer codes were developed for data manipulation, application of the three linearizing numerical methods, and curve fitting. Separate software is required for graphing purposes. All three linearizing techniques facilitate an eruption window based on a data envelope according to the linear least-squares fit, at a specific level of confidence, and an estimated rate at time of failure.

Journal of Volcanology and Geothermal Research↗

USGS2013

Photogrammetric monitoring of lava dome growth during the 2009 eruption of Redoubt Volcano

The 2009 eruption of Redoubt Volcano, Alaska, began with a phreatic explosion on 15 March followed by a series of at least 19 explosive events and growth and destruction of at least two, and likely three, lava domes between 22 March and 4 April. On 4 April explosive activity gave way to continuous lava effusion within the summit crater. We present an analysis of post-4 April lava dome growth using an oblique photogrammetry approach that provides a safe, rapid, and accurate means of measuring dome growth. Photogrammetric analyses of oblique digital images acquired during helicopter observation flights and fixed-wing volcanic gas surveys produced a series of digital elevation models (DEMs) of the lava dome from 16 April to 23 September. The DEMs were used to calculate estimates of volume and time-averaged extrusion rates and to quantify morphological changes during dome growth. Effusion rates ranged from a maximum of 35 m 3 s − 1 during the initial two weeks to a low of 2.2 m 3 s − 1 in early summer 2009. The average effusion rate from April to July was 9.5 m 3 s − 1 . Early, rapid dome growth was characterized by extrusion of blocky lava that spread laterally within the summit crater. In mid-to-late April the volume of the dome had reached 36 × 10 6 m 3 , roughly half of the total volume, and dome growth within the summit crater began to be limited by confining crater walls to the south, east, and west. Once the dome reached the steep, north-sloping gorge that breaches the crater, growth decreased to the south, but the dome continued to inflate and extend northward down the gorge. Effusion slowed during 16 April–1 May, but in early May the rate increased again. This rate increase was accompanied by a transition to exogenous dome growth. From mid-May to July the effusion rate consistently declined. The decrease is consistent with observations of reduced seismicity, gas emission, and thermal anomalies, as well as declining rates of geodetic deflation or inflation. These trends suggest dome growth ceased by July 2009. The volume of the dome at the end of the 2009 eruption was about 72 × 10 6 m 3 , more than twice the estimated volume of the largest dome extruded during the 1989–1990 eruption. In total, the 2009 dome extends over 400 m down the glacial gorge on the north end of the crater, with a total length of 1 km, width of 500 m and an average thickness of 200 m.

Alaska↗

USGS1965

Effects of heptachlor-contaminated earthworms on woodcocks

The effects on woodcocks ( Philohela minor ) of eating heptachlor-contaminated earthworms were studied experimentally in a series of feeding trials in Louisiana in the winter of 1960-61. Six of 12 woodcocks fed worms which had been contaminated at an average of 2.86 ppm of heptachlor epoxide died within 35 days; 4 more had died by the fifty-third day, when the other 2 were killed for analysis. Worms from areas in Louisiana treated with 2 pounds of heptachlor per acre often contain more than 3 ppm of heptachlor epoxide. Eleven of 12 woodcocks fed worms contaminated at an average of 0.65 ppm survived the full 60 days of the experiment; one died on the forty-fifth day, apparently from other causes. All 11 untreated birds survived. Survivors were kept on one-quarter rations of untreated worms for 11 days. Two woodcocks, untreated previously, died during this starvation period. Five previously treated died; two were observed in spasms at death, and these contained 5.9 and 7.2 ppm heptachlor epoxide in their tissues, suggesting that the previous contaminated diet may have influenced mortality, even though the difference between two of nine dying and five of nine dying is not statistically significant. Surviving starved birds given an unrestricted supply of treated or untreated worms for 5 days survived and gained weight. Residues accumulated in their tissues in this time approached levels in birds that died of heptachlor poisoning. Residues in tissues of birds with different histories suggested residue loss at a rate of approximately 2.8 percent per day. Toxicant absorption was estimated to be in the approximate range of 16-20 percent. Residues in birds fed worms containing 0.65 ppm heptachlor epoxide were in the same general magnitude as those in field-caught birds, suggesting a similar average contamination of food supply. Weights and weight changes did not differ significantly between untreated birds and those receiving the lower level of toxicant. Among birds on one-quarter rations, the percentage of weight that could be lost without danger seemed to be near 20 percent. Woodcocks ate 18-208 grams of worms per day (average, 121 grams), representing 11-143 percent (average, 77 percent) of their body weights; birds ate contaminated and uncontaminated food in essentially equivalent amounts. Symptoms of heptachlor poisoning differed considerably between birds.

Journal of Wildlife Management↗

USGS2012

Cretaceous-Eocene (Laramide) landscape development and Oligocene- Pliocene drainage reorganization of transition zone and Colorado Plateau, Arizona

Landscape development of central and northern Arizona can no longer be ascribed mainly to events of Miocene and Pliocene age. New information on the age and distribution of older Cenozoic deposits has led to the recognition of a regional Cretaceous-Paleocene(?) surface of erosion that conforms to major elements of the present topography and to the recognition that a formerly thick deposit of gravel accumulated on this regional surface of erosion. These relations cast new light on the history of evolution of the landscape and indicate a much greater age for the main landscape elements and a more complicated and prolonged history of erosion and deposition than has been previously supposed. The timing of events postulated for development of drainage on the Colorado Plateau can now be compared and partly reconciled with events recognized in the adjacent closely related Mountain Region (Transition Zone) of central Arizona. As a consequence of Late Cretaceous-Paleocene (Laramide) compression, central and northern Arizona underwent at least 1200 m of uplift, documented by paleochannels cut into erosionally truncated Paleozoic strata on the Hualapai Plateau of the southwestern Colorado Plateau. During this time, a highly irregular erosion surface was developed on Proterozoic rocks across the Transition Zone south of the Mogollon Rim, the scarp of the Mogollon Rim was eroded to its present height (600–900 m), and an extensive stripped surface was developed on resistant upper Paleozoic strata north of the rim. Deposition of several hundred meters of Paleocene-Eocene “Rim gravels” derived from highlands south and west of the region followed, covering much of the Cretaceous-Paleocene erosion surface. Nearly complete burial of the rim is suggested by the distribution of remnants of the Rim gravels across the erosional scarps and on high plateau areas north of the rim. A second increment of uplift, apparently occurring in late Eocene time and apparently recorded by a series of fission track cooling ages from the Marble and Grand canyons, is inferred to have been responsible for ending deposition of the Rim gravels, for initiating differential uplift of contemporaneous deposits (Canaan Peak and Claron formations) to their positions in the high plateaus of central Utah, and for causing the drainage reorganization required to explain the extensive removal of Rim gravels from much of the region. A southerly flowing ancestral Verde River related to the drainage reorganization removed much of the older gravel cover from the Transition Zone of central Arizona, resulting in a younger regional erosion surface having 600–900 m of relief, a surface closely approximating the Cretaceous-Paleocene erosion surface. Late Oligocene and early Miocene rocks locally rest unconformably on remnants of Rim gravels in the Transition Zone, indicating that the second episode of regional erosion had been completed by late Oligocene time. North of the Mogollon Rim, a west flowing(?) ancestral Colorado River is inferred to have become established on the Rim gravels, draining the interior parts of the Colorado Plateau and transporting detritus off the plateau. Exhumation of the Mogollon Rim and development of 600–900 m of topographic relief in the Transition Zone by an ancestral Verde River system suggests the potential for a comparable, coeval entrenchment of an ancestral Colorado River in Paleozoic strata north of the Mogollon Rim. Regional extension and volcanic activity ensued in late Oligocene to Pliocene time. The Oligocene erosion surface in the extensional basins of central Arizona became largely concealed by Miocene and Pliocene deposits as the Neogene climate became drier. In late Miocene and Pliocene time, perennial streams appear to have been lacking, transport of detritus appears to have been principally by flash flooding, little or no detritus appears to have been removed from the region, and much of the precipitation presumably moved by groundwater flow through the deposits of aggradation. A coeval episode of aggradation in the Grand Canyon is suggested by deposits that appear to have once choked much of the canyon. If this event parallels the episode of late Miocene and Pliocene aggradation recorded east, south, and west of the Grand Canyon, the Colorado River could have been incised to its present level by late Miocene time. A return to wetter conditions in late Pliocene time presumably was responsible for renewed erosion and reexcavation of older drainages and basins. An understanding of this Tertiary structural, erosional, and depositional history can be important for the geological analysis of geophysical transects across the region.

Journal of Geophysical Research↗

USGS2020

Combining ash analyses with remote sensing to identify juvenile magma involvement and fragmentation mechanisms during the 2018/19 small eruption of Peteroa volcano (Southern Andes)

The Planchón Peteroa Volcanic Complex (PPVC) is located on the border of Chile and Argentina, and is one of the most active volcanic systems in the Andes. Holocene activity has included magma-water interaction with an evolving series of crater lakes, mainly sourced from Peteroa volcano. This study examines data from the 2018/19 eruption, together with the volcanic history of the PPVC, to elucidate the complex interplay between magmatic activity and summit water and ice. From February 2016 to mid-2019, three seismic swarms occurred in the PPVC, preceding the explosive eruption from September 2018 to April 2019. The activity originated from a small vent nested within the easternmost crater, the most active portion of the complex (Peteroa). The explosions interacted with a crater lake, producing ash plumes up to 2 km above the crater and building a small tephra cone. To investigate the eruption mechanisms, we performed remote sensing analysis of plume dispersal, thermal anomalies and ground deformation, and characterized the volcanic products, including grain size, componentry, morphology, internal textures, composition and mineralogy. Our results suggest that the precursory seismicity beginning in 2016 was related to the intrusion of a new magma batch that reached the surface during the 2018/19 eruption. The eruption was also preceded by thermal anomalies, geomorphic changes and increased hydrothermal activity at the surface, though without any ground deformation recognized through radar interferometry (InSAR). The eruption initially produced predominantly recycled ash (phreatic activity), then evolved to increasing proportions of juvenile magma (phreatomagmatic) by April 2019. The juvenile clasts had a trachyandesite composition (~59 wt% SiO 2 ), with vesicular and dense scoria containing plagioclase and pyroxene. The ash surfaces show external quenching cracks and step fractures consistent with phreatomagmatic fragmentation within the active crater lake. Textural characteristics also point to a slowly ascending batch of magma that was relatively viscous by the time it interacted with water in the crater lake. Notably, these juvenile particles are distinctive from the pre-2018 products. Ash erupted from 2010/11 did not contain recognizable juvenile material, and is inferred to have been a mainly phreatic eruption. Our findings suggest that the interplay between phreatic and phreatomagmatic eruptions fed by small magma batches intruding at shallow levels characterize much of the eruptive behavior of the PPVC during the last three decades. Multi-parametric assessment is a powerful tool to discriminate between phreatic and phreatomagmatic eruptions.

Southern Andes↗

USGS2025

Analysis of a human-mediated microbioinvasion: The global spread of the benthic foraminifer Trochammina hadai Uchio, 1962

A non-indigenous species (NIS) of benthic foraminifera was first identified in a core collected in 1993 in San Francisco Bay, California, USA, and subsequently identified as Trochammina hadai Uchio, 1962. Archived samples and literature reviews were used to determine that the species, which is native to Asia, arrived in San Francisco Bay between the early 1960s and 1983. Through molecular analyses of specimens, archived samples and literature reviews from 1930–1983, and site surveys of harbors and estuaries along the western North American seaboard in 1994–2024, in total more than 2500 samples, we documented the presence of T. hadai at 73 locations in the USA and four in Canada. Trochammina hadai has also been recovered at nine sites in Sweden, two in France, three in Brazil, and two locations at one site in Australia. The rapid temporal and geographic spread of the NIS T. hadai in a non-native location is illustrated by a time series from 1930 to 2024 in San Francisco Bay. Between 1980 and 1986, the species' range expanded from low abundance (1.5 %) at a single site to cover nearly the entire South Bay with > 70 % abundance at some locations. By 1995 and continuing into 2010, the species expanded its range into the central and northern portions of San Francisco Bay, commonly with abundances of > 30 % and sometimes exceeding 70 %. This expansion may predate 1995, but a lack of samples makes it difficult to be more precise. Unfortunately, two Pb-210 and Cs-137-dated cores (BC01 and BC02) recovered from northern South Bay and Central Bay did not clarify this point, but additional cores may. Trochammina hadai is an infaunal opportunist that thrives in polluted locations. We surmise the species was introduced along the west coast of the USA in Puget Sound between 1902 and the 1920s, with cultivated oysters and oyster larvae and associated plant matter and residual sediment. This probably also happened in some areas of France, Sweden, and Brazil, where Japanese oysters were introduced in 1966, 1970, and 1975, respectively. After World War II, commercial shipping expanded dramatically and, with it, the release of ballast water and sediment in receiving ports, which introduced NIS worldwide. This primary vector of introduction occurred in large industrial harbors in several countries, sometimes followed by secondary introductions in small industrial centers and marinas by mud attached to the anchors and anchor chains of smaller boats.

Journal of Micropalaeontology↗

USGS2022

Guidelines for volcano-observatory operations during crises: Recommendations from the 2019 Volcano Observatory Best Practices meeting

In November 2019, the fourth meeting on Volcano Observatory Best Practices workshop was held in Mexico City as a series of talks, discussions, and panels. Volcanologists from around the world offered suggestions for ways to optimize volcano-observatory crisis operations. By crisis, we mean unrest that may or may not lead to eruption, the eruption itself, or its aftermath, all of which require analysis and communications by the observatory. During a crisis, the priority of the observatory should be to acquire, process, analyze, and interpret data in a timely manner. A primary goal is to communicate effectively with the authorities in charge of civil protection. Crisis operations should rely upon exhaustive planning in the years prior to any actual unrest or eruptions. Ideally, nearly everything that observatories do during a crisis should be envisioned, prepared, and practiced prior to the actual event. Pre-existing agreements and exercises with academic and government collaborators will minimize confusion about roles and responsibilities. In the situation where planning is unfinished, observatories should prioritize close ties and communications with the land and civil-defense authorities near the most threatening volcanoes. To a large extent, volcanic crises become social crises, and any volcano observatory should have a communication strategy, a lead communicator, regular status updates, and a network of colleagues outside the observatory who can provide similar messaging to a public that desires consistent and authoritative information. Checklists permit tired observatory staff to fulfill their duties without forgetting key communications, data streams, or protocols that need regular fulfilment (Bretton et al. 2018; Newhall et al. 2020). Observatory leaders need to manage staff workload to prevent exhaustion and ensure that expertise is available as needed. Event trees and regular group discussions encourage multi-disciplinary thinking, consideration of disparate viewpoints, and documentation of all group decisions and consensus. Though regulations, roles and responsibilities differ around the world, scientists can justify their actions in the wake of an eruption if they document their work, are thoughtful and conscientious in their deliberations, and carry out protocols and procedures developed prior to volcanic unrest. This paper also contains six case studies of volcanic eruptions or observatory actions that illustrate some of the topics discussed herein. Specifically, we discuss Ambae (Vanuatu) in 2017–2018, Kīlauea (USA) in 2018, Etna (Italy) in 2018, Bárðarbunga (Iceland) in 2014, Cotopaxi (Ecuador) in 2015, and global data sharing to prepare for eruptions at Nyiragongo (Democratic Republic of Congo).

Hawaii↗

USGS2019

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation↗

USGS2014

The Lepanto Cu–Au deposit, Philippines: A fossil hyperacidic volcanic lake complex

Hyperacidic lakes and associated solfatara in active volcanoes are the expression of magmatic gas expansion from source to surface. Here we show for the first time, that the vein system that comprises the ~ 2 Ma high-sulfidation, Lepanto copper–gold deposit in the Mankayan district (Philippines) was associated with a contemporary hyperacidic volcanic lake complex—possibly the first such lake recognized in the geological record. A 15–20‰ difference in sulfur isotopic composition between barite and sulfides and sulfosalts in the vent fumarole encrustations supports the interpretation that SO 2 -rich volcanic gas vented into the base of the lake and marginal to it and ties the mineralization directly to magmatic gas expansion, fracture propagation, and mineralization that occurred through a series of decompression steps within the feeder fracture network. These data confirm that crater lake environments such as Kawah Ijen (Java, Indonesia) provide modern day analogs of the Lepanto and other high sulfidation Cu–Au depositing environments. We also provide extensive analysis of sulfosalt–sulfide reactions during vein formation within the hyperacidic lake complex. Pyrite ± silica deposited first at high temperature followed by enargite that preserves the vapor–solid diffusion of, for example, antimony, tin, and tellurium into the vapor from the crystallizing solid. Subsolidus, intra-crystalline diffusion continued as temperature declined. Pyrite and enargite are replaced by Fe-tennantite in the lodes which initially has low Sb/(Sb + As) atomic ratios around 13.5% close to the ideal tennantite formula, but evolves to higher ratios as crystallization proceeds. Fumarole encrustation clasts and sulfosalts in the lake sediment are more highly evolved with a larger range of trace element substitutions, including antimony. Substitution of especially Zn, Te, Ag, and Sn into tennantite records metal and semi-metal fractionation between the expanding magmatic gas and deposited sulfide sublimates provides a rare insight into the fate of metals and semi-metals in the shallower parts of fracture arrays that feed modern hyperacidic lakes. These data support a growing understanding of the formation of high-sulfidation gold deposits as the consequence of single-phase expansion of gas from magmatic-gas reservoirs beneath the surface of active volcanoes without the intervention of a later aqueous fluid including groundwater. Aggressive sulfide–sulfosalt reactions, including pitting and the almost complete dissolution of earlier minerals, are persistent characteristics of the vein assemblages and precious metals typically occur late in pits or along brittle fractures. These characteristics support a hypothesis of mineral deposition at temperatures of the order of 600 °C in contrast to available fluid inclusion data from enargite that record temperatures following phase transitions in the sulfosalt during the retrograde devolution of the deposit in the presence of groundwater.

Benguet↗