土壤有机碳和全氮含量及其与海拔、植被和气候要素的关系——以祁连山中段北坡为研究对象

调查分析祁连山北坡土壤有机碳和全氮含量随海拔变化的趋势及其与气候要素和植被的关系。结果显示:就土壤有机碳含量,在0～5cm、5～1 5cm和1 5～30cm土层,低海拔(2200m)和高海拔(3600m)处较低,中间海拔(3000～3500m)处较高;灌丛草甸森林高寒草甸干旱草原荒漠草原(P0.05)。就土壤全氮含量,0～5cm和5～1 5cm土层,3400m和3500m处较高,2200m和2800m处较低;1 5～35cm土层,3400m处较高、2200m和2800m处较低;0～5cm土层,灌丛草甸高寒草甸森林干旱草原荒漠草原;5～1 5和1 5～35cm土层,灌丛草甸森林高寒草甸干旱草原荒漠草原(P0.05)。就土壤碳氮比,在0～5cm、5～1 5cm和15～30 cm土层,3000m和3200m处最高;2300m和2800m处最低;森林灌丛草甸高寒草甸干旱草原或荒漠草原。土壤有机碳和全氮含量及碳氮比总体上都随年均温度增加而降低,随年降水量增加而增加。不考虑海拔差异,0～5cm、5～1 5cm和1 5～30cm土层土壤有机碳和全氯含量相关系数较高;考虑海拔差异,在不同土层和海拔的差异较大。     

Comparison of available soil water capacity estimated from topography and soil series information

We present a simple and generalized method to predict Available Soil Water Capacity (ASWC-TOP) for a given area using a topographic index, defined as ln(/tan), where is the upslope area draining past a certain point per unit width of slope, and is the local surface slope angle. The estimated results (ASWC-TOP) were then compared with the available soil water capacity calculated from soil series information provided by Soil Conservation Service, U.S. Department of Agriculture (ASWC-SCS). The model implementation was tested with three study cases: the Seeley-Swan valley, Montana, with pixel resolutions of 100 m and 1 km, respectively; and the state of Montana, U.S.A., with a pixel resolution of 1 km. A linear relationship exists between ASWC-SCS and ln(/tan). Standard errors between ASWC-TOP and ASWC-SCS were about 4.4 cm in the Seeley-Swan valley and 5.5 cm in the state. The number of pixels with absolute residuals 4 cm between ASWC-TOP and ASWC-SCS accounted for 68.2, 64.4, and 51.9% for the valley 100 m, valley 1 km, and the state respectively. Some of the mismatches between ASWC-TOP and ASWC-SCS may indicate an improvement using this method compared to existing data because the topographic method reflects the higher spatial variation of the inputs. The increasing availability of digital elevation data, at various resolutions, may provide an alternative to soil series for estimating ASWC. The accuracy of ASWC-TOP depends on estimation of mean and maximum ASWC for a study area.     

Spatial variability of soil metabolic rate along a dryland elevation gradient

A general framework of ecosystem hotspots suggests variation in soil metabolic activity can be understood through the relative distribution and intensity of patches of disproportionately high ecosystem process rates. To better understand the causes of soil metabolic spatial variability and the variation in ecosystem hotspots we quantified soil respiration (R) spatial heterogeneity across a network of seven sites spanning a 2,489 m elevation gradient in the Santa Rosa Mountains of Southern California. At each site, soil samples were collected from 0–5 and 5–15 cm soil depths at 2 m intervals along three 100 m transects. Each soil sample was analyzed for organic matter content (SOM) and was incubated at 40% water holding capacity for 20 days. R was measured at days 5 and 20. Strong contrasts were observed between the relationships of soil physical variables and R at scales of individual landscapes and the whole region. Notably, the relationship between SOM and R was positive within individual landscapes and negative across the entire region. Plant canopy microenvironments were associated with elevated SOM and R relative to the interspaces. This microenvironment effect on R was reduced by elevation, incubation interval, and soil depth. Geostatistical analyses conducted individually for each site identified an increasing range of autocorrelation from 2 to 10 m and a decreasing proportion of variation that was included in this range with elevation. These results suggest hotspots increase in size but decrease in intensity with elevation thereby creating a maximum hotspot effect at middle elevations.     

