SEARCH · USGS Science

Results for “One Earth”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,189 records · Page 66Linked to original sources

USGS1979

Lunar magnetic anomalies detected by the Apollo substatellite magnetometers

Properties of lunar crustal magnetization thus far deduced from Apollo subsatellite magnetometer data are reviewed using two of the most accurate presently available magnetic anomaly maps - one covering a portion of the lunar near side and the other a part of the far side. The largest single anomaly found within the region of coverage on the near-side map correlates exactly with a conspicuous, light-colored marking in western Oceanus Procellarum called Reiner Gamma. This feature is interpreted as an unusual deposit of ejecta from secondary craters of the large nearby primary impact crater Cavalerius. An age for Cavalerius (and, by implication, for Reiner Gamma) of 3.2 ?? 0.2 ?? 109 y is estimated. The main (30 ?? 60 km) Reiner Gamma deposit is nearly uniformly magnetized in a single direction, with a minimum mean magnetization intensity of ???7 ?? 10-2 G cm3/g (assuming a density of 3 g/cm3), or about 700 times the stable magnetization component of the most magnetic returned samples. Additional medium-amplitude anomalies exist over the Fra Mauro Formation (Imbrium basin ejecta emplaced ???3.9 ?? 109 y ago) where it has not been flooded by mare basalt flows, but are nearly absent over the maria and over the craters Copernicus, Kepler, and Reiner and their encircling ejecta mantles. The mean altitude of the far-side anomaly gap is much higher than that of the near-side map and the surface geology is more complex, so individual anomaly sources have not yet been identified. However, it is clear that a concentration of especially strong sources exists in the vicinity of the craters Van de Graaff and Aitken. Numerical modeling of the associated fields reveals that the source locations do not correspond with the larger primary impact craters of the region and, by analogy with Reiner Gamma, may be less conspicuous secondary crater ejecta deposits. The reason for a special concentration of strong sources in the Van de Graaff-Aitken region is unknown, but may be indirectly related to the existence of strongly modified crustal terrain which also occurs in the same region. The inferred directions of magnetization for the several sources of the largest anomalies are highly inclined with respect to one another, but are generally depleted in the north-south direction. The north-south depletion of magnetization intensity appears to continue across the far-side within the region of coverage. The mechanism of magnetization and the origin of the magnetizing field remain unresolved, but the uniformity with which the Reiner Gamma deposit is apparently magnetized, and the north-south depletion of magnetization intensity across a substantial portion of the far side, seem to require the existence of an ambient field, perhaps of global or larger extent. The very different inferred directions of magnetization possessed by nearly adjacent sources of the Van de Graaff-Aitken anomalies, and the depletion in their north-south component of magnetization, do not favor an internally generated dipolar field oriented parallel to the present spin axis. A variably oriented interplanetary magnetizing field that was intrinsically strong or locally amplified by unknown surface processes is least inconsistent with the data. ?? 1979.

Physics of the Earth and Planetary Interiors↗

USGS2020

Spatial and temporal variations in SO2 and PM2.5 levels around Kīlauea volcano, Hawai'i during 2007–2018

