Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth in Space”

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 757 records · Page 42Linked to original sources
USGS2013

Topological and canonical kriging for design flood prediction in ungauged catchments: an improvement over a traditional regional regression approach?

In the United States, estimation of flood frequency quantiles at ungauged locations has been largely based on regional regression techniques that relate measurable catchment descriptors to flood quantiles. More recently, spatial interpolation techniques of point data have been shown to be effective for predicting streamflow statistics (i.e., flood flows and low-flow indices) in ungauged catchments. Literature reports successful applications of two techniques, canonical kriging, CK (or physiographical-space-based interpolation, PSBI), and topological kriging, TK (or top-kriging). CK performs the spatial interpolation of the streamflow statistic of interest in the two-dimensional space of catchment descriptors. TK predicts the streamflow statistic along river networks taking both the catchment area and nested nature of catchments into account. It is of interest to understand how these spatial interpolation methods compare with generalized least squares (GLS) regression, one of the most common approaches to estimate flood quantiles at ungauged locations. By means of a leave-one-out cross-validation procedure, the performance of CK and TK was compared to GLS regression equations developed for the prediction of 10, 50, 100 and 500 yr floods for 61 streamgauges in the southeast United States. TK substantially outperforms GLS and CK for the study area, particularly for large catchments. The performance of TK over GLS highlights an important distinction between the treatments of spatial correlation when using regression-based or spatial interpolation methods to estimate flood quantiles at ungauged locations. The analysis also shows that coupling TK with CK slightly improves the performance of TK; however, the improvement is marginal when compared to the improvement in performance over GLS.

Hydrology and Earth System Sciences↗
USGS2017

Analysis of the variability in ground-motion synthesis and inversion

In almost all past inversions of large-earthquake ground motions for rupture behavior, the goal of the inversion is to find the “best fitting” rupture model that predicts ground motions which optimize some function of the difference between predicted and observed ground motions. This type of inversion was pioneered in the linear-inverse sense by Olson and Apsel (1982), who minimized the square of the difference between observed and simulated motions (“least squares”) while simultaneously minimizing the rupture-model norm (by setting the null-space component of the rupture model to zero), and has been extended in many ways, one of which is the use of nonlinear inversion schemes such as simulated annealing algorithms that optimize some other misfit function. For example, the simulated annealing algorithm of Piatanesi and others (2007) finds the rupture model that minimizes a “cost” function which combines a least-squares and a waveform-correlation measure of misfit. All such inversions that look for a unique “best” model have at least three problems. (1) They have removed the null-space component of the rupture model—that is, an infinite family of rupture models that all fit the data equally well have been narrowed down to a single model. Some property of interest in the rupture model might have been discarded in this winnowing process. (2) Smoothing constraints are commonly used to yield a unique “best” model, in which case spatially rough rupture models will have been discarded, even if they provide a good fit to the data. (3) No estimate of confidence in the resulting rupture models can be given because the effects of unknown errors in the Green’s functions (“theory errors”) have not been assessed. In inversion for rupture behavior, these theory errors are generally larger than the data errors caused by ground noise and instrumental limitations, and so overfitting of the data is probably ubiquitous for such inversions. Recently, attention has turned to the inclusion of theory errors in the inversion process. Yagi and Fukahata (2011) made an important contribution by presenting a method to estimate the uncertainties in predicted large-earthquake ground motions due to uncertainties in the Green’s functions. Here we derive their result and compare it with the results of other recent studies that look at theory errors in a Bayesian inversion context particularly those by Bodin and others (2012), Duputel and others (2012), Dettmer and others (2014), and Minson and others (2014). Notably, in all these studies, the estimates of theory error were obtained from theoretical considerations alone; none of the investigators actually measured Green’s function errors. Large earthquakes typically have aftershocks, which, if their rupture surfaces are physically small enough, can be considered point evaluations of the real Green’s functions of the Earth. Here we simulate smallaftershock ground motions with (erroneous) theoretical Green’s functions. Taking differences between aftershock ground motions and simulated motions to be the “theory error,” we derive a statistical model of the sources of discrepancies between the theoretical and real Green’s functions. We use this model with an extended frequency-domain version of the time-domain theory of Yagi and Fukahata (2011) to determine the expected variance 2 τ caused by Green’s function error in ground motions from a larger (nonpoint) earthquake that we seek to model. We also differ from the above-mentioned Bayesian inversions in our handling of the nonuniqueness problem of seismic inversion. We follow the philosophy of Segall and Du (1993), who, instead of looking for a best-fitting model, looked for slip models that answered specific questions about the earthquakes they studied. In their Bayesian inversions, they inductively derived a posterior probability-density function (PDF) for every model parameter. We instead seek to find two extremal rupture models whose ground motions fit the data within the error bounds given by 2 τ , as quantified by using a chi-squared test described below. So, we can ask questions such as, “What are the rupture models with the highest and lowest average rupture speed consistent with the theory errors?” Having found those models, we can then say with confidence that the true rupture speed is somewhere between those values. Although the Bayesian approach gives a complete solution to the inverse problem, it is computationally demanding: Minson and others (2014) needed 1010 forward kinematic simulations to derive their posterior probability distribution. In our approach, only about107 simulations are needed. Moreover, in practical application, only a small set of rupture models may be needed to answer the relevant questions—for example, determining the maximum likelihood solution (achievable through standard inversion techniques) and the two rupture models bounding some property of interest. The specific property that we wish to investigate is the correlation between various rupturemodel parameters, such as peak slip velocity and rupture velocity, in models of real earthquakes. In some simulations of ground motions for hypothetical large earthquakes, such as those by Aagaard and others (2010) and the Southern California Earthquake Center Broadband Simulation Platform (Graves and Pitarka, 2015), rupture speed is assumed to correlate locally with peak slip, although there is evidence that rupture speed should correlate better with peak slip speed, owing to its dependence on local stress drop. We may be able to determine ways to modify Piatanesi and others’s (2007) inversion’s “cost” function to find rupture models with either high or low degrees of correlation between pairs of rupture parameters. We propose a cost function designed to find these two extremal models.

