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 775 records · Page 43Linked to original sources
USGS2010

Inelastic off-fault response and three-dimensional dynamics of earthquake rupture on a strike-slip fault

Large dynamic stress off the fault incurs an inelastic response and energy loss, which contributes to the fracture energy, limiting the rupture and slip velocity. Using an explicit finite element method, we model three-dimensional dynamic ruptures on a vertical strike-slip fault in a homogeneous half-space. The material is subjected to a pressure-dependent Drucker-Prager yield criterion. Initial stresses in the medium increase linearly with depth. Our simulations show that the inelastic response is confined narrowly to the fault at depth. There the inelastic strain is induced by large dynamic stresses associated with the rupture front that overcome the effect of the high confining pressure. The inelastic zone increases in size as it nears the surface. For material with low cohesion (~5 MPa) the inelastic zone broadens dramatically near the surface, forming a "flowerlike" structure. The near-surface inelastic strain occurs in both the extensional and the compressional regimes of the fault, induced by seismic waves ahead of the rupture front under a low confining pressure. When cohesion is large (~10 MPa), the inelastic strain is significantly reduced near the surface and confined mostly to depth. Cohesion, however, affects the inelastic zone at depth less significantly. The induced shear microcracks show diverse orientations near the surface, owing to the low confining pressure, but exhibit mostly horizontal slip at depth. The inferred rupture-induced anisotropy at depth has the fast wave direction along the direction of the maximum compressive stress.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Principal component analysis of geodetically measured deformation in Long Valley caldera, eastern California, 1983-1987

Typical geodetic measurements of deformation consist of repeated surveys of a particular geodetic network. Such deformation data can be interpreted as a consequence of one or more self-coherent sources by means of principal component analysis. A self-coherent source is defined as any source that produces deformation that is time and space separable. Principal component analysis then gives the time and space factors that characterize the deformation attributed to each self-coherent source. Geodetic measurements of deformation at Long Valley caldera provide two examples of the application of principal component analysis. A 40-line trilateration network surrounding the caldera was surveyed in midsummer 1983, 1984, 1985, 1986, and 1987. Principal component analysis indicates that the observed deformation can be represented by a single coherent source. The time dependence for that source displays a rapid rate of deformation in 1983–1984 followed by less rapid but uniform rate in the 1984–1987 interval. The spatial factor seems consistent with expansion of a magma chamber beneath the caldera plus some shallow right-lateral slip on a vertical fault in the south moat of the caldera. An independent principal component analysis of the 1982, 1983, 1984, 1985, 1986, and 1987 leveling across the caldera requires two self-coherent sources to explain the deformation. The deformation pattern produced by the larger of these two sources appears to be roughly consistent with that found from the trilateration data. The deformation due to the second source is a nearly uniform tilt in the uplift profile. Presumably, that tilt is simply an artifact of systematic error in the leveling.

Journal of Geophysical Research Solid Earth↗
USGS2016

Subsidence rates at the southern Salton Sea consistent with reservoir depletion

Space geodetic measurements from the Envisat satellite between 2003 and 2010 show that subsidence rates near the southeastern shoreline of the Salton Sea in Southern California are up to 52mmyr −1 greater than the far-field background rate. By comparing these measurements with model predictions, we find that this subsidence appears to be dominated by poroelastic contraction associated with ongoing geothermal fluid production, rather than the purely fault-related subsidence proposed previously. Using a simple point source model, we suggest that the source of this proposed volumetric strain is at depths between 1.0 km and 2.4 km (95% confidence interval), comparable to generalized boundaries of the Salton Sea geothermal reservoir. We find that fault slip on two previously imaged tectonic structures, which are part of a larger system of faults in the Brawley Seismic Zone, is not an adequate predictor of surface velocity fields because the magnitudes of the best fitting slip rates are often greater than the full plate boundary rate and at least 2 times greater than characteristic sedimentation rates in this region. Large-scale residual velocity anomalies indicate that spatial patterns predicted by fault slip are incompatible with the observations.

California↗
USGS2020

