wpeC.jpg (1715 bytes)  wpe3F.jpg (9195 bytes)        

Table Of Contents


2. Synthesis Activities

A. Field Experiments

Consequences of patterns in food webs on stability:
Data collected from the CPER has been used in conjunction with data from other sites as part of an across site comparison in search of patterns in food webs and the effects that disturbance has on these patterns and system stability. The work has resulted in 6 publications. The more significant results of the project have been reached by coupling the models traditionally used by community ecologists with those used by ecosystem ecologists. The studies listed below have provided insight into the roles of productivity and dynamics in structuring ecosystems, the calculation of interactions strength using field data, and common features of agricultural practices on the structure and function of soil communities.

We previously demonstrated that the resilience of simple food chains was directly related to the level of primary productivity and the rates of detritus inputs from external sources. The key to this study was in the development of Lotka­ Volterra based models that were expressed in terms of ecological energetics. This result would seem to support the "dynamics hypothesis" that proposes that food chain length was a function of the limitations that increased length placed on the likelihood of the system recovering from a minor disturbance. A second important result was the finding that the feasibility (ability of the system to maintain positive population densities at steady state) of the food chains was a function of the level of productivity and detritus inputs. Higher levels of productivity were able to sustain longer food chains. This finding had been hypothesized G.E. Hutchinson (1959 Am. Nat. 93:145)and subsequently termed the "energetics hypothesis." The conclusion of our work was that system dynamics and energetics were inextricably inter­related. Food chain length was a function of both energetics and dynamics.

Considerable progress has been made in the development of procedures used to estimate the interaction strength among organisms with a food web. Paine (1992 Nature 355:73) proposed measuring interaction strengths among species through a series of species removal experiments. As an alternative, we estimated interaction strength with models. The population dynamics of species within the community were described after Moore et al. (1993) and DeRuiter et al. (1994), and the field average population densities, were substituted for the theoretical equilibria. The elements of the community matrix (interaction strengths) are estimated by sampling physiological parameters, consumption coefficients and population densities, from ranges established by microcosm and field data.

An interesting pattern emerged when we estimated the interaction strengths (as above) and the trophic position of each consumer­prey interaction within food webs in North America and Northern Europe. The negative effects of consumers on prey decreases with increased trophic position and the converse was true for the positive effects of prey on consumer. De Ruiter et al. (Science, in review)demonstrated that communities that possess this asymmetric patterning of interaction strength with increased trophic position are more likely to be stable than matrices depicting the same trophic structure but no pattern in the distribution of interaction strength.

The conclusion of our work with food chain length and interaction strength is that the energetic organization of communities forms the basis of ecosystem stability. In 1995­1996 the work will be directed towards empirical tests of these findings. (see Figure webs)

References: Moore et. al 1993, Moore et. al 1993, Moore and DeRuiter 1993, Moore 1994, Moore et. al 1995 (accepted), DeRuiter, et. al 1995 (in review)

 

A quantitative assessment of the effects of grazing on vegetation and soils over a global range of environments:
Multiple regression analyses were performed on a world­wide data set compiled from studies that compared species composition, aboveground primary production (ANPP), root biomass, and soil nutrients of grazed versus protected, ungrazed sites.

Changes in species composition with grazing were primarily a function of ANPP and evolutionary history of grazing of the site, with sensitivity to level of consumption third in importance. Changes in species composition increased with increasing productivity and with longer, more intense evolutionary histories of grazing. These three variables together with latitude explained 50% of the variance in composition of grasslands, and precipitation, evolutionary history, and years of protection explained 72% of the variance for shrublands, even though methods of measurement and grazing systems varied between studies. The response of shrublands was different from that of grasslands, with a change in evolutionary history of grazing of shrublands resulting in less effect on change in species composition than in grasslands.

Similar variables entered models of change in the dominant species with grazing, but explanatory power was less than for species composition. Increases of the dominant species were more likely to occur in shrublands than grasslands, and decreases more likely among bunch grasses than other life­forms and among perennials than annuals.