Open-File Report↗
USGS1996

Channels and valleys on Mars: Cold climate features formed as a result of a thickening cryosphere

Large flood channels, valley networks, and a variety of features attributed to the action of ground ice indicate that Mars emerged from heavy bombardment around 3.8Gyr ago, with an inventory of water at the surface equivalent to at least a few hundred meters spread over the whole planet, as compared with 3 km for the Earth. The surface water resided primarily in a porous, kilometers thick, megaregolith created by the high impact rates. At the end of heavy bombardment a rapid decline in erosion rates by a factor of 1000 suggests a major change in the global climate. It is proposed that at this time the climate became similar to today's and that this climate has been maintained throughout the rest of Mars' history. The various drainage features represent an adjustment of the distribution of water to the surface relief inherited from the period of heavy bombardment and to a thickening of the cryosphere as the heat flow declined. The valley networks formed mostly at the end of heavy bombardment when erosion rates were high and climatic conditions permitted an active water cycle. They continued to form after heavy bombardment when the cryosphere started to form by a combination of episodic flooding and mass-wasting aided by the presence of liquid water at shallow depths. As the cryosphere thickened with declining heat flow, water could no longer easily access the surface and the rate of valley formation declined. Hydrostatic pressures built below the cryosphere. Eruptions of groundwater became more catastrophic and massive floods resulted, mainly in upper Hesperian time. Flood sources were preferentially located in low-lying, low-latitude areas where the cryosphere was thin, or near volcanoes where a thinner than typical cryosphere is also expected. Floods caused a drawdown in the global water table so that few formed in the second half of Mars' history. The floodwaters pooled in low-lying areas, mostly in the northern plains. Some of the water may still be present as thick ice deposits, some has been lost to space, particularly during periods of high obliquity.

Planetary and Space Science↗
USGS2017

Quantifying the heterogeneity of the tectonic stress field using borehole data

The heterogeneity of the tectonic stress field is a fundamental property which influences earthquake dynamics and subsurface engineering. Self-similar scaling of stress heterogeneities is frequently assumed to explain characteristics of earthquakes such as the magnitude-frequency relation. However, observational evidence for such scaling of the stress field heterogeneity is scarce. We analyze the local stress orientations using image logs of two closely spaced boreholes in the Coso Geothermal Field with sub-vertical and deviated trajectories, respectively, each spanning about 2 km in depth. Both the mean and the standard deviation of stress orientation indicators (borehole breakouts, drilling-induced fractures and petal-centerline fractures) determined from each borehole agree to the limit of the resolution of our method although measurements at specific depths may not. We find that the standard deviation in these boreholes strongly depends on the interval length analyzed, generally increasing up to a wellbore log length of about 600 m and constant for longer intervals. We find the same behavior in global data from the World Stress Map. This suggests that the standard deviation of stress indicators characterizes the heterogeneity of the tectonic stress field rather than the quality of the stress measurement. A large standard deviation of a stress measurement might be an expression of strong crustal heterogeneity rather than of an unreliable stress determination. Robust characterization of stress heterogeneity requires logs that sample stress indicators along a representative sample volume of at least 1 km.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