Peak ground velocity spatial variability revealed by dense seismic array in southern California

Understanding and modeling variability of ground motion is essential for building accurate and precise ground motion prediction equations, which can net site‐specific characterization and reduced hazard levels. Here, we explore the spatial variability in peak ground velocity (PGV) at Sage Brush Flats along the San Jacinto Fault in southern California. We use data from a dense array (0.6 x 0.6 km 2 , 1,108 geophones, station spacings 10‐30 m) deployed in 2014 for ~1‐month. These data offer an opportunity to study small scale variability in this region. We examine 38 2≤M L ≤4.2 earthquakes within 200 km of the array. Fault strands and a small basin impact the ground motions, producing PGV variations up to 22% of the mean and a 40% reduction in P and S wave near‐surface velocities. We find along‐fault rupture directivity, source, and path effects can increase PGVs by 167%. Surface PGV measurements exceed the co‐located borehole station (depth 148 m) PGV by factors of 3‐10, confirming the impact on PGV from near surface fault structures, basins, topography, and amplifications from soft sediments. Consistently we find high PGVs within the basin structure. A pair of ~co‐located M L 2.6 events produce repeatable PGV values with similar spatial patterns. The average corner frequencies of these two events are 11‐16 Hz and viable measurements of stress drop can differ by 6.45MPa. Within this small array, the PGV values are variable implying spatial extrapolation of PGV to regions of known faults and basins, even across a small area, should be done with caution.

California↗
USGS1982

Post-depositional alteration of titanomagnetite in a Miocene sandstone, south Texas (U.S.A.)

Petrographic and geochemical studies have yielded information on the time-space relationships of the post-depositional alteration of detrital titanomagnetite (Ti-mt) in fine- to medium-grained sandstone from unoriented core samples (taken below the water table at depths of 30–45 m) of the Miocene Catahoula Sandstone, south Texas. Aqueous sulfide introduced from sour gas reservoirs along a growth fault into part of the Catahoula shortly after deposition resulted in the replacement at the periphery of Ti-mt grains by iron disulfide (FeS 2 ) minerals. Remnants of Ti-mt in cores of the partly sulfidized grains show no evidence of earlier hematitic oxidation. After sulfidization, part of the sandstone body was invaded by oxygenated groundwaters flowing down a shallowly inclined (1°) hydrologic gradient. The boundary between oxidized and reduced facies is clearly defined by the distribution of ferric and ferrous iron minerals, and the concentrations of Mo, U, and Se. In oxidized (light-red) strata that had not been previously subjected to sulfidic-reducing conditions but that are correlative with strata containing FeS 2 minerals, Ti-mt has been partly to entirely replaced pseudomorphously by hematite to form martite. The absence of hematitic alteration of Ti-mt in the reduced facies is strong evidence that martite in the oxidized facies formed after deposition.

Texas↗
USGS2017

Using multi-date satellite imagery to monitor invasive grass species distribution in post-wildfire landscapes: An iterative, adaptable approach that employs open-source data and software

Among the most pressing concerns of land managers in post-wildfire landscapes are the establishment and spread of invasive species. Land managers need accurate maps of invasive species cover for targeted management post-disturbance that are easily transferable across space and time. In this study, we sought to develop an iterative, replicable methodology based on limited invasive species occurrence data, freely available remotely sensed data, and open source software to predict the distribution of Bromus tectorum (cheatgrass) in a post-wildfire landscape. We developed four species distribution models using eight spectral indices derived from five months of Landsat 8 Operational Land Imager (OLI) data in 2014. These months corresponded to both cheatgrass growing period and time of field data collection in the study area. The four models were improved using an iterative approach in which a threshold for cover was established, and all models had high sensitivity values when tested on an independent dataset. We also quantified the area at highest risk for invasion in future seasons given 2014 distribution, topographic covariates, and seed dispersal limitations. These models demonstrate the effectiveness of using derived multi-date spectral indices as proxies for species occurrence on the landscape, the importance of selecting thresholds for invasive species cover to evaluate ecological risk in species distribution models, and the applicability of Landsat 8 OLI and the Software for Assisted Habitat Modeling for targeted invasive species management.