Percentage differences in ANPP between grazed and ungrazed sites decreased with increasingly long evolutionary histories of grazing and increased with increasing ANPP. Although most effects of grazing on ANPP were negative, the statistical models predicted increases in ANPP with grazing under conditions of long history, low consumption, low ANPP, and large range in annual temperature for grasslands plus shrublands. The data and the models support the controversial hypothesis that grazing can increase ANPP in some situations.

Counter to the common held view that grazing negatively impacts root systems, there was no relationship between difference in ANPP with grazing and difference in root mass; as many positive as negative differences occurred even though most ANPP differences were negative. Further, there was a weak relationship between change in species composition and change in ANPP, and no relationship with root mass, soil organic matter, or soil nitrogen. All three belowground variables displayed both positive and negative values in response to grazing. Current management of much of the worlds grazinglands based on species composition criteria may lead to erroneous conclusions concerning the long­term ability of a system to sustain productivity. (see Figure 2.20b)

References: Milchunas and Lauenroth 1993

B. Simulation Modeling

STEPPE model:
We are continuing to use models of soil water dynamics (SOILWAT), plant community structure (STEPPE), and production and turnover of soil organic matter (CENTURY) to better understand various aspects of shortgrass steppe structure and function. We are continuing to use these models both independently and as linked, interacting models to aid in the formulation of questions and hypotheses, and in the interpretation of field and laboratory data.

References: Parton et al. 1988, Coffin and Lauenroth 1990

Species recruitment across an environmental gradient:
We evaluated the effects of precipitation and temperature on the recruitment of two dominant perennial C4 grasses along an environmental transect from the Central Plains Experimental Range in northern Colorado to southern New Mexico. Bouteloua gracilis ranges from Canada to Mexico, and dominates the shortgrass steppe. B. eriopoda is limited to the southern United States and Mexico, and dominates the desert grasslands. Successful sexual reproduction is rare in these species, but may be important in explaining the species geographic ranges. We evaluated the probabilities of recruitment for these two species at 16 sites in Colorado and New Mexico, including three LTER sites, the Central Plains Experimental Range, the Sevilleta, and the Jornada, using a daily time step, multi­ layer soil water model (SOILWAT). We used 30 years of climatic data from weather stations to produce 5000 years of weather data for each site using a Markov weather generator. Both the environmental factors of most importance to the two species and probability of recruitment along the transect differed. Precipitation was more important to the recruitment of B. gracilis while temperature was more important to the recruitment of B. eriopoda. Probability of recruitment decreased from north to south for B.gracilis while increasing from north to south for B. eriopoda. These results suggest that recruitment is important in the distribution and dominance of these two species.

References: Minnick, T.J. and D.P. Coffin 1995

The importance of soil water in the recruitment of Bouteloua gracilis in the shortgrass steppe:
In the shortgrass steppe region of North America there is a controversy about the ability of the dominant species to recruit from seedlings. The prevailing view is that B. gracilis is incapable of recruitment form seedling in areas receiving <380 mm of annual precipitation. A common explanation for this situation is that environmental conditions permitting seedling establishment are infrequent. To assess the frequency of environmental conditions appropriate for the recruitment of B. gracilis we used a soil water simulation model and long­term climatic data in conjunction with detailed information about the ecophysiological requirements for seed germination and growth of seminal and adventitious roots.

We found that recruitment events occur as frequently as every 30­ 50 yr. on silty clay, silty clay loam, and silty loam soils, but less that once in 5000 yr. on sandy soils. Simulated frequencies of recruitment were sufficient to account for the observed abundance of B. gracilis in 7 or 11 soil textures evaluated. The differences in silt content and available water occurrence of recruitment events. Therefore, soil texture variability may explain the spatial pattern of recruitment and of population recovery after disturbance that occur at the soil type and microsite scales.

Annual precipitation explained a large fraction of the temporal variability in recruitment. On average, recruitment occurred in years when precipitation was above the mean. The occurrence of recruitment events in some dry years (precipitation < mean), and their absence during some wet years (precipitation > mean), emphasizes the importance of the intraseasonal distribution of precipitation.

