Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Maps”

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,027 records · Page 57Linked to original sources
USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗
USGS2017

Revised tephra volumes for Cascade Range volcanoes

Isopach maps from tephra eruptions from Mount St. Helens were reported in Carey et al. (1995) and for tephra eruptions from Glacier Peak in Gardner et al. (1998). For exponential thinning, the isopach data only define a single slope on a log thickness versus square root of area plot. Carey et al. (1995) proposed a model that was used to estimate a second slope, and volumes were presented in both studies using this model. A study by Sulpizio (2005) for estimating the second slope and square root of area where the lines intersect involves a systematic analysis of many eruptions to provide correlation equations. The purpose of this paper is to recalculate the volumes of Cascades eruptions and compare results from the two methods. In order to gain some perspective on the methods for estimating the second slope, we use data for thickness versus distance beyond the last isopach that are available for some of the larger eruptions in the Cascades. The thickness versus square root of area method is extended to thickness versus distance by developing an approximate relation between the two assuming elliptical isopachs with the source at one of the foci. Based on the comparisons made between the Carey et al. (1995) and Sulpizio (2005) methods, it is felt that the later method provides a better estimate of the second slope. For Mount St. Helens, the estimates of total volume using the Sulpizio (2005) method are generally smaller than those using the Carey et al. (1995) method. For the volume estimates of Carey et al. (1995), the volume of the May 18, 1980, eruption of Mount St. Helens is smaller than six of the eight previous eruptions. With the new volumes using the Sulpizio (2005) method, the 1980 eruption is smaller in volume than the upper end of the range for only three of the layers (Wn, Ye, and Yn) and is the same size as layer We. Thus the 1980 eruption becomes representative of the mid-range of volumes rather than being in the lower range.

Cascade Range volcanoes↗
USGS2025

Greater sage-grouse seasonal habitat associations: A review and considerations for interpretation and management applications

Habitat features needed by wildlife can change in composition throughout the year, particularly in temperate ecosystems, leading to distinct seasonal spatial-use patterns. Studies of species-habitat associations therefore often focus on understanding relationships within discrete seasonal periods with common goals of prediction (e.g., habitat mapping) and inference (e.g., interpreting model coefficients). Across the range of the greater sage-grouse ( Centrocercus urophasianus ) of western North America, the increasing use of high-frequency tracking devices has led to a surge in habitat association studies covering multiple temporal periods and spatial extents. We reviewed the literature for seasonal habitat association studies corresponding to the second and third orders of selection (Johnson 1980). Our objectives were to summarize the methodological approaches used to estimate habitat associations to aid understanding in cross-study comparisons and identify common habitat features reported as selected or avoided within different seasonal periods. We reviewed 71 second- and third-order studies published from 2007–2023 that evaluated covariates collected in a geographic information system (GIS) and modeled probability of selection or intensity of use. The most common study design evaluated a single level of selection (third-order) and was multi-scale (i.e., covariates were measured at varying spatial grains). The most common model approach estimated habitat associations using resource selection functions (RSFs) fit with logistic regression. Studies mostly focused on the breeding periods and winter, but all seasons throughout the annual cycle were covered. There was clear support for selection of sagebrush and avoidance of trees and rugged terrain across seasons, and strong selection of mesic conditions in summer. However, habitat associations for most covariates were mixed, with proportionally equivalent selection and avoidance reported, even within the same seasons. Different factors hampered cross-study comparisons, including variation in study design, but additional contributors likely included important context-dependent habitat associations, such as functional responses to changing habitat availability. We suggest collaborative studies leveraging multiple datasets can help improve seasonal habitat inference by removing the effects of variable study designs.

Journal of Wildlife Management↗
USGS2008

Response of Late Cretaceous migrating deltaic facies systems to sea level, tectonics, and sediment supply changes, New Jersey Coastal Plain, U.S.A.

