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 937 records · Page 52Linked to original sources

Diet of the Antillean manatee (Trichechus manatus manatus) in Belize, Central America

Belize contains important habitat for Antillean manatees (Trichechus manatus manatus) and provides refuge for the highest known population density of this subspecies. As these animals face impending threats, knowledge of their dietary habits can be used to interpret resource utilization. The contents of 13 mouth, 6 digestive tract (stomach, duodenum and colon), and 124 fecal samples were microscopically examined using a modified point technique detection protocol to identify key plant species consumed by manatees at two important aggregation sites in Belize: Southern Lagoon and the Drowned Cayes. Overall, 15 different items were identified in samples from manatees in Belize. Five species of seagrasses (Halodule wrightii, Thalassia testudinum, Ruppia maritima, Syringodium filiforme, and Halophila sp.) made up the highest percentage of items. The red mangrove (Rhizophora mangle), was also identified as an important food item. Algae (Ulva sp., Chara sp., Lyngbya sp.) and invertebrates (sponges and diatoms) were also consumed. Variation in the percentage of seagrasses, other vascular plants, and algae consumption was analyzed as a 4-factor analysis of variance (ANOVA) with main effects and interactions for locality, sex, size classification, and season. While sex and season did not influence diet composition, differences for locality and size classification were observed. These results suggest that analysis of diet composition of Antillean manatees may help to determine critical habitat and use of associated food resources which, in turn can be used to aid conservation efforts in Belize.

Journal of the Marine Biological Association of th↗

Modeling consolidation and dewatering near the toe of the northern Barbados accretionary complex

At the toe of the northern Barbados accretionary complex, temperature and pore water chemistry data indicate that fluid flow is channeled along the décollement and other shallow thrust faults. We examine mechanisms that may prevent consolidation and maintain high permeability over large sections of the décollement. High-resolution bulk density data from five boreholes show that the décollement is well consolidated at some sites while other sites remain underconsolidated. Underconsolidated décollement behavior is associated with kilometer-scale negative-polarity seismic reflections from the décollement plane that have been interpreted to be fluid conduits. We use a coupled fluid flow/consolidation model to simulate the loading response of a 10-km-long by 680-m-thick slice of sediment as it enters the accretionary complex. The simulations capture 185 ka (5 km) of subduction, with a load function representing the estimated effective stress of the overriding accretionary prism (3.8° taper angle). Simulation results of bulk density in the décollement 3.2 km arcward of the deformation front are compared with observations. The results show that persistent high pore pressures at the arcward edge of the simulation domain can explain underconsolidated behavior. The scenario is consistent with previous modeling results showing that high pore pressures can propagate intermittently along the décollement from deeper in the complex. Simulated seaward fluxes in the décollement (1–14 cm yr −1 ) lie between previous estimates from modeling studies of steady state (<1 cm yr −1 ) and transient (>1 m yr −1 ) flow. Maximum simulated instantaneous fluid sources (2.5×10 −13 s −1 ) are comparable to previous estimates. The simulations show minor swelling of incoming sediments (fluid sources ∼−3×10 −15 s −1 ) up to 3 km before subduction that may help to explain small-scale shearing and normal faulting proximal to the protodécollement.

Atlantic Ocean, Northern Barbados Accretionary Com↗

Isotopic evidence for early Proterozoic age of the Idono Complex, west-central Alaska