Climatic and topographic controls on patterns of fire in the southern and central Appalachian Mountains,USA

Climate and topography are two important controls on spatial patterns of fire disturbance in forests globally, via their influence on fuel moisture and fuel production. To assess the influences of climate and topography on fire disturbance patterns in a temperate forest region, we analyzed the mapped perimeters of fires that burned during 1930–2003 in two national parks in the eastern United States. These were Great Smoky Mountains National Park (GSMNP) in the southern Appalachian Mountains and Shenandoah National Park (SNP) in the central Appalachian Mountains. We conducted GIS analyses to assess trends in area burned under differing climatic conditions and across topographic gradients (elevation, slope position, and aspect). We developed a Classification and Regression Tree model in order to further explore the interactions between topography, climate, and fire. The results demonstrate that climate is a strong driver of both spatial and temporal patterns of wildfire. Fire was most prevalent in the drier SNP than the wetter GSMNP, and during drought years in both parks. Topography also influenced fire occurrence, with relatively dry south-facing aspects, ridges, and lower elevations burning most frequently. However, the strength of topographic trends varied according to the climatic context. Weaker topographic trends emerged in the drier SNP than GSMNP, and during low-PDSI (dry) years than high-PDSI (wet) years in both parks. The apparent influence of climate on the spatial patterning of fire suggests a more general concept, that disturbance-prone landscapes exhibit weaker fine-scale spatial patterning of disturbance than do less disturbance-prone landscapes.     

Post-fire aspen seedling recruitment across the Yellowstone (USA) Landscape

Landscape patterns of quaking aspen (Populus tremuloides) seedling occurrence and abundance were studied after a rare recruitment event following the 1988 fires in Yellowstone National Park, Wyoming, USA. Belt transects (1 to 17 km in length, 4 m width) along 18 foot trails were surveyed for aspen seedlings on the subalpine plateau of the Park, along gradients of elevation and geologic substrate, during the summer of 1996. Aspen seedling presence and density were characterized as a function of elevation, geologic substrate, slope, aspect, vegetation/cover type, presence of burned forest, and distance to nearest adult aspen stand. Presence of aspen seedlings was best predicted by the incidence of burned forest and proximity to adult aspen; aspen seedlings were only found in burned forest and were more likely to occur closer to adult aspen clones. When tested against independent data collected in 1997, the logistic regression model for aspen seedling presence performed well (overall accuracy = 73%, Tau_p = 0.41). When present, variation in aspen seedling density at local scales (≤ 200 m) was largely explained by elevation, with higher densities observed at lower elevations. At broad scales (> 1 km), seedling density was a function of cover type, elevation, aspect, slope, and burn severity, with greater seedling density in more severely burned forested habitats on southerly, shallow slopes at lower elevations. Aspen seedling densities ranged from 0 to 46,000 seedlings/ha with a median density of 2,000/ha on sites where they occurred. Aspen seedlings were most abundant in the south central and southwest central regions of the park, approximately an order of magnitude less abundant in the southeast region, and nearly absent in the north central area. Establishment of new aspen stands on Yellowstone's subalpine plateau would represent a substantial change in the landscape. However, the long-term fate of these postfire aspen seedlings is not known. This revised version was published online in August 2006 with corrections to the Cover Date.     

Response of high mountain landscape to topographic variables: Central pyrenees