Paleogeographic, isopach, and deltaic lithofacies mapping of thirteen depositional sequences establish a 35 myr high resolution (> 1 Myr) record of Late Cretaceous wave- and tide-influenced deltaic sedimentation. We integrate sequences defined on the basis of lithologic, biostratigraphic, and Sr-isotope stratigraphy from cores with geophysical log data from 28 wells to further develop and extend methods and calibrations of well-log recognition of sequences and facies variations. This study reveals the northeastward migration of depocenters from the Cenomanian (ca. 98 Ma) through the earliest Danian (ca. 64 Ma) and documents five primary phases of paleodeltaic evolution in response to long-term eustatic changes, variations in sediment supply, the location of two long-lived fluvial axes, and thermoflexural basement subsidence: (1) Cenomanian-early Turonian deltaic facies exhibit marine and nonmarine facies and are concentrated in the central coastal plain; (2) high sediment rates, low sea level, and high accommodation rates in the northern coastal plain resulted in thick, marginal to nonmarine mixed-influenced deltaic facies during the Turonign-Coniacian; (3) comparatively low sediment rates and high long-term sea level in the Santonian resulted in a sediment-starved margin with low deltaic influence; (4) well-developed Campanian deltaic sequences expand to the north and exhibit wave reworking and longshore transport of sands, and (5) low sedimentation rates and high long-term sea level during the Maastrichtian resulted in the deposition of a sediment-starved glauconitic shelf. Our study illustrates the widely known variability of mixed-influence deltaic systems, but also documents the relative stability of deltaic facies systems on the 106-107 yr scale, with long periods of cyclically repeating systems tracts controlled by eustasy. Results from the Late Cretaceous further show that although eustasy provides the template for sequences globally, regional tectonics (rates of subsidence and accommodation), changes in sediment supply, proximity to sediment input, and flexural subsidence from depocenter loading determines the regional to local preservation and facies expression of sequences. Copyright ?? 2008, SEPM (Society for Sedimentary Geology).

Journal of Sedimentary Research↗
USGS2017

Spectral matching techniques (SMTs) and automated cropland classification algorithms (ACCAs) for mapping croplands of Australia using MODIS 250-m time-series (2000–2015) data

Mapping croplands, including fallow areas, are an important measure to determine the quantity of food that is produced, where they are produced, and when they are produced (e.g. seasonality). Furthermore, croplands are known as water guzzlers by consuming anywhere between 70% and 90% of all human water use globally. Given these facts and the increase in global population to nearly 10 billion by the year 2050, the need for routine, rapid, and automated cropland mapping year-after-year and/or season-after-season is of great importance. The overarching goal of this study was to generate standard and routine cropland products, year-after-year, over very large areas through the use of two novel methods: (a) quantitative spectral matching techniques (QSMTs) applied at continental level and (b) rule-based Automated Cropland Classification Algorithm (ACCA) with the ability to hind-cast, now-cast, and future-cast. Australia was chosen for the study given its extensive croplands, rich history of agriculture, and yet nonexistent routine yearly generated cropland products using multi-temporal remote sensing. This research produced three distinct cropland products using Moderate Resolution Imaging Spectroradiometer (MODIS) 250-m normalized difference vegetation index 16-day composite time-series data for 16 years: 2000 through 2015. The products consisted of: (1) cropland extent/areas versus cropland fallow areas, (2) irrigated versus rainfed croplands, and (3) cropping intensities: single, double, and continuous cropping. An accurate reference cropland product (RCP) for the year 2014 (RCP2014) produced using QSMT was used as a knowledge base to train and develop the ACCA algorithm that was then applied to the MODIS time-series data for the years 2000–2015. A comparison between the ACCA-derived cropland products (ACPs) for the year 2014 (ACP2014) versus RCP2014 provided an overall agreement of 89.4% (kappa = 0.814) with six classes: (a) producer’s accuracies varying between 72% and 90% and (b) user’s accuracies varying between 79% and 90%. ACPs for the individual years 2000–2013 and 2015 (ACP2000–ACP2013, ACP2015) showed very strong similarities with several other studies. The extent and vigor of the Australian croplands versus cropland fallows were accurately captured by the ACCA algorithm for the years 2000–2015, thus highlighting the value of the study in food security analysis. The ACCA algorithm and the cropland products are released through http://croplands.org/app/map and http://geography.wr.usgs.gov/science/croplands/algorithms/australia_250m.html

International Journal of Digital Earth↗
USGS2004

Shear wave velocity, seismic attenuation, and thermal structure of the continental upper mantle

