Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “One Earth”

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

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

At least 1,207 records · Page 67Linked to original sources
USGS2012

Hydrogeology and water quality of the Floridan aquifer system and effect of Lower Floridan aquifer pumping on the Upper Floridan aquifer, Pooler, Chatham County, Georgia, 2011–2012

Two test wells were completed in Pooler, Georgia, in 2011 to investigate the potential of using the Lower Floridan aquifer as a source of water for municipal use. One well was completed in the Lower Floridan aquifer at a depth of 1,120 feet (ft) below land surface; the other well was completed in the Upper Floridan aquifer at a depth of 486 ft below land surface. At the Pooler test site, the U.S. Geological Survey performed flowmeter surveys, packer-isolated slug tests within the Lower Floridan confining unit, slug tests of the entire Floridan aquifer system, and aquifer tests of the Upper and Lower Floridan aquifers. Drill cuttings, geophysical logs, and borehole flowmeter surveys indicate that the Upper Floridan aquifer extends 333 –515 ft below land surface, the Lower Floridan confining unit extends 515–702 ft below land surface, and the Lower Floridan aquifer extends 702–1,040 ft below land surface. Flowmeter surveys indicate that the Upper Floridan aquifer contains two water-bearing zones at depth intervals of 339 –350 and 375–515 ft; the Lower Floridan confining unit contains one zone at a depth interval of 550–620 ft; and the Lower Floridan aquifer contains five zones at depth intervals of 702–745, 745–925, 925–984, 984–1,015, and 1,015–1,040 ft. Flowmeter testing of the test borehole open to the entire Floridan aquifer system indicated that the Upper Floridan aquifer contributed 92.4 percent of the total flow rate of 708 gallons per minute; the Lower Floridan confining unit contributed 3.0 percent; and the Lower Floridan aquifer contributed 4.6 percent. Horizontal hydraulic conductivity of the Lower Floridan confining unit derived from slug tests within three packer-isolated intervals ranged from 0.5 to 10 feet per day (ft/d). Aquifer-test analyses yielded values of transmissivity for the Upper Floridan aquifer, Lower Floridan confining unit, and the Lower Floridan aquifer of 46,000, 700, and 4,000 feet squared per day (ft 2 /d), respectively. Horizontal hydraulic conductivity of 4 ft/d for the Lower Floridan confining unit, derived from aquifer-test analyses, is near the midrange for values derived from packer-isolated slug tests. The transmissivity of the entire Floridan aquifer system derived from aquifer-test analyses totals about 51,000 ft 2 /d, similar to the value of 58,000 ft 2 /d derived from open slug tests on the entire Floridan aquifer system. Water-level data for each aquifer test were filtered for external influences such as barometric pressure, earth-tide effects, and long-term trends to enable detection of small (less than 1 foot) water-level responses to aquifer-test pumping. During the 72-hour aquifer test of pumping the Lower Floridan aquifer, a drawdown response of 51.7 ft was observed in the Lower Floridan pumped well and a drawdown response of 0.9 foot was observed in the Upper Floridan observation well located 85 ft from the pumped well.

Georgia↗
USGS2009

A critical evaluation of crustal dehydration as the cause of an overpressured and weak San Andreas Fault

Many plate boundary faults, including the San Andreas Fault, appear to slip at unexpectedly low shear stress. One long-standing explanation for a "weak" San Andreas Fault is that fluid release by dehydration reactions during regional metamorphism generates elevated fluid pressures that are localized within the fault, reducing the effective normal stress. We evaluate this hypothesis by calculating realistic fluid production rates for the San Andreas Fault system, and incorporating them into 2-D fluid flow models. Our results show that for a wide range of permeability distributions, fluid sources from crustal dehydration are too small and short-lived to generate, sustain, or localize fluid pressures in the fault sufficient to explain its apparent mechanical weakness. This suggests that alternative mechanisms, possibly acting locally within the fault zone, such as shear compaction or thermal pressurization, may be necessary to explain a weak San Andreas Fault. More generally, our results demonstrate the difficulty of localizing large fluid pressures generated by regional processes within near-vertical fault zones. ?? 2009 Elsevier B.V.

Earth and Planetary Science Letters↗
USGS2024

Subduction zone geometry modulates the megathrust earthquake cycle: Magnitude, recurrence, and variability