The Idono Complex of west-central Alaska is a fault-bounded, fragment of Early Proterozoic continental crust surrounded by overlap assemblages and younger terranes accreted in Mesozoic time. It is composed of granitic to dioritic orthogneiss, amphibolite, and metasedimentary rocks. Trace element compositions of the granitoids and amphibolite suggest rock formation in a subduction-related volcanic arc terrain. Nine zircon fractions from three samples of granitoid orthogneiss define a U-Pb discordia line intersecting concordia at and . The upper and lower intercepts are interpreted, respectively, as approximations of the time of granitoid crystallization and major episodic Pb-loss. The lower intercept age is similar to that indicated by some hornblende and biotite K-Ar ages. Other biotite, hornblende, and white mica K-Ar ages record later isotopic disturbance in Early Cretaceous time. Nd isotopic compositions at 2.06 Ga for tonalite orthogneiss ( ) and amphibolite ( and +2.3) indicate that Early Proterozoic crust formation in the Idono Complex involved significant additions of mantle-derived magma. These rocks yield (depleted mantle) model ages of 2.0 to 2.1 Ga, similar to the U-Pb zircon upper intercept age. Involvement of Archean crust in the formation of at least some rocks of the Idono Complex is indicated by a granitic orthogneiss ( ), which yields a model age of 2.5 Ga. The ~2.06-Ga age of the Idono Complex is similar to that indicated for provenance(s) of widely distributed sedimentary rocks in western North America and may represent a fragment of the North American craton displaced by northward movement of the Pacific plate. At such, rocks of the Idono Complex may provide important insights into both Early Proterozoic evolution along the craton margin, and subsequent displacements.

Journal of Geology↗

Basis for paleoenvironmental interpretation of magnetic properties of sediment from Upper Klamath Lake (Oregon): Effects of weathering and mineralogical sorting

Studies of magnetic properties enable reconstruction of environmental conditions that affected magnetic minerals incorporated in sediments from Upper Klamath Lake. Analyses of stream sediment samples from throughout the catchment of Upper Klamath Lake show that alteration of Fe-oxide minerals during subaerial chemical weathering of basic volcanic rocks has significantly changed magnetic properties of surficial deposits. Titanomagnetite, which is abundant both as phenocrysts and as microcrystals in fresh volcanic rocks, is progressively destroyed during weathering. Because fine-grained magnetite is readily altered due to large surface-to-volume ratios, weathering causes an increase in average magnetic grain size as well as reduction in the quantity of titanomagnetite both absolutely and relative to hematite. Hydrodynamic mineralogical sorting also produces differences in magnetic properties among rock and mineral grains of differing sizes. Importantly, removal of coarse silicate and Fe-oxide grains by sorting concentrated extremely fine-grained magnetite in the resulting sediment. The effects of weathering and sorting of minerals cannot be completely separated. These processes combine to produce the magnetic properties of a non-glacial lithic component of Upper Klamath Lake sediments, which is characterized by relatively low magnetite content and coarse magnetic grain size. Hydrodynamic sorting alone causes significant differences between the magnetic properties of glacial flour in lake sediments and of fresh volcanic rocks in the catchment. In comparison to source volcanic rocks, glacial flour in the lake sediment is highly enriched in extremely fine-grained magnetite.

Journal of Paleolimnology↗

Seismic velocities within the sedimentary succession of the Canada Basin and southern Alpha-Mendeleev Ridge, Arctic Ocean: evidence for accelerated porosity reduction?

The Canada Basin and the southern Alpha-Mendeleev ridge complex underlie a significant proportion of the Arctic Ocean, but the geology of this undrilled and mostly ice-covered frontier is poorly known. New information is encoded in seismic wide-angle reflections and refractions recorded with expendable sonobuoys between 2007 and 2011. Velocity–depth samples within the sedimentary succession are extracted from published analyses for 142 of these records obtained at irregularly spaced stations across an area of 1.9E + 06 km 2 . The samples are modelled at regional, subregional and station-specific scales using an exponential function of inverse velocity versus depth with regionally representative parameters determined through numerical regression. With this approach, smooth, non-oscillatory velocity–depth profiles can be generated for any desired location in the study area, even where the measurement density is low. Practical application is demonstrated with a map of sedimentary thickness, derived from seismic reflection horizons interpreted in the time domain and depth converted using the velocity–depth profiles for each seismic trace. A thickness of 12–13 km is present beneath both the upper Mackenzie fan and the middle slope off of Alaska, but the sedimentary prism thins more gradually outboard of the latter region. Mapping of the observed-to-predicted velocities reveals coherent geospatial trends associated with five subregions: the Mackenzie fan; the continental slopes beyond the Mackenzie fan; the abyssal plain; the southwestern Canada Basin; and, the Alpha-Mendeleev magnetic domain. Comparison of the subregional velocity–depth models with published borehole data, and interpretation of the station-specific best-fitting model parameters, suggests that sandstone is not a predominant lithology in any of the five subregions. However, the bulk sand-to-shale ratio likely increases towards the Mackenzie fan, and the model for this subregion compares favourably with borehole data for Miocene turbidites in the eastern Gulf of Mexico. The station-specific results also indicate that Quaternary sediments coarsen towards the Beaufort-Mackenzie and Banks Island margins in a manner that is consistent with the variable history of Laurentide Ice Sheet advance documented for these margins. Lithological factors do not fully account for the elevated velocity–depth trends that are associated with the southwestern Canada Basin and the Alpha-Mendeleev magnetic domain. Accelerated porosity reduction due to elevated palaeo-heat flow is inferred for these regions, which may be related to the underlying crustal types or possibly volcanic intrusion of the sedimentary succession. Beyond exploring the variation of an important physical property in the Arctic Ocean basin, this study provides comparative reference for global studies of seismic velocity, burial history, sedimentary compaction, seismic inversion and overpressure prediction, particularly in mudrock-dominated successions.