Seismic velocity and attenuation anomalies in the mantle are commonly interpreted in terms of temperature variations on the basis of laboratory studies of elastic and anelastic properties of rocks. In order to evaluate the relative contributions of thermal and non-thermal effects on anomalies of attenuation of seismic shear waves, Q −1 s , and seismic velocity, V s , we compare global maps of the thermal structure of the continental upper mantle with global Q −1 s and V s maps as determined from Rayleigh waves at periods between 40 and 150 s. We limit the comparison to three continental mantle depths (50, 100 and 150 km), where model resolution is relatively high. The available data set does not indicate that, at a global scale, seismic anomalies in the upper mantle are controlled solely by temperature variations. Continental maps have correlation coefficients of <0.56 between V s and T and of <0.47 between Q s and T at any depth. Such low correlation coefficients can partially be attributed to modelling artefacts; however, they also suggest that not all of the V s and Q s anomalies in the continental upper mantle can be explained by T variations. Global maps show that, by the sign of the anomaly, V s and Q s usually inversely correlate with lithospheric temperatures: most cratonic regions show high V s and Q s and low T , while most active regions have seismic and thermal anomalies of the opposite sign. The strongest inverse correlation is found at a depth of 100 km, where the attenuation model is best resolved. Significantly, at this depth, the contours of near-zero Q s anomalies approximately correspond to the 1000 °C isotherm, in agreement with laboratory measurements that show a pronounced increase in seismic attenuation in upper mantle rocks at 1000–1100 °C. East–west profiles of V s , Q s and T where continental data coverage is best (50°N latitude for North America and 60°N latitude for Eurasia) further demonstrate that temperature plays a dominant, but non-unique, role in determining the value of lithospheric V s and Q s . At 100 km depth, where the resolution of seismic models is the highest, we compare observed seismic V s and Q s with theoretical V T s and Q T s values, respectively, that are calculated solely from temperature anomalies and constrained by experimental data on temperature dependencies of velocity and attenuation. This comparison shows that temperature variations alone are sufficient to explain seismic V s and Q s in ca 50 per cent of continental regions. We hypothesize that compositional anomalies resulting from Fe depletion can explain the misfit between seismic and theoretical V s in cratonic lithosphere. In regions of active tectonics, temperature effects alone cannot explain seismic V s and Q s in the lithosphere. It is likely that partial melts and/or fluids may affect seismic parameters in these regions. This study demonstrates that lithospheric temperature plays the dominant role in controlling V s and Q s anomalies, but other physical parameters, such as compositional variations, fluids, partial melting and scattering, may also play a significant role in determining V s and Q s variations in the continental mantle.

Geophysical Journal International↗
USGS2000

Style and age of late Oligocene-early Miocene deformation in the southern Stillwater Range, west central Nevada: Paleomagnetism, geochronology, and field relations

Paleomagnetic and geochronologic data combined with geologic mapping tightly restrict the timing and character of a late Oligocene to early Miocene episode of large magnitude extension in the southern Stillwater Range and adjacent regions of west central Nevada. The southern Stillwater Range was the site of an Oligocene to early Miocene volcanic center comprising (1) 28.3 to 24.3 Ma intracaldera ash flow tuffs, lava flows, and subjacent plutons associated with three calderas, (2) 24.8 to 20.7 Ma postcaldera silicic dikes and domes, and (3) unconformably overlying 15.3 to 13.0 Ma dacite to basalt lava flows, plugs, and dikes. The caldera-related tuffs, lava flows, and plutons were tilted 60°-70° either west or east during the initial period of Cenozoic deformation that accommodated over 100% extension. Directions of remanent magnetization obtained from these extrusive and intrusive, caldera-related rocks are strongly deflected from an expected Miocene direction in senses appropriate for their tilt. A mean direction for these rocks after tilt correction, however, suggests that they were also affected by a moderate (33.4° ± 11.8°) component of counterclockwise vertical axis rotation. Paleomagnetic data indicate that the episode of large tilting occurred during emplacement of 24.8 to 20.7 Ma postcaldera dikes and domes. In detail, an apparent decrease in rotation with decreasing age of individual, isotopically dated bodies of the postcaldera group indicates that most tilting occurred between 24.4 and 24.2 Ma. The onset of tilting immediately following after the final caldera eruptions suggests that the magmatism and deformation were linked. Deformation was not driven by magma buoyancy, however, because tilting equally affected the caldera systems of different ages, including their plutonic roots. It is more likely that regional extension was focused in the southern Stillwater Range due to magmatic warming and reduction of tensile strength of the brittle crust. Faults that accommodated deformation in the southern Stillwater Range initially dipped steeply and cut deeply to expose more than 9 km of crustal section. The exposed crustal sections are probably rotated blocks above an unexposed basal detachment that lay near the early Miocene brittle-ductile transition.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Strain accumulation rates in the San Francisco Bay area, 1972–1989