An objective method for inductively modelling the distribution of mountain land units using GIS managed topographic variables is presented. The landscape of a small high mountain catchment in the Spanish Pyrenees, covered with grassland, was classified in ten land units by hierarchical agglomerative clustering, using a sample of 194 random plots, in which classes of vegetation, soils and landforms were defined. Additionally, seven layers of topographic variables (altitude, slope angle, aspect, solar radiation, topographic wetness index, specific catchment area, and regolith thickness) were created from a Digital Elevation Model. The affinity of each land unit to the topographic variables was calculated using Binary Discriminant Analysis (BDA), after dichotomising the latter around their mean values. Then, the distribution of each land unit was predicted by boolean operations combining step by step distributions for the seven topographic variables ordered, for each unit, after the absolute values of the Haberman’s residuals in BDA. The predicted distributions were tested (χ²) against that of the observed sampling plots. From the original ten land units, the distributions of eight of them were successfully predicted (four are related to the slope sequence, two reflect the water accumulation in the soil, and two respond to geomorphic processes) while the remaining two had to be rejected. Part of the catchment (39%) was not assigned to any land unit, probably because more distributed variables accounting for snow distribution are necessary.     

Relating tree vigour to the soil and landscape characteristics of an olive orchard in a marly area of southern Spain

The relationships between tree vigour and soil and landscape properties of a rainfed olive orchard located on a marly slope of southern Spain were studied. The different properties and thickness of three soil layers (viz. the secondary carbonate-free layer, the solum, and the layer penetrated by roots), and the properties of the 0–30 cm layer, were used as soil variables. Olive tree vigour and soil properties were found to be strongly correlated. Much of the variability in tree vigour was related to soil morphology. The thickness of the root zone and the average CaCO₃ content in the secondary carbonate-free layer and the upper 30 cm of soil accounted for ∼58 and ∼68%, respectively, of the variance in canopy volume. Tree vigour was related to some properties of the root zone, such as clay content and cation exchange capacity, expressed on a surface basis. Topographic properties such as distance to summit, slope gradient and profile curvature were also found to be significantly correlated with tree vigour, but had less influence than intrinsic soil properties.     

Spatial dynamics of a gypsy moth defoliation outbreak and dependence on habitat characteristics

Forest insects cause defoliation disturbances with complex spatial dynamics. These are difficult to measure but critical for models of disturbance risk that inform forest management. Understanding of spatial dynamics has lagged behind other disturbance processes because traditional defoliation sketch map data often suffered from inadequate precision or spatial resolution. We sought to clarify the influence of underlying habitat characteristics on outbreak patterns by combining forest plots, GIS data and defoliation intensity maps modeled from Landsat imagery. We quantified dependence of defoliation on spatial patterns of host abundance, phenology, topography, and pesticide spray for a recent gypsy moth outbreak (2000–2001), in a mixed deciduous forest in western Maryland, USA. We used semivariograms and hierarchical partitioning to quantify spatial patterns and variable importance. Habitat characteristics from plot data explained 21 % of defoliation variance in 2000 from tree density, phenological asynchrony, pesticide spray status, and landform index and 34 % of the variance in 2001 from previous-year defoliation, relative abundance of non-host species, phenological asynchrony, pesticide spray status, and relative slope position. Spatial autocorrelation in residual defoliation ranged over distances of 788 m in 2000 and 461 m in 2001, corresponding well with gypsy moth larval dispersal distances (100 m to 1 km). Un-measured processes such as predation, virus and pathogen occurrence likely contribute to unexplained variance. Because the spatial dynamics of these factors are largely unknown, our results support modeling gypsy moth defoliation as a function of dependence on significant exogenous characteristics and residual spatial pattern matching.     

Simulating Spatial Nitrogen Dynamics in a Forested Reference Watershed, Hubbard Brook Watershed 6, New Hampshire, USA

