Effects of simulated daily precipitation patterns on annual plant populations depend on life stage and climatic region
© Köchy; licensee BioMed Central Ltd. 2008
Received: 08 December 2006
Accepted: 27 March 2008
Published: 27 March 2008
To improve the understanding of consequences of climate change for annual plant communities, I used a detailed, grid-based model that simulates the effect of daily rainfall variability on individual plants in five climatic regions on a gradient from 100 to 800 mm mean annual precipitation (MAP). The model explicitly considers moisture storage in the soil. I manipulated daily rainfall variability by changing the daily mean rain (DMR, rain volume on rainy days averaged across years for each day of the year) by ± 20%. At the same time I adjusted intervals appropriately between rainy days for keeping the mean annual volume constant. In factorial combination with changing DMR I also changed MAP by ± 20%.
Increasing MAP generally increased water availability, establishment, and peak shoot biomass. Increasing DMR increased the time that water was continuously available to plants in the upper 15 to 30 cm of the soil (longest wet period, LWP). The effect of DMR diminished with increasing humidity of the climate. An interaction between water availability and density-dependent germination increased the establishment of seedlings in the arid region, but in the more humid regions the establishment of seedlings decreased with increasing DMR. As plants matured, competition among individuals and their productivity increased, but the size of these effects decreased with the humidity of the regions. Therefore, peak shoot biomass generally increased with increasing DMR but the effect size diminished from the semiarid to the mesic Mediterranean region. Increasing DMR reduced via LWP the annual variability of biomass in the semiarid and dry Mediterranean regions.
More rainstorms (greater DMR) increased the recharge of soil water reservoirs in more arid sites with consequences for germination, establishment, productivity, and population persistence. The order of magnitudes of DMR and MAP overlapped partially so that their combined effect is important for projections of climate change effects on annual vegetation.
Drylands are characterized not only by low water availability but also by great variability of water supply, which play important roles in structuring ecosystems and maintaining biodiversity [1–3]. However, our knowledge about the effects of climate change on vegetation with regard to changes in daily precipitation patterns is still limited . In this paper I use a simulation approach to systematically examine the effect of daily rainfall variability on the growth of annuals in the Middle East along a climatic gradient from arid to mesic Mediterranean.
Precipitation in dry climates has a high variability among and within years . The distribution of annual and daily rain amounts have already changed in the 20th century and are predicted to continue changing . In many Mediterranean regions the number of days in a year with heavy precipitation is increasing, whereas mean annual precipitation is decreasing . For the years 2071–2100, global climate models project a regionally varying shift of both mean annual precipitation and distribution of daily rainfall intensity. For subtropical/warm-temperate regions (25–40° latitude), which include the Middle East, mean annual precipitation is projected to increase by up to 100 mm this century . At the same time mean surface temperatures are projected to increase by c. 1.5°C in tropical and subtropical latitudes (0–35°, ). Compared to the change in water supply, the increase in temperature seems less important for plant growth in the arid and semi-arid climates within this range of latitude .
The effect of the temporal distribution of water availability on the growth or survival of herbaceous plants has been studied in field and pot experiments in temperate to tropical grasslands [9–13], deserts [14–16], and crops [17, 18]. In most experiments the treatments consisted of applying a fixed amount of water in either small amounts with short intervals or in large amounts with long intervals. Generally, the amounts and intervals were uniform during the duration of the experiment. Only the RaMPs project [9, 19] varied naturally occurring intervals, although by a fixed percentage. In contrast, climate models and observations suggest that changes will be stronger in the distribution of rainstorms than in the distribution of light rainfalls [7, 20]. Experiments that have examined water variability in mesic environments and most studies of crop irrigation methods found that more frequent watering increased plant growth or survival [9, 10, 14, 16, 17]. In contrast, studies [13, 16, 18, 21] and simulations  related to arid environments show that longer intervals between water pulses can have positive effects when rain pattern and the water holding capacity of the soil interact to produce a longer-lasting soil water reservoir [11, 18, 22]. The interaction between rainfall variability, soil texture, and water recharge has been described qualitatively by Noy-Meir as the 'inverse texture effect' , but the effect of different rain patterns on plant growth in arid and mesic regions has not been studied and quantified yet.
Process-based models represent a good way for studying the presumed causes of the effects of precipitation patterns in greater detail, resolving apparently contradictory effects, and projecting consequences of climate change . Models can be used to vary systematically and in combination characteristics of rain patterns, soil properties, and plant functional types, which would otherwise be hard to achieve logistically or technically in experiments. For example, simulations for North-American deserts have indicated that soil texture, the context of rain events regarding antecedent soil moisture, and the clustering of rain events are important factors controlling plant growth .
Standard climate parameters used in the model.
dry Mediterranean, DM
typical Mediterranean, TM
mesic Mediterranean, MM
Precipitation and soil moisture
Increasing the amplitude of a bell-shaped, seasonal distribution of DMR by 20%, while keeping the mean annual precipitation constant, decreased the average number of days with >5 mm of rain and increased the number of days with >15 mm and > 25 mm. Correspondingly, the median interval between days with >5 mm rain decreased, but it increased between days with >15 and >25 mm rain. Thus, the manipulation of rain pattern shifted the distribution of daily rain volumes to more frequent heavy rainfalls at the cost of light rains as observed in the Mediterranean .
The annual mean longest wet period (LWP) was assessed as an integrative measure of water availability. LWP was defined as the maximum number of consecutive days where the water potential (Ψ) of at least one of the four mineral ("A") soil layers is higher than the permanent wilting point (-3 MPa) representative of the herbaceous plants along the aridity gradient.
ANCOVAs of LWP, seedling density, and peek shoot mass.
Peak shoot mass
regions × MAP
regions × DMR
MAP × DMR
regions × MAP × DMR
The results presented here are based on LWP calculated for a permanent wilting point of -3 MPa. Wilting points down to -32 MPa, much lower than the conventional -1.5 MPa for cultivated plants , are characteristic for dryland annual species . A wilting point of -3 MPa seems representative for a range of annual grass species from dry climates . This includes Brachypodium distachyon, a common grass along the precipitation gradient in Israel. The discussion pertains nonetheless to species with low, intermediate, and high permanent wilting points as similar results were obtained for LWP calculated for different wilting points of -1.5 and -5 MPa.
Seedling density was analyzed as an integrative measure of plant responses to DMR via water availability during the early establishment phase. Results reported here, in the preceding and the following section are based on simulations where seed bank density was reset at the beginning of each simulated vegetation year to a value characteristic for each region. The resetting allowed the analysis of rain pattern effects without the influence of inter-annual carry-over effects caused by sequences of wet or dry years. Simulations without seed bank resetting are reported in the section 'Population persistence'.
Seedling density varied with DMR depending on interactions with regions and regions × MAP (Tab. 2). Change of MAP, and its interaction with DMR had no significant effects. The average seedling density of each region increased along the humidity gradient in a hyperbolic fashion (Fig. 2C). Within the arid and semi-arid regions average seedling density increased with MAP, but decreased with increasing MAP in the two most humid regions (Fig. 2D). Increasing DMR at constant nominal MAP increased seedling density in the arid region, whereas increasing DMR generally caused density to decrease in the other regions (Fig. 2D) and did not vary with MAP in the dry Mediterranean region. These contrasting findings suggest that density-dependent germination became more important as water supply increased. The three-way interaction resulted from no effect of MAP on the increase of seedling density with DMR in the arid region, a negative effect of raising MAP on the decrease of density with increasing DMR in the semiarid region, and an on average positive effect of raising MAP on the decrease of density with increasing DMR in the dry to mesic Mediterranean regions.
Peak shoot mass
The effects of rain manipulations on peak shoot mass were similar to that on LWP. Peak shoot mass depended on region, DMR, MAP, and all of their interactions (Tab. 2). The greatest absolute differences among levels were among regions. The average peak shoot mass of each region increased slightly sigmoidally along the humidity gradient. Within each region, mass increased with MAP (Fig. 2E). The 'absolute' rate of increase diminished from the semi-arid to the mesic Mediterranean region (Fig. 2F), whereas the 'relative' rate of increase was greatest in the dry-Mediterranean region (Additional file 4: RelativeChanges.pdf). Further, peak shoot mass (averaged across sites) increased with increasing DMR, but the rate of increase was smaller than that caused by a similar change in MAP. The rate of increase with DMR (averaged across levels of MAP) was greatest in the semiarid region (Fig. 2F). Within the arid region, the effect of increasing DMR was positive and increased with increasing MAP (Fig. 2F). In the other four regions, however, the positive effect of DMR decreased with increasing MAP so that it became insignificant, or sometimes negative, in more mesic regions. In summary, increasing daily mean rain volume and, more strongly, mean annual precipitation both increased peak shoot mass, but their relative effects generally decreased from the semiarid or dry Mediterranean to more arid or more mesic regions.
ANCOVAs of the coefficient of variation of rain, LWP, and peak shoot mass.
Peak shoot Mass
regions × MAP
regions × DMR
MAP × DMR
regions × MAP × DMR
Correlation between coefficients of variation of annual rain, LWP, and peak shoot mass.
CVLWP × CVrain
CVmass × CVrain
CVmass × CVLWP
CVmass.LWP × CVrain
I simulated the change of stochastic daily rain patterns by increasing the seasonal amplitude of daily mean rain volume (DMR) matched by a decrease in the seasonal occurrence of rainy days so that the mean annual volume remained unchanged. As a result of the increase of DMR, heavy rainfalls contributed more and light rainfalls less to annual precipitation than under current conditions.
The shift of rain pattern had clear effects on the average longest wet period (LWP; Fig. 2A, B). The effect of changing mean annual precipitation (MAP) was stronger than that of changing DMR and keeping MAP constant but the orders of magnitude were similar. The increase of LWP based on absolute change of MAP became flatter along the humidity gradient (Fig. 2B) because additional rain was less important for maintaining the soil moisture of an already moist soil in the more humid regions. In addition, the soil became saturated and additional water drained to lower soil layers or ran off at the surface. The increase of LWP with relative change of MAP was unimodal and greatest in the dry Mediterranean region (Additional file 4: RelativeChanges.pdf) because absolute changes of MAP were smaller in the more arid regions and the effect of increasing MAP diminished with the humidity of the region.
Changing DMR produced a complex response of LWP. In almost all instances, a greater amplitude of DMR, i.e., more days with rainstorms, prolonged the wet period but the effect size decreased with the humidity of the region. The effect of DMR occurred to a small degree because a higher likelihood of rainstorms advanced and delayed the average dates of the first and last rainstorms, correspondingly, prolonging the wet period. More importantly, rainstorms filled up the soil immediately and to a greater depth than did light rains, especially in the more arid regions with sandier soils. The depth of recharge considered in the model is within the soil's A layer (extending to 16–31 cm, depending on region), where annuals have most of their roots  and from where grassland plants take up most of their water . The water in deeper parts of the A layer evaporated more slowly than it would do near the surface, because the soil dried out from the surface downwards and dry soil was less permeable than wet soil. This response of LWP to a change in DMR was due to the texture-dependent, non-linear relation between soil water potential and soil water content as shown by the high ranks of retention-curve parameters in the sensitivity analysis (see Sensitivity analyses, below). It corroborates other simulations and experimental results [15, 22] that sandier soils have deeper soil water recharge from rainstorms than clayey soils, whose water availability is increased by more regular rainfalls (inverse-texture effect ). The latter effect did not appear in my simulations because the clayey soils were restricted to the more humid climates, where the mean interval between rainfalls > 2 mm (3 d) was sufficient to keep the soil's A layer above the LWP threshold of -3 MPa. The relative effect of soil texture on peak shoot mass is comparable to that of ± 20% change of DMR in simulations . Therefore, the parameters of the soil module should be adapted when simulations are carried out for other target sites. The relation between soil texture and water holding capacity also affects nutrient availability , which was not considered in the simulations, and plant growth [22, 25, 40], which will be discussed farther down.
In summary, the effects of changing MAP and DMR on soil water supply to plants were additive, depended on soil texture, and increased with the aridity of the region. Changes in DMR became important only when they were strong and the change in MAP was small.
Seedling density increased among regions along the humidity gradient (Fig. 2C, D) because the number of seeds was a fixed input parameter that increased from the arid to the mesic Mediterranean region and germination rates did not differ enough among sites to change this pattern. The increase of seedling density with MAP within regions was positive in the arid and semiarid regions but negative in the Mediterranean regions (Fig. 2D). This was the outcome of combining the concept of hydrothermal time for germination  with density-dependent germination fractions  in the model. Evidently, additional rain enhanced the number of days when conditions for germination were met and increased establishment in the arid and semiarid regions. In contrast, more rain in the Mediterranean regions did not greatly improve the conditions for germination. More rain, however, raised the germination fraction (Detailed model description – Seed bank module, eq. 2) so that competition among germinating seeds intensified and relative seedling density decreased with increasing MAP in the Mediterranean regions. Test simulations without density regulation showed that heavy rainfalls would cause mass germination in the mesic regions and later strong competition among adult plants. Therefore, density-dependent germination is an evolutionary stable strategy . Empirical evidence for this strategy in natural communities is accumulating [42–47].
Independent of changes to MAP, increasing DMR increased seedling density in the arid region, but decreased density in the other regions. In all regions, germination fractions increased with more regular rainfalls according to the applied hydrothermal time concept, but seedlings had access only to water in the upper 5 cm of the soil, which dried out quickly. Therefore, successful establishment of seedlings in the arid region was mostly controlled by access of seedling roots to water provided by rainstorms (high DMR). The increase of seedling density with less frequent watering was also observed in dry temperate grassland, but the effect was attributed to the breaking of seed dormancy . In the less arid regions more regular rainfalls increased the germination fraction and rainstorms were less important for providing water for seedlings. In summary, seedling density was determined mostly by seed availability, density-dependent germination, number of rainstorms under arid conditions, and regularity of rain under more mesic conditions.
Peak shoot mass
Generally, the effect of rain manipulations on peak shoot mass (Fig. 2E) followed a similar pattern as that on LWP (Fig. 2A), showing the strong control that water as a limiting resource exerted on growth. The change of peak shoot mass with relative change of MAP was greatest in the dry Mediterranean region. There, more seeds could germinate and more rain improved soil moisture conditions during establishment than in the more arid regions. On the other hand, individuals grew more sparsely so that competitive effects were weaker than in the typical and mesic Mediterranean regions. In general, the importance of rainstorms (high DMR) for growth decreased via LWP, whereas the importance of competition and the density of individuals increased along the humidity gradient. This change of ranking also showed in the sensitivity analysis (see Sensitivity analyses, below).
The decrease of seedling density with MAP in the non-arid regions did not cause a similar decrease in peak shoot mass. Instead, competition for water among established plants reduced their density. Since the effect of competition decreased with increasing LWP, the final density of mature individuals increased in the semiarid to typical Mediterranean regions but still decreased with LWP in the mesic Mediterranean region. In this region the decrease of density, however, was compensated by higher production per individual so that peak community shoot mass still increased with MAP.
To summarize, the sigmoid increase of community mass with MAP resulted from a similar increase of LWP with relative change of MAP along the humidity gradient, a larger effect of additional rain on establishment in arid than in non-arid regions, and increasing importance of competition for water. The increase of community mass with DMR was the sum of a positive effect of rainstorms on establishment in more arid regions, a negative effect on germination in non-arid regions, and a positive effect on LWP resulting in more intense competition but also higher production with increasing MAP.
Comparison with experiments
The result that the effect of daily rain distribution on the growth of annuals depends on the climatic humidity of the region was also found in experiments where the temporal distribution of water applications was varied independent of the total amount. In more mesic regions, plant growth either remains unaffected or increases with more evenly distributed water applications (lower DMR). For example, in a temperate grassland, water applications every second day vs. monthly more than doubled plant biomass . Similarly, in tall-grass prairie, doubling the intervals between rainfalls (higher DMR) reduced the standing crop by about 10% . Tropical savanna grasses, however, did not respond to watering variance, neither in the field nor in pot experiments (unclipped treatment) . In contrast, plant growth in arid regions tends to be enhanced by stronger pulses with longer intervals (higher DMR). Irrigation every two weeks instead of daily increased the yield of onions by about 40% in sub-Sahelian Burkina Faso . Greater water volumes applied with longer intervals increased the density of annuals and growth of perennials in an arid pasture . The significance of climate was directly shown experimentally . Increasing DMR by increasing intervals between water applications reduced community productivity in a mesic tall-grass prairie but increased it in the drier mixed- and short-grass prairie . The climate effect was also shown in a pot experiment using species of Mediterranean and desert origin . The two desert species grown with low total seasonal water supply had highest survival rates when pulse intervals after establishment were long (11 or 23 d). In contrast, the two Mediterranean species grown with higher total seasonal water supply survived better when time between pulses was 23 d rather than 3 d .
Obviously, intervals between rains cannot be extended indefinitely without having a negative effect on plant production. The intervals used in this model (3 d [mesic Mediterranean] - 7 d [semi-arid], 18 d [arid]) are similar to those found in nature , where longer intervals between rainstorms are interrupted by light rainfalls. In contrast, 'long interval' treatments in experiments are often regularly spaced and much longer than either natural or projected frequencies for a given location; a factor which may have contributed to the observation of negative correlations between productivity and interval length in mesic climates [9, 10, 14, 16].
In the second set of simulations, using a dynamic seed bank that included carry-over effects between years, communities in semi-arid to mesic Mediterranean regions persisted for the full simulation time (30 yr) in almost all combinations of change to MAP and DMR. In these regions, the communities produced enough seeds to compensate for low seed production in drought years. The storage effect of the seed bank to buffer environmental variability  has also been observed in the field in arid and semiarid ecosystems . Since annual precipitation was not autocorrelated, sequences of drought years were rare and community mass recovered within a few years. Most communities failed to recover in the arid region. Average seed production was not high enough to compensate for drought years. In natural sites in arid regions, populations persist for more than 30 yr because shrubs, crevices, and shade of rocks facilitate the establishment of annuals by indirectly increasing soil moisture . Furthermore, seed dispersal by animals, wind, and surface runoff from these protected sites would also increase community persistence . This continuous replenishment of the seed bank from external sources corresponds to the constant seed bank scenario and underlines the importance of external seed input.
Although populations in the Mediterranean regions persisted for at least 30 years, rainfall variability affected community dynamics across years. The annual variability of peak shoot mass decreased along with the variability of annual precipitation and LWP as MAP increased (Fig. 4). Generally, annual rainfall is less variable in more humid climates  and variability of grassland production increases with the variability of precipitation [53–55]. Increasing DMR increased the variability of annual rain volume but decreased the variability of LWP and peak shoot mass in drier regions and had neutral or positive effects in the more mesic regions. The negative correlation of DMR with the variability of community mass via LWP on the decadal time scale was in line with the beneficial effect of DMR on community mass at the annual scale, despite the increase of annual rainfall variability with DMR. This supports Reynolds et al.  that in terms of the pulse-reserve hypothesis  "pulse" ([sensu 15] must refer directly to the resource as it is available to the plant and not to the resource as it is supplied. In this case, the pulse must be sought in the soil moisture reservoir and not in the rain.
How does precipitation change projected by global climate models translate into community change along the gradient? The RegCM3 climate model generally predicts a 10 to 20% increase of DMR for the Middle East, but change of MAP varies between -20 and +10% [56–58]. Although the community response to an increase of DMR would have an opposite trend to the response to a decrease in MAP, the latter is much stronger. Consequently, the net effect is negative unless the decrease in MAP is small and the increase in DMR is strong. At sites with soils that cause a negative response of community mass to the increase of DMR, however, the effect of decreasing MAP would be amplified.
Model extensions and suggestions for further research
The present model concentrates on simulating annuals as one functional group. Apparently, the general results also hold for perennial herbs . Simulating this growth form would require an adaptation of parameters but only small changes to model structure. Vegetative reproduction would have to be introduced and the possibility that one individual may occupy more space than one grid cell. Although herbs take up water mostly from upper soil layers , potential uptake from deeper layers could be explored by increasing the parameter for layer thickness or by adding additional layers, requiring just small changes to the soil and plant modules. The model already includes structures for simulating competition for light and nitrogen and interactions with shrubs, although they have not yet been validated. This would allow studying competitive relationships with woody plants in relation to soil moisture distribution . The model can be parameterized for different annual species. Simulations for characteristic species from each climatic region showed that all species responded in a similar way to rain pattern, but some were more sensitive than others . Extending the model to simulate individuals of several species or functional types simultaneously would allow studying effects of rain pattern on, e.g., community composition or invisibility. This would improve our understanding of how different functional types coexist along a gradient of aridity . Specifically, one could test mechanisms for coexistence proposed by a one-dimensional model for plant functional types in deserts , especially seasonal differences in water use, differences in preparedness for water flushes, and differences in spatial uptake of water from upper soil layers. The model can also be used in a hierarchical modelling approach  to specify the reactions of the herbaceous vegetation for a larger area and coarser scale. For this end I carried out simulations for specific seed bank and rain volume scenarios to obtain transition probabilities of biomass production that were used in models of landscape productivity [63, 64].
Among the populations along the humidity gradient from arid to mesic Mediterranean those in the semi-arid and dry Mediterranean regions were the most sensitive to relative changes in annual precipitation. Changes to daily rain patterns had marked effects on community biomass in more arid regions through the amount of water stored in the soil. The same amount of annual precipitation distributed as sparser, heavier rainfall events recharged the soil water storage in the upper 15–30 cm more effectively than frequent, light rainfall events. Water recharge was increased by a high sand content that increased infiltration rates and reduced hydraulic conductivity in dried out surface layers. In mesic regions changes to daily rain patterns had small or no effects, because the soil remained moist enough between rain events to sustain growth. The effects of daily rain pattern were smaller but of the same order of magnitude as changes to mean annual precipitation. Therefore, these effects should be included in assessments of climate change.
The standard model parameters refer to four experimental sites along a 245-km gradient of humidity in Israel: an arid (AR) site near Sedé Boqér (N 30°52' E 34°46', 470 m a.s.l., 90 mm MAP), a semi-arid (SA) site near Lahav (N 31°23' E 34°54', 590 m a.s.l., 300 mm MAP), a "typical" Mediterranean (TM) site near Matta' (N 31°42' E 35°3', 620 m a.s.l., 540 mm MAP), and a mesic Mediterranean (MM) site near 'En Ya'aqov (N 33°0' E 35°14', 500 m a.s.l., 780 mm MAP). Details of rain fall frequency and distribution at nearby climate stations can be found in Additional file 5: ClimateStations.pdf. The soil types range from sandy loam at the AR site to loam at SA and clay at the MM and TM sites (Sarah Pariente, pers. comm.). Mean annual temperatures vary between 19°C and 21°C. Shrub cover increases along the gradient from <5% (AR) to about 80% (MM). Space between shrubs is covered mostly by annuals, while the percentage of bare ground decreases from about 95% to about 5% along the gradient. The vegetation and experimental setup at the sites are described in greater detail by Holzapfel et al. .
The four experimental sites are spread more or less evenly along the humidity gradient with a gap between the SA and TM site. Therefore, I used a hypothetical "dry Mediterranean" site (DM) with interpolated site characteristics in the simulations.
The model (version 2.6.0; Additional file 6: AnnualsModelSourceCode.zip) simulates the effect of rainfall variability on soil moisture and annual vegetation. The model is spatially explicit  because it includes the competitive interactions among individual neighbouring annual plants and uses principles of cellular automata models. The model's spatial extent is a soil surface of 25 cm × 25 cm per simulated site, divided into 1-cm2 cells representing soil columns (Fig. 1A). Each soil column contains a seed bank and space for one plant to reach maturity. The model consists of four modules describing the dynamics of climate, soil, seed bank, and plants (Fig. 1B).
I used five parameter sets to represent the different climates, soils, and plant communities in the five modelled climatic regions (Tab. 1; Additional file 1: ParametersSoil.pdf, Additional file 2: ParametersSeedBank.pdf, Additional file 3: ParametersPlants.pdf). The parameters were determined from experiments at each field site, greenhouse experiments, and published data or by choosing values that reproduced field data (for details see Detailed model description, below). All modelled processes are defined by physical or logical rules for time steps of one day, except for seed dispersal and seed bank mortality, which are evaluated in annual time steps.
The climate module (for details see Detailed model description – Climate module, below) determines daily temperature and precipitation, which are equally distributed in the grid. Temperature is calculated deterministically by a cosine function. Daily values do not vary among years. In contrast, rainfall is generated stochastically by the ReGen algorithm  (see below: Precipitation scenarios).
The soil module (for details see Detailed model description – Soil module, below) simulates water infiltration along a slope. A part of the surface water infiltrates into a soil column that consists of one organic (O) layer and four mineral (A) layers. In the simulations excess rain (runoff) was removed from the grid because I was interested in the direct effect of rainfall volume on each cell and similar moisture conditions across the grid. Water is lost to the atmosphere by evapotranspiration. Soil characteristics like hydraulic conductivity, field capacity, residual (minimum) water content, and the parameters of the soil moisture retention curve  are kept constant during simulations and are the same for each layer. Infiltration coefficients decrease with soil moisture and increase with vegetation cover.
The seed bank module (for details see Detailed model description – Seed bank module, below) simulates granivory, seed bank mortality, and germination. Before dispersed seeds enter the seed bank, a portion of the seeds in the seed bank dies because of decay, burial, or other causes. A portion of dispersed seeds is eaten before they enter the seed bank if their density is higher than the giving-up density of granivores. I implemented seasonal seed dormancy as an inverse U-shaped function. Germination of germinable seeds is determined by a hydrothermal time approach [41, 67] that incorporates the average soil moisture and temperature before germination. The actual germination fraction is further controlled by the local density of seeds [42, 43].
The plant module (for details see Detailed model description – Plant module, below) simulates the growth of individual annual plants. Plant growth is sigmoid, defined by maximum mass and maximum growth rate. The actual growth rate depends on temperature, water availability, the species' permanent wilting point (PWP), and net size-asymmetric competition  with neighbours within a radius of up to 2.5 cm. A fixed portion of produced biomass is allocated to reproduction and converted to seeds when the plant dies. Plant mortality is caused by lack of water. Seeds are dispersed with a negative exponential distribution, i.e., the likelihood to disperse farther from the mother plant decreases with distance.
The simulation of soil moisture dynamics and vegetation production by the model was validated by comparison with independent data that were not used to parameterize the model (for details see Additional file 7: Validation.pdf). The sensitivity of the model to changes of the model parameters was tested separately for the abiotic and biotic modules (see Additional file 8: SensitivityAnalyses.pdf).
My model does not take into account increasing annual mean temperatures, CO2 concentrations or nitrogen deposition from the atmosphere. A multi-factorial field experiment including these factors and precipitation in an annual grassland showed that interactions among these factors become rare over time and nitrogen deposition produced the strongest effect . Nonetheless, in the long-term, aboveground production of herbaceous arid communities is most strongly correlated with annual precipitation , which is also shown by the tight correlation of community biomass with annual precipitation along the climate gradient (see Validation, below).
I varied the distribution of volume and frequency of daily rainfall for five sets of nominal mean annual precipitation (100, 300, 450, 600, and 800 mm – corresponding roughly to the mean annual precipitation of the field sites) to investigate the effect of the daily rain pattern on the longest wet period (LWP), seedling density, and peak shoot mass of populations of annual plants. The change of daily rain pattern was achieved by increasing the seasonal amplitude of daily mean rain volume (DMR) by -20%, -10%, 0%, 10%, and 20% in the ReGen time series generator . Increasing DMR prolongs intervals between rains, reduces the contribution of days with light rains and increases the contribution of rainstorms to the annual rain volume . In order to compare the magnitude of this treatment, I also varied the mean annual precipitation volume (MAP) by -20%, -10%, 0%, 10%, and 20% in factorial combination with DMR. I used relative instead of absolute changes of precipitation parameters because regional climate models indicate that the changes will be tightly correlated with current values and that relative changes are remarkably constant along the precipitation gradient (M. Köchy, unpublished results). For example, the RegCM3 circulation model projects for the pessimistic A2 climate change scenario an average increase of DMR by 3% and a decrease of MAP by 13% for 2070–2100 for Israel North of the Negev desert, whereas the optimistic B2 scenario results in an average increase of DMR by 27% and an increase of MAP by 12% between 31° and 33° N (calculations based on [56–58]).
Seed bank scenarios
For the assessment of rain pattern effects on LWP, seedling density, and peak biomass, seed bank density in each region was reset to a constant value at the beginning of each vegetation year (AR: 2000, SA: 16000, DM: 17000, TM: 18000, MM: 20000 seeds/m2 corresponding to typical values observed in the field). This constant-seed bank scenario facilitates comparisons among years because carry-over effects of seed bank density due to wet and dry years and mortality are excluded. In a second set of simulations I examined the effect of change of daily rain pattern on population persistence. For population persistence inter-annual seed bank dynamics carry-over effects are important and were included (dynamic-seed bank scenario). To quantify community persistence I recorded the year in which the seed bank had dropped to zero or whether the population survived for the total simulated period.
I simulated 30 years of stochastic rainfall for each combination of levels of DMR and MAP in each climatic region and both seed bank scenarios. This was repeated five times for each combination. I used different random spatial distributions of soil surface heterogeneity in each simulation run. The great detail simulated by the model prevented me from storing daily data and only annual data across the whole lattice was recorded. Therefore, I sometimes phrase interpretations of results cautiously where omniscience of causality in the model could be expected. Among the many vegetation and soil variables simulated by the model, I concentrated on formally analyzing the effect of rain pattern on annual longest wet period (LWP), seedling density, and peak shoot mass (maximum annual aboveground biomass) because each is an integrative measure of, respectively, available soil moisture, establishment, and community performance. The effect of daily rain pattern on LWP, seedling density, and peak shoot mass for the constant-seed bank scenario was examined with analyses of covariance (ANCOVA) across regions. The full-factorial ANCOVA used change of DMR and change of MAP as covariates and regions as a nominal factor. As the data were very heteroscedastic, I used the means across years for each combination of daily pattern × mean annual precipitation × climatic region × repetition. Slopes and means were compared post-hoc using 95% confidence intervals. In addition, I calculated the change of peak shoot mass per 10% change of DMR and per 10% change of MAP to facilitate the comparison of effect sizes.
In the dynamic-seed bank scenario few populations persisted for the total simulation time (30 yr) in the arid region in contrast to the more humid regions. Therefore, an ANCOVA of persistence comprising all sites as for the constant seed bank scenario was not meaningful and I analyzed the arid region only. This turned the ANCOVA for the arid region into a multiple regression. For the other four regions, where time of persistence was almost always the full 30 years of the simulation runs, I compared the coefficient of variation of annual rain volume, longest wet period, and peak shoot mass using an ANCOVA to test how changes in rain pattern affected environmental and production variability. In addition I tested whether the CVs were correlated.
Detailed model description
The climate module generates temperature and precipitation for each day d. Mean daily temperature Td is a function of a seasonal index (S): Td = Ta + r·S, where Ta is the mean and r half the range of monthly mean temperatures (Tab. 1). Index S is defined as a cosine function with maximum in mid July, minimum in mid January, and starting on August 1 when rainfall probability is lowest along the climatic gradient. Thus, S = cos([d + dm]·2π/365.25) with dm = 15.625.
Stochastic time series were determined by the ReGen algorithm , a so-called two-part model . It uses a regular bell-shaped seasonal function that determines the daily mean likelihood averaged across years that a given day of the year is a rainy day. A flattened bell-shaped seasonal function determines the daily mean rain volume (DMR), i.e., the mean rain volume for a given day of the year averaged for rainy days across years. The DMR for a given rainy day of the year is used as the mean of a random negative exponential distribution from which the rain volume on that day is drawn. The parameters of each of the two seasonal functions were determined from monthly aggregated historic data from climate stations with 100 to 800 mm MAP in the interior of Israel. The time series do not include autocorrelation between annual rain volumes because autocorrelation in the historic data was not significant. ReGen allows the manipulation of DMR by ± 30%. In order to keep MAP at the specified value, the change of DMR is compensated by decreasing the occurrence of rainy days
Rainfall distribution at this spatial scale is assumed to be equal, thus all soil columns receive the same amount of rain. I assume that surface unevenness like small depressions or individual stones randomly change water reaching the soil by a fraction of ± H%. This fraction H is constant for each cell across years in each simulation run. If, however, the rainfall is >10 mm/day, I assume that this unevenness has no effect. Water that has reached the soil surface partly infiltrates the soil, the rest (runoff) is removed from the cell. Runoff in the model was removed in order to exclude effects on soil water due to relative location on the simulated grid. The amount of infiltrating water is calculated as volumein = [ICdry - v·ΔIC·min(θA1, θA2)]·volumesurface, where ICdry is a soil-specific infiltration coefficient (infiltrated volume/daily rain volume), v is a fitted parameter accounting for higher soil infiltration in vegetated (v = 0.9) than in bare soil (1.0) , ΔIC is the difference to the infiltration rate in a saturated soil, and θ is the volumetric water content of the soil. Thus, the infiltration coefficient decreases as the upper two A layers become wet. I obtained ICdry by fitting the model output for annual runoff to the mean annual value measured in 0.15 m2 runoff plots at each site in October 2002–April 2003 (Sarah Pariente & H. Lavee, unpublished data). Parameter ΔIC was calculated from runoff experiments along an aridity gradient near Jerusalem .
Infiltrating water fills up the soil downwards layer by layer. This tipping-bucket approach is sufficient for simulating water distribution in the soil after one day . Each layer cannot contain more water corresponding to a soil water potential Ψ of -0.01 MPa, that is, between field capacity (-0.036 MPa) and saturation (0 MPa). This value was suggested by measurements of the A1 layer at the typical and mesic Mediterranean field sites. The change of soil properties between two layers can reduce the permeability P of the soil. In the model this is implemented by letting only a fraction of water (Pc) flow from the A4 into the C layer. This becomes important only in the mesic Mediterranean region where the soil of the field site is compacted (H. Lavee, pers. comm.). There, I assume a permeability Pc of 70%. At the other sites, the C layer is more porous and I assume 100% permeability in these regions. Processes in the C layer are not considered in the model.
If the water remaining in the A4 layer is greater than its capacity, excess water will back up and fill the upper layers to their capacity and, if there is still excess water, contribute to surface runoff. Thus, runoff in the model can be generated in two different ways, infiltration excess (Hortonian runoff) and saturation excess. The model does not provide a mechanism for interflow (water movement in the soil parallel to the surface), because I assume that lateral differences in soil moisture are much smaller than vertical ones.
Daily change of simulated soil water potential by transpiration depending on aboveground mass.
Total shoot mass (mt,d)
no effect of transpiration
15 cm – C
>16 ... >30 cm
no effect of transpiration
Beyond the quick process of infiltration driven by gravitation potential, I also consider the slower process of water movement driven by differences in water potential in the O and A layers. This water flux is limited by the moisture- and texture-dependent hydraulic conductivity of the soil . In the model, I approximate this non-linear process by first calculating the amount of water (ΔV) necessary to change the potential of the thinner of two layers to their geometric mean weighted by layer thickness. This amount is reduced by a conductivity factor KA that varies with soil texture as . This semi-mechanistic approach prevents a drier thinner layer from becoming saturated or from drying out too much before a moister thicker layer's potential has balanced. The change of water content is calculated layer by layer downwards. As a result of balancing the potential, evaporation sucks water upwards, water reaches previously dry, lower layers, and the potentials of adjacent layers become more similar over time. The mechanism also reproduces the positive correlation between soil conductivity and water content , resulting in a lower evaporation from the top layer when the soil is dry than when it is wet. Water does not flow between the A4 and C layer by this mechanism.
where θ is the actual soil moisture (cm3/cm3), θS soil moisture at saturation, θR residual moisture, and α and β are parameters describing the curvature of the function. In the model I use a different set of parameters for each region (Additional file 1: ParametersSoil.pdf). α and β were fit to measurements of moisture and water potential on similar soils along a gradient of aridity (Sarah Pariente, unpubl. data). For θS I used a value slightly higher than the highest soil moisture actually observed at the experimental sites. θR was estimated based on values for standard soils . Soil water potential was limited to range between 0 (saturation) and -150 MPa (Ψmin), i.e., I disregard situations of ponding and extreme desiccation.
Seed bank module
where Ψb50 is the mean Ψb of the seeds with standard variation σ(Ψb). Mean minimum water potential Ψb50 has been found to vary among species , populations , and sites , whereas Tb is assumed mostly constant within species . Germination fraction of annuals in semi-arid regions is strongly density dependent [42, 43]. This is incorporated in the model by calculating the maximum seedling density (Dmax) within a distance of two cells around a target cell (18 neighbours) on a hexagonal lattice as Dmax = N(a ± Δa), where N is the mean number of germinable seeds per square metre in the neighbourhood, a is the density dependence (≈0.5...1.0 ) with random uniform variation Δa. Which of the germinable seeds in a neighbourhood will germinate is decided by chance with the probability adjusted to match Dmax. Germinated seeds prevent the germination of other seeds in the same cell so that only one plant can establish per square centimetre at maximum density. This conforms to plant densities observed in the field (pers. obs.) and generates a permanent seed bank.
Competition reduces resource uptake und thus the relative growth rate (eq. 3a). I use a simplified field-of-neighbourhood approach [87, 88] to simulate net size-asymmetric competition . Each individual is assumed to acquire above- and belowground resources from its 'home' cell, and the 18 surrounding neighbour cells in a hexagonal lattice if the individual's mass is greater than an upper critical mass (mC2); else, the individual can only access the six adjacent cells. If the mass is smaller than a lower critical mass (mC1), the individual's uptake is restricted to its home cell. The maximum area of 19 cm2 (19 cells) is 2–3 times greater than the average area per mature individual (7–9 cm2; Hadas Parag, unpublished data) in the SA, TM, and MM field site in spring 2003. (The average area per mature individual in the AR site was 255 cm2.) Although the radius of this area of potential interaction is smaller than the average lateral spread of roots  or shoots of annuals, it will likely include the zone of strongest interaction. For simplicity of programming I assume that each neighbour cell has the same resource availability as the home cell. Resource uptake from the home cell is exclusive to the resident individual but proportional to the individual's competitiveness (C) for belowground resources, ranging from 0 to 100%. This means that a plant is most competitive if it takes up resources only from the home cell, where it cannot be suppressed by neighbours. Resource uptake from each unoccupied neighbouring cell is relative to the resident individual's mass compared to the total mass of plants sharing access to the neighbour cell. If competition reduces relative resource uptake below 40% (density-dependent mortality, MD), the plant dies. Competitive ability is made to have a size-asymmetric effect on growth rate by using its square in eq. 3a.
For a mature plant with mt,d > mmature, a fixed portion of the daily produced aboveground biomass (Δm) is allocated to flowers and seed production without contributing to new vegetative growth. mmature is increasing along the humidity gradient because individuals in arid regions are smaller and often flower and reproduce at smaller size than in more humid regions . Plants are assumed to wilt when their daily production is < 1‰. This may be caused by drought, competition, or attaining maximum mass mmax. The reproductive biomass of wilted plants is converted to number of seeds by dividing by fruit mass. Here, 'fruit mass' is the average mass of reproductive structures per seed in a species.
The number of seeds dispersed from a plant generally decreases with distance. Most annuals in the study region lack a specialized mode of dispersal, their seeds just fall off the plant. In the model, I assume that seeds are dispersed into neighbouring cells with random direction and a mean distance of 5 cm with a negative exponential distribution according to consensus estimates for dominant species by several field experts (C. Holzapfel, M. Sternberg, M. Petrů, pers. comm.).
I thank Irit Konsens, Hadas Parag, Marcelo Sternberg, Cristina Cardiel, Martina Petrů, Eyal Sachs, Claus Holzapfel, Sara Pariente, Hanoch Lavee, Jaime Kigel, and Katja Tielbörger for measuring, collecting, processing, and sharing field data. They also discussed with us the model outputs and helped refining model rules. I appreciate critical comments on the manuscript by the members of the Plant Ecology and Nature Conservation section at Potsdam University, especially Florian Jeltsch, Monika Schwager, Frank Schurr, Britta Tietjen, and Richard Walters, and by four anonymous reviewers. This study is part of the GLOWA Jordan River project financed by the German Federal Ministry of Science and Education (BMBF), contract 01LW0306(A). The author alone is responsible for the content of this publication.
- Noy-Meir I: Desert Ecosystems: environment and producers. Ann Rev Ecol Syst. 1973, 4: 25-51. 10.1146/annurev.es.04.110173.000325.View ArticleGoogle Scholar
- Wilson SD: Heterogeneity, diversity and scale in plant communities. Ecological consequences of habitat heterogeneity. Edited by: Hutchings M, John EA, Stewart AJ. 2000, Oxford (U.K.): Blackwell Science, 53-69.Google Scholar
- Hansson L: Why ecology fails at application: should we consider variability more than regularity?. Oikos. 2003, 100: 624-627. 10.1034/j.1600-0706.2003.12479.x.View ArticleGoogle Scholar
- Weltzin JF, Loik ME, Schwinning S, Williams DG, Fay PA, Haddad BM: Assessing the response of terrestrial ecosystems to potential changes in precipitation. BioScience. 2003, 53: 941-952. 10.1641/0006-3568(2003)053[0941:ATROTE]2.0.CO;2.View ArticleGoogle Scholar
- Davidowitz G: Does precipitation variability increase from mesic to xeric biomes?. Glob Ecol Biogeogr. 2002, 11: 143-154. 10.1046/j.1466-822X.2002.00271.x.View ArticleGoogle Scholar
- Easterling DR, Meehl GA, Parmesan C, Changnon SA, Karl T, Mearns LO: Climate extremes: observations, modelling, and impacts. Science. 2000, 289: 2068-2074. 10.1126/science.289.5487.2068.View ArticlePubMedGoogle Scholar
- Alpert P, Ben-Gai T, Baharad A, Benjamini Y, Yekutielli D, Colacino M: The paradoxical increase of Mediterranean extreme daily rainfall in spite of decrease in total values. Geophys Res Let. 2002, 29: 31.1-31.4. 10.1029/2001GL013554.View ArticleGoogle Scholar
- Govindasamy B, Duffy PB, Coquard J: High resolution simulations of global climate, part 2: effects of increased greenhouse gases. Clim Dyn. 2003, 21: 391-404. 10.1007/s00382-003-0340-6.View ArticleGoogle Scholar
- Fay PA, Carlisle JD, Knapp AK, Blair JM, Collins SL: Productivity responses to altered rainfall patterns in a C4-dominated grassland. Oecologia. 2003, 137: 245-251. 10.1007/s00442-003-1331-3.View ArticlePubMedGoogle Scholar
- Laporte MF, Duchesne LC, Wetzel W: Effect of rainfall patterns on soil surface CO2 efflux, soil moisture, soil temperature and plant growth in a grassland ecosystem of northern Ontario, Canada: implications for climate change. BMC Ecol. 2002, 2: 10-10.1186/1472-6785-2-10.PubMed CentralView ArticlePubMedGoogle Scholar
- Potts DL, Huxman TE, Enquist BJ, Weltzin JF, Williams DG: Resilience and resistance of ecosystem functional response to a precipitation pulse in a semi-arid grassland. J Ecol. 2006, 94: 23-30. 10.1111/j.1365-2745.2005.01060.x.View ArticleGoogle Scholar
- Weltzin JF, McPherson GR: Implications of precipitation redistribution for shifts in temperate savanna. Ecology. 2000, 81: 1902-1913. 10.1890/0012-9658.View ArticleGoogle Scholar
- Williams KJ, Wilsey BJ, McNaughton SJ, Banyikwa FF: Temporally variable rainfall does not limit yields of Serengeti grasses. Oikos. 1998, 81: 463-470. 10.2307/3546768.View ArticleGoogle Scholar
- Novoplansky A, Goldberg DE: Effects of water pulsing on individual performance and competitive hierarchies in plants. J Veg Sci. 2001, 12: 199-208. 10.2307/3236604.View ArticleGoogle Scholar
- Reynolds JF, Kemp PR, Ogle K, Fernandez RJ: Modifying the 'pulse-reserve' paradigm for deserts of North America: precipitation pulses, soil water, and plant responses. Oecologia. 2004, 141: 194-210. 10.1007/s00442-004-1524-4.View ArticlePubMedGoogle Scholar
- Sher AA, Goldberg DE, Novoplansky A: The effect of mean and variance in resource supply on survival of annuals from Mediterranean and desert environments. Oecologia. 2004, 141: 353-362. 10.1007/s00442-003-1435-9.View ArticlePubMedGoogle Scholar
- Clothier BE, Green SR: Rootzone processes and the efficient use of irrigation water. Agric Water Manage. 1994, 25: 1-12. 10.1016/0378-3774(94)90048-5.View ArticleGoogle Scholar
- Mermoud A, Tamini TD, Yacouba H: Impacts of different irrigation schedules on the water balance components of an onion crop in a semi-arid zone. Agric Water Manage. 2005, 77: 282-295. 10.1016/j.agwat.2004.09.033.View ArticleGoogle Scholar
- Harper CW, Blair JM, Fay PA, Knapp AK, Carlisle JD: Increased rainfall variability and reduced rainfall amount decreases soil CO2 flux in a grassland ecosystem. Glob Change Biol. 2005, 11: 322-334. 10.1111/j.1365-2486.2005.00899.x.View ArticleGoogle Scholar
- Köchy M: Stochastic time series of daily precipitation for the interior of Israel. Israel Journal of Earth Sciences. 2006, 55: 103-109. 10.1560/IJES_55_2_103.View ArticleGoogle Scholar
- Jankju-Borzelabad M: Responses of arid land plants to simulated changes in rain frequency and intensity. 90th Annual Meeting of the Ecological Society of America held jointly with the IX International Congress of Ecology, 7.-12. August 2005, Montréal (Canada). 2005, Washington, D.C. (U.S.A.): Ecological Society of America, 304-Google Scholar
- English NB, Weltzin JF, Fravolini A, Thomas L, Williams DG: The influence of soil texture and vegetation on soil moisture under rainout shelters in a semi-desert grassland. J Arid Environ. 2005, 63: 324-343. 10.1016/j.jaridenv.2005.03.013.View ArticleGoogle Scholar
- Lavorel S, Canadell J, Rambal S, Terrades J: Mediterranean terrestrial ecosystems: research priorities on global change effects. Global Ecology and Biogeography Letters. 1998, 7: 157-166. 10.2307/2997371.View ArticleGoogle Scholar
- Mouillot F, Rambal S, Joffre R: Simulating climate change impacts on fire frequency and vegetation dynamics in a Mediterranean-type ecosystem. Glob Change Biol. 2002, 8: 423-437. 10.1046/j.1365-2486.2002.00494.x.View ArticleGoogle Scholar
- Hoff C, Rambal S: An examination of the interaction between climate, soil and leaf area index in a Quercus ilex ecosystem. Ann For Sci. 2003, 60: 153-161. 10.1051/forest:2003008.View ArticleGoogle Scholar
- De Lillis M, Manes F, Tufaro F: Root plasticity in Mediterranean herbaceous species. Plant Biosystems. 2005, 139: 214-221. 10.1080/11263500500160542.View ArticleGoogle Scholar
- Schenk HJ, Jackson RB: Rooting depth, lateral root spreads and below-ground/above-ground allometries of plants in water-limited ecosystems. J Ecol. 2002, 90: 480-494. 10.1046/j.1365-2745.2002.00682.x.View ArticleGoogle Scholar
- Holzapfel C, Tielbörger K, Parag HA, Kigel J, Sternberg M: Annual plant-shrub interactions along an aridity gradient. Basic Appl Ecol. 2006, 7: 268-279. 10.1016/j.baae.2005.08.003.View ArticleGoogle Scholar
- Petrů M, Tielbörger K, Belkin R, Sternberg M, Jeltsch F: Life history variation in an annual plant under two opposing environmental constraints along an aridity gradient. Ecography. 2006, 29: 66-74. 10.1111/j.2005.0906-7590.04310.x.View ArticleGoogle Scholar
- Schiffers K, Tielbörger K: Ontogenetic shifts in interactions among annual plants. J Ecol. 2006, 94: 336-341. 10.1111/j.1365-2745.2006.01097.x.View ArticleGoogle Scholar
- Dunne JA, Saleska SR, Fischer ML, Harte J: Integrating experimental and gradient methods in ecological climate change research. Ecology. 2004, 85: 904-916. 10.1890/03-8003.View ArticleGoogle Scholar
- Larcher W: Ökologie der Pflanzen auf physiologischer Grundlage. (UTB 232). 1984, Stuttgart: UlmerGoogle Scholar
- Jenny M, Smettan U, Facklam-Moniak M: Soil-vegetation relationship at several arid microsites in the Wadi Araba (Jordan). Vegetatio. 1990, 89: 149-164. 10.1007/BF00032167.View ArticleGoogle Scholar
- Holmes TH, Rice KJ: Patterns of growth and soil-water utilization in some exotic annuals and native perennial bunchgrasses of California. Ann Bot. 1996, 78: 233-243. 10.1006/anbo.1996.0117.View ArticleGoogle Scholar
- Tuljapurkar SD, Orzack SH: Population dynamics in variable environments. I. Long-run growth rates and extinction. Theoretical Population Biology. 1980, 18: 314-342. 10.1016/0040-5809(80)90057-X.View ArticleGoogle Scholar
- Drake JM, Lodge DM: Effects of environmental variation on extinction and establishment. Ecol Lett. 2004, 7: 26-30. 10.1046/j.1461-0248.2003.00546.x.View ArticleGoogle Scholar
- Nippert JB, Knapp A: Linking water uptake with rooting patterns in grassland species. Oecologia. 2007, 153: 261-272. 10.1007/s00442-007-0745-8.View ArticlePubMedGoogle Scholar
- Köchy M, Jeltsch F: Standing crop of an abundant grass species is not affected by rainfall intensity. 35 Jahrestagung der Gesellschaft für Ökologie, 19.-23. September 2005, Regensburg (Germany). Verh Ges Ökol. 2005, 36: 250-Google Scholar
- Austin AT, Yahdjian L, Stark JM, Belnap J, Porporato A, Norton U, Ravetta D, Schaeffer SM: Water pulses and biogeochemical cycles in arid and semiarid ecosystems. Oecologia. 2004, 141: 221-235. 10.1007/s00442-004-1519-1.View ArticlePubMedGoogle Scholar
- Sala OE, Parton WJ, Joyce LA, Lauenroth WK: Primary production of the central grassland region of the United States. Ecology. 1988, 69: 40-45. 10.2307/1943158.View ArticleGoogle Scholar
- Köchy M, Tielbörger K: Hydrothermal time model of germination: parameters for 36 Mediterranean annual species based on a simplified approach. Basic Appl Ecol. 2007, 8: 171-182. 10.1016/j.baae.2006.04.002.View ArticleGoogle Scholar
- Goldberg DE, Turkington R, Olsvig-Whittaker L, Dyer AR: Density dependence in an annual plant community: variation among live history stages. Ecol Monogr. 2001, 71: 423-446. 10.1890/0012-9615.View ArticleGoogle Scholar
- Tielbörger K, Valleriani A: Can seeds predict their future? Germination strategies of density-regulated desert annuals. Oikos. 2005, 111: 235-244. 10.1111/j.0030-1299.2005.14041.x.View ArticleGoogle Scholar
- Bergelson J, Perry R: Interspecific competition between seeds: relative planting date and density affect seedling emergence. Ecology. 1989, 70: 1639-1644. 10.2307/1938097.View ArticleGoogle Scholar
- Dyer AR, Fenech A, Rice KJ: Accelerated seedling emergence in interspecific competitive neighborhoods. Ecol Lett. 2000, 3: 523-529. 10.1046/j.1461-0248.2000.00187.x.View ArticleGoogle Scholar
- Inouye RS: Density-dependent germination response by seeds of desert annuals. Oecologia. 1980, 46: 235-238. 10.1007/BF00540131.View ArticleGoogle Scholar
- Murray BR: Density dependent germination and the role of seed leachate. Aust J Ecol. 1998, 23: 411-418. 10.1111/j.1442-9993.1998.tb00748.x.View ArticleGoogle Scholar
- Lundholm JT, Larson DW: Experimental separation of resource quantity from temporal variability: seedling responses to water pulses. Oecologia. 2004, 141: 346-352. 10.1007/s00442-003-1454-6.View ArticlePubMedGoogle Scholar
- Heisler JL, Knapp AK, Kelly GM, Blair JM: Climate change in the Great Plains: More extreme precipitation patterns impact ecosystem function differentally across the central US. Ecological Society of America/Society for Ecological Restoration Joint Meeting, San Jose, California, 5.-10. 2007, AugustGoogle Scholar
- Higgins SI, Pickett STA, Bond WJ: Predicting extinction risks for plants: environmental stochasticity can save declining populations. Trends Ecol Evol. 2000, 15: 516-520. 10.1016/S0169-5347(00)01993-5.View ArticlePubMedGoogle Scholar
- Holmgren M, Stapp P, Dickman C, Gracia C, Graham S, Gutiérrez JR, Hice C, Jaksic FM, Kelt DA, Letnic M, Lima M, López BC, Meserve PL, Milstead WB, Polis GA, Previtali MA, Richter M, Sabaté S, Squeo FA: Extreme climatic events shape arid and semiarid ecosystems. Front Ecol Environ. 2006, 4: 87-95. 10.1890/1540-9295(2006)004[0087:ECESAA]2.0.CO;2.View ArticleGoogle Scholar
- Shmida A, Wilson M: Biological determinants of species diversity. J Biogeogr. 1985, 12: 1-20. 10.2307/2845026.View ArticleGoogle Scholar
- Le Houérou HN, Bingham RL, Skerbek W: Relationship between the variability of primary production and the variability of annual precipitation in world arid lands. J Arid Environ. 1988, 5: 1-18.Google Scholar
- Fang J, Piao S, Tang Z, Peng C, Ji W: Interannual variability in net primary production and precipitation. Science. 2001, 293: 1723a-10.1126/science.293.5536.1723a.View ArticleGoogle Scholar
- Jobbágy EG, Sala OE, Paruelo JM: Patterns and controls of primary production in the Patagonian steppe: a remote sensing approach. Ecology. 2002, 83: 307-319. 10.1890/0012-9658.Google Scholar
- Giorgi F, Bi X, Pal JS: Mean, interannual variability and trends in a regional climate change experiment over Europe. I. Present-day climate (1961–1990). Clim Dyn. 2004, 22: 733-756. 10.1007/s00382-004-0409-x.View ArticleGoogle Scholar
- Giorgi F, Bi X, Pal J: Mean, interannual variability and trends in a regional climate change experiment over Europe. II: climate change scenarios (2071–2100). Clim Dyn. 2004, 23: 839-858. 10.1007/s00382-004-0467-0.View ArticleGoogle Scholar
- Alpert P, Krichak SO, Dayan M, Osetinsky I: Regional climate modeling over the E. Mediterranean. Scientific Assembly of the International Association of Meterology and Atmospheric Science, Beijing, China, 2.-11. 2005, C51-AugustGoogle Scholar
- Nippert JB, Knapp AK: Soil water partitioning contributes to species coexistence in tallgrass prairie. Oikos. 2007, 116: 1017-1029. 10.1111/j.0030-1299.2007.15630.x.View ArticleGoogle Scholar
- Köchy M, Jeltsch F: Effect of rainfall characteristics on performance of ten annual species. 89. Annual Meeting of the Ecological Society of America, Portland, Oregon (U.S.A.), 1.-6. 2004, AugustGoogle Scholar
- Ogle K, Reynolds JF: Plant responses to precipitation in desert ecosystems: integrating functional types, pulses, thresholds, and delays. Oecologia. 2004, 141: 282-294. 10.1007/s00442-004-1507-5.View ArticlePubMedGoogle Scholar
- Jeltsch F, Moloney K, Schurr F, Köchy M, Schwager M: The state of plant population modelling in light of environmental change. Persp Plant Ecol Evol Syst. 2008, 9: 171-189. 10.1016/j.ppees.2007.11.004.View ArticleGoogle Scholar
- Köchy M, Mathaj M, Jeltsch F, Malkinson D: Resilience of stocking capacity to changing climate in arid to Mediterranean landscapes. Reg Environ Change. 10.1007/s10113-008-0048-6.Google Scholar
- Mathaj M: Modellierung von Vegetation und Erosion entlang eines Klimagradienten von mediterran bis semiarid. Diplomarbeit (Msc thesis). 2007, Potsdam (Germany): Universität PotsdamGoogle Scholar
- Jeltsch F, Moloney KA: Spatially explicit vegetation models: what have we learned?. Progress in Botany. Edited by: Esser K, Lüttge U, Beyschlag W, Hellwig F. 2002, Berlin, Germany: Springer, 326-343.View ArticleGoogle Scholar
- van Genuchten MT: A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980, 44: 892-898.View ArticleGoogle Scholar
- Allen P: When and how many? Hydrothermal models and the prediction of seed germination. New Phytol. 2003, 158: 1-3. 10.1046/j.1469-8137.2003.00729.x.View ArticleGoogle Scholar
- Bauer S, Wyszomirski T, Berger U, Hildenbrandt H, Grimm V: Asymmetric competition as a natural outcome of neighbour interactions among plants: results from the field-of-neighbourhood modelling approach. Plant Ecol. 2004, 170: 135-145. 10.1023/B:VEGE.0000019041.42440.ea.View ArticleGoogle Scholar
- Dukes JS, Chiariello NR, Cleland EE, Moore LA, Shaw MR, Thayer S, Tobeck T, Mooney HA, Field C: Responses of grassland production to single and multiple global environmental changes. PLoS Biol. 2005, 3: e319-10.1371/journal.pbio.0030319.PubMed CentralView ArticlePubMedGoogle Scholar
- Lauenroth WK, Sala OE: Long-term forage production of North American shortgrass steppe. Ecol Appl. 1992, 2: 397-403. 10.2307/1941874.View ArticleGoogle Scholar
- Srikanthan R, McMahon TA: Stochastic generation of annual, monthly and daily climate data: A review. Hydrol Earth System Sci. 2001, 5: 653-670. [http://www.hydrol-earth-syst-sci.net/5/653/2001/hess-5-653-2001.pdf]View ArticleGoogle Scholar
- Durán Zuazo VH, Francia Martínez JR, Martínez Raya A: Impact of vegetative cover on runoff and soil erosion at hillslope scale in Lanjaron, Spain. Environmentalist. 2004, 24: 39-48. 10.1023/B:ENVR.0000046345.44569.35.View ArticleGoogle Scholar
- Cerdà A, Sarah P: The effect of climate on overland flow: laboratory rainfall simulation experiments. Revista Cuaternario y Geomorfología. 2000, 14: 63-75. [http://tierra.rediris.es/CuaternarioyGeomorfologia/images/vol14_1/cuaternario14(1-2)_05_.pdf]Google Scholar
- Baier W, Robertson GW: A new versatile soil moisture budget. Can J Plant Sci. 1966, 46: 299-315.View ArticleGoogle Scholar
- Meidan E: The effects of soil water potential on seed germination of four winter annuals in the Negev desert. J Arid Environ. 1990, 19: 77-83.Google Scholar
- Kigel J, Perevolotsky A, Roundy BA, Allen P, Mayer SE: Predicting seed-bank germination in semiarid rangelands under grazing. International Arid Lands Consortium Project 98R-27, Research Grant 191-1-610-8 Final report. 2002Google Scholar
- Maidment DR, (ed): Handbook of Hydrology. 1993, New York, U.S.A.: MacGraw HillGoogle Scholar
- Louda SM: Predation in the dynamics of seed regeneration. Ecology of seed banks. Edited by: Leck MA, Parker VT, Simpson RL. 1989, San Diego, California, U.S.A.: Academic Press, 25-51.Google Scholar
- Cabin RJ, Marshall DL, Mitchell RJ: The demographic role of soil seed banks. II. Investigations of the fate of experimental seeds of the desert mustard Lesquerella fendleri. J Ecol. 2000, 88: 293-302. 10.1046/j.1365-2745.2000.00444.x.View ArticleGoogle Scholar
- Lortie CJ, Ganey DT, Kotler BP: The effects of gerbil foraging on the natural seedbank and consequences on the annual plant community. Oikos. 2000, 90: 399-407. 10.1034/j.1600-0706.2000.900221.x.View ArticleGoogle Scholar
- Detrain C, Tasse O, Versaen M, Pasteels JM: A field assessment of optimal foraging in ants: trail patterns and seed retrieval by the European harvester ant Messor barbarus. Insectes sociaux. 2000, 47: 56-62. 10.1007/s000400050009.View ArticleGoogle Scholar
- Kunin WE: Density-dependent foraging in the harvester ant Messor ebeninus – 2 experiments. Oecologia. 1994, 98: 328-335. 10.1007/BF00324221.View ArticleGoogle Scholar
- Allen PS, Meyer SE, Khan MA: Hydrothermal time as a tool in comparative germination studies. Seed Biology: Advances and Applications. Edited by: Black M, Bradford KJ, Vázquez-Ramos J. 2000, Wallingford (U.K.): CAB International, 401-410.Google Scholar
- Bauer MC, Meyer SE, Allen PS: A simulation model to predict seed dormancy loss in the field for Bromus tectorum L. J Exp Bot. 1998, 49: 1235-1244. 10.1093/jexbot/49.324.1235.Google Scholar
- Meyer SE, Debaene-Gill SB, Allen PS: Using hydrothermal time concepts to model seed germination response to temperature, dormancy loss, and priming effects in Elymus elymoides. Seed Sci Res. 2000, 10: 213-223. 10.1017/S0960258500000246.View ArticleGoogle Scholar
- Majerus ME: Response of root and shoot growth of three grass species to decreases of soil water potential. J Range Manage. 1975, 28: 473-476. 10.2307/3897225.View ArticleGoogle Scholar
- Weiner J: A neighbourhood model of annual-plant interference. Ecology. 1982, 63: 1237-1241. 10.2307/1938849.View ArticleGoogle Scholar
- Berger U, Hildebrandt H, Grimm V: Towards a standard for the individual-based modeling of plant populations: self-thinning and the field of neighborhood approach. Natural Resource Modelling. 2002, 15: 39-51. 10.1216/nrm/1030539096.View ArticleGoogle Scholar
- Aronson JA, Kigel J, Shmida A, Klein J: Adaptive phenology of desert and mediterranean populations of annual plants grown with and without water stress. Oecologia. 1992, 89: 17-26. 10.1007/BF00319010.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.