USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Interpretation”

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 883 records · Page 49Linked to original sources

Voluminous lava-like precursor to a major ash-flow tuff: Low-column pyroclastic eruption of the Pagosa Peak Dacite, San Juan volcanic field, Colorado

The Pagosa Peak Dacite is an unusual pyroclastic deposit that immediately predated eruption of the enormous Fish Canyon Tuff (~5000 km3) from the La Garita caldera at 28 Ma. The Pagosa Peak Dacite is thick (to 1 km), voluminous (>200 km3), and has a high aspect ratio (1:50) similar to those of silicic lava flows. It contains a high proportion (40-60%) of juvenile clasts (to 3-4 m) emplaced as viscous magma that was less vesiculated than typical pumice. Accidental lithic fragments are absent above the basal 5-10% of the unit. Thick densely welded proximal deposits flowed rheomorphically due to gravitational spreading, despite the very high viscosity of the crystal-rich magma, resulting in a macroscopic appearance similar to flow-layered silicic lava. Although it is a separate depositional unit, the Pagosa Peak Dacite is indistinguishable from the overlying Fish Canyon Tuff in bulk-rock chemistry, phenocryst compositions, and 40Ar/39Ar age. The unusual characteristics of this deposit are interpreted as consequences of eruption by low-column pyroclastic fountaining and lateral transport as dense, poorly inflated pyroclastic flows. The inferred eruptive style may be in part related to synchronous disruption of the southern margin of the Fish Canyon magma chamber by block faulting. The Pagosa Peak eruptive sources are apparently buried in the southern La Garita caldera, where northerly extensions of observed syneruptive faults served as fissure vents. Cumulative vent cross-sections were large, leading to relatively low emission velocities for a given discharge rate. Many successive pyroclastic flows accumulated sufficiently rapidly to weld densely as a cooling unit up to 1000 m thick and to retain heat adequately to permit rheomorphic flow. Explosive potential of the magma may have been reduced by degassing during ascent through fissure conduits, leading to fracture-dominated magma fragmentation at low vesicularity. Subsequent collapse of the 75 x 35 km2 La Garita caldera and eruption of the Fish Canyon Tuff were probably triggered by destabilization of the chamber roof as magma was withdrawn during the Pagosa Peak eruption. (C) 2000 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Proportions, timing, and re-equilibration progress during the 1959 Summit Eruption of Kīlauea: An example of magma mixing processes operating during OIB petrogenesis

Petrographic and chemical analysis of scoria samples collected during the 1959 Kīlauea summit eruption illustrates the progress of thermal and chemical homogenization of the melts, and the gradual growth and/or re-equilibration of olivine phenocrysts, over the course of the eruption. Glass compositions show that thermal equilibration was largely complete within the span of the eruption, whereas chemical homogenization was a work in progress. The olivine phenocryst population, known to contain conspicuous antecrystic components, is also hybrid within the euhedral population. The bulk of the olivine reached the level of the erupting magma on November 18–19, 1959. Zoning patterns in olivine phenocrysts show that initially unzoned grains developed normal zoning by the end of the eruption. Reverse zoning in relatively Fe-rich olivine phenocrysts (interpreted as cognate to the stored magma) was progressively eliminated from November 21 to December 19, 1959, by diffusive re-equilibration between crystals and melt. Toward the end of the eruption, the only olivine composition in direct contact with the melt was Fo 84–86 , with the original rim compositional heterogeneity gone in 4–5 weeks’ time. Activity in December 1959 differed from that in November, as high fountaining events were more closely spaced and almost all samples were picritic, with bulk MgO ≥16·5 wt%. Three different levels were in play during the 1959 eruption: a deep source for high-MgO melts and forsteritic (Fo 87–89 ) olivines, an intermediate source for the bulk of the stored magma, and a shallower source for the most differentiated magma. This model is consistent with geophysical, petrological and chemical observations. Comparison of the 1959 eruption with results from older explosive deposits suggests that stored and recharge melts and olivine from the deeper parts of Kīlauea’s plumbing are similar in composition to those observed or inferred in the 1959 eruption, so they behave similarly during extrusive and explosive periods alike.