We demonstrate that available information on spatial heterogeneity in biotic, topographic, and climatic variables within a forested watershed, Hubbard Brook Experimental Forest (HBEF) Watershed 6, New Hampshire, USA, was sufficient to reproduce the observed elevational pattern in stream NO₃ concentration during the 1982–1992 period. Five gridded maps (N mineralization factor, N uptake factor, precipitation, elevation, and soil depth factor) were created from spatial datasets and successively added to the spatially explicit model SINIC-S as spatially varying input parameters. Adding more spatial information generally improved model predictions, with the exception of the soil depth factor. Ninety percent of the variation in the observed stream NO₃ concentration was explained by the combination of the spatial variation of the N mineralization and N uptake factors. Simulated streamflow NO₃ flux at the outlet point was improved slightly by introducing spatial variability in the model parameters. The model exhibited substantial cell-to-cell variation in soil N dynamics and NO₃ loss within the watershed during the simulation period. The simulation results suggest that the spatial distributions of forest floor organic matter and standing biomass are most responsible for creating the elevational pattern in stream NO₃ concentration within this watershed.     

Utilizing GIS Technologies in Selection of Suitable Vineyard Sites

Abstract

Geographic information system (GIS) technologies and a weighted linear indexing model were used for suitability analysis of potential vineyard sites in Illinois. The model included a macroscale climate variables layer (40 points), a mesoscale climate variables layer (40 points), a soil properties layer (10 points) and a current land use variables layer (10 points) for a possible 100 points. Macroscale climate variables, growing degree day summation for a 33-year period (1969-2002) and occurrence of -26°C were interpolated using thin plate smoothing splines over the Illinois terrain using 100 m2 resolution digital elevation models (DEM). Using the same DEMs, absolute elevation, slope, and aspect were reclassified using surface analysis of the terrain to model the effects of mesoscale climate variables in Jackson and Union Counties in Illinois (study area). Locations in the study area above 259 m in absolute elevation (above sea level), with gently rolling slopes (5 to 10%), and facing North, East, or Northeast received the most points within the mesoscale climate layer. Soils that drained well or moderately well with moderate organic matter content (2 to 3%) received the most points in the soil properties layer. Lands that promoted ease of conversion to vineyards were assigned the most points in the current land use layer. The resulting four layers were arithmetically summed and suitability maps with 760 m2 resolution were developed. In the study area, over 18,155 hectares of land highly suitable or suitable for viticulture were identified. A portion of the existing orchard and vineyard acreage (223 ha) in the study area were surveyed with a global positioning system. Of the surveyed vineyard acreage, over 81% lies within the suitable or better ranking according to the model. Of the surveyed orchard acreage, over 50% lies within the suitable or better ranking.     

Scale-dependent determinants of heterogeneity in fire frequency in a coniferous boreal forest of eastern Canada

Despite the recognized importance of fire in North American boreal forests, the relative importance of stochastic and determinist portions of intra-regional spatial variability in fire frequency is still poorly understood. The first objective of this study is to identify sources of spatial variability in fire frequency in a landscape of eastern Quebec’s coniferous boreal forest. Broad-scale environmental factors considered included latitude, longitude, human activities and belonging to a given bioclimatic domain, whereas fine-scale factors included slope, position on the slope, aspect, elevation, surficial deposit and drainage. The average distance to waterbodies was also considered as a potential intermediate-scale source of variability in fire frequency. In order to assess these environmental factors’ potential influence, they were incorporated into a proportional hazard model, a semi-parametric form of survival analysis. We also used a digital elevation model in order to evaluate the dominant aspect within neighborhoods of varying sizes and successively incorporated these covariates into the proportional hazard model. We found that longitude significantly affects fire frequency, suggesting a maritime influence on fire frequency in this coastal landscape. We also found that position on the slope was related to fire frequency since hilltops and upperslopes were subject to a lower fire frequency. Dominant aspect was also related to fire frequency, but only when characterized within a neighborhood delimited by 4,000 to 10,000-m radii (5,027–31,416 ha). A 2–6-fold variation in fire frequency can be induced by geographic and topographic contexts, suggesting a substantial intra-regional heterogeneity in disturbance regime with potential consequences on forest dynamics and biodiversity patterns. Implications for forest management are also briefly discussed. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Postcards from the past: charting the landscape-scale conversion of tropical Australian savanna to closed forest during the 20th century