The added value of time-variable microgravimetry to the understanding of how volcanoes work

During the past few decades, time-variable volcano gravimetry has shown great potential for imaging subsurface processes at active volcanoes (including some processes that might otherwise remain “hidden”), especially when combined with other methods (e.g., ground deformation, seismicity, and gas emissions). By supplying information on changes in the distribution of bulk mass over time, gravimetry can provide information regarding processes such as magma accumulation in void space, gas segregation at shallow depths, and mechanisms driving volcanic uplift and subsidence. Despite its potential, time-variable volcano gravimetry is an underexploited method, not widely adopted by volcano researchers or observatories. The cost of instrumentation and the difficulty in using it under harsh environmental conditions is a significant impediment to the exploitation of gravimetry at many volcanoes. In addition, retrieving useful information from gravity changes in noisy volcanic environments is a major challenge. While these difficulties are not trivial, neither are they insurmountable; indeed, creative efforts in a variety of volcanic settings highlight the value of time-variable gravimetry for understanding hazards as well as revealing fundamental insights into how volcanoes work. Building on previous work, we provide a comprehensive review of time-variable volcano gravimetry, including discussions of instrumentation, modeling and analysis techniques, and case studies that emphasize what can be learned from campaign, continuous, and hybrid gravity observations. We are hopeful that this exploration of time-variable volcano gravimetry will excite more scientists about the potential of the method, spurring further application, development, and innovation.

Earth-Science Reviews↗
USGS2022

Barkley Canyon gas hydrates: A synthesis based on two decades of seafloor observation and remote sensing

Barkley Canyon is one of the few known sites worldwide with the occurrence of thermogenic gas seepage and formation of structure-II and structure-H gas hydrate mounds on the seafloor. This site is the location of continuous seafloor monitoring as part of the Ocean Networks Canada (ONC) cabled observatory off the west coast off Vancouver Island, British Columbia, Canada. We combine repeat remotely operated vehicle (ROV) seafloor video observations, mapping with an autonomous underwater vehicle (AUV), ship-, ROV-, and AUV-based identification of gas flares, as well as seismic and Chirp data to investigate the distribution of fluid migration pathways. Geologically, the site with the prominent gas hydrate mounds and associated fluid seepage is covering an area of ∼0.15 km 2 and is situated on a remnant of a rotated fault block that had slipped off the steep flanks of the north-east facing canyon wall. The gas hydrate mounds, nearly constant in dimension over the entire observation period, are associated with gas and oil seepage and surrounded by debris of chemosynthetic communities and authigenic carbonate. The formation of gas hydrate at and near the seafloor requires additional accommodation space created by forming blisters at the seafloor that displace the regular sediments. An additional zone located centrally on the rotated fault block with more diffuse seepage (∼0.02 km 2 in extent) has been identified with no visible mounds, but with bacterial mats, small carbonate concretions, and clam beds. Gas venting is seen acoustically in the water column up to a depth of ∼300 m. However, acoustic water-column imaging during coring and ROV dives showed rising gas bubbles to much shallower depth, even <50 m, likely a result of degassing of rising oil droplets, which themselves cannot be seen acoustically. Combining all observations, the location of the gas hydrate mounds is controlled by a combination of fault-focused fluid migration from a deeper reservoir and fluid seepage along more permeable strata within the rotated slope block. Fluids must be provided continuously to allow the sustained presence of the gas hydrate mounds at the seafloor.

Barkley Canyon↗
USGS2015

A 2-D process-based model for suspended sediment dynamics: A first step towards ecological modeling

In estuaries suspended sediment concentration (SSC) is one of the most important contributors to turbidity, which influences habitat conditions and ecological functions of the system. Sediment dynamics differs depending on sediment supply and hydrodynamic forcing conditions that vary over space and over time. A robust sediment transport model is a first step in developing a chain of models enabling simulations of contaminants, phytoplankton and habitat conditions. This works aims to determine turbidity levels in the complex-geometry delta of the San Francisco estuary using a process-based approach (Delft3D Flexible Mesh software). Our approach includes a detailed calibration against measured SSC levels, a sensitivity analysis on model parameters and the determination of a yearly sediment budget as well as an assessment of model results in terms of turbidity levels for a single year, water year (WY) 2011. Model results show that our process-based approach is a valuable tool in assessing sediment dynamics and their related ecological parameters over a range of spatial and temporal scales. The model may act as the base model for a chain of ecological models assessing the impact of climate change and management scenarios. Here we present a modeling approach that, with limited data, produces reliable predictions and can be useful for estuaries without a large amount of processes data.

Hydrology and Earth System Sciences↗
USGS2023