Hawaii↗

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↗

Point spread functions for earthquake source imaging: An interpretation based on seismic interferometry

Recently, various methods have been proposed and applied for earthquake source imaging, and theoretical relationships among the methods have been studied. In this study, we make a follow-up theoretical study to better understand the meanings of earthquake source imaging. For imaging problems, the point spread function (PSF) is used to describe the degree of blurring and degradation in an obtained image of a target object as a response of an imaging system. In this study, we formulate PSFs for earthquake source imaging. By calculating the PSFs, we find that waveform source inversion methods remove the effect of the PSF and are free from artifacts. However, the other source imaging methods are affected by the PSF and suffer from the effect of blurring and degradation due to the restricted distribution of receivers. Consequently, careful treatment of the effect is necessary when using the source imaging methods other than waveform inversions. Moreover, the PSF for source imaging is found to have a link with seismic interferometry with the help of the source-receiver reciprocity of Green’s functions. In particular, the PSF can be related to Green’s function for cases in which receivers are distributed so as to completely surround the sources. Furthermore, the PSF acts as a low-pass filter. Given these considerations, the PSF is quite useful for understanding the physical meaning of earthquake source imaging.

Geophysical Journal International↗

Occam's shadow: levels of analysis in evolutionary ecology - where to next?

Evolutionary ecology is the study of evolutionary processes, and the ecological conditions that influence them. A fundamental paradigm underlying the study of evolution is natural selection. Although there are a variety of operational definitions for natural selection in the literature, perhaps the most general one is that which characterizes selection as the process whereby heritable variation in fitness associated with variation in one or more phenotypic traits leads to intergenerational change in the frequency distribution of those traits. The past 20 years have witnessed a marked increase in the precision and reliability of our ability to estimate one or more components of fitness and characterize natural selection in wild populations, owing particularly to significant advances in methods for analysis of data from marked individuals. In this paper, we focus on several issues that we believe are important considerations for the application and development of these methods in the context of addressing questions in evolutionary ecology. First, our traditional approach to estimation often rests upon analysis of aggregates of individuals, which in the wild may reflect increasingly non-random (selected) samples with respect to the trait(s) of interest. In some cases, analysis at the aggregate level, rather than the individual level, may obscure important patterns. While there are a growing number of analytical tools available to estimate parameters at the individual level, and which can cope (to varying degrees) with progressive selection of the sample, the advent of new methods does not reduce the need to consider carefully the appropriate level of analysis in the first place. Estimation should be motivated a priori by strong theoretical analysis. Doing so provides clear guidance, in terms of both (i) assisting in the identification of realistic and meaningful models to include in the candidate model set, and (ii) providing the appropriate context under which the results are interpreted. Second, while it is true that selection (as defined) operates at the level of the individual, the selection gradient is often (if not generally) conditional on the abundance of the population. As such, it may be important to consider estimating transition rates conditional on both the parameter values of the other individuals in the population (or at least their distribution), and population abundance. This will undoubtedly pose a considerable challenge, for both single- and multi-strata applications. It will also require renewed consideration of the estimation of abundance, especially for open populations. Thirdly, selection typically operates on dynamic, individually varying traits. Such estimation may require characterizing fitness in terms of individual plasticity in one or more state variables, constituting analysis of the norms of reaction of individuals to variable environments. This can be quite complex, especially for traits that are under facultative control. Recent work has indicated that the pattern of selection on such traits is conditional on the relative rates of movement among and frequency of spatially heterogeneous habitats, suggesting analyses of evolution of life histories in open populations can be misleading in some cases.

Journal of Applied Statistics↗

Tidal triggering of low frequency earthquakes near Parkfield, California: Implications for fault mechanics within the brittle-ductile transition