Repeated sequences of digitised and geo-referenced historical aerial photography provide a powerful means of understanding landscape change. We use this method to demonstrate a landscape wide expansion of closed forest (42% increase in total coverage) in the Australian monsoon tropics over the past five decades. Retrospective habitat suitability models (HSI) of closed forest derived using four landscape measures (drainage distance, slope angle, aspect and elevation) for imagery taken in 1947 correctly forecast the subsequent spatial distribution of the expansion, with topographic fire protection primarily determining the closed-forest distribution. The dynamics of the closed forest-savanna boundary were predicted accurately by generalised linear models, with closed-forest expansion in fire-protected sites along forest edges and regression in the more fire-prone areas. Two factors may plausibly explain the expansion of closed forests. First, eco-ethnographic records stress the skilful use of fire by Aboriginal people in protecting isolated and locally resource-rich closed-forest patches. Second, the recent global increase in atmospheric CO₂ may be changing the competitive balance between savanna and forest by enabling C₃ trees to grow fast enough to escape the fire trap presented by flammable C₄ grasses.     

Quantifying the agricultural landscape and assessing spatio-temporal patterns of precipitation and groundwater use

Quantitative agricultural landscape indices are useful to describe functional relationships among climatic conditions, groundwater dynamics, soil properties and agricultural land use for mathematical models. We applied methods of regression statistics, variance component estimation and a Geographical Information System (GIS) to construct indices describing crops and soils and to establish functional relationships among these variables. This paper describes the development of indices and the partitioning of the spatial and temporal variation in groundwater models using the data from Tulare County, California, which was selected as the study area. Indices of ground surface elevation, total crop water demand, soil water infiltration rate, and soil production index explain 91% of the variation in average spring groundwater level. After relating spatial patterns of groundwater use to indices of crop and soil properties, we found that mean groundwater use is positively related to total crop water demand and soil water infiltration rate while the variation in groundwater use was negatively correlated with the crop water demand and soil water infiltration rate and positively related to soil water holding capacity. The spatial variation in groundwater use was largely influenced by crops and soil types while the temporal variation was not. We also found that groundwater use increased exponentially with decreasing annual precipitation for most townships. Based on these associations, groundwater use in each township can be forecast from relative precipitation under current methods of agricultural production. Although groundwater table depth is strongly affected by topography, the statistically significant indices observed in the model clearly show that agricultural land use influences groundwater table depth. These simple relationships can be used by agronomists to make water management decisions and to design alternative cropping systems to sustain agricultural production during periods of surface water shortages.     

The influence of sampling scheme and interpolation method on the power to detect spatial effects of forest birds in Ontario (Canada)

Spatial ecology is becoming an increasingly important component of resource management, and the general monitoring of how human activities affect the distribution and abundance of wildlife. Yet most work on the reliability of sampling strategies is based on a non-spatial analysis of variance paradigm, and little work has been done assessing the power of alternative spatial methods for creating reliable maps of animal abundance. Such a map forms a critical response variable for multiple scale studies relating landscape structure to biotic function. The power to reconstruct patterns of distribution and abundance is influenced by sample placement strategy and density, the nature of spatial auto-correlation among points, and by the technique used to extrapolate points into an animal abundance map. Faced with uncertainty concerning the influence of these factors, we chose to first synthesize a model reference system of known properties and then evaluate the relative performance of alternative sampling and mapping procedures using it. We used published habitat associations of tree nesting boreal neo-tropical birds, a classified habitat map from the Manitou Lakes area of northwestern Ontario, and point count means and variances determined from field studies in boreal Canada to create 4 simulated models of avian abundance to function as reference maps. Four point sampling strategies were evaluated by 4 spatial mapping methods. We found mixed-cluster sampling to be an effective point sampling strategy, particularly when high habitat fragmentation was avoided by restricting samples to habitat patches >10 ha in size. We also found that of the 4 mapping methods, only stratified ordinary point kriging (OPK) was able to generate maps that reproduced an embedded landscape-scale spatial effect that reduced nesting bird abundance in areas of higher forest age-class fragmentation. Global OPK was effective only for detecting broader, regional-scale differences. This revised version was published online in July 2006 with corrections to the Cover Date.     