Geophysical Journal International↗

Effects of crustal stresses on fluid transport in fractured rock: Case studies from northeastern and southwestern USA

The link between stress and hydrologic properties was examined at two sites that are distinguished by different rock types and different stress states. This investigation is based upon the analysis and interpretation of geophysical logs obtained in water wells at the two locations. At the northeast site (Newark Basin), the hydrologic characteristics of sedimentary rocks are dependent upon the relationship to the current regional stress field of two primary types of orthogonal features that serve as preferential pathways for fluid flow. Subhorizontal bedding-plane partings are highly transmissive near the surface and delineate transversely isotropic fluid flow at shallow depths. With increasing depth, the subhorizontal planes become less dominant and steeply dipping fractures become more influential hydrologically. These high-angle features define anisotropic flow pathways that are preferentially oriented along strike. At the southwest site (west Texas), extrusive rocks are subjected to topographically modified tectonic and gravitational stresses that vary spatially within a valley setting. The attendant changes in stress invariants cause fracture connectivity within the rock mass to systematically increase with depth along the valley flanks, but to remain relatively low in the central valley. The degree of fracture connectivity predicted within this valley configuration is consistent with variations in transmissivity determined at several well locations. In each of these cases, the idealized understanding of the hydrologic system is enhanced by considering the effects of regional and local stresses that act upon the fractured-rock aquifer.

Hydrogeology Journal↗

238U–230Th–226Ra–210Pb–210Po disequilibria constraints on magma generation, ascent, and degassing during the ongoing eruption of Kīlauea

The timescales of magma genesis, ascent, storage and degassing at Kīlauea volcano, Hawai‘i are addressed by measuring 238 U-series radionuclide abundances in lava and tephra erupted between 1982 and 2008. Most analyzed samples represent lavas erupted by steady effusion from Pu‘u ‘Ō‘ō and Kūpahianaha from 1983 to 2008. Also included are samples erupted at the summit in April 1982 and March 2008, along the East Rift Zone at the onset of the ongoing eruption in January 1983, and during vent shifting episodes 54 and 56, at Nāpau crater in January 1997, and Kane Nui O Hamo in June 2007. In general, samples have small (∼4%) excesses of ( 230 Th) over ( 238 U) and ∼3 to ∼17% excesses of ( 226 Ra) over ( 230 Th), consistent with melting of a garnet peridotite source at melting rates between 1 × 10 –3 and 5 × 10 –3 kg m –3 a –1 , and melting region porosity between ∼2 and ∼10%, in agreement with previous studies of the ongoing eruption and historical eruptions. A small subset of samples has near-equilibrium ( 230 Th/ 238 U) values, and thus were generated at higher melting rates. Based on U–Th–Ra disequilibria and Th isotopic data from this and earlier studies, melting processes and sources have been relatively stable over at least the past two centuries or more, including during the ongoing unusually long (>30 years) and voluminous (4 km 3 ) eruption. Lavas recently erupted from the East Rift Zone have average initial ( 210 Pb/ 226 Ra) values of 0·80 ± 0·11 (1σ), which we interpret to be the result of partitioning of 222 Rn into a persistently generated CO 2 -rich gas phase over a minimum of 8 years. This ( 210 Pb) deficit implies an average magma ascent rate of ≤3·7 km a –1 from ∼30 km depth to the surface. Spatter and lava associated with vent-opening episodes erupt with variable ( 210 Pb) deficits ranging from 0·7 to near-equilibrium values in some samples. The samples with near-equilibrium ( 210 Pb/ 226 Ra) are typically more differentiated, suggesting decadal timescales of magma storage in shallow conduits or reservoirs that were not degassing. Lava and spatter samples erupted in the East Rift Zone and at the summit had ( 210 Po) ∼0 at the time of eruption, which results from efficient partitioning of Po into the CO 2 - and SO 2 -rich gas phases during and prior to eruption. Summit ash and Pele’s hair samples from 2008 differ from lava and lapilli samples in that they have elevated initial ( 210 Po), ( 210 Pb/ 226 Ra), and Pb concentrations because of Po condensation on tephra particles, and incorporation of fumarolic Po and Pb into erupted tephra fragments during quenching.