Studies of nonvolcanic tremor (NVT) have established the significant impact of small stress perturbations on NVT generation. Here we analyze the influence of the solid earth and ocean tides on a catalog of ∼550,000 low frequency earthquakes (LFEs) distributed along a 150 km section of the San Andreas Fault centered at Parkfield. LFE families are identified in the NVT data on the basis of waveform similarity and are thought to represent small, effectively co-located earthquakes occurring on brittle asperities on an otherwise aseismic fault at depths of 16 to 30 km. We calculate the sensitivity of each of these 88 LFE families to the tidally induced right-lateral shear stress (RLSS), fault-normal stress (FNS), and their time derivatives and use the hypocentral locations of each family to map the spatial variability of this sensitivity. LFE occurrence is most strongly modulated by fluctuations in shear stress, with the majority of families demonstrating a correlation with RLSS at the 99% confidence level or above. Producing the observed LFE rate modulation in response to shear stress perturbations requires low effective stress in the LFE source region. There are substantial lateral and vertical variations in tidal shear stress sensitivity, which we interpret to reflect spatial variation in source region properties, such as friction and pore fluid pressure. Additionally, we find that highly episodic, shallow LFE families are generally less correlated with tidal stresses than their deeper, continuously active counterparts. The majority of families have weaker or insignificant correlation with positive (tensile) FNS. Two groups of families demonstrate a stronger correlation with fault-normal tension to the north and with compression to the south of Parkfield. The families that correlate with fault-normal clamping coincide with a releasing right bend in the surface fault trace and the LFE locations, suggesting that the San Andreas remains localized and contiguous down to near the base of the crust. The deep families that have high sensitivity to both shear and tensile normal stress perturbations may be indicative of an increase in effective fault contact area with depth. Synthesizing our observations with those of other LFE-hosting localities will help to develop a comprehensive understanding of transient fault slip below the “seismogenic zone” by providing constraints on parameters in physical models of slow slip and LFEs.

California↗

Interplay of late Cenozoic siliciclastic supply and carbonate response on the southeast Florida platform

High-resolution seismic-reflection data collected along the length of the Caloosahatchee River in southwestern Florida have been correlated to nannofossil biostratigraphy and strontium-isotope chemostratigraphy at six continuously cored boreholes. These data are interpreted to show a major Late Miocene(?) to Early Pliocene fluvial-deltaic depositional system that prograded southward across the carbonate Florida Platform , interrupting nearly continuous carbonate deposition since early in the Cretaceous. Connection of the platform top to a continental source of siliciclastics and significant paleotopography combined to focus accumulation of an immense supply of siliciclastics on the southeastern part of the Florida Platform . The remarkably thick (> 100 m), sand-rich depositional system, which is characterized by clinoformal progradation, filled in deep accommodation, while antecedent paleotopography directed deltaic progradation southward within the middle of the present-day Florida Peninsula. The deltaic depositional system may have prograded about 200 km southward to the middle and upper Florida Keys, where Late Miocene to Pliocene siliciclastic form the foundation of the Quaternary carbonate shelf and shelf margin of the Florida Keys. These far-traveled siliciclastic deposits filled accommodation on the southeastern part of the Florida Platform so that paleobathymetry was sufficiently shallow to allow Quaternary recovery of carbonate sedimentation in the area of southern peninsular Florida and the Florida Keys.

Florida↗

Habitat use and movement of the endangered Arroyo Toad ( Anaxyrus californicus ) in coastal southern California