Maps of the strain accumulation rate in the San Francisco Bay area have been constructed from trilateration observations extending from about 1972 until the Loma Prieta earthquake in late 1989. The observations were corrected to remove offsets imposed by shallow fault creep and by four M ∼6 earthquakes that occurred in the Bay area during that time interval. The Bay area was divided into 32 contiguous polygons, and the uniform (in both space and time) strain rates that best explain the changes in the corrected (earthquake and shallow fault creep offsets removed) distances within each polygon were calculated. In a coordinate system with the 1 axis directed N58°E and the 2 axis N32°W (perpendicular and parallel to the local tangent to the small circle drawn about the Pacific-Sierra Nevada pole of rotation) the averages of these 32 strain rates (each weighted by the area of the polygon) are ε11=9.2±7.4, ε12=−160.7±4.6, and ε22=8.2±6.2 nanostrain/yr, where extension is reckoned positive and quoted uncertainties are standard deviations. As expected from the Pacific-Sierra Nevada relative plate motion, the overall strain rate is predominantly right-lateral shear across a vertical plane striking N32°W. The net increase in the 12,225 km2 area of the trilateration network is only 212±110 m2/yr, which arises from almost equal extensions in the N32°W and N58°E directions. Within the network the strain rates vary from polygon to polygon. Significant areal dilatation rates are observed in almost ⅓ of the individual polygons and the N32°W extension rates tend to be negative to the west of the Hayward-Rodgers Creek fault trend and positive east of it. The N58E extension rate is positive in 22 out of 32 polygons, a proportion that is significantly larger than would be expected by chance if the N58°E extension rate were zero or negative. The pre-1989 strain accumulation across the eventual site of the Loma Prieta rupture involves fault normal contraction as well as right-lateral shear, consistent with the rupture mechanism.

California↗
USGS2017

Confirmation of the southwest continuation of the Cat Square terrane, southern Appalachian Inner Piedmont, with implications for middle Paleozoic collisional orogenesis

Detailed geologic mapping, U-Pb zircon geochronology and whole-rock geochemical analyses were conducted to test the hypothesis that the southwestern extent of the Cat Square terrane continues from the northern Inner Piedmont (western Carolinas) into central Georgia. Geologic mapping revealed the Jackson Lake fault, a ∼15 m-thick, steeply dipping sillimanite-grade fault zone that truncates lithologically distinct granitoids and metasedimentary units, and roughly corresponds with a prominent aeromagnetic lineament hypothesized to represent the southern continuation of the terrane-bounding Brindle Creek fault. Results of U-Pb SHRIMP geochronology indicate Late Ordovician to Silurian granitoids (444–439 Ma) occur exclusively northwest of the fault, whereas Devonian (404–371 Ma) granitoids only occur southeast of the fault. The relatively undeformed Indian Springs granodiorite (three individual bodies dated 317–298 Ma) crosscuts the fault and occurs on both sides, which indicates the Jackson Lake fault is a pre-Alleghanian structure. However, detrital zircon signatures from samples southeast of the Jackson Lake fault reveal dominant Grenville provenance, in contrast to Cat Square terrane detrital zircon samples from the northern Inner Piedmont, which include peri-Gondwanan (600–500 Ma) and a prominent Ordovician-Silurian (∼430 Ma) signature. We interpret the rocks southeast of the Jackson Lake fault to represent the southwestern extension of the Cat Square terrane primarily based on the partitioning of granitoid ages and lithologic distinctions similar to the northern Inner Piedmont. Data suggest Cat Square terrane metasedimentary rocks were initially deposited in a remnant ocean basin setting and developed into an accretionary prism in front of the approaching Carolina superterrane, ultimately overridden by it in Late Devonian to Early Mississippian time. Burial to >20 km resulted in migmatization of lower plate rocks, forming an infrastructure beneath the Carolina superterrane suprastructure. Provenance patterns support ∼250 km of Devonian dextral translation of the composite Inner Piedmont, which places the northern portion of the Inner Piedmont adjacent to a suite of ∼430 Ma plutons in the Virginia Blue Ridge during deposition. The megascopic thrust-nappe structural style of the northern Inner Piedmont, combined with southwest-directed lateral extrusion at mid-crustal depths, may reconcile differences in timing of metamorphism between the Carolina and central Georgia Inner Piedmont and structural contrasts between the Brindle Creek and Jackson Lake faults.

American Journal of Science↗
USGS2003

The crustal structure from the Altai Mountains to the Altyn Tagh fault, northwest China