Among the hazards posed by volcanoes are the emissions of gases and particles that can affect air quality and damage agriculture and infrastructure. A recent intense episode of volcanic degassing associated with severe impacts on air quality accompanied the 2018 lower East Rift Zone (LERZ) eruption of Kīlauea volcano, Hawai'i. This resulted in a major increase in gas emission rates with respect to usual emission values for this volcano, along with a shift in the source of the dominant plume to a populated area on the lower flank of the volcano. This led to reduced air quality in downwind communities. We analyse open-access data from the permanent air quality monitoring networks operated by the Hawai'i Department of Health (HDOH) and National Park Service (NPS), and report on measurements of atmospheric sulfur dioxide (SO 2 ) between 2007 and 2018 and PM 2.5 (aerosol particulate matter with diameter <2.5 μm) between 2010 and 2018. Additional air quality data were collected through a community-operated network of low-cost PM 2.5 sensors during the 2018 LERZ eruption. From 2007 to 2018 the two most significant escalations in Kīlauea's volcanic emissions were: the summit eruption that began in 2008 (Kīlauea emissions averaged 5–6 kt/day SO 2 from 2008 until summit activity decreased in May 2018) and the LERZ eruption in 2018 when SO 2 emission rates reached a monthly average of 200 kt/day during June. In this paper we focus on characterizing the airborne pollutants arising from the 2018 LERZ eruption and the spatial distribution and severity of volcanic air pollution events across the Island of Hawai'i. The LERZ eruption caused the most frequent and severe exceedances of the Environmental Protection Agency (EPA) PM 2.5 air quality threshold (35 μg/m 3 as a daily average) in Hawai'i in the period 2010–2018. In Kona, for example, the maximum 24-h-mean mass concentration of PM 2.5 was recorded as 59 μg/m 3 on the twenty-ninth of May 2018, which was one of eight recorded exceedances of the EPA air quality threshold during the 2018 LERZ eruption, where there had been no exceedances in the previous 8 years as measured by the HDOH and NPS networks. SO 2 air pollution during the LERZ eruption was most severe in communities in the south and west of the island, as measured by selected HDOH and NPS stations in this study, with a maximum 24-h-mean mass concentration of 728 μg/m 3 recorded in Ocean View (100 km west of the LERZ emission source) in May 2018. Data from the low-cost sensor network correlated well with data from the HDOH PM 2.5 instruments, confirming that these low-cost sensors provide a robust means to augment reference-grade instrument networks.

Hawaii↗

USGS2023

ENSO and NAO linkages to interannual salinity variability in north central Gulf of Mexico estuaries through teleconnections with precipitation

Though the importance of Earth's internal climate modes such as the El Niño-Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO) to regional-scale climate variability is well recognized, the degree to which these oscillations are reflected by spatio-temporal salinity variability over interannual timescales in estuaries is less understood. Here an 11-year continuous salinity monitoring dataset spanning 223 stations across Louisiana's coastal wetlands along the northern Gulf of Mexico is examined with empirical orthogonal function (EOF) analysis to identify dominant modes of interannual variability in the salinity field. The first EOF mode accounts for 72% of the variance in the salinity field and captures a domain-wide pattern where salinities vary in-phase through space in response to local precipitation anomalies occurring in the vicinity of the study area. This local precipitation anomaly is positively correlated with ENSO (Nino3.4 index), consistent with the El Niño – wet (La Niña – dry) precipitation teleconnection that is prevalent throughout the northern Gulf of Mexico coast. The second EOF mode, which accounts for 13% of the variance in the salinity field, is expressed primarily in the marshes across the lower reaches of the Mississippi River deltaic plain (MRDP). EOF2 is anticorrelated with annual Mississippi River discharge anomaly such that salinities in the lower MRDP decrease as discharge increases, pointing to enhanced advection of fresh river plume waters over the shelf into the estuary via estuary-ocean exchange during years of anomalously high river discharge. Mississippi River discharge anomaly is positively correlated with the NAO at a one-year time lag, through a teleconnection with precipitation throughout much of the central region of the Mississippi River drainage basin. Together, these findings indicate that most of the interannual salinity variability across Louisiana's coastal wetlands can be linked to climate variability through teleconnections with precipitation. Incorporating these dynamics into restoration planning, monitoring, and adaptive management efforts may help constrain background environmental variation and better isolate restoration effects.

Louisiana↗

USGS2023

Debris-flow process controls on steepland morphology in the San Gabriel Mountains, California

Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris-flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady-state fluvial channel longitudinal profiles are concave-up and exhibit a power-law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope-area signature in log-log space. Here, we use a one-dimensional landform evolution model that incorporates debris-flow erosion to reproduce the relationship between this curved slope-area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris-flow erosion, valley-bottom slope, and debris-flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris-flow incision law. Finally, modeled fluvial incision outpaces debris-flow erosion at drainage areas less than those at which valleys morphologically transition from near-invariant slopes to concave profiles. This result emphasizes the critical role of debris-flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

California↗

USGS2003

Paleomagnetism and geochronology of an Early Proterozoic quartz diorite in the southern Wind River Range, Wyoming, USA