Hawaii↗

Climate change implications for tropical islands: Interpolating and interpreting statistically downscaled GCM projections for management and planning

The potential ecological and economic effects of climate change for tropical islands were studied using output from 12 statistically downscaled general circulation models (GCMs) taking Puerto Rico as a test case. Two model selection/model averaging strategies were used: the average of all available GCMs and the average of the models that are able to reproduce the observed large-scale dynamics that control precipitation over the Caribbean. Five island-wide and multidecadal averages of daily precipitation and temperature were estimated by way of a climatology-informed interpolation of the site-specific downscaled climate model output. Annual cooling degree-days (CDD) were calculated as a proxy index for air-conditioning energy demand, and two measures of annual no-rainfall days were used as drought indices. Holdridge life zone classification was used to map the possible ecological effects of climate change. Precipitation is predicted to decline in both model ensembles, but the decrease was more severe in the “regionally consistent” models. The precipitation declines cause gradual and linear increases in drought intensity and extremes. The warming from the 1960–90 period to the 2071–99 period was 4.6°–9°C depending on the global emission scenarios and location. This warming may cause increases in CDD, and consequently increasing energy demands. Life zones may shift from wetter to drier zones with the possibility of losing most, if not all, of the subtropical rain forests and extinction risks to rain forest specialists or obligates.

Puerto Rico↗

Geometry of the décollement below eastern Bangladesh and implications for seismic hazard

Eastern Bangladesh sits on the seismically active Chittagong-Myanmar fold and thrust belt (CMFB), a north-trending accretionary wedge on the eastern side of the India-Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time-velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east-central Bangladesh. The décollement has an area of 7.25 × 10 4 km 2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual-failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior.

Journal of Geophysical Research: Solid Earth↗

Mountains on Io: High-resolution Galileo observations, initial interpretations, and formation models

During three close flybys in late 1999 and early 2000 the Galileo spacecraft acquired new observations of the mountains that tower above Io's surface. These images have revealed surprising variety in the mountains' morphologies. They range from jagged peaks several kilometers high to lower, rounded structures. Some are very smooth, others are covered by numerous parallel ridges. Many mountains have margins that are collapsing outward in large landslides or series of slump blocks, but a few have steep, scalloped scarps. From these observations we can gain insight into the structure and material properties of Io's crust as well as into the erosional processes acting on Io. We have also investigated formation mechanisms proposed for these structures using finite-element analysis. Mountain formation might be initiated by global compression due to the high rate of global subsidence associated with Io's high resurfacing rate; however, our models demonstrate that this hypothesis lacks a mechanism for isolating the mountains. The large fraction (∼40%) of mountains that are associated with paterae suggests that in some cases these features are tectonically related. Therefore we have also simulated the stresses induced in Io's crust by a combination of a thermal upwelling in the mantle with global lithospheric compression and have shown that this can focus compressional stresses. If this mechanism is responsible for some of Io's mountains, it could also explain the common association of mountains with paterae.

Journal of Geophysical Research E: Planets↗