Geochronology and mapping constraints on the time-space evolution of the igneous and hydrothermal systems in the Taurus Cu-Mo district, eastern Alaska

The Taurus porphyry Cu-Mo district contains four mineralized porphyry centers in the eastern interior of Alaska. All four centers were emplaced during a magmatic episode that spanned from ca. 72 to 67 Ma, with seven distinct igneous suites. Each igneous suite resulted in hydrothermal alteration and mineralization, with younger pulses overprinting older pulses. Each magmatic-hydrothermal system is not present at all four mineralized centers. Apart from the Dennison occurrence, each mineralized center records pulses of repeated intermediate-silicic magmatism and associated alteration and mineralization. Laser ablation-inductively coupled plasma-mass spectrometry U-Pb zircon crystallization ages indicate that an early quartz porphyry dike swarm ranges in age from ca. 71 to 70 Ma and is associated with potassic, sericitic, and propylitic alteration. Quartz latite intrusions were emplaced at ca. 69 Ma and exhibit early sodiccalcic alteration overprinted by potassic, sericitic, and propylitic alteration. The Taurus monzonite suite is cut by quartz latite but yielded an ca. 70 Ma emplacement age and exhibits the largest footprint of potassic and sericitic alteration. Feldspar porphyry dikes were emplaced ca. 69 Ma and have significant tourmaline-bearing potassic and sericitic alteration. This suite was followed by development of an igneous breccia with a monzonitic igneous matrix. Sodic-calcic alteration was associated with the igneous brecciation. A small stock of monzonite was emplaced at ca. 68 Ma causing locally pervasive sericite-tourmaline-pyrite alteration. The youngest suite of magmatism dated in the district is a series of granodiorite porphyry dikes with weak sodic-calcic and propylitic alteration that truncates earlier alteration assemblages. Mineralization in the district consists of chalcopyrite and molybdenite associated with sugary quartz veins with potassium feldspar and biotite alteration envelopes (A veins). Less common banded quartz-molybdenite veins (B veins) occur with potassium feldspar envelopes. Gold occurs throughout the district and is strongly correlated with copper grade. Sericitic alteration contains lower copper contents and is predominantly associated with quartz-pyrite veins with sericite envelopes (D veins). Pyrrhotite and local arsenopyrite are present in sericitic assemblages. Pyrrhotite also occurs as inclusions in pyrite within D veins. Magmas across the district exhibit oxidized characteristics, evidenced by the presence of abundant magnetite, rare titanite, and elevated Eu/Eu* and Ce/Ce* in zircon. Zircon Th/U and Yb/Gd compositions suggest a fractionation path controlled by apatite, titanite, and hornblende. Zircon rare earth element ratios and trace element data indicate two distinct batches of magma evolved from mafic parental compositions to monzonite and granodioritic compositions via fractional crystallization. In the early pulse of magma (ca. 72–69 Ma), fractional crystallization was key to ore formation. Earlier, better mineralized suites evolve to less negative Eu anomalies (Eu/Eu* > 0.7), indicating more oxidized and higher-water-pressure conditions evidenced by the suppression of plagioclase crystallization, compared to later, more poorly mineralized suites. The temporal and spatial evolution of the district was determined from mapping and U-Pb and Re-Os geochronology. Mapping of igneous and hydrothermal assemblages indicates that the locus of the intrusive suites and hydrothermal systems shifted spatially over time, based on the presence of high-temperature (K-silicate–dominant) alteration, which is coincident with the highest Cu and Au grades. The earliest hydrothermal system was centered at Bluff and East Taurus and transitioned to West Taurus during emplacement of the second magmatic suite. Emplacement of the third magmatic suite was centered back at East Taurus, and the fourth and fifth suites were centered at West Taurus. The latest suites were widespread without a core of high-temperature alteration marking a central locus. East Taurus contains the overlap of six of the seven magmatic and hydrothermal suites and has the highest intersected grades and tonnages in the district. The Bluff and Dennison occurrences exhibit fewer igneous suites and hydrothermal assemblages with weak mineralization. Sodic-calcic alteration, common on the deep and distal flanks of porphyry systems, is only present at West Taurus and is indicative of a localized source of high-salinity nonmagmatic fluids.

Alaska, Yukon Territory↗
USGS2024

Wildfire, extreme precipitation and debris flows, oh my! Channel response to compounding disturbances in a mountain stream in the Upper Colorado Basin, USA