We present geochronologic and paleomagnetic data from a north-trending quartz diorite intrusion that cuts Archean metasedimentary and metaigneous rocks of the South Pass Greenstone Belt of the Wyoming craton. The quartz diorite was previously thought to be either Archean or Early Proterozoic (?) in age and is cut by north and northeast-trending Proterozoic diabase dikes of uncertain age, for which we also report paleomagnetic data. New U-Pb analyses of baddeleyite and zircon from the quartz diorite yield a concordia upper intercept age of 2170 ?? 8 Ma (95% confidence). An 40Ar/39Ar amphibole date from the same sample yields a similar apparent age of about 2124 ?? 30 Ma (2??), thus confirming that the intrusion is Early Proterozoic in age and that it has probably not been thermally disturbed since emplacement. A magmatic event at ca. 2.17 Ga has not previously been documented in the Wyoming craton. The quartz diorite and one of the crosscutting diabase dikes yield essentially identical, well-defined characteristic remanent magnetizations. Results from eight sites in the quartz diorite yield an in situ mean direction of north declination and moderate to steep positive inclination (Dec.=355??, Inc.=65??, k=145, ??95=5??) with a paleomagnetic pole at 84??N, 215??E (??m=6??, ??p=7??). Data from other diabase dike sites are inconsistent with the quartz diorite results, but the importance of these results is uncertain because the age of the dikes is not well known. Interpretation of the quartz diorite remanent magnetization is problematic. The in situ direction is similar to expected directions for magnetizations of Late Cretaceous/early Tertiary age. However, there is no compelling evidence to suggest that these rocks were remagnetized during the late Mesozoic or Cenozoic. Assuming this magnetization to be primary, then the in situ paleomagnetic pole is strongly discordant with poles of 2167, 2214, and 2217 Ma from the Canadian Shield, and is consistent with proposed separation of the Wyoming Craton and Laurentia prior to about 1.8 Ga. Correcting the quartz diorite pole for the possible effects of Laramide-age tilting of the Wind River Range, based on the attitude of nearby overlying Cambrian Flathead Sandstone (dip=20??, N20??E), gives a tilt corrected pole of 75??N, 58??E (??m=4??, ??p=6??), which is also discordant with respect to time-equivalent poles from the Superior Province. Reconstruction of the Superior and Wyoming Province using a rotation similar to that proposed by Roscoe and Card [Can. J. Earth Sci. 46(1993)2475] is problematic, but reconstruction of the Superior and Wyoming Provinces based on restoring them to their correct paleolatitude and orientation using a closest approach fit indicates that the two cratons could have been adjacent at about 2.17 Ga prior to rifting at about 2.15 Ga. The paleomagnetic data presented are consistent with the hypothesis that the Huronian and Snowy Pass Supergroups could have evolved as part of a single epicratonic sedimentary basin during the Early Proterozoic. ?? 2002 Elsevier Science B.V. All rights reserved.

Tectonophysics↗

USGS2017

Delayed seismicity rate changes controlled by static stress transfer

On 15 June 2010, a M w 5.7 earthquake occurred near Ocotillo, California, in the Yuha Desert. This event was the largest aftershock of the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake in this region. The EMC mainshock and subsequent Ocotillo aftershock provide an opportunity to test the Coulomb failure hypothesis (CFS). We explore the spatiotemporal correlation between seismicity rate changes and regions of positive and negative CFS change imparted by the Ocotillo event. Based on simple CFS calculations we divide the Yuha Desert into three subregions, one triggering zone and two stress shadow zones. We find the nominal triggering zone displays immediate triggering, one stress shadowed region experiences immediate quiescence, and the other nominal stress shadow undergoes an immediate rate increase followed by a delayed shutdown. We quantitatively model the spatiotemporal variation of earthquake rates by combining calculations of CFS change with the rate-state earthquake rate formulation of Dieterich (1994), assuming that each subregion contains a mixture of nucleation sources that experienced a CFS change of differing signs. Our modeling reproduces the observations, including the observed delay in the stress shadow effect in the third region following the Ocotillo aftershock. The delayed shadow effect occurs because of intrinsic differences in the amplitude of the rate response to positive and negative stress changes and the time constants for return to background rates for the two populations. We find that rate-state models of time-dependent earthquake rates are in good agreement with the observed rates and thus explain the complex spatiotemporal patterns of seismicity.

California↗

USGS2009

Environmental influences on speleothem growth in southwestern Oregon during the last 380, 000 years

