SEARCH · Search USGS

Results for “Fractals”

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 55 records · Page 3Linked to original sources

USGS1996

Imaging surface contacts: Power law contact distributions and contact stresses in quartz, calcite, glass and acrylic plastic

A procedure has been developed to obtain microscope images of regions of contact between roughened surfaces of transparent materials, while the surfaces are subjected to static loads or undergoing frictional slip. Static loading experiments with quartz, calcite, soda-lime glass and acrylic plastic at normal stresses to 30 MPa yield power law distributions of contact areas from the smallest contacts that can be resolved (3.5 μm 2 ) up to a limiting size that correlates with the grain size of the abrasive grit used to roughen the surfaces. In each material, increasing normal stress results in a roughly linear increase of the real area of contact. Mechanisms of contact area increase are by growth of existing contacts, coalescence of contacts and appearance of new contacts. Mean contacts stresses are consistent with the indentation strength of each material. Contact size distributions are insensitive to normal stress indicating that the increase of contact area is approximately self-similar. The contact images and contact distributions are modeled using simulations of surfaces with random fractal topographies. The contact process for model fractal surfaces is represented by the simple expedient of removing material at regions where surface irregularities overlap. Synthetic contact images created by this approach reproduce observed characteristics of the contacts and demonstrate that the exponent in the power law distributions depends on the scaling exponent used to generate the surface topography.

Tectonophysics↗

USGS1993

Geomorphic observations of rivers in the Oregon Coast Range from a regional reconnaissance perspective

Changes in long profile, gradient, gradient index, pseudo-hypsometric integral, valley incision, and sinuosity fractal dimension for rivers in western Oregon were studied to determine their usefulness in assessing an hypothesis of differential uplift within the Coast Range. All data were gathered from topographic quadrangles and geologic maps, and so the results of this paper are limited to the description of river forms with only limited interpretations. Rivers were naturally divided into three groups, those in the western Coast Range, the eastern Coast Range, and the Klamath and Cascade Ranges, and differences in river morphometry were generalized. Rivers in the western Coast Range had more divergent characteristics, whereas those in the other groups were more similar within each group. For example, long profiles of western Coast Range rivers had diverse forms, but in the other groups, long profiles were similar within the group. Pseudo-hypsometric integrals had a wide scatter in the western Coast Range, but in the other ranges, the pseudo-hypsometric integral was more narrowly defined. Valley entrenchment and high sinuosity were common in the western Coast Range, and were less visible in the eastern Coast Range. These regional generalizations do not hold for the central Coast Range near 44.5°N where the Yaquina and Marys Rivers are located. In the eastern Coast Range, the Marys River is unique in that it has a convex bulge in its long profile, does not have an increasing valley-floor width to valley height ratio downstream, and has a high sinuosity fractal dimension. In the western Coast Range, the Yaquina River is unique in that it has a smoothly decreasing long profile, has a very low pseudo-hypsometric integral, and does not have a decreasing valley-floor width to valley height ratio downstream. There may be many explanations for these observations, but these observations are also consistent with tectonic interpretations that the central western Coast Range is the locus of synclinal tilting and that the entire Oregon coast is experiencing landward tilting and uplift.

Geomorphology↗

USGS2012

Spatial patterns of aquatic habitat richness in the Upper Mississippi River floodplain, USA