Compounding changes in climate and human activities stand to increase sediment input to rivers in many landscapes, including via discrete perturbations such as post-fire debris flows. Because sediment supply is a dominant control on river morphology, understanding mountain river responses to sediment regime perturbations is critical to predicting and addressing downstream effects to infrastructure, water security and aquatic habitat. A growing body of literature explores the causes, likelihood, size and composition of post-fire debris flows, but the channel response to these disturbances remains poorly studied. This study used repeat field surveys, time-lapse photographs and pre- and post-disturbance remote sensing datasets to document and analyse space- and time-varying channel response to post-fire debris flows along a steep mountain stream in the Upper Colorado River Basin, USA. Specifically, we evaluated channel morphology and bed composition changes, correlations between channel changes and valley and channel attributes, and the relative importance of spring snowmelt versus summer monsoon events. Several cross-sectional channel change types were observed from lidar a month after post-fire debris-flow events, including channelized and braided incision into deposits, incision into the pre-fire channel bed, bank erosion and no change. Channel changes were most correlated with pre-fire channel width, valley width and unit stream power, and these relationships could be tested in other burned locations to evaluate their transferability. Repeat channel surveys before and after snowmelt indicate rapid recovery and channel narrowing following major sediment disturbances, although sediment deposits remained in the channel margins. Together, these results highlight the importance of field and remote sensing-based channel surveys to improve understanding of, and potential to predict, mountain channel response to compounding climate disturbances.

Earth Surface Processes and Landforms↗
USGS2014

Monitoring conterminous United States (CONUS) land cover change with Web-Enabled Landsat Data (WELD)