Information on the habitat use and movement patterns of Arroyo Toads ( Anaxyrus californicus ) is limited. The temporal and spatial characteristics of terrestrial habitat use, especially as it relates to upland use in coastal areas of the species' range, are poorly understood. We present analyses of radiotelemetry data from 40 individual adult toads tracked at a single site in coastal southern California from March through November of 2004. We quantify adult Arroyo Toad habitat use and movements and interpret results in the context of their life history. We show concentrated activity by both male and female toads along stream terraces during and after breeding, and, although our fall sample size is low, the continued presence of adult toads in the floodplain through the late fall. Adult toads used open sandy flats with sparse vegetation. Home-range size and movement frequency varied as a function of body mass. Observed spatial patterns of movement and habitat use both during and outside of the breeding period as well as available climatological data suggest that overwintering of toads in floodplain habitats of near-coastal areas of southern California may be more common than previously considered. If adult toads are not migrating out of the floodplain at the close of the breeding season but instead overwinter on stream terraces in near-coastal areas, then current management practices that assume toad absence from floodplain habitats may be leaving adult toads over-wintering on stream terraces vulnerable to human disturbance during a time of year when Arroyo Toad mortality is potentially highest.

California↗

Uplift and subsidence associated with the great Aceh-Andaman earthquake of 2004

Rupture of the Sunda megathrust on 26 December 2004 produced broad regions of uplift and subsidence. We define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great Mw ??? 9.2 earthquake. In the region of the Andaman and Nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model. At the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls. Uplift extends from the middle of Simeulue Island, Sumatra, at ??? 2.5??N, to Preparis Island, Myanmar (Burma), at ??? 14.9??N. Thus the rupture is ??? 1600 km long. The distance from the pivot line to the trench varies appreciably. The northern and western Andaman Islands rose, whereas the southern and eastern portion of the islands subsided. The Nicobar Islands and the west coast of Aceh province, Sumatra, subsided. Tilt at the southern end of the rupture is steep; the distance from 1.5 m of uplift to the pivot line is just 60 km. Our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Field Reconnaissance Geologic Mapping of the Columbia Hills, Mars: Results from MER Spirit and MRO HiRISE Observations

Chemical, mineralogic, and lithologic ground truth was acquired for the first time on Mars in terrain units mapped using orbital Mars Reconnaissance Orbiter's High Resolution Imaging Science Experiment (MRO HiRISE) image data. Examination of several dozen outcrops shows that Mars is geologically complex at meter length scales, the record of its geologic history is well exposed, stratigraphic units may be identified and correlated across significant areas on the ground, and outcrops and geologic relationships between materials may be analyzed with techniques commonly employed in terrestrial field geology. Despite their burial during the course of Martian geologic time by widespread epiclastic materials, mobile fines, and fall deposits, the selective exhumation of deep and well-preserved geologic units has exposed undisturbed outcrops, stratigraphic sections, and structural information much as they are preserved and exposed on Earth. A rich geologic record awaits skilled future field investigators on Mars. The correlation of ground observations and orbital images enables construction of a corresponding geologic reconnaissance map. Most of the outcrops visited are interpreted to be pyroclastic, impactite, and epiclastic deposits overlying an unexposed substrate, probably related to a modified Gusev crater central peak. Fluids have altered chemistry and mineralogy of these protoliths in degrees that vary substantially within the same map unit. Examination of the rocks exposed above and below the major unconformity between the plains lavas and the Columbia Hills directly confirms the general conclusion from remote sensing in previous studies over past years that the early history of Mars was a time of more intense deposition and modification of the surface. Although the availability of fluids and the chemical and mineral activity declined from this early period, significant later volcanism and fluid convection enabled additional, if localized, chemical activity.

Journal of Geophysical Research E: Planets↗

Northern East Pacific Rise: Magnetic anomaly and bathymetric framework