The growth of carbonate formations in caves (speleothems) is sensitive to changes in environmental conditions at the surface (temperature, precipitation and vegetation) and can provide useful paleoclimatic and paleoenvironmental information. We use 73 230 Th dates from speleothems collected from a cave in southwestern Oregon (USA) to constrain speleothem growth for the past 380 000 years. Most speleothem growth occurred during interglacial periods, whereas little growth occurred during glacial intervals. To evaluate potential environmental controls on speleothem growth we use two new modeling approaches: i) a one-dimensional thermal advection–diffusion model to estimate cave temperatures during the last glacial cycle, and ii) a regional climate model simulation for the Last Glacial Maximum (21 000 years before present) that assesses a range of potential controls on speleothem growth under peak glacial conditions. The two models are mutually consistent in indicating that permafrost formation did not influence speleothem growth during glacial periods. Instead, the regional climate model simulation combined with proxy data suggest that the influence of the Laurentide and Cordilleran ice sheets on atmospheric circulation induced substantial changes in water balance in the Pacific Northwest and affected speleothem growth at our location. The overall drier conditions during glacial intervals and associated periods of frozen topsoil at times of maximum surface runoff likely induced drastic changes in cave recharge and limited speleothem growth. This mechanism could have affected speleothem growth in other mid-latitude caves without requiring the presence of permafrost.

Oregon↗

USGS1997

Morphology and distribution of seamounts surrounding Easter Island

We investigate the morphology and distribution of a seamount population on a section of seafloor influenced by both superfast seafloor spreading and hotspot volcanism. The population under investigation is part of a broad chain of seamounts extending eastward from the East Pacific Rise, near Easter Island. In order to define the morphological variability of the seamounts, basal shape, cross-sectional area, volume, flatness, and flank slope are plotted against height for 383 seamounts with heights greater than 200 m, based on bathymetry data collected by GLORI-B and SeaBeam 2000, during three cruises onboard the R/V Melville in the spring of 1993. Nearly complete swath mapping coverage of the seamounts is available for the analysis of size and shape distribution. We quantitatively describe the seamount population of this active region, in which seamounts cover ???27% of the seafloor, and account for ???4.2% of the total crustal volume. Over 50% of the total volume (61,000 km3) of seamounts used in this study is made up by the 14 largest seamounts, and the remaining volume is made up by the 369 smaller seamounts (>200 m in height). Our analysis indicates there are at least two seamount populations in the Easter Island-Salas y Gomez Island (25??-29??S, 113??-104??W) study area. One population of seamounts is composed of short seamounts (<1200 m in height) with variable flatness from pointy cones to flattened domes (flatness from 0.01 to 0.57) and predominantly steep flanks (slopes from 5?? to 32??). A second population is of massive (>1200 m), shield-like, pointy cones (flatness < 0.2) and gentle slopes (from 5?? to 15??). An exponential maximum likelihood distribution is fit to the binned raw frequency of height and gives a characteristic height of the seamount population of 308??12m and an expected number of seamounts per 1000 km2 of 2.7??0.15. Many seamounts that have different slope and flatness relationships with height are formed next to each other. We speculate that the larger volcanoes (>???1200 m) originate exclusively from a hotspot source, but only a portion of the smaller volcanoes (<???1200 m) are formed from a hotspot source. The remainder would be presumably formed by a normal mantle or mixed source.

Journal of Geophysical Research B: Solid Earth↗

USGS1992

Plutonism in the central part of the Sierra Nevada Batholith, California