Water-soluble material on aerosols collected within volcanic eruption clouds

In February and March of 1978, filter samplers mounted on an aircraft were used to collect the aerosol fraction of the eruption clouds from three active Guatemalan volcanoes (Fuego, Pacaya, and Santiaguito). The samples were collected on Teflon (Fluoropore) filters with a nominal pore diameter of 0.5μm. The mass of air sampled by the filters ranged from 0.15 to 6.6 kg. The particulate material collected consisted of fragments of angular silicate ash and droplets of what is interpreted as dilute H 2 SO 4 and HCl. After collection of the samples, each filter was rinsed with 60 ml of distilled-deionized water. Splits of each extract were centrifuged to remove particles greater than or equal to 0.1 μm in diameter, acidified, and analyzed for B, Ba, Be, Ca, Cd, Co, Cu, Fe, Li, Mg, Mn, Mo, Na, Pb, Si, Sr, V, and Zn by inductively coupled plasma—optical emission spectroscopy. Separate splits were analyzed for F and Cl by specific-ion-electrode methods and for U by a fission track technique. The elements dissolved in the aqueous extracts represent components of water-soluble material either formed directly in the eruption cloud or derived from interaction of ash particles and aerosol components of the plume. Calculations of enrichment factors, based upon concentration ratios, showed the elements most enriched in the extracts relative to bulk ash composition were Cd, Cu, V, F, Cl, Zn, and Pb. These elements represent a subset (with the addition of Cl and F) of elements previously reported enriched in atmospheric aerosols in remote regions as well as in volcanic areas. This suggests that some of the enriched elements were widely dispersed as volatile halides emitted from a volcanic source.

Journal of Geophysical Research - Oceans↗

Structure and Paleozoic stratigraphy of a complex of thrust plates in the Fish Creek Reservoir area, south-central Idaho

Permian, Pennsylvanian, Mississippian, Devonian, and Silurian marine rocks of diverse facies are brought together in a complex of six thrust sheets in the Fish Creek Reservoir area on the north edge of the Snake River Plain, Idaho. The lowest structural element, the parautochthon, is made of more than 450 m (1,500 ft) of folded and faulted Devonian miogeosynclinal carbonate rocks present in a 6.5-km 2 (2.5-mi 2 ) window. Along the east margin of the window, a sliver of continental margin transitional carbonate rocks of Early Devonian and Late Silurian age assigned to the Roberts Mountains Formation is thrust over the miogeosynclinal rocks. The window of middle Paleozoic rocks is overridden along the Fish Creek thrust fault by the flysch facies of the Copper Basin Formation, a turbidite-submarine-fan sequence more than 1,000 m (3,300 ft) thick, of Mississippian age. About 4.8 km (3 mi) southwest of the window, about 100 m (300 ft) of deepwater siliceous oceanic facies clastic rocks are exposed, which are assigned with question to the Milligen(?) Formation of Devonian age. These clastic rocks are interpreted to be thrust over the Copper Basin Formation. The highest structural elements are sequences more than 610 m (2,000 ft) thick of interbedded sandy and conglomeratic limestones, quartzites, and conglomerates and interbedded siltstones and argillites of the Wood River Formation of Middle Pennsylvanian to Early Permian age. The Wood River Formation is in thrust contact with the Milligen(?) Formation in the southwest part of the mapped area and with Copper Basin Formation along the west side of Fish Creek Reservoir. All the thrust sheets have moved eastward. The minimum distance moved is estimated from sedimentation models and facies reconstructions to range from perhaps several kilometres for the allochthon of the Roberts Mountains Formation to 48 km (30 mi) for the Milligen(?) Formation allochthon. The principal period of thrusting was post-Early Permian (post-Wood River Formation) and preEocene (pre-Challis Volcanics) and is of probable Sevier age. Middle Paleozoic rocks of the Milligen and Roberts Mountains Formations, however, also may have been involved in an earlier period of thrusting of latest Devonian to earliest Mississippian age related to the Antler orogeny. The thrust sheets were deformed into a northwest-trending dome in late Mesozoic time and were broken by basin-range faults during the Tertiary.

Idaho↗