Forest cover loss and bare ground gain from 2006 to 2010 for the conterminous United States (CONUS) were quantified at a 30 m spatial resolution using Web-Enabled Landsat Data available from the USGS Center for Earth Resources Observation and Science (EROS) (http://landsat.usgs.gov/WELD.php). The approach related multi-temporal WELD metrics and expert-derived training data for forest cover loss and bare ground gain through a decision tree classification algorithm. Forest cover loss was reported at state and ecoregional scales, and the identification of core forests' absent of change was made and verified using LiDAR data from the GLAS (Geoscience Laser Altimetry System) instrument. Bare ground gain correlated with population change for large metropolitan statistical areas (MSAs) outside of desert or semi-desert environments. GoogleEarth™ time-series images were used to validate the products. Mapped forest cover loss totaled 53,084 km2 and was found to be depicted conservatively, with a user's accuracy of 78% and a producer's accuracy of 68%. Excluding errors of adjacency, user's and producer's accuracies rose to 93% and 89%, respectively. Mapped bare ground gain equaled 5974 km2 and nearly matched the estimated area from the reference (GoogleEarth™) classification; however, user's (42%) and producer's (49%) accuracies were much less than those of the forest cover loss product. Excluding errors of adjacency, user's and producer's accuracies rose to 62% and 75%, respectively. Compared to recent 2001–2006 USGS National Land Cover Database validation data for forest loss (82% and 30% for respective user's and producer's accuracies) and urban gain (72% and 18% for respective user's and producer's accuracies), results using a single CONUS-scale model with WELD data are promising and point to the potential for national-scale operational mapping of key land cover transitions. However, validation results highlighted limitations, some of which can be addressed by improving training data, creating a more robust image feature space, adding contemporaneous Landsat 5 data to the inputs, and modifying definition sets to account for differences in temporal and spatial observational scales. The presented land cover extent and change data are available via the official WELD website (ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/ftp://weldftp.cr.usgs.gov/CONUS_5Y_LandCover/).

Remote Sensing of Environment↗
USGS1985

Mineralization of breccia pipes in northern Arizona

The Paleozoic sedimentary rocks on the Colorado Plateau of northern Arizona are host to hundreds ofbreccia pipes. The uranium and copper deposits in these breccia pipes transgress formation boundaries from the Mississippian Redwall Limestone to the Triassic Chinle Formation. They are not classic breccia pipes in that there is no volcanic rock associated with them in time or space. They are the result of solution-collapse within the Redwall Limestone and stoping of the overlying strata. The karst development in the Redwall Limestone began in the Mississippian and apparently either continued to the Triassic or was at least once again active during that time. The mineralization apparently occurred shortly thereafter, sometime during the Mesozoic. Mining activity in breccia pipes of the Grand Canyon region began during the nineteenth century and continues today with the operation of the Hack I, II, and III mines, although the exploited commodity has changed from Cu to U. Although small in size, these pipes contain samples with up to 55 percent U 3 O 8 and can yield ore averaging between 0.30 and 0.60 percent U 3 O 8 .Mineralization at the surface commonly occurs within nodules and concretions associated with pyrite and goethite and along fractures, while the primary ore of the unoxidized zones is commonly within a comminuted sandstone matrix surrounding breccia fragments of overlying formations. The ore mineral is uraninite, although associated with it are sphalerite, galena, chalcopyrite, tennantite, millerite, siegenite, and/molybdenite. Some of the surface nodules are encrusted with malachite and are exceptionally enriched in Ag. Pyrite is abundant, and the organic carbon content of some rocks is high enough to suggest that it, along with the pyrite, may be a reductant for uranium. In contrast, it is possible, if uranium were transported as a bicarbonate or carbonate complex, that only a conduit of brecciated rock was necessary to release CO 2 , thus disrupting the equilibrium and allowing uraninite to precipitate. An extensive suite of elements is significantly enriched in the mineralized rock: Ag, As, Ba, Cd, Co, Cr, Cs, Cu, Hg, Mo, Ni, Pb, Sb, Se, Sr, U, V, Zn, and the rare earth elements. Of these, Cu, Pb, Zn, Ag, and particularly As appear to be the best geochemical indicators of mineralized pipes. At present the origin of the mineralizing fluids is not known. The lack of extensive silification within the breccia, along with the 80 degrees to 173 degrees C fluid inclusion-filling temperatures on sphalerite, dolomite, quartz, and calcite, suggests relatively low-temperature mineralizing fluids, although heated in excess of what would be expected from the normal geothermal gradient on the Colorado Plateau. With the exception of the U-mineralized rock, the mineral assemblage and geochemistry is similar to Mississippi Valley-type deposits.

Economic Geology↗
USGS2012

Stratigraphic, lithologic, and major element geochemical constraints on magmatic evolution at Lassen volcanic center, California

The Lassen volcanic center is the most recent of several long-lived volcanic centers in the southernmost Cascade Range. These centers have erupted products ranging from basaltic andesite to rhyolite and are superimposed on a background of regional basaltic to andesitic volcanism. The evolution of the Lassen volcanic center is described in three stages. Stages I and II comprise the Brokeoff volcano, and 80 km 3 andesitic stratocone, active from 600 to 400 ka. Brokeoff volcano is compositionally equivalent to the regional basaltic andesite to andesite volcanism in the Lassen region and is the result of structurally controlled focusing of the diffuse regional mafic magmatism. Stage III comprises a silicic dome field and adjacent area of hybrid andesites and has a total volume of about 100 km 3 . Volcanism during stage III was episodic and is subdivided into four sequences of lithologically and temporally distinct lavas. Stage III began at 400 ka with a rhyolitic, caldera-forming pyroclastic eruption and chemically related lavas. Additional sequences of dacite erupted between 250–200 ka and 100–0 ka. Hybrid andesites erupted adjacent to the silicic dome field between 300 and 0 ka. Porphyritic andesite and dacite with high Al 2 O 3 , low TiO 2 , medium K 2 O and FeO/MgO ratios of 1.5–2.0 are the most abundant rock types in the Lassen volcanic center. However, the single most voluminous unit is sparsely phyric rhyolite pumice. In general, the lavas of Lassen volcanic center form a single coherent trend on major element variation diagrams and in pseudo-quaternary phase space, consistant with an origin either by fractional crystallization or magma mixing. In detail, however, the lack of systematic temporal change in silica and subtly crossing trends indicate a complex origin. A variety of statistically successful fractional crystallization models can be constructed that derive Brokeoff andesites from regional magmas. An important conclusion of the modeling is that if fractional crystallization is the process responsible for generation of Brokeoff andesite, then the parent magma must be low to medium K in geochemical affinity in order to explain the variation in K 2 O. However, although major element variation can be modeled by fractional crystallization, petrographic and stratigraphic evidence indicates that magma mixing is an important but subtle process in Brokeoff lavas and suggests that lavas evolved in small independent batches. Lavas erupted during stage III, while predominantly silicic, range from 53 to 75% SiO 2 . Disequilibrium mineral assemblages in the stage III lavas indicate that they are not directly derived from Brokeoff andesite by fractional crystallization. Mixing of silicic magma with regional mafic magma and disaggregation of andesitic quenched magmatic inclusions play dominant roles in the compositional diversity of stage III lavas.

Journal of Geophysical Research Solid Earth↗
USGS2019

Agricultural cropland extent and areas of South Asia derived using Landsat satellite 30-m time-series big-data using random forest machine learning algorithms on the Google Earth Engine cloud

The South Asia (India, Pakistan, Bangladesh, Nepal, Sri Lanka and Bhutan) has a staggering 900 million people (~43% of the population) who face food insecurity or severe food insecurity as per United Nations, Food and Agriculture Organization’s (FAO) the Food Insecurity Experience Scale (FIES). The existing coarse-resolution (>250-m) cropland maps lack precision in geo-location of individual farms and have low map accuracies. This also results in uncertainties in cropland areas calculated from such products. Thereby, the overarching goal of this study was to develop high spatial resolution (30-m or better) baseline cropland extent product of South Asia for the year 2015 using Landsat satellite time-series big-data and machine learning algorithms (MLAs) on the Google Earth Engine (GEE) cloud computing platform. To eliminate the impact of clouds, ten time-composited Landsat bands (blue, green, red, NIR, SWIR1, SWIR2, Thermal, EVI, NDVI, NDWI) were derived for each of the 3 time-periods over 12 months (monsoon: Julian days 151-300; winter: Julian days 301-365 plus 1-60; and summer: Julian days 61-150), taking the every 8-day data from Landsat-8 and 7 for the years 2013-2015, for a total of 30-bands plus global digital elevation model (GDEM) derived slope band. This 31-band mega-file big data-cube was composed for each of the 5 agro-ecological zones (AEZ’s) of South Asia and formed a baseline data for image classification and analysis. Knowledge-base for the Random Forest (RF) MLAs were developed using spatially well spread-out reference training data (N=2179) in 5 AEZs. Classification was performed on GEE for each of the 5 AEZs using well-established knowledge-based and RF MLAs on the cloud. Map accuracies were measured using independent validation data (N=1185). The survey showed that the South Asia cropland product had a producer’s accuracy of 89.9% (errors of omissions of 10.1%), user’s accuracy of 95.3% (errors of commission of 4.7%) and an overall accuracy of 88.7%. The National and sub-national (districts) areas computed from this cropland extent product explained 80-96% variability when compared with the National statistics of the South Asian Countries. The full resolution imagery can be viewed at full-resolution, by zooming-in to any location in South Asia or the world, at www.croplands.org and the cropland products of South Asia downloaded from The Land Processes Distributed Active Archive Center (LP DAAC) of National Aeronautics and Space Administration (NASA) and the United States Geological Survey (USGS): https://lpdaac.usgs.gov/products/gfsad30saafgircev001/

GIScience and Remote Sensing↗
USGS2012

EAARL coastal topography--Alligator Point, Louisiana, 2010

This project provides highly detailed and accurate datasets of a portion of Alligator Point, Louisiana, acquired on March 5 and 6, 2010. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the National Aeronautics and Space Administration (NASA) Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multispectral color-infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine aircraft, but the instrument was deployed on a Pilatus PC-6. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for presurvey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the "bare earth" under vegetation from a point cloud of last return elevations.

Louisiana↗
USGS1949

Populations and home range relationships of the box turtle, Terrapene carolina (Linnaeus)

A population study of Terrapene carolina (Linnaeus) was made at the Patuxent Research Refuge, Maryland, from 1944 to 1947. A thirty acre area in bottomland forest was selected for intensive study. Turtles were marked by filing notches in marginal scutes according to a code. Turtles make extensive use of brushy shelter during the day as well.as at night. Gully banks and woods openings are used for sunning. Nights are usually spent in a 'form,' constructed by the turtle in leaves, debris, or earth. A form may be used once or it may be used repeatedly by the same or different turtles. Weather conditions most favorable to turtle activity are high humidity, warm sunny days, and frequent rains. Periods of activity are alternated with periods of quiet, even in favorable weather. There is no evidence for territorialism. Ranges of turtles of all ages and both sexes overlap grossly. Turtles are frequently found near each other but no antagonistic behavior has been observed. Adult turtles occupy specific home ranges which they maintain from year to year. Turtles retained their ranges even though a flood that completely covered the study area. Maximum home range diameters were determined by measurements of the mapped ranges of individual turtles. There was no significant difference between sizes of male and female ranges: males 33O+ 26 feet, females 37O+29 feet. A trail-laying device was used in following travel routes for 456 turtle days. Normal movements within the home range are characterized by (1) turns, doublings, detours, and criss-crossing paths, (2) interspersion of fairly direct traverses of the home range, (3) frequently repeated travels over certain routes. Maximum limits of the home range are ordinarily reached within a few days or weeks, although some turtles cover only one portion of the range at a time. Some turtles have two home ranges. One of these turtles was followed with a trailer for 161 days in 1946 and 1947. Trips outside the home range are made by females for egg laying. Trips of unexplained nature are made by both males and females. Turtles from other areas are occasional transients through the study plot. Standardized census trips provided data for estimating the size of the population. Census data were taken after females had returned from egg laying. The samples were spaced at intervals of a week or more to allow free movement of turtles over their ranges and assure more nearly equal availability of all turtles. These standard samples were compared by a collection ratio. Assumptions involved in the use of this ratio are discussed. Correction factors were applied to make allowance for turtles whose ranges.were partly inside and partly outside the study area, and for transient turtles. A second estimate, on the basis of the entire season's collecting, gave closely comparable results. The population of the study area is estimated to be between four and five turtles per acre, with juveniles constituting less than ten per cent of the total.

Dissertation Abstracts International↗
USGS2013

EAARL coastal topography and imagery–Western Louisiana, post-Hurricane Rita, 2005: First surface

These remotely sensed, geographically referenced color-infrared (CIR) imagery and elevation measurements of lidar-derived first-surface (FS) topography datasets were produced by the U.S. Geological Survey (USGS), St. Petersburg Coastal and Marine Science Center, St. Petersburg, Florida, and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, Virginia. This project provides highly detailed and accurate datasets of a portion of the Louisiana coastline beachface, acquired post-Hurricane Rita on September 27-28 and October 2, 2005. The datasets are made available for use as a management tool to research scientists and natural-resource managers. An innovative airborne lidar instrument originally developed at the National Aeronautics and Space Administration (NASA) Wallops Flight Facility, and known as the Experimental Advanced Airborne Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532-nanometer) lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive lidar, a down-looking red-green-blue (RGB) digital camera, a high-resolution multispectral color-infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system, and the resulting data were then processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of lidar data in an interactive or batch mode. Modules for presurvey flight-line definition, flight-path plotting, lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is used routinely to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the "bare earth" under vegetation from a point cloud of last return elevations. For more information about similar projects, please visit the Lidar for Science and Resource Management Website.