[1] We present a new crustal section across northwest China based on a seismic refraction profile and geologic mapping. The 1100‐km‐long section crosses the southern margin of the Chinese Altai Mountains, Junggar Accretional Belt and eastern Junggar basin, easternmost Tianshan Mountains, and easternmost Tarim basin. The crustal velocity structure and Poisson's ratio (σ), which provide a constraint on crustal composition, were determined from P and S wave data. Despite the complex geology, the crustal thickness along the entire profile is nearly uniform at 50 km. The thickest crust (56 km) occurs at the northern end of the profile beneath the Altai Mountains and the thinnest (46 km) crust is beneath the Junggar basin. Beneath surficial sediments, the crust is found to have three layers with P wave velocities ( V p ) of 6.0–6.3, 6.3–6.6, and 6.9–7.0 km/s, respectively. The southern half of the profile, including the eastern Tianshan Mountains and eastern margin of the Tarim basin, shows low P wave velocities and σ = 0.25 to a depth of 30 km, which suggests a quartz‐rich, granitic upper crustal composition. The northern half of the profile below the Altai Mountains and Junggar Accretional Belt has a higher Poisson's ratio of σ = 0.26–0.27 to a depth of 30 km, indicative of an intermediate crustal composition. The entire 1100‐km‐long profile is underlain by a 15–30 km thick high velocity (6.9–7.0 km/s; σ = 0.26–0.28) lower‐crustal layer that we interpret to have a bulk composition of mafic granulite. At the southern end of the profile, a 5‐km‐thick midcrustal low‐velocity layer ( V p = 5.9 km/s, σ = 0.25) underlies the Tianshan and the region to the south, and may be indicative of a near‐horizontal detachment interface. P n velocities are ∼7.7–7.8 km/s between the Tianshan and the Junggar basin, and ∼7.9–8.0 km/s below the Altai Mountains and eastern margin of the Tarim basin. We interpret the consistent three‐layer stratification of the crust to indicate that the crust has undergone partial melting and differentiation after Paleozoic terrane accretion. The thickness (50 km) of the crust appears to be related to compression resulting from the Indo‐Asian collision.

Altai Mountains↗
USGS2005

Repeated surveys by acoustic Doppler current profiler for flow and sediment dynamics in a tidal river

A strategy of repeated surveys by acoustic Doppler current profiler (ADCP) was applied in a tidal river to map velocity vectors and suspended-sediment indicators. The Sacramento River at the junction with the Delta Cross Channel at Walnut Grove, California, was surveyed over several tidal cycles in the Fall of 2000 and 2001 with a vessel-mounted ADCP. Velocity profiles were recorded along flow-defining survey paths, with surveys repeated every 27 min through a diurnal tidal cycle. Velocity vectors along each survey path were interpolated to a three-dimensional Cartesian grid that conformed to local bathymetry. A separate array of vectors was interpolated onto a grid from each survey. By displaying interpolated vector grids sequentially with computer animation, flow dynamics of the reach could be studied in three-dimensions as flow responded to the tidal cycle. Velocity streamtraces in the grid showed the upwelling of flow from the bottom of the Sacramento River channel into the Delta Cross Channel. The sequential display of vector grids showed that water in the canal briefly returned into the Sacramento River after peak flood tides, which had not been known previously. In addition to velocity vectors, ADCP data were processed to derive channel bathymetry and a spatial indicator for suspended-sediment concentration. Individual beam distances to bed, recorded by the ADCP, were transformed to yield bathymetry accurate enough to resolve small bedforms within the study reach. While recording velocity, ADCPs also record the intensity of acoustic backscatter from particles suspended in the flow. Sequential surveys of backscatter intensity were interpolated to grids and animated to indicate the spatial movement of suspended sediment through the study reach. Calculation of backscatter flux through cross-sectional grids provided a first step for computation of suspended-sediment discharge, the second step being a calibrated relation between backscatter intensity and sediment concentration. Spatial analyses of ADCP data showed that a strategy of repeated surveys and flow-field interpolation has the potential to simplify computation of flow and sediment discharge through complex waterways. The use of trade, product, industry, or firm names in this report is for descriptive purposes only and does not constitute endorsement of products by the US Government. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2002

Plume-driven plumbing and crustal formation in Iceland

Through combination of surface wave and body wave constraints we derive a three-dimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT- SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0–15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (−7%) are found at 5–10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15–20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snæfellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last ∼20 Myr, the Snæfellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15–20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Journal of Geophysical Research B: Solid Earth↗
USGS1998

Paleoseismic investigations in the Santa Cruz mountains, California: Implications for recurrence of large-magnitude earthquakes on the San Andreas fault