Recalculated probability of M ≥ 7 earthquakes beneath the Sea of Marmara, Turkey

New earthquake probability calculations are made for the Sea of Marmara region and the city of Istanbul, providing a revised forecast and an evaluation of time-dependent interaction techniques. Calculations incorporate newly obtained bathymetric images of the North Anatolian fault beneath the Sea of Marmara [Le Pichon et al., 2001; Armijo et al., 2002]. Newly interpreted fault segmentation enables an improved regional A.D. 1500-2000 earthquake catalog and interevent model, which form the basis for time-dependent probability estimates. Calculations presented here also employ detailed models of coseismic and postseismic slip associated with the 17 August 1999 M = 7.4 Izmit earthquake to investigate effects of stress transfer on seismic hazard. Probability changes caused by the 1999 shock depend on Marmara Sea fault-stressing rates, which are calculated with a new finite element model. The combined 2004-2034 regional Poisson probability of M≥7 earthquakes is ~38%, the regional time-dependent probability is 44 ± 18%, and incorporation of stress transfer raises it to 53 ± 18%. The most important effect of adding time dependence and stress transfer to the calculations is an increase in the 30 year probability of a M ??? 7 earthquake affecting Istanbul. The 30 year Poisson probability at Istanbul is 21%, and the addition of time dependence and stress transfer raises it to 41 ± 14%. The ranges given on probability values are sensitivities of the calculations to input parameters determined by Monte Carlo analysis; 1000 calculations are made using parameters drawn at random from distributions. Sensitivities are large relative to mean probability values and enhancements caused by stress transfer, reflecting a poor understanding of large-earthquake aperiodicity.

Sea Of Marmara↗

Discrimination of carbon and nitrogen isotopes from milk to serum and vibrissae in Alaska Steller sea lions (Eumetopias jubatus)

Knowledge of diet-tissue stable isotope discrimination is required to properly interpret stable isotope values and to identify possible diet shifts, such as might be expected from nursing through weaning. This study compared ??13C and ??15N of paired serum and vibrissal roots with those of ingested milk (n = 52) from free-ranging Steller sea lion (Eumetopias jubatus (Schreber, 1776)) pups (1-11 months) and juveniles (14-27 months) to estimate diet-tissue discrimination. Mean 15N enrichment from ingested milk to serum was 2.1??? ?? 0.6%??? and ??15N at the root of the vibrissae (representing current growth) were not significantly different from serum values. Milk was enriched for mean 13C by 5.0??? ?? 1.0%??? and 7.3??? ?? 1.2??? relative to serum and vibrissal roots, respectively, which was due to the presence of 13C-depleted lipids in milk. This was confirmed by lipid extraction from a subset of milk and serum samples, resulting in a 5.8??? ?? 1.0??? change only in milk. This study established that vibrissal roots and serum are reflective of a milk diet with approximately 2.0??? 15N enrichment, and vibrissal roots reflect serum and lipid-extracted milk values with approximately 2.0??? 13C enrichment. These discrimination factors are important to establish for stable isotope studies assessing diet shifts. ?? 2008 NRC.

Canadian Journal of Zoology↗

Interpretation of hydraulic conductivity in a fractured-rock aquifer over increasingly larger length dimensions

A comparison of the hydraulic conductivity over increasingly larger volumes of crystalline rock was conducted in the Piedmont physiographic region near Bethesda, Maryland, USA. Fluid-injection tests were conducted on intervals of boreholes isolating closely spaced fractures. Single-hole tests were conducted by pumping in open boreholes for approximately 30 min, and an interference test was conducted by pumping a single borehole over 3 days while monitoring nearby boreholes. An estimate of the hydraulic conductivity of the rock over hundreds of meters was inferred from simulating groundwater inflow into a kilometer-long section of a Washington Metropolitan Area Transit Authority tunnel in the study area, and a groundwater modeling investigation over the Rock Creek watershed provided an estimate of the hydraulic conductivity over kilometers. The majority of groundwater flow is confined to relatively few fractures at a given location. Boreholes installed to depths of approximately 50 m have one or two highly transmissive fractures; the transmissivity of the remaining fractures ranges over five orders of magnitude. Estimates of hydraulic conductivity over increasingly larger rock volumes varied by less than half an order of magnitude. While many investigations point to increasing hydraulic conductivity as a function of the measurement scale, a comparison with selected investigations shows that the effective hydraulic conductivity estimated over larger volumes of rock can either increase, decrease, or remain stable as a function of the measurement scale. Caution needs to be exhibited in characterizing effective hydraulic properties in fractured rock for the purposes of groundwater management.