International Journal of Applied Earth Observation↗
USGS2023

Tracing magmatic genesis and evolution through single zircon crystals from successive supereruptions from the Socorro Caldera Complex, USA

Large volume rhyolitic ignimbrite volcanism is a significant contributor to the evolving crust. The introduction of high-silica material into the upper crust, differentiation within the middle crust, and partial melting in the lower crust contributes to geochemical and isotopic evolution of the crust. Developing accurate models for the genetic evolution of these events is dependent upon geochronology to determine rates of magmatic processes as model constraints. We present new zircon high-precision CA-ID-TIMS U-Pb geochronology and MC-ICPMS Hf isotope geochemistry for four ignimbrites from the nested caldera complex near Socorro, New Mexico (USA), within the Mogollon-Datil volcanic field. In agreement with past 40 Ar- 39 Ar data, interpretations of new U-Pb data indicate eruptions from the Socorro caldera cluster were pulsed. These pulses were intermittently spaced, and a volcanic hiatus following the Hells Mesa Tuff at 33.442 ± 0.015 Ma was interrupted by four successive eruptions, beginning with the La Jencia Tuff at 29.158 ± 0.025 Ma and finishing with the South Canyon Tuff at 28.066 ± 0.021 Ma. Zircon age spectra became more protracted with each eruption, exhibiting age dispersions ranging from 0.347 Myr in the Hells Mesa Tuff to 4.502 Myr in the South Canyon Tuff. The increased dispersion is paralleled by an increase in the proportion of normally discordant grains, indicative of xenocryst incorporation. These protracted age spectra are not necessarily a function of thermal maturation in the middle to upper crust due to long-lived magma chambers. Rather, they are likely the result of increased melting of zircon-bearing lower crust due to deep thermal maturation from repeated juvenile magma injections based on the incorporation of zircon material at the melt source. In contrast, the Hf isotope record is volumetrically dominated by autocrystic zircon domains and becomes more radiogenic through time, recording juvenile replenishment of the lower crust during progressive melting. Together, these data record the protracted evolution of the lower crust sampled by ignimbrites, lend insight into that evolution, and emphasize the need for detailed interpretation of high-precision datasets to advance volcanic models.

Arizona, New Mexico↗
USGS2010

Coherence of Mach fronts during heterogeneous supershear earthquake rupture propagation: Simulations and comparison with observations

We study how heterogeneous rupture propagation affects the coherence of shear and Rayleigh Mach wavefronts radiated by supershear earthquakes. We address this question using numerical simulations of ruptures on a planar, vertical strike-slip fault embedded in a three-dimensional, homogeneous, linear elastic half-space. Ruptures propagate spontaneously in accordance with a linear slip-weakening friction law through both homogeneous and heterogeneous initial shear stress fields. In the 3-D homogeneous case, rupture fronts are curved owing to interactions with the free surface and the finite fault width; however, this curvature does not greatly diminish the coherence of Mach fronts relative to cases in which the rupture front is constrained to be straight, as studied by Dunham and Bhat (2008a). Introducing heterogeneity in the initial shear stress distribution causes ruptures to propagate at speeds that locally fluctuate above and below the shear wave speed. Calculations of the Fourier amplitude spectra (FAS) of ground velocity time histories corroborate the kinematic results of Bizzarri and Spudich (2008a): (1) The ground motion of a supershear rupture is richer in high frequency with respect to a subshear one. (2) When a Mach pulse is present, its high frequency content overwhelms that arising from stress heterogeneity. Present numerical experiments indicate that a Mach pulse causes approximately an ω −1.7 high frequency falloff in the FAS of ground displacement. Moreover, within the context of the employed representation of heterogeneities and over the range of parameter space that is accessible with current computational resources, our simulations suggest that while heterogeneities reduce peak ground velocity and diminish the coherence of the Mach fronts, ground motion at stations experiencing Mach pulses should be richer in high frequencies compared to stations without Mach pulses. In contrast to the foregoing theoretical results, we find no average elevation of 5%-damped absolute response spectral accelerations (SA) in the period band 0.05–0.4 s observed at stations that presumably experienced Mach pulses during the 1979 Imperial Valley, 1999 Kocaeli, and 2002 Denali Fault earthquakes compared to SA observed at non-Mach pulse stations in the same earthquakes. A 20% amplification of short period SA is seen only at a few of the Imperial Valley stations closest to the fault. This lack of elevated SA suggests that either Mach pulses in real earthquakes are even more incoherent that in our simulations or that Mach pulses are vulnerable to attenuation through nonlinear soil response. In any case, this result might imply that current engineering models of high frequency earthquake ground motions do not need to be modified by more than 20% close to the fault to account for Mach pulses, provided that the existing data are adequately representative of ground motions from supershear earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Decaying lava extrusion rate at El Reventador Volcano, Ecuador measured using high-resolution satellite radar

