Search USGS⌕ Search

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.

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

Temporal scaling of hydraulic head and river base flow and its implication for groundwater recharge

Spectral analyses were conducted for hourly hydraulic head ( h ) data observed over a 4‐year period at seven monitoring wells in the Walnut Creek watershed, Iowa. The log power spectral density of the hydraulic head fluctuations versus log frequency ( f ) at all seven wells is shown to have a distinct slope or fractal dimension ( D ), indicating temporal scaling in the time series of water level fluctuations. The fractal dimension of the time series varies from well to well, and the spectrum for the average h over all seven wells has a fractal dimension of 1.46 and Hurst coefficient of 0.54. The log power spectral density of estimated base flow in the Walnut Creek and four other watersheds versus log f is shown to have two distinct slopes with a break in scaling at about 30 days. It is shown that the groundwater recharge process in a basin can be estimated from a head spectrum based on existing theoretical results. Hydraulic head in an aquifer may fluctuate as a fractal in time in response to either a white noise or fractal recharge process, depending on physical parameters (i.e., transmissivity and specific yield) of the aquifer. The recharge process at the Walnut Creek watershed is shown to have a white noise spectrum based on the observed head spectrum.

Water Resources Research↗

Comparison of methods used to estimate conventional undiscovered petroleum resources: World examples

Various methods for assessing undiscovered oil, natural gas, and natural gas liquid resources were compared in support of the USGS World Petroleum Assessment 2000. Discovery process, linear fractal, parabolic fractal, engineering estimates, PETRIMES, Delphi, and the USGS 2000 methods were compared. Three comparisons of these methods were made in: (1) the Neuquen Basin province, Argentina (different assessors, same input data); (2) provinces in North Africa, Oman, and Yemen (same assessors, different methods); and (3) the Arabian Peninsula, Arabian (Persian) Gulf, and North Sea (different assessors, different methods). A fourth comparison (same assessors, same assessment methods but different geologic models), between results from structural and stratigraphic assessment units in the North Sea used only the USGS 2000 method, and hence compared the type of assessment unit rather than the method. In comparing methods, differences arise from inherent differences in assumptions regarding: (1) the underlying distribution of the parent field population (all fields, discovered and undiscovered), (2) the population of fields being estimated; that is, the entire parent distribution or the undiscovered resource distribution, (3) inclusion or exclusion of large outlier fields; (4) inclusion or exclusion of field (reserve) growth, (5) deterministic or probabilistic models, (6) data requirements, and (7) scale and time frame of the assessment. Discovery process, Delphi subjective consensus, and the USGS 2000 method yield comparable results because similar procedures are employed. In mature areas such as the Neuquen Basin province in Argentina, the linear and parabolic fractal and engineering methods were conservative compared to the other five methods and relative to new reserve additions there since 1995. The PETRIMES method gave the most optimistic estimates in the Neuquen Basin. In less mature areas, the linear fractal method yielded larger estimates relative to other methods. A geologically based model, such as one using the total petroleum system approach, is preferred in that it combines the elements of petroleum source, reservoir, trap and seal with the tectono-stratigraphic history of basin evolution with petroleum resource potential. Care must be taken to demonstrate that homogeneous populations in terms of geology, geologic risk, exploration, and discovery processes are used in the assessment process. The USGS 2000 method (7th Approximation Model, EMC computational program) is robust; that is, it can be used in both mature and immature areas, and provides comparable results when using different geologic models (e.g. stratigraphic or structural) with differing amounts of subdivisions, assessment units, within the total petroleum system. ?? 2005 International Association for Mathematical Geology.

Conference Paper↗

Physical modeling of the formation and evolution of seismically active fault zones

