Effects of simulated daily precipitation patterns on annual plant populations depend on life stage and climatic region
BMC Ecology volume 8, Article number: 4 (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 .
I use a process-based model (Fig. 1) to quantify the effect of different daily rainfall distributions and compare it to the effect of a change in mean annual amount on vegetation. The model simulates explicitly the response of soil moisture and individual plants to rainfall variability. In addition, I consider the effects on different life stages and spatial interactions (competition, dispersal). The model simulates vegetation in the Mediterranean region, situated between water-stressed subtropical and mesic temperate regions. Therefore, the region may be particularly sensitive to climate change effects . Earlier models addressing precipitation effects on Mediterranean vegetation focused on shrubs and trees in the western Mediterranean region [24, 25] but did not systematically consider the effect of daily within-year rainfall variability. Annual plants are the dominant or co-dominant life form at lower elevations in the eastern Mediterranean. Furthermore, annuals may respond more directly to water variability than perennials given the annuals' shallower and less plastic root system [26, 27]. My model is the first to systematically examine the effect of projected changes in daily rain pattern on natural vegetation along a climatic gradient from arid to mesic. Hence, the sensitivity of different climatic regions to daily rain patterns can be compared explicitly. The gradient is also used by complementary field studies [28–30], which provided the opportunity to combine simulations with field experiments for parameterization and validation of the model . The present study provides a comprehensive examination of net precipitation effects on representative annual communities by investigating soil processes, plant germination, growth, and seed bank dynamics in five climatic regions.
I assessed the importance of intra-annual rainfall distribution on soil moisture, seedling density, peak shoot mass, and population persistence by varying independently the amplitude of daily mean rain volume (DMR) and the mean annual precipitation (MAP) and comparing their relative effects. DMR and MAP were varied from -20% to +20% relative to current conditions. Daily and thus annual rain amounts were simulated as stochastic time series with specified means. I used relative changes in order to make the effects of DMR and MAP comparable and because my analysis of projections of regional climate models indicated that these factors change relative to the present value. Therefore, results are preferably expressed with regard to the relative change of DMR and MAP. The comparison of effects was conducted for one site in each of five regions (AR, SA, DM, TM, MM) along a climatic gradient in the interior of Israel (Tab. 1). One soil type and one artificial species with representative characteristics to match overall community structure per site were used in the simulations (Additional file 1: ParametersSoil.pdf, Additional file 2: ParametersSeedBank.pdf, Additional file 3: ParametersPlants.pdf). Artificial species differed only in few characteristics (soil moisture required for germination, permanent wilting point, minimum mass for reproduction) to facilitate comparisons. Each site was represented by a grid of 25 × 25 1-cm2 cells (Fig. 1A).
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.
LWP varied with DMR depending on interactions with regions and regions × MAP (Tab. 2). In addition, LWP varied with MAP depending on interactions with regions. The average LWP of each region increased in a slightly sigmoid way along the humidity gradient (Fig. 2A). Within sites, LWP increased with MAP. The 'absolute' rate of increase (ΔLWD/ΔMAP) diminished from the arid to the mesic Mediterranean region (Fig. 2B), whereas the 'relative' rate of increase (ΔLWD/10% MAP), indicating the effect of expected relative changes, peaked in the dry Mediterranean region (Additional file 4: RelativeChanges.pdf). LWP increased with DMR when MAP was <600 mm. Beyond this point, LWP did not change significantly with DMR. The highest relative rate of increase occurred around 200 mm MAP. Thus, the relative rate of increase rose with MAP in the arid region, decreased in the semiarid and dry Mediterranean regions, and did not change with MAP in the two most mesic sites.
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.
Simulations with a seed bank that carries over to the following year was conducted to assess the consequences of change of rain pattern on population persistence. In the arid region only 25 of 125 populations persisted until the end of the simulations (30 yr, Fig. 3A). In contrast, in the semiarid region 107 populations (Fig. 3B), in the dry Mediterranean region all but one population (combination -20% MAP, -10% DMR), and in the typical and mesic Mediterranean regions all populations persisted for 30 years. In the arid region, populations persisted with full growth cycles only in the combination +20% MAP, +20% DMR. In the other instances, populations survived up to 30 years via a dormant seed bank without germination. Thus, persistence decreased with increasing MAP (P = 0.04) and increased with the interaction between DMR and MAP (P = 0.001) because persistence was on average higher when both change of DMR and MAP were either both positive or both negative (Fig. 3A). DMR had no significant main effect. In the semiarid region, population persistence increased with DMR when MAP was reduced by 20% (Fig. 3B, P < 0.001). Persistence in other combinations and regions was not assessed statistically because the pattern was obvious.
I compared the effects of DMR and MAP on the coefficients of variation of annual rain volume (CVR), longest wet period (CVLWP), and peak shoot mass (CVmass) because greater environmental and population variability is theoretically related to shorter persistence [35, 36]. I did not include the arid region in this analysis because its CVs were based on less than half the number of years and CVmass was strongly influenced by a decline in mean mass over time. CVR varied among regions, decreased with change of MAP, and increased strongly with change of DMR (Tab. 3; each P < 0.0001). These three factors did not interact. CVR could be predicted by the regression equation CVR = 19.8 + [3.9 | SA, 1.5 | DM, -1.5 | TM, -3.8 | MM] - 0.084·ΔMAP + 0.128·ΔDMR (adjusted R2 = 0.64). CVLWP varied strongly among regions (Tab. 3; SA: 53, DM: 31, TM: 17, MM: 12) and decreased overall with increasing MAP. The rate of decrease with MAP within regions became smaller along the humidity gradient. The effect of DMR depended on region. CVLWP decreased with increasing DMR across the semiarid region, did not change with DMR across the dry Mediterranean region, increased with DMR across the typical Mediterranean region, and did not change with DMR across the mesic Mediterranean region. CVLWP increased with the interaction of MAP and DMR in the semiarid region but in the other regions the interaction effect was not significant. This indicates that water storage in the soil reduced the variability of the rain. The CV of peak shoot mass (CVmass) behaved in a similar way as CVLWP (Fig. 4; Tab. 3) except that all interactions were significant. CVmass decreased along the humidity gradient (Fig. 4A, Tab. 3), with increasing MAP and with increasing DMR. The rate of decrease with MAP became smaller along the humidity gradient. The confidence intervals of the rate of decrease included zero in the mesic Mediterranean region. The average effect of DMR differed among regions because its direction switched from negative to neutral along the humidity gradient (Fig. 4B). The interaction effect on CVmass by MAP and DMR was positive on average. In addition, the negative effect of DMR diminished with increasing MAP in the semiarid region only, causing a significant three-way interaction.
The CVs of all three variables were positively correlated among each other (Tab. 4). CVmass was generally correlated more strongly with CVLWP than with CVR. The partial correlation between CVmass and CVR, accounting for the correlation with CVLWP, was negative in the semiarid, dry, and typical Mediterranean regions and neutral in the mesic Mediterranean region. The negative partial correlation indicated that soil water storage (via LWP) also reduced the inter-annual effect of rain variability on the variability of shoot mass.
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
The soil module determines the water availability for plant growth by simulating water infiltration, surface runoff, and evapotranspiration of individual soil columns with a surface of 1 cm2 (one cell; Figs. 1, 5) with different standard parameters for each region (Additional file 1: ParametersSoil.pdf). A soil column consists of one organic ([O layer] 0 – 1 cm) and 4 mineral soil layers (A1: 1 – 5 cm, A2: 5 – 10 cm, A3: 10 – 15 cm, A4: 15 cm down to the weathered mineral layer [C layer] at 16 to 30 cm depending on region).
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.
Analyses of daily soil moisture in the O and A1 layer after heavy rainfalls followed by a dry spell at each of the experimental sites (Sarah Pariente & H. Lavee, unpublished data) and measurements of similar soils [75, Taylor et al. in 76] indicated that evaporation from the A1 layer can be represented by a linear decrease of soil water potential Ψ with time as Ψd+1 = EA1·Ψd, with daily evaporation factor EA1 ranging between 2 and 3, depending on site. Evaporation from the O layer was also linear but its evaporation factor EO (≈3) increased linearly with mean daily temperature Td (Ψd+1 = EO·Td·Ψd). Transpiration by plants is simulated by transpiration factors (Tab. 5) that increase the evaporation factors, differ among soil layers, and vary with aboveground mass to reflect the influence of root distribution. I assume that transpiration factors do not vary among regions because grasses were the dominant annual growth form at each field site.
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.
To calculate evaporation and water movement, volumetric water content in each soil layer is converted to water potential using a soil moisture retention curve according to van Genuchten :
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
The seed bank module simulates the dynamics of the seed bank by germination and seed mortality in annual time steps each summer (Fig. 6) with different standard parameters for each region (Additional file 2: ParametersSeedBank.pdf). Before dispersed seeds enter the modelled seed bank, the seed bank is decreased by 10% to reflect seed mortality, downward movement of seeds in the soil, and other processes that reduce seed viability. Most newly dispersed seeds in the field are eaten by ants or rodents [78, 79]. Rodents stop searching for seeds when their density is below a giving up density . I also assume a giving-up density for ants although evidence is inconsistent [81, 82]. I assume in the model that 25 – 35% of the seeds survive granivory and that granivores leave a cell when the seed density is <20'000/m2. These values were tuned so that simulated populations with unmodified regional parameters persisted for several decades. All seeds in the seed bank are considered viable but dormant at the beginning of the vegetation year (August). In the research area, germination events are observed only from October to March. Germinability of any species is highest in January (K. Tielbörger, J. Kigel, pers. comm.). In the model, these relations are expressed as an inverse-U shaped seasonal function of germinability G = (-S)0.1 if S < 0, else G = 0 (Fig. 7), where S is the seasonal index of the climate module.
The number of seeds per cell ready to germinate on a specific day (germinable seeds) is determined by their accumulated "hydrothermal time", θHT , a concept related to temperature sum or degree-days. Hydrothermal time is defined as θHT = ∑([ΨO - Ψb]·[Td-Tb]), where ΨO is the soil water potential in the top layer, Ψb is the minimum water potential required for germination, Td is mean daily temperature, and Tb is the species-specific minimum temperature required for germination. Hydrothermal time is accumulated daily if the increment is positive. As a strategy to prevent untimely germination, I assume that hydrothermal time is reset to zero if the increment has not been positive for more than seven days. When θHT is greater than a species-specific threshold, the seed is ready to germinate. If a cell contains more than one seed, the germination fraction (g) is defined as
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.
The plant module simulates the daily growth and seed production of annual plants dependent on water availability and competition (Fig. 8) with different standard parameters for each region (Additional file 3: ParametersPlants.pdf). Growth of germinated seeds starts with seed mass and increases logistically with relative growth rate on day d(RGRd) reduced by a function of temperature (f(Td)), relative resource availability (RAd), and competitiveness (C) as
with mv,d: vegetative aboveground mass, mt: total aboveground mass, RGRmax: maximum relative growth rate and mmax: maximum attainable total mass per individual. RGR is halved for every 10 K when the temperature drops below 25°C, reflecting the effect of temperature on biochemical growth processes. Daily relative resource availability RAd equals water availability (WA) in this study because I consider only patches between shrubs with full sun and assume that water is the only limiting resource in this situation. Effective rooting depth of herbs is assumed restricted to the A layer [27, 59, 76]. Water availability is calculated as a sigmoid function of soil water potential ranging from 0% at the species-specific permanent wilting point (PWP = -1.5 MPa for mesic region plant species, -6.5 for arid region species, based on ) to 100% at field capacity (-0.036 MPa) as WA = 1/(1 + exp(r/PWP·(Ψ - PWP) + s)), with fitted parameters r = 12.5 and s = 5.3 (Fig. 9). I assume that the highest water potential of the four A layers controls the growth rate and that access to deeper soil layers increases with plant mass. These categories were estimated and adjusted to fit final community biomass in 2002/2003.
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.).
Noy-Meir I: Desert Ecosystems: environment and producers. Ann Rev Ecol Syst. 1973, 4: 25-51. 10.1146/annurev.es.04.110173.000325.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Weltzin JF, McPherson GR: Implications of precipitation redistribution for shifts in temperate savanna. Ecology. 2000, 81: 1902-1913. 10.1890/0012-9658.
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.
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.
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.
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.
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.
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.
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.
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.
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-
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.
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.
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.
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.
De Lillis M, Manes F, Tufaro F: Root plasticity in Mediterranean herbaceous species. Plant Biosystems. 2005, 139: 214-221. 10.1080/11263500500160542.
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.
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.
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.
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.
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.
Larcher W: Ökologie der Pflanzen auf physiologischer Grundlage. (UTB 232). 1984, Stuttgart: Ulmer
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.
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.
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.
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.
Nippert JB, Knapp A: Linking water uptake with rooting patterns in grassland species. Oecologia. 2007, 153: 261-272. 10.1007/s00442-007-0745-8.
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-
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.
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.
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.
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.
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.
Bergelson J, Perry R: Interspecific competition between seeds: relative planting date and density affect seedling emergence. Ecology. 1989, 70: 1639-1644. 10.2307/1938097.
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.
Inouye RS: Density-dependent germination response by seeds of desert annuals. Oecologia. 1980, 46: 235-238. 10.1007/BF00540131.
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.
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.
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, August
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.
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.
Shmida A, Wilson M: Biological determinants of species diversity. J Biogeogr. 1985, 12: 1-20. 10.2307/2845026.
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.
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.
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.
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.
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.
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-August
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.
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, August
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.
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.
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.
Mathaj M: Modellierung von Vegetation und Erosion entlang eines Klimagradienten von mediterran bis semiarid. Diplomarbeit (Msc thesis). 2007, Potsdam (Germany): Universität Potsdam
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.
van Genuchten MT: A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980, 44: 892-898.
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.
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.
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.
Lauenroth WK, Sala OE: Long-term forage production of North American shortgrass steppe. Ecol Appl. 1992, 2: 397-403. 10.2307/1941874.
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]
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.
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]
Baier W, Robertson GW: A new versatile soil moisture budget. Can J Plant Sci. 1966, 46: 299-315.
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.
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. 2002
Maidment DR, (ed): Handbook of Hydrology. 1993, New York, U.S.A.: MacGraw Hill
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.
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.
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.
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.
Kunin WE: Density-dependent foraging in the harvester ant Messor ebeninus – 2 experiments. Oecologia. 1994, 98: 328-335. 10.1007/BF00324221.
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.
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.
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.
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.
Weiner J: A neighbourhood model of annual-plant interference. Ecology. 1982, 63: 1237-1241. 10.2307/1938849.
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.
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.
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.
MK designed and wrote the program, planned and evaluated the simulations and scenarios, and wrote the manuscript.
Electronic supplementary material
Additional file 1: Standard parameters used in the soil module. Standard parameters used in the soil module for all simulations. (PDF 45 KB)
Additional file 2: Standard parameters used in the seed bank module. Standard parameters used in the seed bank module for all simulations. (PDF 61 KB)
Additional file 3: Standard parameters used in the plant module. Standard parameters used in the plant module for all simulations. (PDF 84 KB)
Additional file 4: Absolute changes of variables. Alternative view of data shown in Fig. 2 based on relative changes of mean annual precipitation. (PDF 104 KB)
Additional file 5: Precipitation distribution at climate stations. Precipitation distribution at climate stations close to the field sites. (PDF 40 KB)
Additional file 6: Model source code. The archive contains directories with the source code in the c++ language distributed over several text files. The archive also includes a sample input file for running the model in batch mode. (ZIP 89 KB)
Additional file 7: Validation. The file describes how the soil, seed bank, and plant modules of the model were validated. The file contains text and three figures. (PDF 830 KB)
Additional file 8: Sensitivity analyses. The file describes the method and results of testing the sensitivity of runoff, soil moisture, and peak shoot mass to changes in parameters. The file contains text and two figures. (PDF 476 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Köchy, M. Effects of simulated daily precipitation patterns on annual plant populations depend on life stage and climatic region. BMC Ecol 8, 4 (2008). https://doi.org/10.1186/1472-6785-8-4
- Seed Bank
- Arid Region
- Mean Annual Precipitation
- Seedling Density
- Rain Volume