Lava extrusion at erupting volcanoes causes rapid changes in topography and morphology on the order of tens or even hundreds of meters. Satellite radar provides a method for measuring changes in topographic height over a given time period to an accuracy of meters, either by measuring the width of radar shadow cast by steep sided features, or by measuring the difference in radar phase between two sensors separated in space. We measure height changes, and hence estimate extruded lava volume flux, at El Reventador, Ecuador, between 2011 and 2016, using data from the RADARSAT-2 and TanDEM-X satellite missions. We find that 39 new lava flows were extruded between 9 February 2012 and 24 August 2016, with a cumulative volume of 44.8M m 3 dense rock equivalent and a gradually decreasing eruption rate. The average dense rock rate of lava extrusion during this time is 0.31 ± 0.02 m 3 s −1 , which is similar to the long-term average from 1972 to 2016. Apart from a volumetrically small dyke opening event between 9 March and 10 June 2012, lava extrusion at El Reventador is not accompanied by any significant magmatic ground deformation. We use a simple physics-based model to estimate that the volume of the magma reservoir under El Reventador is greater than 3 km 3 . Our lava extrusion data can be equally well fit by models representing a closed reservoir depressurising during the eruption with no magma recharge, or an open reservoir with a time-constant magma recharge rate of up to 0.35 ± 0.01 m 3 s −1 .

El Reventador Volcano↗
USGS2024

A two-dimensional, reach-scale implementation of space-time image velocimetry (STIV) and comparison to particle image velocimetry (PIV)

Image-based algorithms have become a powerful tool for estimating flow velocities in rivers. In this study, we generalize the space-time image velocimetry (STIV) framework for reach-scale application rather than along a cross section. The new algorithm provides information on both the magnitude and orientation of velocity vectors, and we refer to the algorithm as two-dimensional STIV, or 2D-STIV. The workflow involves setting up a grid, using centreline tangent vectors as initial estimates of flow direction, and then extracting space-time images (STIs) along search lines radiating from each grid node. The autocorrelation function is used to infer the inclination of streak lines present in STIs, which represents the advection of water surface features. Information on flow direction is obtained by evaluating various candidate search lines and identifying that which yields the highest velocity. This search can be performed exhaustively or via optimization. We applied the new 2D-STIV algorithm to three test cases, one simulated data set and two natural channels, and compared image-derived velocities to modelled or measured values. We also applied two established particle image velocimetry (PIV) algorithms to the same data sets. 2D-STIV performed as well as the two PIV algorithms for simulated images. For a natural river with distinct water surface features, 2D-STIV was effective for much of the channel but also led to a more patchy, irregular velocity field than the two PIV algorithms. For a site lacking obvious surface features, exhaustive 2D-STIV led to velocity estimates uncorrelated with field data while the optimization-based version produced erratic flow directions. 2D-STIV also required greater image sequence durations, higher frame rates, and generally longer computational run times. Overall, ensemble PIV was the most reliable algorithm.

Earth Surface Processes and Landforms↗
USGS1989

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

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

Journal of African Earth Sciences↗
USGS1968

Water resources inventory of Connecticut Part 3: lower Thames and southeastern coastal river basins