Acoustic emission (AE) in rocks is studied as a model of natural seismicity. A special technique for rock loading has been used to help study the processes that control the development of AE during brittle deformation. This technique allows us to extend to hours fault growth which would normally occur very rapidly. In this way, the period of most intense interaction of acoustic events can be studied in detail. Characteristics of the acoustic regime (AR) include the Gutenberg-Richter b -value, spatial distribution of hypocenters with characteristic fractal (correlation) dimension d , Hurst exponent H , and crack concentration parameter P c . The fractal structure of AR changes with the onset of the drop in differential stress during sample deformation. The change results from the active interaction of microcracks. This transition of the spatial distribution of AE hypocenters is accompanied by a corresponding change in the temporal correlation of events and in the distribution of event amplitudes as signified by a decrease of b -value. The characteristic structure that develops in the low-energy background AE is similar to the sequence of the strongest microfracture events. When the AR fractal structure develops, the variations of d and b are synchronous and d = 3 b . This relation which occurs once the fractal structure is formed only holds for average values of d and b . Time variations of d and b are anticorrelated. The degree of temporal correlation of AR has time variations that are similar to d and b variations. The observed variations in laboratory AE experiments are compared with natural seismicity parameters. The close correspondence between laboratory-scale observations and naturally occurring seismicity suggests a possible new approach for understanding the evolution of complex seismicity patterns in nature.

Tectonophysics↗

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↗

Spatial scaling of core and dominant forest cover in the Upper Mississippi and Illinois River floodplains, USA

Different organisms respond to spatial structure in different terms and across different spatial scales. As a consequence, efforts to reverse habitat loss and fragmentation through strategic habitat restoration ought to account for the different habitat density and scale requirements of various taxonomic groups. Here, we estimated the local density of floodplain forest surrounding each of ~20 million 10-m forested pixels of the Upper Mississippi and Illinois River floodplains by using moving windows of multiple sizes (1–100 ha). We further identified forest pixels that met two local density thresholds: ‘core’ forest pixels were nested in a 100% (unfragmented) forested window and ‘dominant’ forest pixels were those nested in a >60% forested window. Finally, we fit two scaling functions to declines in the proportion of forest cover meeting these criteria with increasing window length for 107 management-relevant focal areas: a power function (i.e. self-similar, fractal-like scaling) and an exponential decay function (fractal dimension depends on scale). The exponential decay function consistently explained more variation in changes to the proportion of forest meeting both the ‘core’ and ‘dominant’ criteria with increasing window length than did the power function, suggesting that elevation, soil type, hydrology, and human land use constrain these forest types to a limited range of scales. To examine these scales, we transformed the decay constants to measures of the distance at which the probability of forest meeting the ‘core’ and ‘dominant’ criteria was cut in half ( S 1/2 , m). S 1/2 for core forest was typically between ~55 and ~95 m depending on location along the river, indicating that core forest cover is restricted to extremely fine scales. In contrast, half of all dominant forest cover was lost at scales that were typically between ~525 and 750 m, but S 1/2 was as long as 1,800 m. S 1/2 is a simple measure that (1) condenses information derived from multi-scale analyses, (2) allows for comparisons of the amount of forest habitat available to species with different habitat density and scale requirements, and (3) can be used as an index of the spatial continuity of habitat types that do not scale fractally.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

Constancy of the relation between floc size and density in San Francisco Bay

The size and density of fine-sediment aggregates, or flocs, govern their transport and depositional properties. While the mass and volume concentrations of flocs can be measured directly or by optical methods, they must be determined simultaneously to gain an accurate density measurement. Results are presented from a tidal cycle study in San Francisco Bay, where mass concentration was determined directly, and volume concentration was measured in 32 logarithmically spaced size bins by laser-diffraction methods. The relation between floc size and density is investigated assuming a constant primary particle size and fractal floc dimension. This relation is validated with measurements from several sites throughout San Francisco Bay. The constancy of this relation implies a uniform primary particle size throughout the Bay, as well as uniform aggregation/disaggregation mechanisms (which modify fractal dimension). The exception to the relation is identified during near-bed measurements, when adverted flocs mix with recently resuspended flocs from the bed, which typically have a higher fractal dimension than suspended flocs. The constant relation for suspended flocs simplifies monitoring and numerical modeling of suspended sediment in San Francisco Bay.

California↗

Nature and origin of upper crustal seismic velocity fluctuations and associated scaling properties: Combined stochastic analyses of KTB velocity and lithology logs