Megathrust geometric properties exhibit some of the strongest correlations with maximum earthquake magnitude in global surveys of large subduction zone earthquakes, but the mechanisms through which fault geometry influences subduction earthquake cycle dynamics remain unresolved. Here, we develop 39 models of sequences of earthquakes and aseismic slip (SEAS) on variably-dipping planar and variably-curved nonplanar megathrusts using the volumetric, high-order accurate code tandem to account for fault curvature. We vary the dip, downdip curvature and width of the seismogenic zone to examine how slab geometry mechanically influences megathrust seismic cycles, including the size, variability, and interevent timing of earthquakes. Dip and curvature control characteristic slip styles primarily through their influence on seismogenic zone width: wider seismogenic zones allow shallowly-dipping megathrusts to host larger earthquakes than steeply-dipping ones. Under elevated pore pressure and less strongly velocity-weakening friction, all modeled fault geometries host uniform periodic ruptures. In contrast, shallowly-dipping and sharply-curved megathrusts host multi-period supercycles of slow-to-fast, small-to-large slip events under higher effective stresses and more strongly velocity-weakening friction. We discuss how subduction zones' maximum earthquake magnitudes may be primarily controlled by the dip and dimensions of the seismogenic zone, while second-order effects from structurally-derived mechanical heterogeneity modulate the recurrence frequency and timing of these events. Our results suggest that enhanced co- and interseismic strength and stress variability along the megathrust, such as induced near areas of high or heterogeneous fault curvature, limits how frequently large ruptures occur and may explain curved faults' tendency to host more frequent, smaller earthquakes than flat faults.

Journal of Geophysical Research: Solid Earth↗
USGS2014

Post-Hurricane Sandy coastal oblique aerial photographs collected from Cape Lookout, North Carolina, to Montauk, New York, November 4-6, 2012

The U.S. Geological Survey (USGS) conducts baseline and storm response photography missions to document and understand the changes in vulnerability of the Nation's coasts to extreme storms. On November 4-6, 2012, approximately one week after the landfall of Hurricane Sandy, the USGS conducted an oblique aerial photographic survey from Cape Lookout, N.C., to Montauk, N.Y., aboard a Piper Navajo Chieftain (aircraft) at an altitude of 500 feet (ft) and approximately 1,000 ft offshore. This mission was flown to collect post-Hurricane Sandy data for assessing incremental changes in the beach and nearshore area since the last survey in 2009. The data can be used in the assessment of future coastal change. The photographs provided here are Joint Photographic Experts Group (JPEG) images. The photograph locations are an estimate of the position of the aircraft and do not indicate the location of the feature in the images. These photos document the configuration of the barrier islands and other coastal features at the time of the survey. Exiftool was used to add the following to the header of each photo: time of collection, Global Positioning System (GPS) latitude, GPS longitude, keywords, credit, artist (photographer), caption, copyright, and contact information. Photographs can be opened directly with any JPEG-compatible image viewer by clicking on a thumbnail on the contact sheet. Table 1 provides detailed information about the GPS location, image name, date, and time each of the 9,481 photographs were taken, along with links to each photograph. The photographs are organized in segments, also referred to as contact sheets, and represent approximately 5 minutes of flight time. In addition to the photographs, a Google Earth Keyhole Markup Language (KML) file is provided and can be used to view the images by clicking on the marker and then clicking on either the thumbnail or the link above the thumbnail. The KML files were created using the photographic navigation files.

Delaware, Maryland, New Jersey, New York, North Ca↗
USGS1991

A comparison of Sr-Nd-Pb isotopes in young and old continental lithospheric mantle: Patagonia and eastern China

It is commonly accepted that beneath the continental crust lies a keel of lithospheric mantle, which extends 50–200 kilometres downward to a transition zone into the asthenosphere. The chemical and physical properties of this reservoir are best known through studies of the basalts and xenoliths that provide samples of the subcrustal mantle. Although sharing many characteristics with oceanic island basalts, some continental basalts become increasingly distinct isotopically as crustal age increases, strongly supporting a permanent association between crust and mantle. Consequently, the distinctive trace element and isotope composition of the lithospheric mantle is able to give important clues to its origin and evolution. The mantle under newly‐created crust is typified by a radiogenic isotope variability that emphasizes the materials from which the continental lithosphere is assembled. Old lithospheric mantle, on the other hand, exhibits more evolved isotopic patterns that attest to the existence of long‐lived, chemically complex systems. A comparison of the Pb, Sr and Nd isotopes in alkalic to sub‐alkalic basalt derived from Phanerozoic (Patagonia) and Middle Archaean to Early Proterozoic (eastern China) subcrustal mantle is useful for identifying ‘end‐member’ components of the lithosphere. One component, having an isotopic composition close to PREMA, either continues to evolve virtually unchanged after incorporation into the lithosphere or is, itself, a relatively new addition even to old lithosphere. Another component, beginning with the isotopic composition of BSE, undergoes significant reduction in U/Pb and Sm/Nd (but not Rb/Sr) upon incorporation into the lithosphere and, with time, shows an increasingly retarded evolution of 206 Pb/ 204 Pb and negative ε Nd‐values approaching the isotopic composition of EMI. Five models are discussed that relate the isotopic composition of the continental lithospheric mantle to that of other parts of the terrestrial system, which may be involved in its origin and evolution. The potential locations of the contributing components and the mechanisms and timing of their assembly into lithosphere are considered. Current knowledge, however, does not allow us to distinguish unequivocally among the various scenarios for the creation and evolution of this reservoir.