The Sierra Nevada batholith comprises the plutonic rocks of Mesozoic age that underlie most of the Sierra Nevada, a magnificent mountain range that originated in the Cenozoic by the westward tilting of a huge block of the Earth's crust. Scattered intrusions west of the batholith in the western metamorphic belt of the Sierra Nevada and east of the Sierra Nevada in the Benton Range and the White and Inyo Mountains are satellitic to but not strictly parts of the Sierra Nevada batholith. Nevertheless, all the plutonic rocks are related in origin. The batholith lies along the west edge of the Paleozoic North American craton, and Paleozoic and early Mesozoic oceanic crust underlies its western margin. It was emplaced in strongly deformed but weakly metamorphosed strata ranging in age from Proterozoic to Cretaceous. Sedimentary rocks of Proterozoic and Paleozoic age crop out east of the batholith in the White and Inyo Mountains, and metamorphosed sedimentary and volcanic rocks of Paleozoic and Mesozoic age crop out west of the batholith in the western metamorphic belt. A few large and many small, generally elongate remnants of metamorphic rocks lie within the batholith. Sparse fossils from metasedimentary rocks and isotopic ages for metavolcanic rocks indicate that the metamorphic rocks in the remnants range in age from Early Cambrian to Early Cretaceous. Within the map area (the Mariposa 1 0 by 2 0 quadrangle), the bedding, cleavage, and axial surfaces of folds generally trend about N. 35 0 W., parallel to the long axis of the Sierra Nevada. The country rocks comprise strongly deformed but generally coherent sequences; however, some units in the western metamorphic belt may partly consist of melanges. Most sequences are in contact with other sequences, at least for short distances, but some sequences within the batholith are bounded on one or more sides by plutonic rocks. Proterozoic and Paleozoic sedimentary strata east of the Sierra Nevada and Paleozoic strata in remnants of country rocks within the eastern part of the batholith, although strongly deformed, are autochthonous or have been displaced only short distances, whereas some Mesozoic strata in the western metamorphic belt may be allochthonous. Probably the strata in the western metamorphic belt were deposited in marginal basins and island arcs, but the possibility that they were transported from distant places has not been disproved. All the country rocks have been strongly deformed, most of them more than once. Tectonic disturbances occurred during the Devonian and Mississippian (Antler? orogeny), the Permian and (or) Early Triassic (Sonoman? orogeny), the Late Jurassic (Nevadan orogeny), and at various other times during emplacement of the batholith and uplift that accompanied and followed its emplacement. The strata in the western metamorphic belt probably were deformed in an early Mesozoic subduction complex. The plutonic rocks range in composition from gabbro to leucogranite, but tonalite, granodiorite, and granite are the most common rock types. Most are medium to coarse grained, but some small rock masses are fine grained. Most have hypidiomorphic-granular textures and are equigranular, but some having compositions close to the boundary between granite and granodiorite contain large crystals of alkali feldspar. Serpentinized ultramafic rocks are present locally in the western metamorphic belt within and adjacent to the Melones fault zone. Except for serpentinized ultramafic rocks, trondhjemite, and most granites, all the plutonic rocks contain significant amounts of hornblende. Most of the granitoids are metaluminous or weakly peraluminous; strongly peraluminous granites are present only in the White Mountains. Most of the granitoids are assigned to units of lithodemic rank, and most of these units are assigned to intrusive suites. Plutons assigned to the same lithodeme are composed of rock of similar composition, fabric, and age and are presumed to h

Professional Paper↗

USGS2023

A one-dimensional volcanic plume model for predicting ash aggregation

During explosive volcanic eruptions, volcanic ash is ejected into the atmosphere, impacting aircraft safety and downwind communities. These volcanic clouds tend to be dominated by fine ash (<63 μm in diameter), permitting transport over hundreds to thousands of kilometers. However, field observations show that much of this fine ash aggregates into clusters or pellets with faster settling velocities than individual particles. Models of ash transport and deposition require an understanding of aggregation processes, which depend on factors like moisture content and local particle collision rates. In this study, we develop a Plume Model for Aggregate Prediction, a one-dimensional (1D) volcanic plume model that predicts the plume rise height, concentration of water phases, and size distribution of resulting ash aggregates from a set of eruption source parameters. The plume model uses a control volume approach to solve mass, momentum, and energy equations along the direction of the plume axis. The aggregation equation is solved using a fixed pivot technique and incorporates a sticking efficiency model developed from analog laboratory experiments of particle aggregation within a novel turbulence tower. When applied to the 2009 eruption of Redoubt Volcano, Alaska, the 1D model predicts that the majority of the plume is over-saturated with water, leading to a high rate of aggregation. Although the mean grain size of the computed Redoubt aggregates is larger than the measured deposits, with a peak at 1 mm rather than 500 μm, the present results provide a quantitative estimate for the magnitude of aggregation in an eruption.

JGR Solid Earth↗

USGS2022

Possible effects of multiphase methane evolution during a glacial cycle on underpressure development in sedimentary basins: An analysis with application to the northeast Michigan Basin