Trenching, microgeomorphic mapping, and tree ring analysis provide information on timing of paleoearthquakes and behavior of the San Andreas fault in the Santa Cruz mountains. At the Grizzly Flat site alluvial units dated at 1640–1659 A.D., 1679–1894 A.D., 1668–1893 A.D., and the present ground surface are displaced by a single event. This was the 1906 surface rupture. Combined trench dates and tree ring analysis suggest that the penultimate event occurred in the mid-1600 s, possibly in an interval as narrow as 1632–1659 A.D. There is no direct evidence in the trenches for the 1838 or 1865 earthquakes, which have been proposed as occurring on this part of the fault zone. In a minimum time of about 340 years only one large surface faulting event (1906) occurred at Grizzly Flat, in contrast to previous recurrence estimates of 95–110 years for the Santa Cruz mountains segment. Comparison with dates of the penultimate San Andreas earthquake at sites north of San Francisco suggests that the San Andreas fault between Point Arena and the Santa Cruz mountains may have failed either as a sequence of closely timed earthquakes on adjacent segments or as a single long rupture similar in length to the 1906 rupture around the mid-1600 s. The 1906 coseismic geodetic slip and the late Holocene geologic slip rate on the San Francisco peninsula and southward are about 50–70% and 70% of their values north of San Francisco, respectively. The slip gradient along the 1906 rupture section of the San Andreas reflects partitioning of plate boundary slip onto the San Gregorio, Sargent, and other faults south of the Golden Gate. If a mid-1600 s event ruptured the same section of the fault that failed in 1906, it supports the concept that long strike-slip faults can contain master rupture segments that repeat in both length and slip distribution. Recognition of a persistent slip rate gradient along the northern San Andreas fault and the concept of a master segment remove the requirement that lower slip sections of large events such as 1906 must fill in on a periodic basis with smaller and more frequent earthquakes.

Journal of Geophysical Research B: Solid Earth↗
USGS2021

OpenET: Filling a critical data gap in water management for the western United States

The lack of consistent, accurate information on evapotranspiration (ET) and consumptive use of water by irrigated agriculture is one of the most important data gaps for water managers in the western United States (U.S.) and other arid agricultural regions globally. The ability to easily access information on ET is central to improving water budgets across the West, advancing the use of data-driven irrigation management strategies, and expanding incentive-driven conservation programs. Recent advances in remote sensing of ET have led to the development of multiple approaches for field-scale ET mapping that have been used for local and regional water resource management applications by U.S. state and federal agencies. The OpenET project is a community-driven effort that is building upon these advances to develop an operational system for generating and distributing ET data at a field scale using an ensemble of six well-established satellite-based approaches for mapping ET. Key objectives of OpenET include: Increasing access to remotely sensed ET data through a web-based data explorer and data services; supporting the use of ET data for a range of water resource management applications; and development of use cases and training resources for agricultural producers and water resource managers. Here we describe the OpenET framework, including the models used in the ensemble, the satellite, meteorological, and ancillary data inputs to the system, and the OpenET data visualization and access tools. We also summarize an extensive intercomparison and accuracy assessment conducted using ground measurements of ET from 139 flux tower sites instrumented with open path eddy covariance systems. Results calculated for 24 cropland sites from Phase I of the intercomparison and accuracy assessment demonstrate strong agreement between the satellite-driven ET models and the flux tower ET data. For the six models that have been evaluated to date (ALEXI/DisALEXI, eeMETRIC, geeSEBAL, PT-JPL, SIMS, and SSEBop) and the ensemble mean, the weighted average mean absolute error (MAE) values across all sites range from 13.6 to 21.6 mm/month at a monthly timestep, and 0.74 to 1.07 mm/day at a daily timestep. At seasonal time scales, for all but one of the models the weighted mean total ET is within ±8% of both the ensemble mean and the weighted mean total ET calculated from the flux tower data. Overall, the ensemble mean performs as well as any individual model across nearly all accuracy statistics for croplands, though some individual models may perform better for specific sites and regions. We conclude with three brief use cases to illustrate current applications and benefits of increased access to ET data, and discuss key lessons learned from the development of OpenET.

western United States↗
USGS2008

Reflectance spectroscopy for the assessment of soil salt content in soils of the yellow river delta of China