Australian Journal of Earth Sciences↗
USGS2006

Submarine landslides in the Santa Barbara Channel as potential tsunami sources

Recent investigations using the Monterey Bay Aquarium Research Institutes (MBARI) Remotely Operated Vehicles (ROVs) "Ventana" and "Tiburon" and interpretation of MBARI's EM 300 30 kHz multibeam bathymetric data show that the northern flank of the Santa Barbara Basin has experienced massive slope failures. Of particular concern is the large (130 km2) Goleta landslide complex located off Coal Oil Point near the town of Goleta, that measures 14.6-km long extending from a depth of 90 m to nearly 574 m deep and is 10.5 km wide. We estimate that approximately 1.75 km3 has been displaced by this slide during the Holocene. This feature is a complex compound submarine landslide that contains both surfical slump blocks and mud flows in three distinct segments. Each segment is composed of a distinct head scarp, down-dropped head block and a slide debris lobe. The debris lobes exhibit hummocky topography in the central areas that appear to result from compression during down slope movement. The toes of the western and eastern lobes are well defined in the multibeam image, whereas the toe of the central lobe is less distinct. Continuous seismic reflection profiles show that many buried slide debris lobes exist and comparison of the deformed reflectors with ODP Drill Site 149, Hole 893 suggest that at least 200 000 years of failure have occurred in the area (Fisher et al., 2005a). Based on our interpretation of the multibeam bathymetry and seismic reflection profiles we modeled the potential tsunami that may have been produced from one of the three surfical lobes of the Goleta slide. This model shows that a 10 m high wave could have run ashore along the cliffs of the Goleta shoreline. Several other smaller (2 km2 and 4 km2) slides are located on the northern flank of the Santa Barbara Basin, both to the west and east of Goleta slide and on the Concepcion fan along the western flank of the basin. One slide, named the Gaviota slide, is 3.8 km2, 2.6 km long and 1.7 km wide. A distinct narrow scar extends from near the eastern head wall of this slide for over 2 km eastward toward the Goleta slide and may represent either an incipient failure or a remnant of a previous failure. Push cores collected within the main head scar of this slide consisted of hydrogen sulfide bearing mud, possibly suggesting active fluid seepage and a vibra-core penetrated ???50 cm of recent sediment overlying colluvium or landslide debris confirming the age of ???300 years as proposed by Lee et al. (2004). However, no seeps or indications of recent movement were observed during our ROV investigation within this narrow head scar indicating that seafloor in the scar is draped with mud. ?? 2006 Author(s). This work is licensed under a Creative Commons License.

Natural Hazards and Earth System Sciences↗
USGS2008

Categorizing natural disaster damage assessment using satellite-based geospatial techniques

Remote sensing of a natural disaster's damage offers an exciting backup and/or alternative to traditional means of on-site damage assessment. Although necessary for complete assessment of damage areas, ground-based damage surveys conducted in the aftermath of natural hazard passage can sometimes be potentially complicated due to on-site difficulties (e.g., interaction with various authorities and emergency services) and hazards (e.g., downed power lines, gas lines, etc.), the need for rapid mobilization (particularly for remote locations), and the increasing cost of rapid physical transportation of manpower and equipment. Satellite image analysis, because of its global ubiquity, its ability for repeated independent analysis, and, as we demonstrate here, its ability to verify on-site damage assessment provides an interesting new perspective and investigative aide to researchers. Using one of the strongest tornado events in US history, the 3 May 1999 Oklahoma City Tornado, as a case example, we digitized the tornado damage path and co-registered the damage path using pre- and post-Landsat Thematic Mapper image data to perform a damage assessment. We employed several geospatial approaches, specifically the Getis index, Geary's C , and two lacunarity approaches to categorize damage characteristics according to the original Fujita tornado damage scale (F-scale). Our results indicate strong relationships between spatial indices computed within a local window and tornado F-scale damage categories identified through the ground survey. Consequently, linear regression models, even incorporating just a single band, appear effective in identifying F-scale damage categories using satellite imagery. This study demonstrates that satellite-based geospatial techniques can effectively add spatial perspectives to natural disaster damages, and in particular for this case study, tornado damages.