Sod establishment and turfgrass growth as affected by urea–formaldehyde resin foam soil amendment

Urea–formaldehyde resin foam (UFRF) has recently been introduced as a soil amendment for turfgrass culture. A field study was performed to evaluate the impact of UFRF on soil physical and chemical properties, sod establishment and turfgrass shoot and root growth. Treatments included a non-amended, sandy loam soil and the same soil amended with UFRF at a rate of 20% (v/v) incorporated into the upper 0.1 m of the substrate. Physical and chemical analyses of the two substrates involved the determination of bulk density, total porosity, air-filled porosity at 0.4 m, moisture release curve, pH and electrical conductivity. Turfgrass growth was determined through several measurements that included shoot growth rate, root growth rate and monitoring of three root architecture parameters namely, total root length, total root area and mean root diameter. Turfgrass establishment was evaluated by measuring the vertical force applied to detach the sod from the substrate.UFRF amendment did not influence pH and EC and provided minimal alterations (non-significant) in soil physical properties by slightly increasing total porosity, easily available water and air-filled porosity, and by reducing bulk density. The impact of UFRF amendment improved cumulative clipping yield in five sampling dates but root growth and sod establishment rate were better than the control only on one sampling date. There appears to be limited benefit in amending soil with UFRF.     

Small-scale effects of historical land use and topography on post-cultural tree species composition in an Alpine valley in southern Switzerland

Investigations of spatial patterns in forest tree species composition are essential in the understanding of landscape dynamics, especially in areas of land-use change. The specific environmental factors controlling the present patterns, however, vary with the scale of observation. In this study we estimated abundance of adult trees and tree regeneration in a Southern Alpine valley in Ticino, Switzerland. We hypothesized that, at the present scale, spatial pattern of post-cultural tree species does not primarily depend on topographic features but responds instead to small-scale variation in historical land use. We used multivariate regression trees to relate species abundances to environmental variables. Species matrices were comprised of single tree species abundance as well as species groups. Groups were formed according to common ecological species requirements with respect to shade tolerance, soil moisture and soil nutrients. Though species variance could only be partially explained, a clear ranking in the relative importance of environmental variables emerged. Tree basal area of formerly cultivated Castanea sativa (Mill.) was the most important factor accounting for up to 50% of species’ variation. Influence of topographic attributes was minor, restricted to profile curvature, and partly contradictory in response. Our results suggest the importance of biotic factors and soil properties for small-scale variation in tree species composition and need for further investigations in the study area on the ecological requirements of tree species in the early growing stage.     

Total soil carbon assessment: linking field,lab, and landscape through VNIR modelling

Context

Point based measurements provide only a limited overview of landscape variation in measured properties. Upscaling of measurements from point to landscape comes with challenges particularly considering error propagation.

Objectives

We investigated the impact of using proximal derived measurements of soil total carbon taken at point locations on upscaling to landscape levels.

Methods

1087 soil samples across Florida, USA were collected, laboratory (LAB) analysed for total carbon (TC), and then measured using visible-/near-infrared (VNIR) spectroscopy. Proximal TC values were generated through chemometric modelling using random forest (RF) and partial least squares (PLS) regression. These three datasets were then upscaled to the State of Florida, USA using ordinary kriging and compared.

Results