The sensitivity of recruitment to soil water availability suggests that climate change, particularly changes that increase or decrease the amount or the effectiveness of soil water, could have important effects on the future of populations of B.gracilis. (see Figure pop_b)

Reference: Lauenroth et al. 1994

Seedling recruitment potential of Bouteloua gracilis in the central grassland region
We evaluated the recruitment potential by seedlings of the perennial C4 grass,Bouteloua gracilis (blue grama) for the central grassland region of the U.S., including the Central Plains Experimental Range, under current climatic conditions and for changes in climate. Four grassland types are found within this large region. B. gracilis dominates the vegetation for one grassland type found at the CPER, the shortgrass steppe, and is common in the other three types (northern and southern mixedgrass, and tallgrass prairies). The geographic distribution of abundance of this species may be explained by its potential for recruitment by seedlings. We used a multi­layer daily time step soil water model to evaluate the probability of recruitment of B. gracilis seedlings for a range of soil textures and a range of current and expected changes in climatic conditions representative of the region. Simulations were conducted using daily precipitation and temperature data for > 200 weather stations. Probability of recruitment increased with increasing temperature and precipitation, and was also positively related to silt content of the soil. Probabilities were lowest in the coolest and driest areas, including parts of the shortgrass steppe. The effects of climate change on recruitment was dependent upon the balance between increases in temperature and decreases in precipitation. Sites most susceptible to change in climate were the shortgrass steppe and northern mixedgrass prairie where B. gracilis is most important to plant community structure and ecosystem processes.

Simulations conducted with the STEPPE model also indicated the importance of soil texture to the recovery of B. gracilis after disturbances of various sizes.  Average probabilities of seedling establishment from the SOILWAT simulations were used in STEPPE based on silt and content. Disturbances on the soil with a large silt content (silt loam) had large B. gracilis biomass through time compared to soils with low silt content (sandy clay loam, loamy sand). Additionally, recovery was dependent upon disturbance size, where small disturbances recover faster than large disturbances. A nonlinear relationship was found between the proportion of plots recovered by B. gracilis in year 500 and silt content of the soil. Disturbance size was less important than soil texture in determining plant recovery. (see Figures 1.15 2.33)

References:

CENTURY model:
We used the Century model to simulate the impact of climatic change and increased atmospheric CO2 on dynamics of the shortgrass steppe. The results showed that increasing air temperature resulted in increased plant production and reduced soil organic matter levels. Increasing atmospheric CO2 levels caused plant production and soil carbon levels to increase. The combined impact of increasing atmospheric CO2 and climatic change is to reduce the decrease in SOM level caused by increasing air temperatures and to further increase the plant production. As part of the ecosystem modeling work, we compared the Century model (monthly time step simplified model) with the Grass model (daily time step process oriented model) and found that the Century and Grass models both did a good job of simulating annual plant production and seasonal patterns of live biomass at CPER. Results from the modeling work are described in five papers (see attached references).

During the last year we have worked on improving the soil water and soil temperature model (SOWAT) and have developed a new model for simulating snow melt, sublimation and bare soil evaporation during the winter. The SOWAT model was tested using data from the PILPS soil water comparison workshop in Australia and compared to 12 other water flow models. A paper describing this work was written and submitted for publication. A paper describing the new winter water budget model has been written, however, it has not yet been submitted for publication. The winter water budget model has been developed using data from the CPER lysimeter and microclimate station. (see Figure 2.34)

References:

Soil­atmosphere exchange of CH4, N2O and CO2: using field­flux measurements for development of process­based gas­flux models for CH4 and N2O:
Using the information collected during the past four years from a wide variety of research sites, wedeveloped a set of first generation, process based models which describe the oxidation of CH4 and production of N2O in the soil. A general model was developed to simulate N2 and N2O fluxes from nitrification and denitrification. N2O flux from nitrification is a function of the soil water filled pore space (WFPS), soil temperature and the maximum nitrification rate (soil specific parameter). The effect of WFPS on nitrification is a function of soil texture with the optimum value for nitrification occurring at higher WFPS for fine textured soils. Total N2 and N2O gasfluxes from denitrification are a function of the soil NO3 level, soil respiration rateat 0.6 WFPS (index of available soil C) and soil WFPS. The maximum N gas flux from denitrification is calculated as the minimum of the soil respiration and soil NO3 functions and reduced by WFPS function. Denitrification rates decrease as WFPS decreases below 0.9, with the rate decreasing most rapidly in the fine textured soils. The ratio of N2:N2O gas fluxes from denitrification is a function of soil respiration rate, soil NO3 level and soil WFPS. The highest N2:N2O ratios occur for soils with lowsoil NO3 levels, and high values of soil respiration and WFPS. The model is used to relate field­measured fluxes to the processes which drive the fluxes.