Natural Hazards and Earth System Sciences↗
USGS2009

Global irrigated area map (GIAM), derived from remote sensing, for the end of the last millennium

A Global Irrigated Area Map (GIAM) has been produced for the end of the last millennium using multiple satellite sensor, secondary, Google Earth and groundtruth data. The data included: (a) Advanced Very High Resolution Radiometer (AVHRR) 3-band and Normalized Difference Vegetation Index (NDVI) 10 km monthly time-series for 1997-1999, (b) Syste me pour l'Observation de la Terre Vegetation (SPOT VGT) NDVI 1 km monthly time series for 1999, (c) East Anglia University Climate Research Unit (CRU) rainfall 50km monthly time series for 1961-2000, (d) Global 30 Arc-Second Elevation Data Set (GTOPO30) 1 km digital elevation data of the World, (e) Japanese Earth Resources Satellite-1 Synthetic Aperture Radar (JERS-1 SAR) data for the rain forests during two seasons in 1996 and (f) University of Maryland Global Tree Cover 1 km data for 1992-1993. A single mega-file data-cube (MFDC) of the World with 159 layers, akin to hyperspectral data, was composed by re-sampling different data types into a common 1 km resolution. The MFDC was segmented based on elevation, temperature and precipitation zones. Classification was performed on the segments. Quantitative spectral matching techniques (SMTs) used in hyperspectral data analysis were adopted to group class spectra derived from unsupervised classification and match them with ideal or target spectra. A rigorous class identification and labelling process involved the use of: (a) space-time spiral curve (ST-SC) plots, (b) brightness-greenness-wetness (BGW) plots, (c) time series NDVI plots, (d) Google Earth very-high-resolution imagery (VHRI) 'zoom-in views' in over 11 000 locations, (e) groundtruth data broadly sourced from the degree confluence project (3 864 sample locations) and from the GIAM project (1 790 sample locations), (f) high-resolution Landsat-ETM+ Geocover 150m mosaic of the World and (g) secondary data (e.g. national and global land use and land cover data). Mixed classes were resolved based on decision tree algorithms and spatial modelling, and when that did not work, the problem class was used to mask and re-classify the MDFC, and the class identification and labelling protocol repeated. The sub-pixel area (SPA) calculations were performed by multiplying full-pixel areas (FPAs) with irrigated area fractions (IAFs) for every class. A 28 class GIAM was produced and the area statistics reported as: (a) annualized irrigated areas (AIAs), which consider intensity of irrigation (i.e. sum of irrigated areas from different seasons in a year plus continuous year-round irrigation or gross irrigated areas), and (b) total area available for irrigation (TAAI), which does not consider intensity of irrigation (i.e. irrigated areas at any given point of time plus the areas left fallow but 'equipped for irrigation' at the same point of time or net irrigated areas). The AIA of the World at the end of the last millennium was 467 million hectares (Mha), which is sum of the non-overlapping areas of: (a) 252 Mha from season one, (b) 174 Mha from season two and (c) 41 Mha from continuous year-round crops. The TAAI at the end of the last millennium was 399 Mha. The distribution of irrigated areas is highly skewed amongst continents and countries. Asia accounts for 79% (370 Mha) of all AIAs, followed by Europe (7%) and North America (7%). Three continents, South America (4%), Africa (2%) and Australia (1%), have a very low proportion of the global irrigation. The GIAM had an accuracy of 79-91%, with errors of omission not exceeding 21%, and the errors of commission not exceeding 23%. The GIAM statistics were also compared with: (a) the United Nations Food and Agricultural Organization (FAO) and University of Frankfurt (UF) derived irrigated areas and (b) national census data for India. The relationships and causes of differences are discussed in detail. The GIAM products are made available through a web portal (http://www.iwmigiam.org). ?? 2009 Taylor & Francis.

International Journal of Remote Sensing↗
USGS2010

Oblique Aerial Photography of the Arctic Coast of Alaska, Cape Sabine to Milne Point, July 16-19, 2009