Interactions among hydrology and geomorphology create shifting mosaics of aquatic habitat patches in large river floodplains (e.g., main and side channels, floodplain lakes, and shallow backwater areas) and the connectivity among these habitat patches underpins high levels of biotic diversity and productivity. However, the diversity and connectivity among the habitats of most floodplain rivers have been negatively impacted by hydrologic and structural modifications that support commercial navigation and control flooding. We therefore tested the hypothesis that the rate of increase in patch richness (# of types) with increasing scale reflects anthropogenic modifications to habitat diversity and connectivity in a large floodplain river, the Upper Mississippi River (UMR). To do this, we calculated the number of aquatic habitat patch types within neighborhoods surrounding each of the ≈19 million 5-m aquatic pixels of the UMR for multiple neighborhood sizes (1–100 ha). For all of the 87 river-reach focal areas we examined, changes in habitat richness ( R ) with increasing neighborhood length ( L , # pixels) were characterized by a fractal-like power function R = L z ( R 2 > 0.92 ( P < 0.05)). The scaling exponent ( z ) measures the rate of increase in habitat richness with neighborhood size and is related to a fractal dimension. Variation in z reflected fundamental changes to spatial patterns of aquatic habitat richness in this river system. With only a few exceptions, z exceeded the river-wide average of 0.18 in focal areas where side channels, contiguous floodplain lakes, and contiguous shallow-water areas exceeded 5%, 5%, and 10% of the floodplain respectively. In contrast, z was always less than 0.18 for focal areas where impounded water exceeded 40% of floodplain area. Our results suggest that rehabilitation efforts that target areas with <5% of the floodplain in side channels, <5% in floodplain lakes, and/or <10% in shallow-water areas could improve habitat diversity across multiple scales in the UMR.

Upper Mississippi River↗

USGS2007

Ore controls in the Charters Towers goldfield, NE Australia: Constraints from geological, geophysical and numerical analyses

The approach taken in this paper, namely synthesising a wealth of previous information with new data and a genetic model, in combination with integrated numerical analyses, led to new insights into the geological controls on the localisation of auriferous veins and residual prospectivity of the Charters Towers goldfield, NE Australia. The method also has implications for the assessment of other "mature" goldfields worldwide. Despite a number of different ore controls having operated within the Charters Towers goldfield, the controlling factors can be linked to a single genetic model for orogenic, granitoid-hosted lode-gold mineralisation in a brittle deformation regime (D4) of NE-SW to NNE-SSW shortening, under conditions of supralithostatic fluid pressure and low stress difference. Spatial autocorrelation results suggest district-scale alignment of the auriferous veins parallel to and overlapping with the ESE-WNW- to E-W-striking Charters Towers-Ravenswood lineament, a major crustal boundary in the basement to the Ravenswood batholith. At the camp-scale, auriferous veins have abundance and proximity relationships with NW-SE-, NNW-SSE-, NE-SW- and ENE-WSW-oriented lineaments, suggesting that structures that controlled gold deposition in one camp did not necessarily control mineralisation in other camps. Fractal dimensions obtained with the box-counting method range from 1.02 to 1.10, whereas veins in the Charters Towers City camp are characterised by a significantly higher fractal dimension of 1.28. This discrepancy may be taken to imply that most or all outcropping and near-surface deposits within the Charters Towers City camp have been found and that new discoveries are more likely to occur at greater levels of depth, or outside the boundaries of this camp. The new understanding has implications for the assessment of the residual prospectivity of the Charters Towers goldfield, where large areas of prospective rock types and structures (e.g., approximately 40% of the Charters Towers-Ravenswood lineament) are hidden under cover. This parameter space was inaccessible to the historic prospectors and has received relatively little attention from recent explorers. The following steps are suggested for the development of a targeting strategy for lode-gold exploration in areas of the goldfield under cover: (1) identify from geological and geophysical data the ENE-WSW (?? 15??) and NNW-SSE (?? 15??) striking structures and geological boundaries within a 20-km-wide corridor parallel to and centred upon the Charters Towers-Ravenswood lineament, the potential control on gold deposit distribution at the regional- to district-scale, (2) interpret from geological and geophysical data the distribution of pre-Middle Devonian granitoids within these areas that are the preferred host rocks of the payable gold deposits, (3) deduce from geophysical data the ENE-WSW (?? 15??) and NNW-SSE (?? 15??) striking structures that cut or bound the intrusions identified in step 2, (4) locate segments along the structures identified in step 3 that deviate most from the geometry of a straight line (e.g., potential bends or splays) and/or intersect other structures or geological contacts, or both, and (5) define and rank potential targets within the prospective areas identified in step 4 and systematically test the best ones. ?? 2006 Elsevier B.V. All rights reserved.

Ore Geology Reviews↗

USGS2002

Description and comparison of geologic maps with FRAGSTATS - A spatial statistics program

FRAGSTATS is a public-domain GIS implementation of a set of spatial statistics that address a fundamental problem in GIS applications, description and comparison of maps. The spatial statistics from the 1:2,500,000-scale United States geologic map of Nevada, the central United States, and the northeastern United States quantify the differences in complexity and variability between these three geologic terranes. Nevada is defined by a large number of patches of small size and low size variability, whereas the Central area has a small number of patches with smaller relative size variability. All three areas have similar map-unit shape complexity with Nevada having the lowest. Based on the density of edges, the areas can be ranked from highest to lowest, as Nevada, Central, and Northeast. The Shannon diversity index ranks the areas from highest to lowest, as Northeast, Nevada, and Central, but the Shannon evenness index ranks them from highest to lowest, as Northeast, Central, and Nevada. These rankings may reflect the influence of folding in the Northeast and Central areas as opposed to basin and range extension in Nevada. The core areas statistic ranks the areas for spatial accuracy from highest to lowest, as Central, Northeast, and Nevada, with Northeast and Nevada being similar. For a scale comparison, the FRAGSTATS statistics quantify the increased complexity and spatial accuracy that is inherent in going from small- to larger-scale maps. For example for 1:2,500,000-1:500,000-scale maps of Nevada, respectively, the area weighted fractal dimension increase from 1.1 to 1.18, and the total core areas index almost doubles from 39.09 to 63.38. In addition, the fractal dimensions discriminate gross lithology and tectonic terranes. ?? 2002 Elsevier Science Ltd. All rights reserved.

Computers & Geosciences↗

USGS2019

Experimental study on the impact of thermal maturity on shale microstructures using hydrous pyrolysis

Hydrous pyrolysis was applied to four low-maturity aliquots from the Utica, Excello, Monterey, and Niobrara Shale Formations in North America to create artificial maturation sequences, which could be used to study the impact of maturation on geochemical and microstructural properties. Modified Rock-Eval pyrolysis, reflectance, organic petrology, and Fourier transform infrared spectroscopy (FTIR) were employed to analyze their geochemical properties, while gas adsorption (CO 2 and N 2 ) was used to characterize their pore structures (pores < 200 nm). Organic petrography using white and blue light (fluorescence) before and after hydrous pyrolysis showed that amorphous organic matter cracked into solid bitumen, oil, and gas during hydrous pyrolysis. A reduction of the CH 2 /CH 3 ratio in hydrous pyrolysis residues was observed from FTIR analysis. Rock-Eval pyrolysis showed that kerogens in the four samples were dissimilar, and hydrous pyrolysis residues showed smaller hydrogen index and Sh2 values than starting materials. Results from CO 2 and N 2 gas adsorption analysis showed that pore structures (micropore volume, micropore surface area, meso-macropore volume, and meso-macropore surface area) changed significantly during hydrous pyrolysis. However, changes in pore structure were dissimilar among the four samples, which was attributed to different activation energies of organic matter. A thermodynamic fractal model showed a decrease in fractal dimensions of Utica, Monterey, and Excello after hydrous pyrolysis, indicating a decrease in surface roughness. The pore size heterogeneity in the Utica sample increased as hydrous pyrolysis temperature increased, whereas the pore size heterogeneity distributions in the Monterey and Excello decreased based on the N 2 adsorption data.

Energy & Fuels↗

USGS2006

Effects of wildfire and permafrost on soil organic matter and soil climate in interior Alaska

The influence of discontinuous permafrost on ground-fuel storage, combustion losses, and postfire soil climates was examined after a wildfire near Delta Junction, AK in July 1999. At this site, we sampled soils from a four-way site comparison of burning (burned and unburned) and permafrost (permafrost and nonpermafrost). Soil organic layers (which comprise ground-fuel storage) were thicker in permafrost than nonpermafrost soils both in burned and unburned sites. While we expected fire severity to be greater in the drier site (without permafrost), combustion losses were not significantly different between the two burned sites. Overall, permafrost and burning had significant effects on physical soil variables. Most notably, unburned permafrost sites with the thickest organic mats consistently had the coldest temperatures and wettest mineral soil, while soils in the burned nonpermafrost sites were warmer and drier than the other soils. For every centimeter of organic mat thickness, temperature at 5cm depth was about 0.5??C cooler during summer months. We propose that organic soil layers determine to a large extent the physical and thermal setting for variations in vegetation, decomposition, and carbon balance across these landscapes. In particular, the deep organic layers maintain the legacies of thermal and nutrient cycling governed by fire and revegetation. We further propose that the thermal influence of deep organic soil layers may be an underlying mechanism responsible for large regional patterns of burning and regrowth, detected in fractal analyses of burn frequency and area. Thus, fractal geometry can potentially be used to analyze changes in state of these fire prone systems. ?? 2006 Blackwell Publishing Ltd.

Global Change Biology↗

USGS2003

A comparison of measures of riverbed form for evaluating distributions of benthic fishes

A method to quantitatively characterize the bed forms of a large river and a preliminary test of the relationship between bed-form characteristics and catch per unit area of benthic fishes is presented. We used analog paper recordings of bathymetric data from the Missouri River and fish data collected from 1996 to 1998 at both the segment (???101-102-km) and macrohabitat (???10-1-100-km) spatial scales. Bed-form traces were transformed to digital data with image analysis software. The slope, mean residual, and SD of the residuals of the regression of depth versus distance along the bottom, as well as mean depth, were estimated for each trace. These four metrics were compared with sinuosity, fractal dimension, critical scale, and maximum mean angle for the same traces. Mean depth and sinuosity differed among segments and macrohabitats. Fractal-based measures of the relative depth of bottom troughs (critical scale) and smoothness (maximum mean angle) differed among segments. Statistics-based measures of the relative depth of bottom troughs (mean residual) and smoothness (SD of the residuals) differed among macrohabitats. Sites with shovelnose sturgeon Scaphirhynchus platorynchus were shallower and smoother than sites without shovelnose sturgeon. When compared with sites without sicklefin chub Macrhybopsis meeki, sites with sicklefin chub were shallower, had shallower troughs, and sloped more out of the flow of the river. Sites with sturgeon chub M. gelida were shallower, had shallower troughs, and were smoother than sites without sturgeon chub. Sites with and without channel catfish Ictalurus punctatus did not differ for any bed-form variables measured. Nonzero shovelnose sturgeon density increased with depth, whereas nonzero sturgeon chub density decreased with depth. Indices of bed-form structure demonstrated potential for describing the distribution and abundance of Missouri River benthic fishes. The observed fish patterns, though limited, provide valuable direction for future research into the habitat preferences of these fishes.

North American Journal of Fisheries Management↗

USGS2004

Analysis of predator movement in prairie landscapes with contrasting grassland composition

Mammalian predation influences waterfowl breeding success in the U.S. northern Great Plains, yet little is known about the influence of the landscape on the ability of predators to find waterfowl nests. We used radiotelemetry to record nightly movements of red foxes ( Vulpes vulpes ) and striped skunks ( Mephitis mephitis ) in two 41.4-km 2 study areas in North Dakota. Study areas contained either 15–20% grassland (low grassland composition) or 45–55% grassland (high grassland composition). Grasslands included planted cover, pastureland, and hayland. We predicted that the type and composition of cover types in the landscape would influence both predator movement across the landscape (as measured by the fractal dimension and displacement ratio) as well as localized movement (as measured by the rate of movement and turning angle between locations) within patches of different cover types. Red fox movements were straighter (lower fractal dimensions and higher displacements) across landscapes with a low grassland composition, indicating directed movement between the more isolated patches of planted cover. Striped skunk movements did not differ between landscape types, illustrating their movement along wetland edges, which had similar compositions in both landscape types. The high variability in turning angles by red fox in planted cover and pastureland in both landscape types is consistent with restricted-area foraging. The high rate of movement by red foxes in planted cover and by striped skunks in wetland edges suggests that spatial memory may influence movement patterns. Understanding the behavior of predators in fragmented prairie landscape is essential for managing breeding habitat for grassland birds and for predicting the spatial and temporal dynamics of predators and their prey.

North Dakota↗

USGS1999

Calculation of broadband time histories of ground motion: Comparison of methods and validation using strong-ground motion from the 1994 Northridge earthquake

This article compares techniques for calculating broadband time histories of ground motion in the near field of a finite fault by comparing synthetics with the strong-motion data set for the 1994 Northridge earthquake. Based on this comparison, a preferred methodology is presented. Ground-motion-simulation techniques are divided into two general methods: kinematic- and composite-fault models. Green's functions of three types are evaluated: stochastic, empirical, and theoretical. A hybrid scheme is found to give the best fit to the Northridge data. Low frequencies (< 1 Hz) are calculated using a kinematic-fault model and a 3D finite-difference code to propagate energy through a realistic 3D velocity structure. High frequencies (> 1 Hz) are calculated using a composite-fault model with a fractal subevent size distribution and stochastic, bandlimited, white-noise Green's functions. At frequencies below 1 Hz, theoretical elastic-wave-propagation synthetics introduce proper seismic-phase arrivals of body waves and surface waves. The 3D velocity structure more accurately reproduces record durations for the deep sedimentary basin structures found in the Los Angeles region. At frequencies above 1 Hz, scattering effects become important and wave propagation is more accurately represented by stochastic Green's functions. A fractal subevent size distribution for the composite fault model ensures an ω −2 spectral shape over the entire frequency band considered (0.1-20 Hz).

California↗

USGS2009

Scaling the effects of moose browsing on forage distribution, from the geometry of plant canopies to landscapes

Landscape heterogeneity influences large herbivores by altering their feeding rates, but as herbivores attempt to maximize feeding rates they also create spatial heterogeneity by altering plant growth. Herbivore feeding rates thus provide a quantitative link between the causes and consequences of spatial heterogeneity in herbivore-dominated ecosystems. The fractal geometry of plant canopies determines both the density and mass of twigs available to foraging herbivores. These properties determine a threshold distance between plants (d * ) that distinguishes the mechanisms regulating herbivore intake rates. When d * is greater than the actual distance between plants (d), intake is regulated by the rate of food processing in the mouth. But when d * < d, intake is regulated by the rate at which the herbivore encounters new plants. Alterations to plant geometry due to past browsing could change the rate at which herbivores encounter and process bites of plant tissue, modify d * relative to d, and thus change intake rates and the distribution of mechanisms regulating it across landscapes. We measured changes in the geometry of aspen (Populus tremuloides) and balsam fir (Abies balsamea) saplings along gradients of moose browsing from 2001 to 2005 at Isle Royale National Park, Michigan, USA. For aspen saplings, fractal dimension of bite density, bite mass, and forage biomass responded quadratically to increasing moose browsing and were greatest at -3-4 g-g.m-2.yr" 1 consumption. For balsam fir, in contrast, these same measures declined steadily with increasing moose browsing. The different responses of plant canopies to increased browsing altered d * around plants. In summer, d * > d for aspen saplings at all prior consumption levels. Food processing therefore regulated summer moose feeding rates across our landscapes. In winter, changes in bite mass due to past browsing were sufficient to cause d * < d for aspen and balsam fir. Therefore, travel velocity and food processing jointly regulated intake rate during winter. Browsing-induced changes in the small-scale geometry of plant canopies can determine intake rate at larger spatial scales by changing d * relative to d and, hence, which mechanisms determine intake rate, essentially altering how herbivores sense the distribution of their food resources. ?? 2009 by the Ecological Society of America.

Ecological Monographs↗

USGS2009

A multiscale analysis of coral reef topographic complexity using lidar-derived bathymetry

Coral reefs represent one of the most irregular substrates in the marine environment. This roughness or topographic complexity is an important structural characteristic of reef habitats that affects a number of ecological and environmental attributes, including species diversity and water circulation. Little is known about the range of topographic complexity exhibited within a reef or between different reef systems. The objective of this study was to quantify topographic complexity for a 5-km x 5-km reefscape along the northern Florida Keys reef tract, over spatial scales ranging from meters to hundreds of meters. The underlying dataset was a 1-m spatial resolution, digital elevation model constructed from lidar measurements. Topographic complexity was quantified using a fractal algorithm, which provided a multi-scale characterization of reef roughness. The computed fractal dimensions ( D ) are a measure of substrate irregularity and are bounded between values of 2 and 3. Spatial patterns in D were positively correlated with known reef zonation in the area. Landward regions of the study site contain relatively smooth ( D ≈ 2.35) flat-topped patch reefs, which give way to rougher ( D ≈ 2.5), deep, knoll-shaped patch reefs. The seaward boundary contains a mixture of substrate features, including discontinuous shelf-edge reefs, and exhibits a corresponding range of roughness values (2.28 ≤ D ≤ 2.61).

Journal of Coastal Research↗

USGS2026

Self-similar vent clustering in the Pinacate volcanic field (Sonora, Mexico and Arizona, United States)

The spatial distribution of vents in a distributed volcanic field may provide clues about the geometry of the volcanic field’s plumbing system. Distributed volcanism is characterized by magma feeders that are independent of one another, scattered over wide areas, and active for relatively short times (typically less than one hundred years). On the other hand, focused volcanism results in large volcanic edifices, forming central or shield volcanoes that are active for tens of thousands of years or longer. Whether distributed or focused, the form of volcanism may depend on several factors: magma fluxes at depth; magma composition, temperature, and viscosity; and crustal stress and strain states. Over time, some volcanic fields may shift from focused to distributed volcanism. Located in northwestern Sonora, Mexico, and southwestern Arizona, United States, the Pinacate volcanic field is made up of lava flows, 8 maars, and more than 400 cinder cones. It has been active for approximately 2 million years. The oldest features in the field are large-volume lava flows that constructed Sierra El Pinacate (known locally and referred to herein as the Santa Clara shield volcano) and formed the base of the field. Later activity took the form of monogenetic cones and maars, produced by high-alkali basalts with an ocean-island basalt chemical signature. The plumbing system of the Pinacate volcanic field, along with its relations with the structural setting of the area, are analyzed here in terms of the spatial distribution of vents, vent morphology, and main structural trends of the morphologic lineaments in the area. The spatial analysis of vents identified three clusters. The largest cluster (cluster 1) is in the north-northwest portion of the volcanic field, whereas the smallest cluster (cluster 2), formed by younger cones, is on the east side of the volcanic field. The southern cluster (cluster 3) is mainly located over the Santa Clara shield volcano. Self-similar clustering analysis was performed on the whole volcanic field, as well as on the two larger clusters (clusters 1 and 3); cluster 2 was not analyzed individually due to the small number of vents. The self-similar clustering analysis indicated that the dataset containing all the vents in the volcanic field and the dataset containing only the vents in cluster 1 (the largest cluster) have magma storage at approximately 17 kilometers (km) depth. The vents of cluster 3 indicate storage at approximately 9 km depth. The fractal exponent characterizing the self-similar clustering of cluster 3 (1.5905) is lower than the fractal exponent of cluster 1 (1.7252), implying that cluster 1 formed after cluster 3. The available absolute ages of volcanic products in the Pinacate volcanic field indicate a middle to early Pleistocene age for the vents in cluster 3, and younger ages (late Pleistocene to Holocene) for cluster 2. The cone morphology is equivocal for defining relative ages of the clusters; nevertheless, the spatial distribution of the best-preserved cones implies that northwest-southeast, north-south, and northeast-southwest trends are linked to the development of the volcanic field. Lineament mapping, as well as the shape, elongation, and orientation of inferred feeder dikes, indicates that the Pinacate volcanic field formed along a northwest-southeast, regional-scale dextral transtensional shear zone and southwestward shallowing of the Mohorovičić discontinuity to 18–20 km.

Arizona, Sonora↗

USGS2016

Complexity Theory

A complex system consists of many interacting parts, generates new collective behavior through self organization, and adaptively evolves through time. Many theories have been developed to study complex systems, including chaos, fractals, cellular automata, self organization, stochastic processes, turbulence, and genetic algorithms.

Book chapter↗

USGS2002

Structure and physical characteristics of pumice from the climactic eruption of Mount Mazama (Crater Lake), Oregon

The vesicularity, permeability, and structure of pumice clasts provide insight into conditions of vesiculation and fragmentation during Plinian fall and pyroclastic flow-producing phases of the ~7,700 cal. year B.P. climactic eruption of Mount Mazama (Crater Lake), Oregon. We show that bulk properties (vesicularity and permeability) can be correlated with internal textures and that the clast structure can be related to inferred changes in eruption conditions. The vesicularity of all pumice clasts is 75–88%, with >90% interconnected pore volume. However, pumice clasts from the Plinian fall deposits exhibit a wider vesicularity range and higher volume percentage of interconnected vesicles than do clasts from pyroclastic-flow deposits. Pumice permeabilities also differ between the two clast types, with pumice from the fall deposit having higher minimum permeabilities (~5×10 –13 m 2 ) and a narrower permeability range (5–50×10 –13 m 2 ) than clasts from pyroclastic-flow deposits (0.2–330×10 –13 m 2 ). The observed permeability can be modeled to estimate average vesicle aperture radii of 1–5 µm for the fall deposit clasts and 0.25–1 µm for clasts from the pyroclastic flows. High vesicle number densities (~10 9 cm –3 ) in all clasts suggest that bubble nucleation occurred rapidly and at high supersaturations. Post-nucleation modifications to bubble populations include both bubble growth and coalescence. A single stage of bubble nucleation and growth can account for 35–60% of the vesicle population in clasts from the fall deposits, and 65–80% in pumice from pyroclastic flows. Large vesicles form a separate population which defines a power law distribution with fractal dimension D=3.3 (range 3.0–3.5). The large D value, coupled with textural evidence, suggests that the large vesicles formed primarily by coalescence. When viewed together, the bulk properties (vesicularity, permeability) and textural characteristics of all clasts indicate rapid bubble nucleation followed by bubble growth, coalescence and permeability development. This sequence of events is best explained by nucleation in response to a downward-propagating decompression wave, followed by rapid bubble growth and coalescence prior to magma disruption by fragmentation. The heterogeneity of vesicle sizes and shapes, and the absence of differential expansion across individual clasts, suggest that post-fragmentation expansion played a limited role in the development of pumice structure. The higher vesicle number densities and lower permeabilities of pyroclastic-flow clasts indicate limited coalescence and suggest that fragmentation occurred shortly after decompression. Either increased eruption velocities or increased depth of fragmentation accompanying caldera collapse could explain compression of the pre-fragmentation vesiculation interval.

Oregon↗

USGS2005

Temporal variations and scaling of streamflow and baseflow and their nitrate-nitrogen concentrations and loads

The patterns of temporal variations of precipitation (P), streamflow (SF) and baseflow (BF) as well as their nitrate-nitrogen (nitrate) concentrations (C) and loads (L) from a long-term record (28 years) in the Raccoon River, Iowa, were analyzed using variogram and spectral analyses. The daily P is random but scaling may exist in the daily SF and BF with a possible break point in the scaling at about 18 days and 45 days, respectively. The nitrate concentrations and loads are shown to have a half-year cycle while daily P, SF, and BF have a one-year cycle. Furthermore, there may be a low-frequency cycle of 6-8 years in C. The power spectra of C and L in both SF and BF exhibit fractal 1/f scaling with two characteristic frequencies of half-year and one-year, and are fitted well with the spectrum of the gamma distribution. The nitrate input to SF and BF at the Raccoon watershed seems likely to be a white noise process superimposed on another process with a half-year and one-year cycle. ?? 2005 Elsevier Ltd. All rights reserved.

Advances in Water Resources↗

USGS2011

Investigating the complex interface where bedrock transforms to regolith

The interface where bedrock transforms to regolith is not planar but rather has a roughness that varies with the scale of observation. The complexity of this surface is manifested in both element-depth and fragment size-depth distributions and may sometimes be related to the longitudinal profiles of watershed streams. The fractal nature of the bedrock-regolith interface means that the interface has a “thickness” which is >20 m in two ridgetop examples from Pennsylvania and Puerto Rico. Such weathering thicknesses, modeled as a function of one-dimensional fluid flow, are affected by the balance between rates of weathering and erosion. One-dimensional models are consistent with weathering advance rates that vary with equilibrium solubility and porefluid velocities (and not reaction kinetics). However, fluid flow is not strictly downward and one-dimensional. Permeability of regolith changes as particle size and bulk density changes with depth. Thus, both downward and lateral flow occurs especially at reaction fronts where reactions change permeability. The rate of weathering advance is, therefore, affected by the 3-dimensional distribution of reaction zones that affect permeability across the watershed. Quantitative models of such phenomena over a range of spatial and temporal scales are needed.

Applied Geochemistry↗

USGS2018

Porosity of the Marcellus Shale: A contrast matching small-angle neutron scattering study

Neutron scattering techniques were used to determine the effect of mineral matter on the accessibility of water and toluene to pores in the Devonian Marcellus Shale. Three Marcellus Shale samples, representing quartz-rich, clay-rich, and carbonate-rich facies, were examined using contrast matching small-angle neutron scattering (CM-SANS) at ambient pressure and temperature. Contrast matching compositions of H 2 O, D 2 O and toluene, deuterated toluene were used to probe open and closed pores of these three shale samples. Results show that although the mean pore radius was approximately the same for all three samples, the fractal dimension of the quartz-rich sample was higher than for the clay-rich and carbonate-rich samples, indicating different pore size distributions among the samples. The number density of pores was highest in the clay-rich sample and lowest in the quartz-rich sample. Contrast matching with water and toluene mixtures shows that the accessibility of pores to water and toluene also varied among the samples. In general, water accessed approximately 70–80% of the larger pores (>80 nm radius) in all three samples. At smaller pore sizes (~5–80 nm radius), the fraction of accessible pores decreases. The lowest accessibility to both fluids is at pore throat size of ~25 nm radii with the quartz-rich sample exhibiting lower accessibility than the clay- and carbonate-rich samples. The mechanism for this behaviour is unclear, but because the mineralogy of the three samples varies, it is likely that the inaccessible pores in this size range are associated with organics and not a specific mineral within the samples. At even smaller pore sizes (~<2.5 nm radius), in all samples, the fraction of accessible pores to water increases again to approximately 70–80%. Accessibility to toluene generally follows that of water; however, in the smallest pores (~<2.5 nm radius), accessibility to toluene decreases, especially in the clay-rich sample which contains about 30% more closed pores than the quartz- and carbonate-rich samples. Results from this study show that mineralogy of producing intervals within a shale reservoir can affect accessibility of pores to water and toluene and these mineralogic differences may affect hydrocarbon storage and production and hydraulic fracturing characteristics

International Journal of Coal Geology↗