In low-permeability systems, groundwater may be accompanied by separate-phase fluids, and measured pore water pressures may deviate from those expected in steady-state, single-phase systems. These same systems may be of interest for storage of nuclear waste in Deep Geologic Repositories. Therefore, it is important to understand the relationship between the presence of a separate phase and anomalous pressure development. At the Bruce site in Southern Ontario, a significant underpressure was observed, and there is evidence for the presence of gas-phase methane in situ. This study used a one-dimensional (vertical) numerical model of the subsurface down to a depth of 844 m beneath the Bruce site to evaluate possible effects of hydromechanical coupling with multiphase flow on pressure evolution during glacial loading and unloading. The simulated pressure conditions were affected strongly by the amount of methane initially present in the system, and the maximum simulated underpressure varied nonmonotonically with increasing initial methane content. When the initial methane content was below the solubility limit, exsolution led to underpressures that briefly exceeded those that formed in the single-phase case. At intermediate initial methane contents (sufficient to produce an immobile gas phase), the gas phase dampened the hydromechanical effects of the glacial cycle. At large initial methane contents (when a mobile gas phase was present), gas migration caused a large decrease in relative liquid permeability, which further contributed to underpressure development in the pore water. Multiple scenarios that spanned a range of initial methane contents yielded underpressures like those observed at the Bruce site.

Journal of Geophysical Research - Solid Earth↗

USGS2021

Extent of impact of deep-sea nodule mining midwater plumes is influenced by sediment loading, turbulence and thresholds

Deep-sea polymetallic nodule mining research activity has substantially increased in recent years, but the expected level of environmental impact is still being established. One environmental concern is the discharge of a sediment plume into the midwater column. We performed a dedicated field study using sediment from the Clarion Clipperton Fracture Zone. The plume was monitored and tracked using both established and novel instrumentation, including acoustic and turbulence measurements. Our field studies reveal that modeling can reliably predict the properties of a midwater plume in the vicinity of the discharge and that sediment aggregation effects are not significant. The plume model is used to drive a numerical simulation of a commercial-scale operation in the Clarion Clipperton Fracture Zone. Key takeaways are that the scale of impact of the plume is notably influenced by the values of environmentally acceptable threshold levels, the quantity of discharged sediment, and the turbulent diffusivity in the Clarion Clipperton Fracture Zone.

Communications Earth & Environment↗

USGS2011

USGS global change science strategy: A framework for understanding and responding to climate and land-use change

This U.S. Geological Survey (USGS) Global Change Science Strategy expands on the Climate Variability and Change science component of the USGS 2007 Science Strategy, “Facing Tomorrow’s Challenges: USGS Science in the Coming Decade” (U.S. Geological Survey, 2007). Here we embrace the broad definition of global change provided in the U.S. Global Change Research Act of 1990 (Public Law 101–606,104 Stat. 3096–3104)—“Changes in the global environment (including alterations in climate, land productivity, oceans or other water resources, atmospheric chemistry, and ecological systems) that may alter the capacity of the Earth to sustain life”—with a focus on climate and land-use change. There are three major characteristics of this science strategy. First, it addresses the science required to broadly inform global change policy, while emphasizing the needs of natural-resource managers and reflecting the role of the USGS as the science provider for the Department of the Interior and other resource-management agencies. Second, the strategy identifies core competencies, noting 10 critical capabilities and strengths the USGS uses to overcome key problem areas. We highlight those areas in which the USGS is a science leader, recognizing the strong partnerships and effective collaboration that are essential to address complex global environmental challenges. Third, it uses a query-based approach listing key research questions that need to be addressed to create an agenda for hypothesis-driven global change science organized under six strategic goals. Overall, the strategy starts from where we are, provides a vision for where we want to go, and then describes high-priority strategic actions, including outcomes, products, and partnerships that can get us there. Global change science is a well-defined research field with strong linkages to the ecosystems, water, energy and minerals, natural hazards, and environmental health components of the USGS Science Strategy (2007). When science strategies that cover these other components are developed, coordinated implementation will be necessary to achieve Bureau-level synergies and optimize capabilities and expertise. In October 2010, USGS realigned its management and budget structure to implement its 2007 Science Strategy. The new organizational structure, in which “Global Change” is one of seven key mission areas, lends itself to the advancement of the established six strategic goals. USGS global change science is formally represented by the “Climate and Land-Use Change” Mission Area in the FY 2012 budget (USGS, 2011). This plan was developed by the USGS Global Change Science Strategy Planning Team (SSPT) appointed by the USGS Director on March 4, 2010 and charged with developing a Global Change Science Strategy for the coming decade (McNutt, 2010). USGS managers and science staff are the main audience for this science strategy. This document is also intended to serve as the foundation for consistent USGS collaboration and communication with partners and stakeholders.