The CH4 uptake model is based on the observation that CH4 uptake rates are limited by diffusion, and can thus be modeled based on Fick's first law of diffusion in which flux is proportional to the concentration gradient with respect to soil depth. Diffusion coefficients are estimated as functions of pressure, temperature, soil bulk density, texture and moisture. A factor relating soil N turnover to CH4 uptake is also included.

Further refinement of these models and linkage to the CENTURY model is planned. Once that phase is accomplished, then the model package will be linked with GIS systems to develop regional flux models.

References:

C. Regional Analysis

Our regional analysis efforts are in conjunction with another NSF­supported project (BSR 9106183) that provides the funding for these efforts. Below we summarize a few of our important findings from that project that are relevant to our site­level research at the CPER.

Above­ground net primary productivity across a precipitation gradient in the Central Grassland region:
Controls on above­ground net primary productivity (ANPP) and vegetation structure at the CPER and at 12 satellite sites were studied across a precipitation gradient through the Central Grassland region. Patterns of productivity across the gradient may be explained by the inverse texture hypothesis which predicts that in arid and semi­arid regions, plant communities on coarse­textured soils should have higher above­ground net primary productivity (ANPP) than those on fine­textured soils, while the reverse should occur in humid regions. The first objective was to test these predictions by measuring ANPP in 1993 and 1994 at 14 natural grassland sites (CPER + 12 satellite sites), equally divided between coarse and fine­textured soils, along a regional precipitation gradient (333 mm/y to 759 mm/y) where mean annual temperature remains constant. Below 400 mm of mean annual precipitation (MAP), coarse­textured sites were significantly more productive than fine­textured sites in 1993 but not in 1994. For all sites, ANPP was significantly correlated with both MAP and precipitation for the 12 month period preceding sampling. Soil texture explained a significant proportion of variability in ANPP of graminoids and forbs but not total ANPP.

The second objective was to quantify aboveground canopy structure of natural grasslands along the precipitation gradient to examine the consequences of increases in the importance of aboveground competition with increased productivity. We found abrupt shifts in maximum canopy height and in PAR transmission at the transitions between vegetation types. Transmission of PAR through the canopy decreased with increasing precipitation. With increasing precipitation, total biomass and LAI increased, but the proportion of biomass and leaf area near the soil surface decreased. Overall, these results revealed discontinuous changes in canopy structure between different vegetation types.

The third objective was to examine long­term precipitation dynamics and long­ term soil water dynamics for two contrasting soil textures to understand potential mechanisms for observed patterns in ANPP and canopy structure across the precipitation gradient. We used a daily time step multi­layer soil water simulation model to examine these dynamics. With increasing total precipitation, the proportion of growing season (April 1 ­ October 1) precipitation accounted for by large events (> 30 mm) increased significantly. Distribution of precipitation interacted with soil texture in determining transpiration and evaporation. Transpiration from the sandy soil was higher than from the clay loam at the four driest sites, while the reverse was true at the two wettest sites. For the clay loam, the deep 30­100 cm soil layer represented a consistently important water resource only for the two wettest sites. Simulation results suggest that across this precipitation gradient the relative success of plant functional types with different rooting distributions and phenologies may vary according to soil texture.

Results presented here demonstrated that effects of abiotic driving variables on ecosystem dynamics may be mediated through vegetation structure. These results indicate the importance of including vegetation structure in regional analyses of ecosystem structure and function.

References:

A regional analysis of N mineralization in coarse and fine textured soils along a precipitation gradient in the Central Grasslands of the U.S:
We estimated field and laboratory N mineralization rates for seven paired sites of fine and course textured soils along a precipitation gradient extending from the semiarid shortgrass steppe of Eastern Colorado through the mesic tallgrass prairie of Nebraska.

No clear trend in field or laboratory rates of N mineralization was evident along this precipitation gradient. Laboratory rates were generally higher than field rates, especially toward the dry end of the gradient. Fine textured soils generally supported higher rates of N mineralization than did coarse textured soils.

References: Barrett et al. 1995

N retention in the northern mixed­grass prairie and shortgrass steppe:
Anthropogenic sources of fixed N contribute to a suite of environmental problems including acid deposition, groundwater pollution and accumulation of greenhouse and photochemically active grasses. Additions of fixed forms of N in excess of plant demand, may stress native terrestrial systems. We are developing an experiment to study the effects of elevated levels of N on native grassland systems along a temperature gradient of the Central Grasslands. We have collected soils form five sites along this gradient, extending form SE Montana, through Wyoming and Colorado, to the Panhandle of Texas. Preliminary analysis will include estimates of gross and net rates of N mineralization, potential C mineralization and SOM. Plans to initiate a fertilization experiment and monitor N trace gas flux are being developed.

References:

Regional patterns of NDVI in North American shrublands and grasslands:
We characterized the seasonal pattern of the Normalized Difference Vegetation Index (NDVI) of 49 grassland and shrubland sites corresponding to areas with low human impact. Each site was characterized as a vector of 21 elements, where each element corresponded to the average NDVI for one date. We used the integral of the NDVI as an estimate of the total production, and the difference between the maximum and minimum NDVI over the year as a measure of the seasonality.

Based on a Principal Component Analysis, the integral of the NDVI accounted for 63% of the variability of the SITEXDATE matrix. The difference between the maximum and minimum NDVI throughout the year accounted for 19%. Mean annual precipitation explained 67% of the variability of the integral of the NDVI curves (p<0.01). The slope of the relationship between precipitation and integrated NDVI was significantly steeper for grasslands than for shrublands. The difference between the extreme values of NDVI increased with mean annual precipitation and the thermal amplitude, and decreased with mean annual temperature (p<0.01). (see Figure 2-20a)

References: Paruelo and Lauenroth 1995

Effects of grazing on ecosystems of the Great Plains:
Grazing by large herbivores has been part of the evolutionary and developmental environment of the Great Plains for the past 10,000 years. There is evidence that past grazing has had an important positive effect on the ability of the ecosystems of the Great Plains to withstand grazing by domestic livestock. Settlement of the region by European agriculturalists in the past 150 years has had major impacts, but there is no conclusive evidence that the change from wild ungulates to properly managed domestic livestock has had detrimental ecological impacts.

The major west­east and south­north climatic gradients provide the best explanations for differences in the responses of ecosystems to grazing. Experiments conducted in the western, driest portion of the region suggests that plant communities in these areas are the least susceptible to detrimental impacts of improper grazing management. These areas also have the lowest productivity both of plants and animals. By contrast, plant communities in the eastern, wettest part of the plains has the capacity to support the largest number of animals per unit area but they are most susceptible to changes in community composition, diversity, and physiognomy under well managed light to moderate grazing. These communities are also potentially the most vulnerable to the detrimental impacts of poor management. Aboveground net primary production generally declines with grazing, although increases in the wettest areas have been attributed to opening dense canopies and reducing litter layers.

Interactions between livestock grazing and wildlife is a complex issue in the Great Plains as it is every where else. The same grazing management system can have positive impacts for one wildlife species and negative for others. Furthermore, it is clear in may cases that the effects of one management system on a single species can be a mixture of positive and negative impacts.

Despite the fact that grazing by domestic livestock does not have simple relationships to the other elements of the structure and function of Great Plains ecosystems, management of grazing animals using current technology appears to be a sustainable resource use. The conversion of rangeland to cropland has many implications for both sustainable resource use and global environmental issues. (see Figure 2.20d)

References: Lauenroth et al. 1994