R² (RPD) values for the PLS and RF chemometric models were 88% (2.96) and 91% (3.23), respectively. All 3 spatial models had an accuracy of 54% on an independent validation dataset, with greater than 70% accuracy if predicted values were considered within the interpolation variance range. When comparing spatial interpolations derived from the proximally measured samples, only 18% of the PLS versus 51% of the RF fell within a range of 0.05 logTC (g kg^?1) of the LAB measured interpolations.

Conclusions

Using proximal sampling and modelling provides comparable output to laboratory measured soil TC measurements at point level, but when upscaled to landscape level the selection of proximal modelling method will impact the spatial interpolations derived. The error propagation within sequential modelling must be considered particularly when one wishes to use sequential modelling to analyse change in environmental properties.

     

Forest gradient response in Sierran landscapes: the physical template

Vegetation pattern on landscapes is the manifestation of physical gradients, biotic response to these gradients, and disturbances. Here we focus on the physical template as it governs the distribution of mixed-conifer forests in California's Sierra Nevada. We extended a forest simulation model to examine montane environmental gradients, emphasizing factors affecting the water balance in these summer-dry landscapes. The model simulates the soil moisture regime in terms of the interaction of water supply and demand: supply depends on precipitation and water storage, while evapotranspirational demand varies with solar radiation and temperature. The forest cover itself can affect the water balance via canopy interception and evapotranspiration. We simulated Sierran forests as slope facets, defined as gridded stands of homogeneous topographic exposure, and verified simulated gradient response against sample quadrats distributed across Sequoia National Park. We then performed a modified sensitivity analysis of abiotic factors governing the physical gradient. Importantly, the model's sensitivity to temperature, precipitation, and soil depth varies considerably over the physical template, particularly relative to elevation. The physical drivers of the water balance have characteristic spatial scales that differ by orders of magnitude. Across large spatial extents, temperature and precipitation as defined by elevation primarily govern the location of the mixed conifer zone. If the analysis is constrained to elevations within the mixed-conifer zone, local topography comes into play as it influences drainage. Soil depth varies considerably at all measured scales, and is especially dominant at fine (within-stand) scales. Physical site variables can influence soil moisture deficit either by affecting water supply or water demand; these effects have qualitatively different implications for forest response. These results have clear implications about purely inferential approaches to gradient analysis, and bear strongly on our ability to use correlative approaches in assessing the potential responses of montane forests to anthropogenic climatic change.     

Spatial analysis of selected soil attributes across an alpine topographic/snow gradient

The impact of the topographic/snow gradient on soil processes in alpinetundra on Niwot Ridge of the Colorado Front Range (Rocky mts, USA) was assessedusinggeostatistical modeling and a fractal approach. The mean snow depth, whichmeasured between 1984 and 2000, exhibited a smooth spatial continuity acrossthestudy grid area (550 × 400 meter). Soil color variables showed a nestedstructure that was attributed to a confounded effect of various soil-formingfactors on catenary processes. The spatial structure of texture classesexhibited no spatial structure and was explained by data sparsity,cryoturbation, and biological processes that mask the expected long-distancevariations (i.e., 550-m) of the catenary processes. Organic C, pH, bulkdensity,and soil moisture content showed various degrees of spatial continuity, but allindicated that the topographic/snow gradient is not the only dominantsoil-forming factor in this alpine ecosystem. The estimated fractal dimensionDfor the grid landform and the mean snow depth varied between 1.2 and 1.4,indicating that they vary smoothly with long-range variation. The estimatedDofthe soil variables ranged between 1.6 and 1.8, showing a noisy appearance withshort-range variations. These results strongly suggest that most small andmicro-scale variations in the alpine soil environs resulted from the combinedeffect of cryoturbation, biological activity, parent-material and eoliandeposition, whereas the large-scale variations originated as a result of thetopographic/snow gradient.This revised version was published online in May 2005 with corrections to the Cover Date.