Open-File Report↗

USGS2000

Palynomorphs of Permian Gondwana coal from borehole GDH-38, Barapukuria Coal Basin, Bangladesh

Thirty-two core samples of Permian Gondwana coal from three coal beds of borehole GDH-38, Barapukuria Coal Basin, Dinajpur, the north-northwestern part of Bangladesh, have been collected for palynological analysis. All samples except one yielded palynomorphs and some samples contain well-preserved and abundant palynomorphs of the gymnospermal and cryptogamic groups that are considered to be useful for future correlation studies. The lower coal bed (331.6–372.5 m) can easily be differentiated from the upper two coal beds by the presence of Alisporites, Cordaitina, Corisaccites, Hamiapollenites, Leuckisporites, Nuskoisporites, Tumoripollenites, Vestgisporites and Vittatina . It is difficult to palynologically differentiate the middle (198.1–208 m) and upper (162.3–172.9 m) coal beds as they contain a very limited number of specimens by which they can be identified. The middle bed is distinguished by the presence of Microbaculispora and Weylandites and the upper bed by the presence of a single taxon Acanthotriletes . Some of the vesiculate or saccate taxa extracted from these coal beds are typical of those occurring in Permian strata of Gondwana in India, South Africa, South America, Russia, Australia and Antarctica. They are thought to be derived from Glossopteris flora, which is characterised by an abundance of Pteridospermic plants of the gymnosperm group.

Dinajpur↗

USGS2000

The enigma of the Arthur's Pass, New Zealand, earthquake: 1. Reconciling a variety of data for an unusual earthquake sequence

The 1994 Arthur's Pass earthquake ( M W 6.7) is the largest in a recent sequence of earthquakes in the central South Island, New Zealand. No surface rupture was observed, the aftershock distribution was complex, and routine methods of obtaining the faulting orientation of this earthquake proved contradictory. We use a range of data and techniques to obtain our preferred solution, which has a centroid depth of 5 km, M 0 =1.3 × 10 19 N m, and a strike, dip, and rake of 221°, 47°, 112°, respectively. Discrepancies between this solution and the Harvard centroid moment tensor, together with the Global Positioning System (GPS) observations and unusual aftershock distribution, suggest that the rupture may not have occurred on a planar fault. A second, strike slip, subevent on a more northerly striking plane is suggested by these data but neither the body wave modeling nor regional broadband recordings show any complexity or late subevents. We relocate the aftershocks using both one-dimensional and three-dimensional velocity inversions. The depth range of the aftershocks (1–10 km) agrees well with the preferred mainshock centroid depth. The aftershocks near the hypocenter suggest a structure dipping toward the NW, which we interpret to be the mainshock fault plane. This structure and the Harper fault, ∼15 km to the south, appear to have acted as boundaries to the extensive aftershock zone trending NNW-SSE. Most of the M L ≥ 5 aftershocks, including the two largest ( M L 6.1 and M L 5.7), clustered near the Harper fault and have strike slip mechanisms consistent with motion on this fault and its conjugates. Forward modeling of the GPS data suggests that a reverse slip mainshock, combined with strike slip aftershock faulting in the south, is able to match the observed displacements. The occurrence of this earthquake sequence implies that the level of seismic hazard in the central South Island is greater than previous estimates.

Arthur's Pass↗

USGS2008

Landslides Mapped from LIDAR Imagery, Kitsap County, Washington