The lower Thames and southeastern coastal river basins have a relatively abundant supply of water of generally good quality which is derived from streams entering the area and precipitation that has fallen on the area. Annual precipitation has ranged from about 32 inches to 65 inches and has averaged about 48 inches over a 30-year period. Approximately 22 inches of water are returned to the atmosphere each year by evaporation and transpiration; the remainder of the annual precipitation either flows overland to streams or percolates downward to the water table and ultimately flows out of the report area through estuaries and coastal streams or as underflow through the deposits beneath. During the autumn and winter months precipitation normally is sufficient to cause a substantial increase in the amount of water stored underground and in surface reservoirs within the report area, whereas in the summer most of the precipitation is lost through evaporation and transpiration, resulting in sharply reduced stream-flow and lowered ground-water levels. The mean monthly storage of water on an average is about 3.8 inches higher in November than it is in June. The amount of water that flows through and out of the report area represents the total amount of water potentially available for use by man. For the 30-year period 1931 through 1960, the annual runoff from the report area has averaged nearly 26 inches (200 billion gallons), from the entire Thames River basin above Norwich about 24 inches (530 billion gallons), and from the Pawcatuck River basin about 26 inches (130 billion gallons). A total average annual runoff of 860 billion gallons is therefore available. Although runoff indicates the total amount of water potentially available, it is usually not economically feasible for man to use all of it. On the other hand, with increased development, it is possible that some water will be reused several times. The water available may be tapped as it flows through the area or is temporarily stored in streams, lakes, and aquifers. The amounts that can be developed vary from place to place and time to time, depending on the amount of precipitation, on the size of drainage area, on the thickness, permeability, and areal extent of aquifers, and on the variations in chemical and physical quality of the water. Differences in streamflow from point to point are due primarily to differences in the proportion of stratified drift in the drainage basin above each point, which affect the timing of streamflow, and to differences in precipitation, which affect the amount of streamflow. Ground water can be obtained from wells almost anywhere in the area, but the amount obtainable at any particular point depends upon the type and water-bearing properties of the aquifers. For practical purposes, the earth materials in the report area comprise three aquifers--stratified drift, bedrock, and till. Stratified drift is the only aquifer generally capable of yielding more than 100 gpm (gallons per minute) to individual wells. It covers about 20 percent of the area and occurs chiefly in lowlands where it overlies till and bedrock. The coefficient of permeability of the coarse-grained unit of stratified drift averages about 1,500 gbd (gallons per day) per sq ft. Drilled, screened wells tapping this unit are known to yield from 4 to 88o gpm and average 146 gpm. Dug wells in coarse-grained stratified drift supply about 2 gpm per foot of drawdown over a period of a few hours. Fine-grained stratified drift has an average coefficient of permeability of about 300 gpd per sq ft and can usually yield supplies sufficient for household use to dug wells. Bedrock and till are widespread in extent but generally provide only small water supplies. Bedrock is tapped chiefly by drilled wells, about 90 percent of which will supply at least 3 gpm. Very few, however, will supply more than 50 gpm. Till is tapped in a few places by dug wells which can yield small supplies of only a few hundred gpd throughout all or most of the year. The coefficient of permeability of till ranges from about 0.2 gpd per sq ft to 120 gpd per sq ft. The amount of ground water potentially available in the report area depends upon the amount of ground-water outflow, the amount of ground water in storage, and the quantity of water available by induced infiltration from streams and lakes. From data on permeability, saturated thickness, recharge, yield from aquifer storage, well performance, and streamflow, preliminary estimates of ground-water availability can be made for any point in the report area. Long-term yields estimated for 18 areas of stratified drift especially favorable for development of large ground-water supplies ranged from 1.3 to 66 mgd. Detailed site studies to determine optimum yields, drawdowns, and spacing of individual wells are needed before major ground-water development is undertaken in these or other areas. The chemical quality of water in the report area is generally good to excellent. Samples of naturally occurring surface water collected at 24 sites contained less than 151 ppm (parts per million) of dissolved solids and less than 63 ppm of hardness. Water from wells is more highly mineralized than naturally occurring water from streams. Even so only 12 percent of the wells sampled yielded water with more than 200 ppm of dissolved solids and only 8 percent yielded water with more than 120 ppm of hardness. Even in major streams, which are used to transport industrial waste, hardness rarely exceeds 60 ppm and the dissolved mineral content is generally less than 200 ppm. At a few places in the town of Montville however, waters may contain dissolved mineral concentrations of 2,000 to 4,000 ppm. Iron and manganese in both ground water and surface water are the only constituents whose concentrations commonly exceed recommended limits for domestic and industrial use. Most wells in the report area yield clear water with little or no iron or manganese, but distributed among them are wells yielding ground water that contains enough of these dissolved constituents to be troublesome for most uses. Iron concentrations in naturally occurring stream water exceed 0.3 ppm under low-flow conditions at 33 percent of the sites sampled. Large concentrations of iron in stream water result from discharge of iron-bearing water from aquifers or from swamps where it is released largely from decaying vegetation. Ground water more than 30 feet below the land surface has a relatively constant temperature, usually between 48°F and 52°F. Water temperature in very shallow wells may fluctuate from about 38°F in February or March to about 55°F in late summer. Water temperature in the larger streams fluctuates much more widely, ranging from 32°F at least for brief periods in winter, to about 85°F occasionally during summer. The quality of suspended sediment transported by streams in the area is negligible. Turbidity in streams is generally not a problem although amounts large enough to be troublesome may occur locally at times. The total amount of water used in the report area for all purposes during 1964 was about 118,260 million gallons, of which 105,600 million gallons was estuarine water used for cooling by industry. The average per capita water use, excluding estuarine, temporary summer residence, and institutional water was equivalent to 186 gpd. Public water systems supplied the domestic needs of nearly tw0-thirds the population of the report area. All of the 19 systems, which were sampled, provided water of better quality than the U.S. Public Health Service suggests for drinking water standards.