The oceanic crust in the eastern Pacific between 7°N and 30°N and east of 127°W contains a fairly complete history of the spreading centers associated with the East Pacific Rise since 25 m.y. B.P. (late Oligocene). In this paper, we have summarized the seafloor spreading magnetic-anomaly data and the bathymetric data that reflect the record of this tectonic history. The well-defined magnetic lineations north of the Clarion fracture zone, in the mouth of the Gulf of California, and on the east flank of the East Pacific Rise (EPR) are carefully examined and used to provide a guide for interpreting the spreading pattern between the Clarion and Clipperton fracture zones, southward of the Rivera fracture zone over the Mathematician Ridge, and over the entire EPR east of the Mathematician Ridge between the Rivera and Siqueiros fracture zones. The bathymetric data provide a trace of the fracture zone pattern in each of the above mentioned areas. The fracture zone bathymetry and the seafloor spreading magnetic lineations on the EPR south of the Rivera fracture zone have a distinctive fanning pattern caused by close poles of rotation and plate boundary reorganizations. All these data provide a good record of the plate reorganizations in the middle Miocene at magnetic anomaly 5 A time (12.5 to 11 m.y. B.P.), in the late Miocene at magnetic anomaly 3′−4 time (6.5 m.y. B.P.), and in the Pliocene at magnetic anomaly 2′−3 time (3.5 m.y. B.P.). Several abandoned spreading centers, including the Mathematician Ridge, were left behind as a result of these reorganizations. The Mathematician Ridge is shown to be a set of ridges and trough whose origin is related to the tectonic activity associated with each of the above mentioned reorganizations since anomaly 5A.

Journal of Geophysical Research B: Solid Earth↗

Evidence for mechanical and chemical alteration of iron-nickel meteorites on Mars: Process insights for Meridiani Planum

The weathering of meteorites found on Mars involves chemical and physical processes that can provide clues to climate conditions at the location of their discovery. Beginning on sol 1961, the Opportunity rover encountered three large iron meteorites within a few hundred meters of each other. In order of discovery, these rocks have been assigned the unofficial names Block Island, Shelter Island, and Mackinac Island. Each rock presents a unique but complimentary set of features that increase our understanding of weathering processes at Meridiani Planum. Significant morphologic characteristics interpretable as weathering features include (1) a large pit in Block Island, lined with delicate iron protrusions suggestive of inclusion removal by corrosive interaction; (2) differentially eroded kamacite and taenite lamellae in Block Island and Shelter Island, providing relative timing through crosscutting relationships with deposition of (3) an iron oxide-rich dark coating; (4) regmaglypted surfaces testifying to regions of minimal surface modification, with other regions in the same meteorites exhibiting (5) large-scale, cavernous weathering (in Shelter Island and Mackinac Island). We conclude that the current size of the rocks is approximate to their original postfall contours. Their morphology thus likely results from a combination of atmospheric interaction and postfall weathering effects. Among our specific findings is evidence supporting (1) at least one possible episode of aqueous acidic exposure for Block Island; (2) ripple migration over portions of the meteorites; (3) a minimum of two separate episodes of wind abrasion; alternating with (4) at least one episode of coating-forming chemical alteration, most likely at subzero temperatures.

Journal of Geophysical Research E: Planets↗

Pre-eruptive outgassing and pressurization, and post-fragmentation bubble nucleation, recorded by vesicles in breadcrust bombs from vulcanian activity at Guagua Pichincha Volcano, Ecuador

Breadcrust bombs formed during Vulcanian eruptions are assumed to originate from the shallow plug or dome. Their rim to core texture reflects the competition between cooling and degassing timescales, which results in a dense crust with isolated vesicles contrasting with a highly vesicular vesicle network in the interior. Due to relatively fast quenching, the crust can shed light on pre- and syn-eruptive conditions prior to or during fragmentation, whereas the interior allows us to explore post-fragmentation vesiculation. Investigation of pre- to post-fragmentation processes in breadcrust bombs from the 1999 Vulcanian activity at Guagua Pichincha, Ecuador, via 2D and 3D textural analysis reveals a complex vesiculation history, with multiple, spatially localized nucleation and growth events. Large vesicles (Type 1), present in low number density in the crust, are interpreted as pre-eruptive bubbles formed by outgassing and collapse of a permeable bubble network during ascent or stalling in the plug. Haloes of small, syn-fragmentation vesicles (Type 2), distributed about large vesicles, are formed by pressurization and enrichment of volatiles in these haloes. The nature of the pressurization process in the plug is discussed in light of seismicity and ground deformation signals, and previous textural and chemical studies. A third population (Type 3) of post-fragmentation small vesicles appears in the interior of the bomb, and growth and coalescence of Type 2 and 3 vesicles causes the transition from isolated to interconnected bubble network in the interior. We model the evolution of viscosity, bubble growth rate, diffusion timescales, bubble radius and porosity during fragmentation and cooling. These models reveal that thermal quenching dominates in the crust whereas the interior undergoes a viscosity quench caused by degassing, and that the transition from crust to interior corresponds to the onset of percolation and development of permeability in the bubble network.