The Arctic Coastal Plain of northern Alaska, an area of strategic economic importance to the United States, is home to remote Native American communities and encompasses unique habitats of global significance. Coastal erosion along the Arctic coast is chronic and widespread; recent evidence suggests that erosion rates are among the highest in the world (as high as ~16 m/yr) and may be accelerating. Coastal erosion adversely impacts energy-related infrastructure, natural shoreline habitats, and Native American communities. Climate change is thought to be a key component of recent environmental changes in the Arctic. Reduced sea-ice cover in the Arctic Ocean is one of the probable mechanisms responsible for increasing coastal exposure to wave attack and the resulting increase in erosion. Extended periods of permafrost melting and associated decreases in bluff cohesion and stability are another possible source of the increase in erosion. Several studies of selected areas on the Alaska coast document past shoreline positions and coastal change, but none have examined the entire North coast systematically. Results from these studies indicate high rates of coastal retreat that vary spatially along the coast. To address the need for a comprehensive and regionally consistent evaluation of shoreline change along the North coast of Alaska, the U.S. Geological Survey (USGS), as part of their Coastal and Marine Geology Program's (CMGP) National Assessment of Shoreline Change Study, is evaluating shoreline change from Peard Bay to the United States/Canadian border, using historical maps and photography and a standardized methodology that is consistent with other shoreline-change studies along the Nation's coastlines (see, for example, http://coastal.er.usgs.gov/shoreline-change/, last accessed February 12, 2010). This is the second in a series of publications containing photographs collected during reconnaissance surveys conducted in support of the National Assessment of Shoreline Change Study. An accompanying ESRI ArcGIS shape file (and plaintext copy) indicates the position of the aircraft and time when each photograph was taken. The USGS-CMGP Field Activity ID for the survey is A-5-09-AK, and more information on the survey and how to view the photographs using Google Earth software is available online at http://walrus.wr.usgs.gov/infobank/a/a509ak/html/a-5-09-ak.photos.kmz (last accessed February 12, 2010). The initial report ?Oblique aerial photography of the Arctic coast of Alaska, Nulavik to Demarcation Point, August 7-10, 2006? is available online at http://pubs.usgs.gov/ds/436/, and the associated Google Earth .kmz file is available at http://walrus.wr.usgs.gov/infobank/a/a106ak/html/a-1-06-ak.photos.kmz (last accessed February 12, 2010).

Data Series↗
USGS1998

Ductile shear zones beneath strike-slip faults: Implications for the thermomechanics of the San Andreas fault zone

We have carried out two-dimensional (2-D) numerical experiments on the bulk flow of a layer of fluid that is driven in a strike-slip sense by constant velocities applied at its boundaries. The fluid has the (linearized) conventional rheology assumed to apply to lower crust/upper mantle rocks. The temperature dependence of the effective viscosity of the fluid and the shear heating that accompanies deformation have been incorporated into the calculations, as has thermal conduction in an overlying crustal layer. Two end-member boundary conditions have been considered, corresponding to a strong upper crust driving a weaker ductile substrate and a strong ductile layer driving a passive, weak crust. In many cases of practical interest, shear heating is concentrated close to the axial plane of the shear zone for either boundary condition. For these cases, the resulting steady state temperature field is well approximated by a cylindrical heat source embedded in a conductive half-space at a depth corresponding to the top of the fluid layer. This approximation, along with the application of a theoretical result for one-dimensional shear zones, permits us to obtain simple analytical approximations to the thermal effects of 2-D ductile shear zones for a range of assumed rheologies and crustal geotherms, making complex numerical calculations unnecessary. Results are compared with observable effects on heat flux near the San Andreas fault using constraints on the slip distribution across the entire fault system. Ductile shearing in the lower crust or upper mantle can explain the observed increase in surface heat flux southeast of the Mendocino triple junction and match the amplitude of the regional heat flux anomaly in the California Coast Ranges. Because ductile dissipation depends only weakly on slip rate, faults moving only a few millimeters per year can be important heat sources, and the superposition of effects of localized ductile shearing on both currently active and now inactive strands of the San Andreas system can explain the breadth of the heat flux anomaly across central California.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Regional deformation of the Sierra Nevada, California, on conjugate microfault sets

Strike slip microfaults are pervasive throughout the granitic rocks of the eastern Sierra Nevada. Offsets typically range from less than a millimeter to several tens of centimeters but exceed 100 m in some places. The spacing between microfaults varies from a few tens of centimeters to a few tens of meters throughout much of the high Sierra Nevada. Many of these microfaults are loci of slickensided, compact fault gouge, and they are commonly mineralized by quartz veinlets with minor amounts of epidote, chlorite, and rare sulphide minerals. The microfaults are oriented in two nearly vertical conjugate sets; a north to northeast striking set showing right lateral offset and an east to northeast striking set showing left lateral offset. Microfaults with left lateral offset are more common than microfaults with right lateral offset. Most lineaments visible on aerial photographs are microfaults. The age of this microfaulting is not precisely known. It developed after consolidation of the youngest granitic plutons in the Sierra (79 m.y. B.P.) and is known to cut a late Miocene volcanic dike in one area. Slickensides along the microfaults are subhorizontal but show a slight (about 3°) westward plunging average inclination, suggesting that much of the deformation occurred prior to the westward tilting of the Sierran block in late Tertiary time. The direction of maximum horizontal extensional strain (determined as the bisector of average microfault trends) changes systematically from north to south (WNW at 38.5°N; NW at 36.5°N). A pure shear constant volume solution based on a detailed study of microfaults at 37°20′N indicates a maximum extension of 2.3% in a N61°W direction. These extension directions are remarkably parallel to late Mesozoic to present‐day tectonic extension directions in the Basin and Range province. The pattern of microfaulting demonstrates that the supposedly monolithic Sierran terrane was also affected by the late Cenozoic and possibly earlier regional extension of western North America and provides an independent criterion for determining extensional strain directions.

California↗
USGS2010

Shear-rate-dependent strength control on the dynamics of rainfall-triggered landslides, Tokushima Prefecture, Japan

A typhoon (Typhoon No. 10) attacked Shikoku Island and the Tyugoku area of Japan in 2004. This typhoon produced a new daily precipitation record of 1317 mm on Shikoku Island and triggered hundreds of landslides in Tokushima Prefecture. One catastrophic landslide was triggered in the Shiraishi area of Kisawa village, and destroyed more than 10 houses while also leaving an unstable block high on the slope. The unstable block kept moving after the event, showing accelerating and decelerating movement during and after rainfall and reaching a displacement of several meters before countermeasures were put into place. To examine the mechanism for this landsliding characteristic, samples (weathered serpentinite) were taken from the field, and their shear behaviours examined using ring shear tests. The test results revealed that the residual shear strength of the samples is positively dependent on the shear rate, which may provide an explanation for the continuous acceleratingdecelerating process of the landsliding. The roughness of the shear surface and the microstructure of the shear zone were measured and observed by laser microscope and SEM techniques in an attempt to clarify the mechanism of shear rate effect on the residual shear strength. Copyright ?? 2010 John Wiley & Sons, Ltd.

Earth Surface Processes and Landforms↗
USGS2005

Macropolygon morphology, development, and classification on North Panamint and Eureka playas, Death Valley National Park CA

Panamint and Eureka playas, both located within Death Valley National Park, exhibit a host of surficial features including fissures, pits, mounds, and plant-covered ridges, representing topographic highs and lows that vary up to 2 m of relief from the playa surface. Aerial photographs reveal that these linear strands often converge to form polygons, ranging in length from several meters to nearly a kilometer. These features stand out in generally dark contrast to the brighter intervening expanse of flat, plant-free, desiccated mud of the typical playa surface. Ground-truth mapping of playa features with differential GPS (Global Positioning System) was conducted in 1999 (North Panamint Valley) and 2002 (Eureka Valley). High-resolution digital maps reveal that both playas possess macropolygons of similar scale and geometry, and that fissures may be categorized into one of two genetic groups: (1) shore-parallel or playa-interior desiccation and shrinkage; and (2) tectonic-induced cracks. Early investigations of these features in Eureka Valley concluded that their origin may have been related to agricultural activity by paleo-Indian communities. Although human artifacts are abundant at each locale, there is no evidence to support the inference that surface features reported on Eureka Playa are anthropogenic in origin. Our assumptions into the genesis of polygons on playas is based on our fortuitous experience of witnessing a fissure in the process of formation on Panamint Playa after a flash flood (May 1999); our observations revealed a paradox that saturation of the upper playa crusts contributes to the establishment of some desiccation features. Follow-up visits to the same feature over 2 yrs' time are a foundation for insight into the evolution and possible longevity of these features. ?? 2005 Elsevier B.V. All rights reserved.

Earth-Science Reviews↗
USGS2023

Accuracy of finite fault slip estimates in subduction zone regions with topographic Green's functions and seafloor geodesy

Until recently, the lack of seafloor geodetic instrumentation and the use of unrealistically simple, half-space based forward models have resulted in poor resolution of near-trench slip in subduction zone settings. Here, we use a synthetic framework to investigate the impact of topography and geodetic data distribution on coseismic slip estimates in various subduction zone settings. We calculate surface displacements in two synthetic topographic domains that have topography similar to that of Chile and Japan, respectively. We then attempt to image target slip distributions by using a Bayesian approach to solve for slip with two sets of Green's functions—one that accounts for topography and one that does not—and five sets of 50 or more observation points selected from the synthetic surface displacements. Three of these sets of observation points are entirely onland, and two include 5–10 seafloor geodetic sites. We find that the use of topographic Green's functions always improves inferred slip models, and with seafloor geodetic data, it enables an almost perfect recovery of a target slip model, even in the near-trench region. Critically, our results demonstrate that it would be impossible for non-topographic Green's functions to properly recover the true slip distribution, particularly in the near-trench region. We also perform a parameter study with approximately 4,000 slip models estimated using a least-square approach, and find that topographic Green's functions yield significantly more accurate slip models in cases where good data (well distributed and reasonably dense) are available, even in the absence of seafloor geodetic sites.

Journal of Geophysical Research: Solid Earth↗
USGS2023

Modular compositional learning improves 1D hydrodynamic lake model performance by merging process-based modeling with deep learning

Hybrid Knowledge-Guided Machine Learning (KGML) models, which are deep learning models that utilize scientific theory and process-based model simulations, have shown improved performance over their process-based counterparts for the simulation of water temperature and hydrodynamics. We highlight the modular compositional learning (MCL) methodology as a novel design choice for the development of hybrid KGML models in which the model is decomposed into modular sub-components that can be process-based models and/or deep learning models. We develop a hybrid MCL model that integrates a deep learning model into a modularized, process-based model. To achieve this, we first train individual deep learning models with the output of the process-based models. In a second step, we fine-tune one deep learning model with observed field data. In this study, we replaced process-based calculations of vertical diffusive transport with deep learning. Finally, this fine-tuned deep learning model is integrated into the process-based model, creating the hybrid MCL model with improved overall projections for water temperature dynamics compared to the original process-based model. We further compare the performance of the hybrid MCL model with the process-based model and two alternative deep learning models and highlight how the hybrid MCL model has the best performance for projecting water temperature, Schmidt stability, buoyancy frequency, and depths of different isotherms. Modular compositional learning can be applied to existing modularized, process-based model structures to make the projections more robust and improve model performance by letting deep learning estimate uncertain process calculations.

Journal of Advances in Modeling Earth Systems↗
USGS1997

Wrinkle-like slip pulse on a fault between different materials

Pulses of slip velocity can propagate on a planar interface governed by a constant coefficient of friction, where the interface separates different elastic materials. Such pulses have been found in two-dimensional plane strain finite difference calculations of slip on a fault between elastic media with wave speeds differing by 20%. The self-sustaining propagation of the slip pulse arises from interaction between normal and tangential deformation that exists only with a material contrast. These calculations confirm the prediction of Weertman [1980] that a dislocation propagating steadily along a material interface has a tensile change of normal traction with the same pulse shape as slip velocity. The self-sustaining pulse is associated with a rapid transition from a head wave traveling along the interface with the S wave speed of the faster material, to an opposite polarity body wave traveling with the slower S speed. Slip occurs during the reversal of normal particle velocity. The pulse can propagate in a region with constant coefficient of friction and an initial stress state below the frictional criterion. Propagation occurs in only one direction, the direction of slip in the more compliant medium, with rupture velocity near the slower S wave speed. Displacement is larger in the softer medium, which is displaced away from the fault during the passage of the slip pulse. Motion is analogous to a propagating wrinkle in a carpet. The amplitude of slip remains approximately constant during propagation, but the pulse width decreases and the amplitudes of slip velocity and stress change increase. The tensile change of normal traction increases until absolute normal traction reaches zero. The pulse can be generated as a secondary effect of a drop of shear stress in an asperity. The pulse shape is unstable, and the initial slip pulse can change during propagation into a collection of sharper pulses. Such a pulse enables slip to occur with little loss of energy to friction, while at the same time increasing irregularity of stress and slip at the source.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Fracture characterization by means of attenuation and generation of tube waves in fractured crystalline rock at Mirror Lake, New Hampshire

Results are presented from experiments carried out in conjunction with the U. S. Geological Survey at the Hubbard Brook Experimental Forest near Mirror Lake, New Hampshire. The study focuses on our ability to obtain orientation and transmissivity estimates of naturally occurring fractures. The collected data set includes a four-offset hydrophone vertical seismic profile, full waveform acoustic logs at 5, 15, and 34 kHz, borehole televiewer, temperature, resistivity, and self-potential logs, and borehole-to-borehole pump test data. Borehole televiewer and other geophysical logs indicate that permeable fractures intersect the Mirror Lake boreholes at numerous depths, but less than half of these fractures appear to have significant permeability beyond the annulus of drilling disturbance on the basis of acoustic waveform log analysis. The vertical seismic profiling (VSP) data indicate a single major permeable fracture near a depth of 44 m, corresponding to one of the most permeable fractures identified in the acoustic waveform log analysis. VSP data also indicate a somewhat less permeable fracture at 220 m and possible fractures at depths of 103 and 135 m; all correspond to major permeable fractures in the acoustic waveform data set. Pump test data confirm the presence of a hydraulic connection between the Mirror Lake boreholes through a shallow dipping zone of permeability at 44 m in depth. Effective fracture apertures calculated from modeled transmissivities correspond to those estimated for the largest fractures indicated on acoustic waveform logs but are over an order of magnitude larger than effective apertures calculated from tube waves in the VSP data set. This discrepancy is attributed to the effect of fracture stiffness. A new model is presented to account for the mechanical strength of asperities in resisting fracture closure during the passage of seismic waves during the generation of VSPs.

New Hampshire↗
USGS2013

Seismic Station Installation Orientation Errors at ANSS and IRIS/USGS Stations

Many seismological studies depend on the published orientations of sensitive axes of seismic instruments relative to north (e.g., Li et al., 2011). For example, studies of the anisotropic structure of the Earth’s mantle through SKS‐splitting measurements (Long et al., 2009), constraints on core–mantle electromagnetic coupling from torsional normal‐mode measurements (Dumberry and Mound, 2008), and models of three‐dimensional (3D) velocity variations from surface waves (Ekström et al., 1997) rely on accurate sensor orientation. Unfortunately, numerous results indicate that this critical parameter is often subject to significant error (Laske, 1995; Laske and Masters, 1996; Yoshizawa et al., 1999; Schulte‐Pelkum et al., 2001; Larson and Ekström, 2002). For the Advanced National Seismic System (ANSS; ANSS Technical Integration Committee, 2002), the Global Seismographic Network (GSN; Butler et al., 2004), and many other networks, sensor orientation is typically determined by a field engineer during installation. Successful emplacement of a seismic instrument requires identifying true north, transferring a reference line, and measuring the orientation of the instrument relative to the reference line. Such an exercise is simple in theory, but there are many complications in practice. There are four commonly used methods for determining true north at the ANSS and GSN stations operated by the USGS Albuquerque Seismological Laboratory (ASL), including gyroscopic, astronomical, Global Positioning System (GPS), and magnetic field techniques. A particular method is selected based on site conditions (above ground, below ground, availability of astronomical observations, and so on) and in the case of gyroscopic methods, export restrictions. Once a north line has been determined, it must be translated to the sensor location. For installations in mines or deep vaults, this step can include tracking angles through the one or more turns in the access tunnel leading to the vault (e.g., GSN station WCI in Wyandotte Cave, Indiana). Finally, the third source of error comes from the ability of field engineers to orient the sensor relative to the reference line. In order to quantify bounds on the errors in each step in the orientation process, we conducted a series of tests at the ASL using twelve GSN and ANSS field engineers. The results from this exercise allow us to estimate upper bounds on the precision of our ability to orient instruments, as well as identify the sources of error in the procedures. We are also able to identify systematic bias of various true‐north‐finding methods relative to one another. Although we are unable to estimate the absolute accuracy of our orientation measurements due to our inability to identify true north without some error, the agreement between independent methods for finding true north provides confidence in the different approaches, assuming no systematic bias. Finally, our study neglects orientation errors that are beyond the control of the field engineer during a station visit. These additional errors can arise from deviations in the sensitive axes of the instruments relative to the case markings, processing errors (Holcomb, 2002) when comparing horizontal orientations relative to other sensors (e.g., borehole installations), and deviations of the sensitive axes of instruments from true orthogonality (e.g., instruments with separate modules such as the Streckeisen STS‐1).

Seismological Research Letters↗