Connecticut↗
USGS2006

Dislocation pileup as a representation of strain accumulation on a strike-slip fault

The conventional model of strain accumulation on a vertical transform fault is a discrete screw dislocation in an elastic half-space with the Burgers vector of the dislocation increasing at the rate of relative plate motion. It would be more realistic to replace that discrete dislocation by a dislocation distribution, presumably a pileup in which the individual dislocations are in equilibrium. The length of the pileup depends upon the applied stress and the amount of slip that has occurred at depth. I argue here that the dislocation pileup (the transition on the fault from no slip to slip at the full plate rate) occupies a substantial portion of the lithosphere thickness. A discrete dislocation at an adjustable depth can reproduce the surface deformation profile predicted by a pileup so closely that it will be difficult to distinguish between the two models. The locking depth (dislocation depth) of that discrete dislocation approximation is substantially (???30%) larger than that (depth to top of the pileup) in the pileup model. Thus, in inverting surface deformation data using the discrete dislocation model, the locking depth in the model should not be interpreted as the true locking depth. Although dislocation pileup models should provide a good explanation of the surface deformation near the fault trace, that explanation may not be adequate at greater distances from the fault trace because approximating the expected horizontally distributed deformation at subcrustal depths by uniform slip concentrated on the fault is not justified.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

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

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

Alaska↗
USGS2018

Interseismic ground deformation and fault slip rates in the greater San Francisco Bay Area from two decades of space geodetic data