There has been growing interest in the use of reflectance spectroscopy as a rapid and inexpensive tool for soil characterization. In this study, we collected 95 soil samples from the Yellow River Delta of China to investigate the level of soil salinity in relation to soil spectra. Sample plots were selected based on a field investigation and the corresponding soil salinity classification map to maximize variations of saline characteristics in the soil. Spectral reflectances of air-dried soil samples were measured using an Analytical Spectral Device (ASD) spectrometer (350-2500 nm) with an artificial light source. In the Yellow River Delta, the dominant chemical in the saline soil was NaCl and MgCl2. Soil spectra were analysed using two-thirds of the available samples, with the remaining one-third withheld for validation purposes. The analysis indicated that with some preprocessing, the reflectance at 1931-2123 nm and 2153-2254 nm was highly correlated with soil salt content (SSC). In the spectral region of 1931-2123 nm, the correlation R ranged from -0.80 to -0.87. In the region of 2153-2254 nm, the SSC was positively correlated with preprocessed reflectance (0.79-0.88). The preprocessing was done by fitting a convex hull to the reflectance curve and dividing the spectral reflectance by the value of the corresponding convex hull band by band. This process is called continuum removal, and the resulting ratio is called continuum removed reflectance (CR reflectance). However, the SSC did not have a high correlation with the unprocessed reflectance, and the correlation was always negative in the entire spectrum (350-2500 nm) with the strongest negative correlation at 1981 nm (R = -0.63). Moreover, we found a strong correlation (R=0.91) between a soil salinity index (SSI: Constructed using CR reflectance at 2052 nm and 2203 nm) and SSC. We estimated SSC as a function of SSI and SSI' (SSI': Constructed using unprocessed reflectance at 2052 nm and 2203 nm) using univariate regression. Validation of the estimation of SSC was conducted by comparing the estimated SSC with the holdout sample points. The comparison produced an estimated root mean squared error (RMSE) of 0.986 (SSC ranging from 0.06 to 12.30 g kg-1) and R2 of 0.873 for SSC with SSI as independent variable and RMSE of 1.248 and R2 of 0.8 for SSC with SSI' as independent variable. This study showed that a soil salinity index developed for CR reflectance at 2052 nm and 2203 nm on the basis of spectral absorption features of saline soil can be used as a quick and inexpensive method for soil salt-content estimation.

International Journal of Remote Sensing↗
USGS1995

Seafloor environments within the Boston Harbor- Massachusetts Bay sedimentary system: A regional synthesis

Modern seafloor sedimentary environments within the glaciated, topographically complex Boston Harbor and Massachusetts Bay area have been interpreted and mapped from an extensive collection of sidescan sonar records and supplemental marine geologic data. Three categories of environments are present that reflect the dominant long-term processes of erosion or nondeposition, deposition, and sediment reworking. (1) Environments of erosion or nondeposition comprise exposures of bedrock, glacial drift, coarse lag deposits, and possibly coastal plain rocks that contain sediments (where present) ranging from boulder fields to gravelly sands and occur in areas of relatively strong currents. (2) Environments of deposition contain fine-grained sediments ranging from muddy sands to muds that have accumulated in areas of predominantly weak bottom currents. (3) Environments of sediment reworking contain patches with textures ranging from sandy gravels to muds that have been produced by a combination of erosion and deposition in areas with variable bottom currents. The distribution of sedimentary environments across the Boston Harbor-Massachusetts Bay area is extremely patchy. Locally, this patchiness is due either to modifications of bottom-current strength (caused by the irregular topography and differences in water depth) or to small-scale changes in the supply of fine-grained sediments. Regional patchiness, however, reflects differences in geologic and oceanographic conditions among the estuarine, inner shelf, and basinal parts of the sedimentary system. The estuarine part of the system (Boston Harbor) is a depositional trap for fine-grained sediments because it is protected from large waves, has generally weak and variable tidal currents, and receives a large supply of fine grained detritus from natural and anthropogenic sources. The inner shelf, on the other hand, is largely an area of erosion or nondeposition due to sediment removal and redistribution during past sea-level changes, to sediment resuspension and winnowing by modern waves and currents, and to an inadequate supply of fine-grained sediments. The basinal part of the system (Stellwagen Basin) is mainly a tranquil depositional environment in which fine-grained sediments from several potential sources settle through the water column and accumulate under weak bottom currents. This study indicates areas within the Boston Harbor-Massachusetts Bay sedimentary system where fine-grained sediments and associated contaminants are likely to be either moved or deposited. It also provides a guide to the locations and variability of benthic habitats.

Massachusetts Bay↗
USGS2017

Plateaus and sinuous ridges as the fingerprints of lava flow inflation in the Eastern Tharsis Plains of Mars