Louisiana↗
USGS2014

Characterizing riverbed sediment using high-frequency acoustics 1: spectral properties of scattering

Bed-sediment classification using high-frequency hydro-acoustic instruments is challenging when sediments are spatially heterogeneous, which is often the case in rivers. The use of acoustic backscatter to classify sediments is an attractive alternative to analysis of topography because it is potentially sensitive to grain-scale roughness. Here, a new method is presented which uses high-frequency acoustic backscatter from multibeam sonar to classify heterogeneous riverbed sediments by type (sand, gravel,rock) continuously in space and at small spatial resolution. In this, the first of a pair of papers that examine the scattering signatures from a heterogeneous riverbed, methods are presented to construct spatially explicit maps of spectral properties from geo-referenced point clouds of geometrically and radiometrically corrected echoes. Backscatter power spectra are computed to produce scale and amplitude metrics that collectively characterize the length scales of stochastic measures of riverbed scattering, termed &lsquo;stochastic geometries&rsquo;. Backscatter aggregated over small spatial scales have spectra that obey a power-law. This apparently self-affine behavior could instead arise from morphological- and grain-scale roughnesses over multiple overlapping scales, or riverbed scattering being transitional between Rayleigh and geometric regimes. Relationships exist between stochastic geometries of backscatter and areas of rough and smooth sediments. However, no one parameter can uniquely characterize a particular substrate, nor definitively separate the relative contributions of roughness and acoustic impedance (hardness). Combinations of spectral quantities do, however, have the potential to delineate riverbed sediment patchiness, in a data-driven approach comparing backscatter with bed-sediment observations (which is the subject of part two of this manuscript).

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Geoelectric hazard maps for the Pacific Northwest