The detailed spatial variations of strain accumulation and creep on major faults in the northern San Francisco Bay Area (North Bay), which are important for seismic potential and evaluation of natural hazards, remain poorly understood. Here we combine interferometric synthetic aperture radar data from the ERS‐1/2 and Envisat satellites between 1992 and 2010 with continuous and campaign GPS data to obtain a high spatial and temporal coverage of ground deformation of the North Bay. The SAR data from both ascending and descending orbits are combined to separate horizontal and vertical components of the deformation. We jointly invert the horizontal component of the mean velocities derived from these data to infer the deep strike‐slip rates on major locked faults. We use the estimated deep rates to simulate the long‐wavelength deformation due to interseismic elastic strain accumulation along these locked faults. After removing the long‐wavelength signal from the InSAR horizontal mean velocity field, we estimate fault‐parallel surface creep rates of up to 2 mm/year along the central section of the Rodgers Creek fault and surface creep rates ranging between 2 and 4 mm/year along the Concord fault. No surface creep is geodetically resolved along the West Napa and Green Valley fault zones. We identified characteristically repeating earthquakes on the Rodgers Creek fault, the West Napa fault, the Green Valley fault, and the Concord fault. Nontectonic deformation in the Geysers geothermal field and in Late Cenozoic basins (Rohnert Park and Sonoma basins) are also observed, likely due to hydrological and sediment‐compaction processes, respectively.

Journal of Geophysical Research B: Solid Earth↗
USGS2011

EAARL Coastal Topography--Cape Canaveral, Florida, 2009: First Surface

These remotely sensed, geographically referenced elevation measurements of lidar-derived first-surface (FS) topography datasets were produced collaboratively by the U.S. Geological Survey (USGS), St. Petersburg Coastal and Marine Science Center, St. Petersburg, FL, and the National Aeronautics and Space Administration (NASA), Kennedy Space Center, FL. This project provides highly detailed and accurate datasets of a portion of the eastern Florida coastline beachface, acquired on May 28, 2009. 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 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.

Data Series↗
USGS2009

EAARL Topography-Vicksburg National Military Park 2007: First Surface

These remotely sensed, geographically referenced elevation measurements of Lidar-derived first-surface (FS) topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL; the National Park Service (NPS), Gulf Coast Network, Lafayette, LA; and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of the Vicksburg National Military Park in Mississippi, acquired on September 12, 2007. 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 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 multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers, and an integrated miniature digital inertial measurement unit, which provide for submeter 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 first surface topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations.

Data Series↗
USGS1967

Water resources inventory of Connecticut Part 2: Shetucket River Basin