The Tharsis Montes rift aprons are composed of outpourings of lava from chaotic terrains to the northeast and southwest flank of each volcano. Sinuous and branching channel networks that are present on the rift aprons suggest the possibility of fluvial processes in their development, or erosion by rapidly emplaced lavas, but the style of lava flow emplacement throughout rift apron development is not clearly understood. To better characterize the style of lava emplacement and role of fluvial processes in rift apron development, we conducted morphological mapping of the Pavonis Mons southwest rift apron and the eastern Tharsis plains using images from the High Resolution Imaging Science Experiment (HiRISE), Mars Orbiter Camera (MOC), Context Camera (CTX), Thermal Emission Imaging System (THEMIS), and High Resolution Stereo Camera (HRSC) along with the Mars Orbiter Laser Altimeter (MOLA) Precision Experiment Data Records (PEDRs) and gridded data. Our approach was to: (1) search for depositional fans at the slope break between the rift apron and adjacent low slope plains; (2) determine if there is evidence that previously formed deposits might have been buried by plains units; (3) characterize the Tharsis plains morphologies east of Pavonis Mons; and (4) assess their relationship to the rift apron units. We have not identified topographically significant depositional fans, nor did we observe evidence to suggest that plains units have buried older rift apron units. Flow features associated with the rift apron are observed to continue across the slope break onto the plains. In this area, the plains are composed of a variety of small fissures and low shield vents around which broad channel-fed and tube-fed flows have been identified. We also find broad, flat-topped plateaus and sinuous ridges mixed among the channels, tubes and vents. Flat-topped plateaus and sinuous ridges are morphologies that are analogous to those observed on the coastal plain of Hawai‘i, where lava flows have advanced from the volcano's several degree flank onto the nearly zero degree coastal plain. When local volumetric flow rates are low, flow fronts tend to spread laterally and often thicken via endogenous growth, or inflation, of the sheet-like flow units. If flow advance is restricted by existing topography into narrow pathways, inflation can be focused into sinuous, elongate ridges. The presence of plateaus and ridges—emplaced from the rift zones, across the plains to the east of Pavonis Mons—and a lack of fan-like features, or evidence for their burial, are consistent with rift apron lavas crossing a slope break with low local volumetric flow rates that led to inflation of sheet-like and tube-fed lava flows.

Journal of Volcanology and Geothermal Research↗
USGS2012

Geology of the Mount St. Helens area: Record of discontinuous volcanic and plutonic activity in the Cascade Arc of southern Washington

The Quaternary edifice of Mount St. Helens volcano was built upon a deeply eroded terrane of gently folded and altered volcanic and plutonic rocks that represent the core of the Tertiary Cascade magmatic arc. These rocks constitute an east dipping homoclinal sequence, several kilometers thick, of subaerially erupted mafic to silicic flows and volcaniclastic strata; K‐Ar ages from this section range from about 28 to 23 Ma (late Oligocene and earliest Miocene), which corresponds to an apparent lull in Cascade volcanism to the north of Mount Rainier. Volcanism was essentially continuous during this period of time, and neither a well‐defined base nor top of the section is exposed within the mapped area. Basalt and basaltic andesite dominate the lower part of the mid‐Tertiary section, whereas andesitic and dacitic rocks comprise most of the upper part. This section was intruded by numerous mafic to silicic dikes, sills, and irregular plutonic bodies, most no more than a few million years younger than their host rocks, and subjected to pervasive burial metamorphism and widespread hydrothermal alteration. Large areas of hornfelsed rock surrounding even relatively small intrusions indicate that the proportion of plutonic rock becomes significantly greater at shallow depth beneath the existing erosion surface. A large granitic pluton intruded the mid‐Tertiary section north of Spirit Lake at about 21 Ma. The Earl porphyry copper deposit occurs within the pluton but appears too young (17 Ma) to be genetically related to it. In contrast to the rather continuous and voluminous Oligocene to early Miocene activity, volcanism since then in the Mount St. Helens area has been localized and Volumetrically minor. Products of three younger eruptive periods have been recognized: a sequence of 15 m.y. old pyroxene andesite flows resting unconformably on mid‐Tertiary strata south of Mount St. Helens, widespread shallow dikes and sills of pyroxene andesite between 10 and 8 m.y. old, and compositionally diverse rocks erupted during the past 3 m.y. The Quaternary lavas are more potassic than the Tertiary lavas and typically contain phenocrysts of hornblende and biotite, which are absent from the older rocks. A number of Tertiary structures define a broad NNE trending zone that may reflect a deep‐seated lithospheric flaw that has controlled the locus of Cascade magmatism in southern Washington for the past 25 m.y. Mount St. Helens lies within this zone at the intersection of the NNW striking St. Helens seismic zone (SHZ) and an ENE trending alignment of Pleistocene silicic plug‐domes. No surface breakage has been detected along the SHZ, which is apparently very young. The linear zone of silicic vents is probably controlled by a fault that has been interpreted from seismic records to occur directly beneath the volcano. This zone parallels the directions of regional maximum horizontal compressive stress and North America/Juan de Fuca plate convergence. Mount St. Helens is an example of a low‐volume tectonically controlled magmatic system in an early stage of development.

Washington↗