ANPP estimates from NDVI for the Central Grassland region of the US:
The lack of a quantitative ecological interpretation of remotely sensed data is a major constraint for the use of satellite imagery in environmental studies. The ability to translate spectral data to biologically meaningful variables is a key step in increasing the use and value of satellite information. In this paper we present an empirical calibration of the relationship between the integral of the Normalized Difference Vegetation Index (NDVI) and aboveground net primary production (ANPP) for the Central Grassland Region of the United States. We calibrated a three year average of the NDVI integral (NDVI­I) (AVHRR/NOAA 1.1 km dataset) using spatially aggregated values of ANPP from rangeland survey data obtained from the Soil Conservation Service. This is the first attempt to calibrate a standard NDVI product for temperate grasslands.

For grassland areas with mean annual precipitation between 280 mm and 1150 mm, and mean annual temperature between 4oC and 20oC, the NDVI­I and ANPPshowed a strong and statistically significant relationship (r=0.93, n=19, p<0.001). NDVI­I was positively related to the ANPP. The relationship found is based on NDVI­I and ANPP data averaged over time and aggregated at the same spatial scale.

References: Paruelo et al. submitted

Climatic controls of the relative abundance of plant functional types in grasslands and shrublands of North America:
It has been widely accepted that temperature is the principal control of the distribution of the number of C3 and C4 species. Our results showed that at a regional scale and based on data of relative abundance, precipitation and its seasonal distribution are important controls, in addition to temperature, of the distribution of C3 and C4 grasses. Latitude and longitude explained a substantial portion of the variability of the distribution of shrubs, C3 grasses, and C4 grasses (55%, 46%, and62% respectively). Along a given longitude, C3 grasses increased with latitude. As one moves westward there is a replacement of C4 grasses by shrubs. The relative abundance of C4 grasses reached a maximum at southern latitudes and eastern longitudes. Succulents showed a marginal decreasing trend with latitude. For forbs, no relationship with geographic variables was detected. C4 grass distribution was positively related with three climatic variables: mean annual precipitation, mean annual temperature, and the proportion of the precipitation falling in summer. These variables accounted for 66% of the total variability of this functional type. C3 grass abundance decreased with mean annual temperature and increased with the proportion of the precipitation falling during winter (r2 = 0.37). Sixty two percent of the relative abundance of shrubs was explained by climate. The abundance of this functional type decreased with mean annual precipitation and increased with the proportion of water falling in winter.

References: Paruelo and Lauenroth 1995

Cross LTER­site comparison of factors that control decomposition: Effects of Bouteloua gracilis and Bouteloua eriopoda and indigenous microbes on rates of decomposition and soil organic matter formation:
We initiated a litter bag root decomposition experiment using transplanted soils from the Jornada and CPER long­term ecological research sites. Roots of the species Bouteloua gracilis and Bouteloua eriopoda were placed in both types of soils. The experiment began in late January 1994. We will harvest litter bags and soils in early July of this year, and calculate mass loss of the root litter.

Generating and testing daily weather data using a regional grassland data set:
Water balance will be a key aspect of ecosystem function in grasslands in response to global change. Most soil water balance models need daily weather data to generate reliable patterns. Those data are scarce for many regions of the world. We developed a procedure that generates daily data from a global data base of monthly precipitation and temperature on a 5o x 5o (lat., long.) scale. We tested the procedure by simulating daily data for the Central Grasslands USA and comparing them to actual data for 300 weather stations collected over a twenty year period. The generated data reproduced the regional and seasonal patterns of climate of the region. Multiple regression analysis showed that monthly averages of temperature, maximum temperature, minimum temperature, precipitation and number of rainy days were associated with geographic variables (lat., long.) for both generated and actual data (p<0.0001) in a similar manner.

References:

 

wpe1C.jpg (1101 bytes)

wpe26.jpg (1525 bytes)                               wpe25.jpg (1298 bytes)                                 wpe16.jpg (1211 bytes)

 

  02/08/01


About - Reports - Research - Data - Publications - Personnel -
News -
Links - What's New - Web Site Map - Home

To contact us, please email: Sallie Sprague  (Sallie.Sprague@colostate.edu)