Search USGS⌕ Search

SEARCH · Search USGS

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 829 records · Page 46Linked to original sources

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 ‘stochastic geometries’. 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↗

Characterization of the spatial variability of channel morphology

The spatial variability of two fundamental morphological variables is investigated for rivers having a wide range of discharge (five orders of magnitude). The variables, water-surface width and average depth, were measured at 58 to 888 equally spaced cross-sections in channel links (river reaches between major tributaries). These measurements provide data to characterize the two-dimensional structure of a channel link which is the fundamental unit of a channel network. The morphological variables have nearly log-normal probability distributions. A general relation was determined which relates the means of the log-transformed variables to the logarithm of discharge similar to previously published downstream hydraulic geometry relations. The spatial variability of the variables is described by two properties: (1) the coefficient of variation which was nearly constant (0.13-0.42) over a wide range of discharge; and (2) the integral length scale in the downstream direction which was approximately equal to one to two mean channel widths. The joint probability distribution of the morphological variables in the downstream direction was modelled as a first-order, bivariate autoregressive process. This model accounted for up to 76 per cent of the total variance. The two-dimensional morphological variables can be scaled such that the channel width-depth process is independent of discharge. The scaling properties will be valuable to modellers of both basin and channel dynamics. Published in 2002 John Wiley and Sons, Ltd.

Earth Surface Processes and Landforms↗

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↗

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↗

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

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

Alaska↗

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↗

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↗

The crystal structure and thermal history of orthopyroxene from lunar anorthosite 15415

A single crystal of untwinned orthopyroxene from lunar anorthosite sample 15415, with composition (Mg 1.14 Fe 0.80 Mn 0.02 Ca 0.04 )(Si 1.97 Al 0.03 )O 6 , has a unit cell in space group Pbca with a = 18.310(15)Å, b = 8.904(10)Å, c = 5.214(7)Å, containing 2 formula units. A set of 742 counter-measured intensity data made with MoK α radiation has been used to refine the crystal structure in isotropic thermal mode to R = 0.116. Anisotropic refinement led to R = 0.092, but thermal parameters are distorted by non-random errors resulting from poor crystal texture. The resulting structure is in close agreement with that obtained by Ghose [9] for a hypersthene from Greenland. A parameter q , which gives ( Mg q Fe 1−q ) for cation site M (1) and (Mg 1.14− q Fe q −0.18 Ca 0.04 ) for site M (2), was included in the least-squares analysis, yielding q = 0.90(1). This orthopyroxene has the high degree of cation order expected of pyroxenes subjected to Apollonian metamorphism at lower than 500–600°C. No evidence exists for a subsequent thermal event of sufficient intensity to disorder the pyroxene. On the basis of previous laboratory studies of argon-release patterns of lunar plagioclase and order-disorder kinetics of terrestrial pyroxenes, we attribute the reported isotopic age (3.9–4.1 AE) to cessation of metamorphism, perhaps caused by impact excavation.

Earth and Planetary Science Letters↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