Hydrogeology Journal↗

Time series of high-resolution images enhances efforts to monitor post-fire condition and recovery, Waldo Canyon fire, Colorado, USA

Interpretations of post-fire condition and rates of vegetation recovery can influence management priorities, actions and perception of latent risks from landslides and floods. In this study, we used the Waldo Canyon fire (2012, Colorado Springs, Colorado, USA) as a case study to explore how a time series (2011–2016) of high-resolution images can be used to delineate burn extent and severity, as well as quantify post-fire vegetation recovery. We applied an object-based approach to map burn severity and vegetation recovery using Worldview-2, Worldview-3 and QuickBird-2 imagery. The burned area was classified as 51% high, 20% moderate and 29% low burn-severity. Across the burn extent, the shrub cover class showed a rapid recovery, resprouting vigorously within 1 year, whereas 4 years post-fire, areas previously dominated by conifers were divided approximately equally between being classified as dominated by quaking aspen saplings with herbaceous species in the understorey or minimally recovered. Relative to using a pixel-based Normalised Difference Vegetation Index (NDVI), our object-based approach showed higher rates of revegetation. High-resolution imagery can provide an effective means to monitor post-fire site conditions and complement more prevalent efforts with moderate- and coarse-resolution sensors.

Colorado↗

Semiempirical model of soil water hysteresis

In order to represent hysteretic soil water retention curves accurately using as few measurements as possible, a new semiempirical model has been developed. It has two postulates related to physical characteristics of the medium, and two parameters, each with a definite physical interpretation, whose values are determined empirically for a given porous medium. One parameter represents the fraction of the pore space that is not subject to Haines jump hysteresis. Its associated postulate is that a single value of this parameter characterizes a given medium, or, equivalently, that the medium is self-similar with respect to the division between hysteretic and nonhysteretic pore space. The second parameter is the effective body-to-neck size ratio of the medium's largest pore. The second postulate specifies a particular relation between the size distributions of pore bodies and of pore necks. Tests of the model show that it provides high-quality optimized fits to measured water content vs. matric pressure wetting curves for a wide variety of media. The parameter values obtained through optimization correlate plausibly with such media properties as uniformity of particles, complexity of structure, and degree of compaction. A practical use of this model is to provide a complete simulated main wetting curve for a medium where only a main drying curve and two points on the wetting curve have been measured. With additional development, it may be possible to do without the two measured wetting points if parameters can be evaluated from other measurements or known properties of the medium.

Soil Science Society of America Journal↗

Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA

An integrated interpretation of field experimental cross-hole radar, tracer, and hydraulic data demonstrates the value of combining time-lapse geophysical monitoring with conventional hydrologic measurements for improved characterization of a fractured-rock aquifer. Time-lapse difference-attenuation radar tomography was conducted during saline tracer experiments at the US Geological Survey Fractured Rock Hydrology Research Site near Mirror Lake, Grafton County, New Hampshire, USA. The presence of electrically conductive saline tracer effectively illuminates permeable fractures or pathways for geophysical imaging. The geophysical results guide the construction of three-dimensional numerical models of ground-water flow and solute transport. In an effort to explore alternative explanations for the tracer and tomographic data, a suite of conceptual models involving heterogeneous hydraulic conductivity fields and rate-limited mass transfer are considered. Calibration data include tracer concentrations, the arrival time of peak concentration at the outlet, and steady-state hydraulic head. Results from the coupled inversion procedure suggest that much of the tracer mass migrated outside the three tomographic image planes, and that solute is likely transported by two pathways through the system. This work provides basic and site-specific insights into the control of permeability heterogeneity on ground-water flow and solute transport in fractured rock.

New Hampshire↗