The main borehole of the German Continental Deep Drilling Program (KTB) extends over 9000 m into a crystalline upper crust consisting primarily of interlayered gneiss and metabasite. We present a joint analysis of the velocity and lithology logs in an effort to extract the lithology component of the velocity log. Covariance analysis of lithology log, approximated as a binary series, indicates that it may originate from the superposition of two Brownian stochastic processes (fractal dimension 1.5) with characteristic scales of ∼2800 m and ∼150 m, respectively. Covariance analysis of the velocity fluctuations provides evidence for the superposition of four stochastic process with distinct characteristic scales. The largest two scales are identical to those derived from the lithology, confirming that these scales of velocity heterogeneity are caused by lithology variations. The third characteristic scale, ∼20 m, also a Brownian process, is probably related to fracturing based on correlation with the resistivity log. The superposition of these three Brownian processes closely mimics the commonly observed 1/ k decay (fractal dimension 2.0) of the velocity power spectrum. The smallest scale process (characteristic scale ∼1.7 m) requires a low fractal dimension, ∼1.0, and accounts for ∼60% of the total rms velocity variation. A comparison of successive logs from 6900–7140 m depth indicates that such variations are not repeatable and thus probably do not represent true velocity variations in the crust. The results of this study resolve disparity between the differing published estimates of seismic heterogeneity based on the KTB sonic logs, and bridge the gap between estimates of crustal heterogeneity from geologic maps and borehole logs.

Journal of Geophysical Research B: Solid Earth↗

Scale-specific metrics for adaptive generalization and geomorphic classification of stream features

The Richardson plot has been used to illustrate fractal dimension of naturally occurring landscape features that are sensitive to changes in scale or resolution, such as coastlines and river channels. The Richardson method estimates the length of a path by traversing (i.e., “walking”) the path with a specific stride length. Fractal dimension is determined as the slope of the Richardson plot, which shows path length over a range of stride lengths graphed on log-log axes. This paper describes a variant of the Richardson plot referred to as the Scale-Specific Sinuosity (S 3 ) plot. S 3 is defined as negative one times the slope of the Richardson plot for a given stride length. A plot of S 3 against stride length offers a frequency distribution whose area under the curve reflects total sinuosity, and whose points mark the amount of sinuosity contributed to the total sinuosity at each stride length. Mathematical relations of S 3 with fractal dimension and sinuosity for linear features are described. The S 3 metric is demonstrated and discussed for several linear stream features distributed over the conterminous United States. The S 3 metric can help guide the preservation of stream feature sinuosity during cartographic generalization and may assist automated geomorphic classification of river systems.

Conference Paper↗

Nonlinear growth dynamics and the origin of fluctuating asymmetry

The nonlinear, complex nature of biosynthesis magnifies the impacts of small, random perturbations on organism growth, leading to distortions in adaptive allometries and, in particular, to fluctuating asymmetry. These distortions can be partly checked by cell-cell and inter-body part feedback during growth and development, though the latter mechanism also may lead to complex patterns in right-left asymmetry. Stress can be expected to increase the degree to which random growth perturbations are magnified and may also result in disruption of the check mechanisms, thus exaggerating fluctuating asymmetry. The processes described not only provide one explanation for the existence of fluctuating asymmetry and its augmentation under stress, but suggest additional effects of stress as well. Specifically, stress is predicted to lead to decreased fractal dimension of bone sutures and branching structures in animals, and in increased dimension of growth trace patterns such as those found in mollusc shells and fish otoliths and scales. A basic yet broad primer on fractals and chaos is provided as background for the theoretical development in this manuscript.

Genetica↗

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↗

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↗

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&ndash;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↗

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↗

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↗

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↗

Simulating and mapping spatial complexity using multi-scale techniques

A central problem in spatial analysis is the mapping of data for complex spatial fields using relatively simple data structures, such as those of a conventional GIS. This complexity can be measured using such indices as multi-scale variance, which reflects spatial autocorrelation, and multi-fractal dimension, which characterizes the values of fields. These indices are computed for three spatial processes: Gaussian noise, a simple mathematical function, and data for a random walk. Fractal analysis is then used to produce a vegetation map of the central region of California based on a satellite image. This analysis suggests that real world data lie on a continuum between the simple and the random, and that a major GIS challenge is the scientific representation and understanding of rapidly changing multi-scale fields.

International Journal of Geographical Information ↗

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↗

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↗