Maps of extreme value, horizontal component geoelectric field amplitude are constructed for the Pacific Northwest United States (and parts of neighboring Canada). Multidecade long geoelectric field time series are calculated by convolving Earth surface impedance tensors from 71 discrete magnetotelluric survey sites across the region with historical 1‐min (2‐min Nyquist) geomagnetic variation time series obtained from two nearby observatories. After fitting statistical models to 1‐min geoelectric amplitudes realized during magnetic storms, extrapolations are made to estimate threshold amplitudes that are only exceeded, on average, once per century. One hundred‐year geoelectric exceedance amplitudes range from 0.06 V/km at a survey site in western Washington State to 9.47 V/km at a site in southeast British Columbia; 100‐year geoelectric exceedance amplitudes equal 7.10 V/km at a site north of Seattle and 2.28 V/km at a site north of Portland. Systematic and random errors are estimated to be less than 20%, much less than site‐to‐site differences in geoelectric amplitude that arise from site‐to‐site differences in surface impedance. Maps of 100‐year exceedance amplitudes are compared with the peak geoelectric amplitudes realized during the March 1989 magnetic superstorm; it is noted that some storms of relatively modest intensity can generate localized geoelectric fields of relatively high amplitude. The geography of geoelectric hazard across the Pacific Northwest is closely related to known geologic and tectonic structures.

Pacific Northwest↗