The Shetucket River basin has a relatively abundant supply of water of generally good quality which is derived from precipitation that has fallen on the basin. Annual precipitation has ranged from about 30 inches to 75 inches and has averaged about 45 inches over a 35-year period. Approximately 20 inches of water are returned to the atmosphere each year by evaporation and transpiration; the remainder of the annual precipitation either flows overland to streams or percolates downward to the water table and ultimately flows out of the basin in the Shetucket River or as underflow through the deposits beneath. During the autumn and winter months precipitation normally is sufficient to cause a substantial increase in the amount of water stored underground and in surface reservoirs within the basins whereas in the summer most of the precipitation is lost through evaporation and transpiration, resulting in sharply reduced streamflow and lowered groundwater levels. The mean monthly storage of water in the basin on an average is 3.5 inches higher in November than it is in June. The amount of water that flows out of the basin in the Shetucket River represents the total amount of water potentlally available for use by man. Annual runoff from the entire basin above the Quinebaug River has ranged from about 13 to 42 inches since 1929, and has averaged about 23 inches (300 billion gallons). Although runoff indicates the total amount of water potentially available, it is usually not economically or legally feasible for man to use all of it. On the other hand, with increased development, It is possible that some water will be reused several times. The water available may be tapped as it flows through the area or is temporarily stored in streams, lakes, and aquifers. The amounts that can be developed vary from place to place and time to time, depending on the amount of precipitation, on the size of drainage area, on the thickness, permeability and areal extent of aquifers, and on the variations in chemical and physical quality of the water. Differences in streamflow from point to point are due primarily to differences in the proportion of stratified drift in the drainage basin above each point, which affect the timing of streamflow, and to differences in precipitation, which affect the amount of streamflow. Information on streamflow from gaging stations may be extended to ungaged sites by accounting for both of these factors ,in calculations. Future floods on the upper Willimantic River or the Shetucket River are unlikely to cause major damage so long as buildings are not constructed below the highest flood elevations to be expected with the present system of reservoirs for flood control. Ground water can be obtained from wells almost anywhere in the Shetucket River basin, but the amount obtainable from individual wells at any particular point depends upon the type and water-bearing properties of the aquifers present. For practical purposes, the earth materials in the basin comprise three aquifers--stratified drift, till, and bedrock, Stratified drift is the only aquifer generally capable of yielding more than 100 gpm to individual wells. This aquifer covers about 18 percent of the basin and occurs chiefly In lowlands where it overlies till or bedrock. Coefficient of permeability of the coarse-grained unit of stratified drift averages about 1,900 gpd per sq ft. Drilled, screened wells tapping this unit, are known to yield from 200 to 675 gpm. Dug wells in coarse-grained stratified drift should supply at least 2 gpm per foot of drawdown over an 8-hour period. Fine-grained stratified drift has an average coefficient of permeability of about 400 gpd per sq ft and can usually yield to dug wells supplies sufficient for household use. Till and bedrock are widespread in extent but can provide only small to moderate water supplies. Till is tapped chiefly by dug wells; permanent supplies of more than 200 gpd can be obtained from dug wells at a majority of sites in areas of till, but there are many sites where the till is too impermeable or too thin to provide this much water throughout the year. The coefficient of permeability of till ranges from about 0.2 gpd per sq ft to 55 gpd per sq it. Bedrock Is tapped chiefly by drilled wells, about 90 percent of which will supply at least 3 gpm. Very few, however, will supply more than 50 gpm. The amount of ground water potentially available In an area depends upon the amount of groundwater outflow, the amount of ground water in storage, and the quantity of water available by Induced infiltration from streams and lakes. From data on permeability, saturated thickness, recharge, yield from aquifer storage, well performance, and streamflow, preliminary estimates of ground-water availability can be made for any point in the basin. Long-term yields estimated for 15 areas especially favorable for development of large ground-water supplies ranged from 1.3 to 61.8 mgd. Detailed site studies to determine optimum yields, drawdowns, and spacing of individual wells are needed before major ground-water development is undertaken In these or other areas. The chemical quality of water in the Shetucket basin Is generally good to excellent. Samples of naturally occurring surface water collected from 32 sites contained less than 61 ppm of dissolved solids and less than 32 ppm of hardness. Water from wells is more highly mineralized than naturally occurring water from streams. Even so only 7 percent of wells sampled yielded water with more than 200 ppm of dissolved sol-ids and only 9 percent yielded water with more than 120 ppm of hardness. Even in the major rivers, which are used to transport industrial waste, the dissolved mineral content is less than 100 ppm and hardness rarely exceeds 40 ppm. One notable exception occurs in the lower reaches of Little River where an exceptional amount of industrial waste is discharged into the river near Versailles. This waste is particularly noticeable during low streamflow. Iron and manganese In both ground water and surface water are the only constituents whose concentrations commonly exceed recommended limits for domestic and industrial use. Most wells in the basin yield clear water with little or no iron or manganese, but distributed among them are wells with ground water that contains enough of these dissolved constituents to be troublesome for most uses. iron concentrations in naturally occurring stream water exceeded 0.3 ppm under tow-flow conditions at 20 percent of the sites sampled. Large concentrations of iron in stream water result from discharge of iron-bearing ground water or from the discharge of water from swamps. In swamps the iron is released largely from decaying vegetation. Ground water more than 30 feet below the land surface has a relatively constant temperature, usually between 48°F and 50°F. Water temperature in very shallow wells may fluctuate from about 38°F in February or March to about 55°F in late summer. Water temperature in the larger streams fluctuates much more widely, ranging from 32°F at least for brief periods in winter, to about 85°F occasionally during The quantity of suspended sediment transported by streams in the basin is negligible, though amounts large enough to be troublesome may occur locally at times. The total amount of water used In the Shetucket Rlver basin for all purposes during 1961 was about 5,810 million gallons~ which is equivalent to 208 gpd per person, Public water systems supplied the domestic needs of nearly half the population of the basin; 10 systems were sampled, all of which provided water of better quality than the U.S. Public Health Service suggests for drinking water standards.

Connecticut↗