Guagua Pichincha Volcano↗

The complex frequencies of long-period seismic events as probes of fluid composition beneath volcanoes

Long-period (LP) events have been widely observed in relation to magmatic and hydrothermal activities in volcanic areas. LP waveforms characterized by their harmonic signature have been interpreted as oscillations of a fluid-filled resonator, and mixtures of liquid and gas in the form of bubbly liquids have been mainly assumed for the fluid. To investigate the characteristic properties of the resonator system, we analyse waveforms of LP events observed at four different volcanoes in Hawaii, Alaska, Colombia and Japan using a newly developed spectral method. This method allows an estimation of the complex frequencies of decaying sinusoids based on an autoregressive model. The results of our analysis show a wide variety of Q factors, ranging from tens to several hundred. We compare these complex frequencies with those predicted by the fluid-filled crack model for various mixtures of liquid, gas and ash. Although the oscillations of LP events with Q smaller than 50 can be explained by various combinations of liquids and gases, we find that ash-laden gases are required to explain long-lasting oscillations with Q larger than 100. The complex frequencies of LP events yield useful information on the types of fluids. Temporal and spatial variations of the complex frequencies can be used as probes of fluid composition beneath volcanoes.

Geophysical Journal International↗

Shrinking channels, growing threats: Habitat degradation from channel narrowing and invasive vegetation in three dryland rivers

Water development and the proliferation of invasive riparian vegetation have led to widespread habitat loss and simplification of rivers in the western United States, contributing to the imperilment of native fishes. Here, we quantify channel narrowing and vegetation encroachment, which are conspicuous indicators of riverine habitat alteration, along ∼400 km of three dryland tributaries of the upper Colorado River. We conducted a comparative analysis of aerial photographs between the 1930s and 2010s/2020s time periods using visual interpretation and used Light Detection and Ranging (LiDAR) data along with Object-Based Image Analysis (OBIA) to quantify canopy cover of woody riparian species. All three rivers underwent substantial channel narrowing, coinciding with a general decrease in spring floods over time. However, the extent of narrowing varied among the rivers (78 %, 73 %, and 29 %) with greater narrowing corresponding to larger reductions in spring flows. In contrast, contemporary woody cover was similarly high among all three rivers (39 %, 41 %, and 36 %), and a woody vegetation analysis we conducted for one river indicated a substantial increase in vegetation along the active channel (4 %–74 %). These findings underscore a common pattern observed in rivers throughout the basin, where river channels often undergo narrowing and encroachment by invasive vegetation following dam construction and/or decreases in flows, ultimately leading to habitat simplification, with negative implications for native fishes and other riparian biota. Our findings also emphasize that, even in the presence of nonnative vegetation establishment, preserving or restoring large magnitude and long duration floods can help conserve diverse habitat in dryland rivers.

Arizona, Colorado, New Mexico, Utah↗

Estimating instream constituent loads using replicate synoptic sampling, Peru Creek, Colorado