Landslides are a recurring problem on hillslopes throughout the Puget Lowland, Washington, but can be difficult to identify in the densely forested terrain. However, digital terrain models of the bare-earth surface derived from LIght Detection And Ranging (LIDAR) data express topographic details sufficiently well to identify landslides. Landslides and escarpments were mapped using LIDAR imagery and field checked (when permissible and accessible) throughout Kitsap County. We relied almost entirely on derivatives of LIDAR data for our mapping, including topographic-contour, slope, and hill-shaded relief maps. Each mapped landslide was assigned a level of 'high' or 'moderate' confidence based on the LIDAR characteristics and on field observations. A total of 231 landslides were identified representing 0.8 percent of the land area of Kitsap County. Shallow debris topples along the coastal bluffs and large (>10,000 m2) landslide complexes are the most common types of landslides. The smallest deposit mapped covers an area of 252 m2, while the largest covers 0.5 km2. Previous mapping efforts that relied solely on field and photogrammetric methods identified only 57 percent of the landslides mapped by LIDAR (61 percent high confidence and 39 percent moderate confidence), although nine landslides previously identified were not mapped during this study. The remaining 43 percent identified using LIDAR have 13 percent high confidence and 87 percent moderate confidence. Coastal areas are especially susceptible to landsliding; 67 percent of the landslide area that we mapped lies within 500 meters of the present coastline. The remaining 33 percent are located along drainages farther inland. The LIDAR data we used for mapping have some limitations including (1) rounding of the interface area between low slope surfaces and vertical faces (that is, along the edges of steep escarpments) which results in scarps being mapped too far headward (one or two meters), (2) incorrect laser-distance measurements resulting in inaccurate elevation values, (3) removal of valid ground elevations, (4) false ground roughness, and (5) faceted surface texture. Several of these limitations are introduced by algorithms in the processing software that are designed to remove non-ground elevations from LIDAR data. Despite these limitations, the algorithm-enhanced LIDAR imagery does effectively 'remove' vegetation that obscures many landslides, and is therefore a valuable tool for landslide inventories and investigations in heavily vegetated regions such as the Puget Lowland.

Open-File Report↗

USGS2016

The Eastern California Shear Zone as the northward extension of the southern San Andreas Fault

Cluster analysis offers an agnostic way to organize and explore features of the current GPS velocity field without reference to geologic information or physical models using information only contained in the velocity field itself. We have used cluster analysis of the Southern California Global Positioning System (GPS) velocity field to determine the partitioning of Pacific-North America relative motion onto major regional faults. Our results indicate the large-scale kinematics of the region is best described with two boundaries of high velocity gradient, one centered on the Coachella section of the San Andreas Fault and the Eastern California Shear Zone and the other defined by the San Jacinto Fault south of Cajon Pass and the San Andreas Fault farther north. The ~120 km long strand of the San Andreas between Cajon Pass and Coachella Valley (often termed the San Bernardino and San Gorgonio sections) is thus currently of secondary importance and carries lesser amounts of slip over most or all of its length. We show these first order results are present in maps of the smoothed GPS velocity field itself. They are also generally consistent with currently available, loosely bounded geologic and geodetic fault slip rate estimates that alone do not provide useful constraints on the large-scale partitioning we show here. Our analysis does not preclude the existence of smaller blocks and more block boundaries in Southern California. However, attempts to identify smaller blocks along and adjacent to the San Gorgonio section were not successful.

Eastern California Shear Zone↗

USGS2000

Synthetic aperture radar interferometry of Okmok volcano, Alaska: radar observations

ERS-1/ERS-2 synthetic aperture radar interferometry was used to study the 1997 eruption of Okmok volcano in Alaska. First, we derived an accurate digital elevation model (DEM) using a tandem ERS-1/ERS-2 image pair and the preexisting DEM. Second, by studying changes in interferometric coherence we found that the newly erupted lava lost radar coherence for 5-17 months after the eruption. This suggests changes in the surface backscattering characteristics and was probably related to cooling and compaction processes. Third, the atmospheric delay anomalies in the deformation interferograms were quantitatively assessed. Atmospheric delay anomalies in some of the interferograms were significant and consistently smaller than one to two fringes in magnitude. For this reason, repeat observations are important to confidently interpret small geophysical signals related to volcanic activities. Finally, using two-pass differential interferometry, we analyzed the preemptive inflation, coeruptive deflation, and posteruptive inflation and confirmed the observations using independent image pairs. We observed more than 140 cm of subsidence associated with the 1997 eruption. This subsidence occurred between 16 months before the eruption and 5 months after the eruption, was preceded by ∼18 cm of uplift between 1992 and 1995 centered in the same location, and was followed by ∼10 cm of uplift between September 1997 and 1998. The best fitting model suggests the magma reservoir resided at 2.7 km depth beneath the center of the caldera, which was ∼5 km from the eruptive vent. We estimated the volume of the erupted material to be 0.055 km 3 and the average thickness of the erupted lava to be ∼7.4 m. Copyright 2000 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