The synoptic mass balance approach is often used to evaluate constituent mass loading in streams affected by mine drainage. Spatial profiles of constituent mass load are used to identify sources of contamination and prioritize sites for remedial action. This paper presents a field scale study in which replicate synoptic sampling campaigns are used to quantify the aggregate uncertainty in constituent load that arises from (1) laboratory analyses of constituent and tracer concentrations, (2) field sampling error, and (3) temporal variation in concentration from diel constituent cycles and/or source variation. Consideration of these factors represents an advance in the application of the synoptic mass balance approach by placing error bars on estimates of constituent load and by allowing all sources of uncertainty to be quantified in aggregate; previous applications of the approach have provided only point estimates of constituent load and considered only a subset of the possible errors. Given estimates of aggregate uncertainty, site specific data and expert judgement may be used to qualitatively assess the contributions of individual factors to uncertainty. This assessment can be used to guide the collection of additional data to reduce uncertainty. Further, error bars provided by the replicate approach can aid the investigator in the interpretation of spatial loading profiles and the subsequent identification of constituent source areas within the watershed. The replicate sampling approach is applied to Peru Creek, a stream receiving acidic, metal-rich effluent from the Pennsylvania Mine. Other sources of acidity and metals within the study reach include a wetland area adjacent to the mine and tributary inflow from Cinnamon Gulch. Analysis of data collected under low-flow conditions indicates that concentrations of Al, Cd, Cu, Fe, Mn, Pb, and Zn in Peru Creek exceed aquatic life standards. Constituent loading within the study reach is dominated by effluent from the Pennsylvania Mine, with over 50% of the Cd, Cu, Fe, Mn, and Zn loads attributable to a collapsed adit near the top of the study reach. These estimates of mass load may underestimate the effect of the Pennsylvania Mine as leakage from underground mine workings may contribute to metal loads that are currently attributed to the wetland area. This potential leakage confounds the evaluation of remedial options and additional research is needed to determine the magnitude and location of the leakage.

Colorado↗

Scaling variables and interpretation of eigenvalues in principal component analysis of geologic data

The dominant feature distinguishing one method of principal components analysis from another is the manner in which the original data are transformed prior to the other computations. The only other distinguishing feature of any importance is whether the eigenvectors of the inner product-moment of the transformed data matrix are taken directly as the Q-mode scores or scaled by the square roots of their associated eigenvalues and called the R-mode loadings. If the eigenvectors are extracted from the product-moment correlation matrix, the variables, in effect, were transformed by column standardization (zero means and unit variances), and the sum of the p-largest eigenvalues divided by the sum of all the eigenvalues indicates the degree to which a model containing p components will account for the total variance in the original data. However, if the data were transformed in any manner other than column standardization, the eigenvalues cannot be used in this manner, but can only be used to determine the degree to which the model will account for the transformed data. Regardless of the type of principal components analysis that is performed-even whether it is R or Q-mode-the goodness-of-fit of the model to the original data is given better by the eigenvalues of the correlation matrix than by those of the matrix that was actually factored. ?? 1980 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗

Evidence for a shallow pluton beneath the Goat Rocks Wilderness, Washington, from gravity and magnetic data

A gravity and aeromagnetic study was conducted in and adjacent to the Goat Rocks Wilderness area, Washington. This work was done in conjunction with geologic and geochemical studies which were critically important to our interpretation. The Bouguer gravity anomaly map displays a large central positive gravity anomaly and three superimposed smaller positive anomalies. The main anomaly, about 26 by 11 km with an amplitude of about 16 mGal seems to delineate a dense large pluton or sheeted dike complex intruding less dense, older volcanic and sedimentary rocks lying directly below the volcano. The pluton or sheeted dike complex may have fed the Goat Rocks volcano throughout its 1–2 m.y. history. Three small feeders, mapped on the surface as andesite intrusions, may be the cause of three smaller gravity anomalies superimposed on the main anomaly. The magnetic anomaly maps also exhibit a large central anomaly whose source is the Goat Rocks pluton. The magnetic data suggest that the three shallow intrusive cupolas have undergone a significant amount of hydrothermal alteration. In addition, to some extent the magnetic data reflect the terrain of the area. Topographic highs produce magnetic highs, and topographic lows produce magnetic lows. The deviations from this pattern delineate subsurface geology. Magnetic lows may be of particular interest in locating zones of mineralization where magnetic minerals have been altered and have become significantly less magnetic.

Washington↗