1. Introduction
Soil moisture plays an integral role within the hydrologic cycle as a critical link between soils, climate, and biogeography (Legates et al. 2011). Soil moisture has been shown to influence soil respiration (Geng et al. 2012), act as a thermal reservoir that impacts cloud formation and wind fields (Ek and Holtslag 2004; Findell and Eltahir 2003; Entekhabi et al. 1996), and directly influence precipitation formation (Koster et al. 2004). As the understanding of the importance of soil moisture within an ecosystem increases, the effect precipitation patterns have on the ability of rainfall to infiltrate into the soil becomes an important area of study. If the relationship of intra-annual precipitation events on soil moisture can be better understood, models of precipitation and soil moisture can be improved.
Climates across the world are undergoing unprecedented changes. While global precipitation has been reported to exceed the 1961–90 average every year since 1995 (Dore 2005), climatic variability at the local scale is shifting toward intra-annual patterns of extreme weather such as longer growing seasons, larger temperature ranges, increased storm intensity, and longer dry periods (Frich et al. 2002; Zhang et al. 2013; Easterling et al. 2000). In the early twenty-first century, grassland regions of the United States have experienced a prolonged warm drought (MacDonald 2010) and a shift to larger, more infrequent storms (Moran et al. 2014). This has raised concern because grasslands have the capacity to respond, through shifts in water use efficiency and biomass production, to variability in precipitation patterns. In fact, there is evidence that grasslands may be the most responsive biome to future climate changes (Ponce-Campos et al. 2013).
Recent studies of the impact of intra-annual rainfall variability on grassland productivity have been based almost exclusively on simulation of increased event size and longer within-season drought periods. Results have shown that a shift toward more intense storms will have a significant impact on grassland soil moisture and, consequently, on aboveground net primary production (ANPP). In semiarid regions, manipulated studies have shown that a shift to fewer but larger events, with no change in total rainfall, led to greater soil water content at depths of 20 cm (Heisler-White et al. 2008). In mesic regions, longer dry intervals between events led to below-average soil water content (Heisler-White et al. 2009; Knapp et al. 2002). These changes in soil moisture were sufficient to increase ANPP in arid grassland by 30% and decrease ANPP in mesic grassland by 18%. Fay et al. (2008) found that the response of grassland ecosystems to extreme rainfall patterns was dependent upon the combination of event size, interval between rainfall events, and total rainfall. It is expected that the sensitivity of grasslands to extreme intra-annual precipitation patterns will be greater during periods of prolonged drought, and particularly in arid regions (Knapp et al. 2002). Cherwin and Knapp (2012) used rainfall manipulations in semiarid grasslands to induce extreme drought and reduce soil moisture at the 20-cm depth; however, they reported that ANPP was not reduced when intra-annual rainfall patterns were characterized by large rain events.
Though these manipulated experiments provide good insights into the mechanisms behind the relation between extreme precipitation patterns and soil moisture, ecosystem responses to precipitation changes are not necessarily a simple combination of the responses of the individual factors (Shaw et al. 2002; Zhang et al. 2013; Knapp et al. 2008). Conditions that affect soil moisture and increase production through 2 years can cause simplifications in the food web after 5 years (Suttle et al. 2007). Chronic intense storms will alter both the mean and temporal variability of soil moisture, resulting in long-term shifts in community composition (Knapp et al. 2002). Even with extended drought, plants show some potential to increase their water use efficiency to buffer the short-term effects of climate change (Ponce-Campos et al. 2013; Knapp and Smith 2001). Long-term field studies provide the opportunity to study how ecosystems are affected by climatic changes.
There remains some controversy about whether the effects of the storm size or frequency will have greater influence on annual and perennial vegetation in grasslands. In a rainfall simulation experiment in an arid region of Iran, Jankju (2008) reported that the effects of storm size were more significant than storm frequency on plant community composition. In contrast, field studies of grasslands in China have shown that precipitation frequency and the timing of dry days in relation to the event intensity had a significant influence on the ability of soils to store water at depths of 20–40 cm (He et al. 2012; Wu et al. 2012). Wu et al. (2012) concluded in their study of worldwide sites that precipitation frequency is an overlooked or underestimated parameter in the relationship between precipitation and soil moisture because of lack of available data. Cherwin and Knapp (2012) reported a strong inverse relation between sensitivity to drought and event size, where larger events resulted in less sensitivity. They explained that when rainfall events were sufficiently large, ANPP was uncoupled from total precipitation amount.
The objective of this research was to study the impact of varying storm intensity and prolonged dry periods on grassland soil moisture during the early twenty-first-century drought at sites across the southern United States. We hypothesized that there exists a fundamental difference in the soil moisture response to extreme precipitation patterns between mesic and semiarid precipitation regimes. This study builds on previously published work, yet it is distinctive in several ways. It is one of only a few studies to include consistent observations over a range of precipitation from semiarid to mesic regimes. It is also one of the first studies to address this objective with long-term observations, in this case extending from 1999 to 2013. Instead of studying the soil moisture at depths 20–30 cm, this study is based on measurements at 5-cm depth for two reasons. First, the bulk of grassland root biomass is reported to be in the top 0–10 cm of the soil, that is, ~50% in mesic grasslands and up to 90% in semiarid systems (Knapp et al. 2002; Cox et al. 1986). Second, the currently orbiting Soil Moisture Ocean Salinity (SMOS) and planned Soil Moisture Active Passive (SMAP) sensors will provide global measurements of soil moisture at this depth (Kerr et al. 2001; Entekhabi et al. 2010).
2. Methods
a. Study sites and data selection
Nine sites were selected across the southern United States (Fig. 1), composed of seven Natural Resources Conservation Service (NRCS) Soil Climate Analysis Network (SCAN) stations, one U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) site, and one DOE AmeriFlux site (Baldocchi et al. 2004). The data were obtained from the International Soil Moisture Network (ISMN), and within this network, all data were subjected to automated quality control to ensure reliability of the dataset (Dorigo et al. 2013). Locations of sites ranged from a small grassland field approximately 150 m × 200 m to a site surrounded by everal square kilometers of grassland or scrubland (Table 1). All locations were classified as either mesic or semiarid according to their mean annual precipitation
Distribution of study sites across the southern United States. Sites with lettering in white are mesic and sites with lettering in black are semiarid.
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
Site locations and descriptions. The EVI site column designates sites where MODIS EVI data were obtained and used for analysis.
The soil moisture and climate datasets were obtained from in situ weather stations that measure daily precipitation and soil moisture at 5-cm depth. At SCAN sites, soil moisture was also measured at 20-cm depth. The soil moisture was measured with different sensors depending on the network being used. SCAN sites utilized an analog or digital HydraProbe sensor, the ARM site utilized a Campbell Scientific Model 229L Matric Potential Sensor, and the AmeriFlux site utilized the Delta-T ThetaProbe ML2. The SCAN HydraProbe has a reported accuracy of ±0.013 m3 m−3 (Seyfried et al. 2005). The ARM Matric Potential Sensor and AmeriFlux ThetaProbe both claim an accuracy of ±0.05 m3 m−3 (Schneider et al. 2003; Delta-T Devices 2014).
Daily precipitation data were compiled from rain gauges at the same location as the soil moisture probes. Rain gauges for SCAN sites consisted of tipping buckets or storage-type rain gauges. SCAN precipitation data were given in cumulative amounts and daily rainfall data were obtained by subtracting the annual cumulative total from the previous day’s annual cumulative total. Daily precipitation data from the ARM and AmeriFlux networks were collected using a tipping-bucket rain gauge.
Even with a quality-controlled dataset, oftentimes sites had periods of missing or unreported data. The threshold for discarding a year of data because of unreported soil moisture or precipitation was set at 31 continuous days. The 31-continuous-day threshold allowed less than 10% of the yearly missing data to be interpolated to create a complete year of representational data. Soil moisture data were interpolated by averaging the soil moisture values from reported values on dates before and after data were missing. Interpolation to fill gaps in data for soil moisture needed to be done in approximately 155 instances; approximately 143 of those times were for 1- or 2-day periods. For the other 12 times, soil moisture data were interpolated for time periods between 3 and 16 days.
If precipitation data were missing for 1 day and there was a change in the accumulated precipitation value from before and after the unreported data, we recorded that the difference in the accumulated values was the amount of precipitation that fell on the unreported date. For the approximately 11 times precipitation data were missing between 2 and 31 days, precipitation data gaps were filled using a two-step process. First, we made a calculation of the precipitation that fell during the period that data that were missing. This was accomplished by subtracting the annual cumulative precipitation value immediately after reliable data were reported from the annual cumulative precipitation immediately before the station began to report errors. If the precipitation that fell during this time period was >1 mm, the total amount of precipitation that fell during the missing gap was then temporally distributed in levels as storms and dry periods based on precipitation records from weather stations within 30 km operated by the National Oceanic and Atmospheric Administration and accessed through the National Climatic Data Center.
b. Satellite data
To study the ecosystem effects of soil moisture on biomass production, the NASA Moderate Resolution Imaging Spectroradiometer (MODIS) enhanced vegetation index (EVI) was used as a cross-site indicator of plant activity (Huete et al. 2002). The dataset for this study was taken from the MODIS product subset MOD13Q1 over the time period from 2001 to 2013. A footprint of 2.25 km × 2.25 km (9 × 9 pixels) was obtained over a homogenous area, where the EVI data could then be averaged to produce a single value representative of the in situ measurements (Table 1). If a site was not surrounded by a homogeneous area large enough for the MODIS footprint, EVI data were not retrieved. MODIS data utilize a pixel-based quality assurance scheme to limit noise from aerosols or other atmospheric interference within the dataset. Software developed by Jönsson and Eklundh (2004) was then used to further smooth the dataset before integrating between the start and end of the growing season to produce a proxy for the ANPP known as the integrated EVI iEVI. The iEVI value has been shown to be a reasonable approximation of ANPP (Zhang et al. 2013; Ponce-Campos et al. 2013).
The SMOS satellite provides global soil moisture measurements once every 3 days with an expected accuracy of 0.04 m3 m−3. One pixel of SMOS data covers approximately 50 km × 50 km on the ground, an area far greater than what is obtained from the in situ soil moisture probes (Barre et al. 2008). SMOS data were obtained in the form of the level 3 (L3) product from the Centre Aval de Traitement des Données SMOS (CATDS), operated for the Centre National d’Études Spatiales (CNES; France) by the Institut français de recherche pour l’exploitation de la mer (IFREMER; Brest, France). The level-3 soil moisture (m3 m−3) global data are a 1-day product and contain filtered data gridded at 25-km spatial resolution. The best estimation of soil moisture and dielectric constant are selected (based on the minimization of the data quality index) for each node when several multiorbit retrievals are available for a given soil moisture User Data Product (UDP). Once data were acquired using the in situ locations, soil moisture values were extracted from the ascending measurements for the 2011, 2012, and 2013 hydrologic years. There were 10 site years where both SMOS and in situ probe measurements were available at each site over the same timeframe that allowed for direct comparison.
c. Data analysis











Summary of high and low SDII and CDD classifications.






3. Results and discussion
a. Extent and depth of annual average soil moisture values
A basic premise of this study was that the results obtained with soil moisture measurements at 5-cm depth would compare well with previously reported results based on measurements at 20-cm depth. At the SCAN sites where mean annual soil moisture at 5- and 20-cm depths were available over the same annual period (n = 40), the soil moisture at the two depths was strongly correlated (coefficient of determination R2 = 0.96, p < 0.001; Fig. 2). The high correlation between mean annual soil moisture at 5- and 20-cm depths provided confidence in the ability of the 5-cm soil moisture to be representative of soil moisture at lower depths at these grassland sites at the annual time scale.
The linear relationship between mean annual soil moisture observed at 5- and 20-cm depths over the dataset when both observations were available (n = 40).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
In addition, our decision to use soil moisture measured at 5-cm depth was based partly on the fact that the currently orbiting SMOS, launched in 2009, and planned SMAP sensors will provide global measurements of soil moisture at this depth. For the SMOS and in situ measurements where both datasets were of good quality and available over the same timeframe (n = 10), we found a reasonable relationship (R2 = 0.48, p = 0.027) between
Linear regression of mean annual soil moisture (cm3 cm−3) measurements at 5-cm depth obtained from the SMOS satellite and from in situ probes within the SCAN network when both SMOS and in situ measurements were available (n = 10, p < 0.05).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
b. Cross-site 
–
and 
–iEVI spatial relationships



Across sites where iEVI was available, we found that
The relationship between iEVI and (a) soil moisture at 5-cm depth and (b) precipitation. Data are given for sites where iEVI was available (n = 30).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
The general within-site temporal relationship between soil moisture and precipitation was studied by using the normalized values of soil moisture and precipitation. For both mesic and semiarid sites, there were weak but significant correlations between
Linear regression of normalized mean annual soil moisture and normalized annual precipitation for (a) mesic and (b) semiarid sites, where each data point represents a site year (n = 32 for each precipitation regime).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
c. SDII influence on linear 
–
relationship


For mesic sites,
Linear regression of normalized mean annual soil moisture and normalized annual precipitation for (a) mesic (n = 32) and (b) semiarid sites (n = 32), where sites were classified into high and low SDII as described in Table 2. (c) The distribution of the difference between field capacity and mean annual soil moisture
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
Linear regression R2 with data combined or split by SDII and CDD. Boldface values indicate significance at p < 0.05.
Linear regression variable coefficient (VC) estimates for Eqs. (8)–(10) that were identified in section 2c. The R2 is reported in the bottom row for the equation of each column. The p values for each variable are presented in parentheses after the estimated coefficient, and boldface values indicate significance at p < 0.05.
For semiarid sites, results were inconsistent. For years with high SDII, the correlation between
Storm intensity is an important factor in the correlation between
As mentioned previously, SDII is a simplification of storm intensity based on precipitation data available at the daily time scale. The suitability of SDII as a proxy for storm intensity is arguable because of a potential greater influence from shifts between summer and winter precipitation events than actual increases in annual average storm intensity. We found that the correlation between
d. CDD influence on linear 
–
relationship


For mesic sites, the relation between
Linear regression of normalized mean annual soil moisture and normalized annual precipitation for (a) mesic and (b) semiarid sites, where sites were classified into high and low CDD, as described in Table 2 (n = 32 for each precipitation regime).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
For semiarid sites, classification of years with high or low CDD improved correlation between
e. SDII and CDD effect on 

For mesic sites, average
Boxplots of the normalized mean annual soil moisture for mesic and semiarid sites split by length of (a) SDII and (b) CDD, as described in Table 2. Boxplots broken up by precipitation regime are shown for (c) SDII and (d) CDD. The boldface line represents the median, the upper and lower limits of the box represent the first and third quartiles, and the dashed lines extend to the max and min values (n = 32 for each precipitation regime within each boxplot).
Citation: Journal of Hydrometeorology 16, 1; 10.1175/JHM-D-14-0057.1
Mean normalized annual soil moisture values broken down by extreme index classification. Boldface values indicate there was a significant difference (p < 0.05) between the mean normalized annual soil moisture of high and low CDD or SDII.
The unique influence of CDD on semiarid sites and the inability of SDII to influence soil moisture without the incorporation of precipitation may be affected by the difference in the length and variation of CDD or SDII that each regime naturally experiences. Even though SDII values between mesic and semiarid sites were significantly different (P < 0.001; Fig. 8c), the SDII influence on soil moisture at mesic sites is codependent with precipitation. On the other hand, semiarid soil moisture was influenced significantly between high and low CDD. This may be a result of semiarid sites experiencing an average CDD of 75 days with a standard deviation of 55.5 days. These values are greater than the mesic sites’ average CDD of 28 days, with a standard deviation of 12.1 days. Even the mesic site maximum CDD value of 62 days would place that year with a semiarid low CDD classification (Fig. 8d). With this difference of CDD values that mesic and semiarid sites experience, mesic sites may not be influenced by CDD because they experience a relatively low CDD and a relatively small variation of CDD.
4. Conclusions
For nine grassland sites across the southern United States, we investigated the differing impact of intra-annual patterns of precipitation on soil moisture in mesic and semiarid precipitation regimes. Soil moisture was strongly related to ANPP, as measured by the iEVI. In fact, we found that soil moisture was a better predictor of iEVI than was total precipitation. We also found that our in situ point measurements of soil moisture were correlated to large-scale remotely sensed soil moisture measurements that will allow for increases in the scale of similar studies.
In this study, we addressed the within-site temporal relation between precipitation and soil moisture by normalizing soil moisture and precipitation to units of standard deviation and analyzing the results for the sets of semiarid and mesic sites. We concluded that the precipitation–soil moisture relation was impacted by variations in storm intensity for mesic sites and by the length of consecutive dry days for semiarid sites. This was explained in relation to general climate patterns and the associated distinctions in soil texture in these two precipitation regimes. First, the soil moisture at mesic sites tends to reside closer to field capacity than soil moisture at semiarid sites. So, precipitation associated with high or low storm intensity events will infiltrate into the drier soils at semiarid sites, whereas for mesic sites, precipitation associated with high-intensity storms will result in greater water loss to runoff and less infiltration into soil, compared to low-intensity storms. Second, this field study attributed the length of CDD as having a significant impact on soil moisture only at semiarid sites. This is related to the fact that the variation in length of CDD was naturally low at mesic sites and not variable from year to year, in contrast to the high variability of CDD at semiarid sites. For semiarid sites, long periods of CDD decreased the mean annual soil moisture regardless of the total precipitation throughout the year.
While many factors influence how soil moisture varies with precipitation, this research offers a consistent hydrometeorological explanation for how mesic and semiarid grasslands will be affected by changing precipitation patterns and a changing climate. Incorporation of intra-annual patterns of storm intensity in mesic grasslands and consecutive dry days in semiarid grasslands will help improve estimates of annual soil moisture and better predict grassland production.
Acknowledgments
The research was funded and supported by the University of Arizona/Arizona NASA Space Grant Consortium, the NASA SMAP Science Definition Team under Agreement 08-SMAPSDT08-0042, and the University of Arizona Graduate College Latin America Summer Research Program. The SMAP Early Adopter Program provided valuable technical support. We acknowledge the U.S. Department of Energy ARM Climate Research Facility and the USDA SCAN network for making these data available. We would also like to thank Dr. Joe Watkins, University of Arizona Department of Mathematics, for his help in the statistical analysis of our results.
REFERENCES
Baldocchi, D. D., Xu L. , and Kiang N. , 2004: How plant functional-type, weather, seasonal drought, and soil physical properties alter water and energy fluxes of an oak–grass savanna and an annual grassland. Agric. For. Meteor., 123, 13–39, doi:10.1016/j.agrformet.2003.11.006.
Barre, H. M. J., Duesmann B. , and Kerr Y. H. , 2008: SMOS: The mission and the system. IEEE Trans. Geosci. Remote Sens.,46, 587–593, doi:10.1109/TGRS.2008.916264.
Cherwin, K., and Knapp A. , 2012: Unexpected patterns of sensitivity to drought in three semi-arid grasslands. Oecologia, 169, 845–852, doi:10.1007/s00442-011-2235-2.
Cox, J. R., Frasier G. W. , and Renard K. G. , 1986: Biomass distribution at grassland and shrubland sites. Rangelands, 8, 67–68.
Delta-T Devices, cited 2014: ML3 ThetaProbe research data sheet. [Available online at www.delta-t.co.uk/product-support-material.asp.]
Dore, M. H. I., 2005: Climate change and changes in global precipitation patterns: What do we know? Environ. Int., 31, 1167–1181, doi:10.1016/j.envint.2005.03.004.
Dorigo, W. A., and Coauthors, 2013: Global automated quality control of in situ soil moisture data from the International Soil Moisture Network. Vadose Zone J.,12, doi:10.2136/vzj2012.0097.
Easterling, D. R., Meehl G. A. , Parmesan C. , Changnon S. A. , Karl T. R. , and Mearns L. O. , 2000: Climate extremes: Observations, modeling, and impacts. Science, 289, 2068–2074, doi:10.1126/science.289.5487.2068.
Ek, M. B., and Holtslag A. A. M. , 2004: Influence of soil moisture on boundary layer cloud development. J. Hydrometeor., 5, 86–99, doi:10.1175/1525-7541(2004)005<0086:IOSMOB>2.0.CO;2.
Entekhabi, D., Rodriguez-Iturbe I. , and Castelli F. , 1996: Mutual interaction of soil moisture state and atmospheric processes. J. Hydrol., 184, 3–17, doi:10.1016/0022-1694(95)02965-6.
Entekhabi, D., and Coauthors, 2010: The Soil Moisture Active Passive (SMAP) Mission. Proc. IEEE, 98, 704–716, doi:10.1109/JPROC.2010.2043918.
Fay, P. A., Kaufman D. M. , Nippert J. B. , Carlisle J. D. , and Harper C. W. , 2008: Changes in grassland ecosystem function due to extreme rainfall events: Implications for responses to climate change. Global Change Biol., 14, 1600–1608, doi:10.1111/j.1365-2486.2008.01605.x.
Findell, K. L., and Eltahir E. A. B. , 2003: Atmospheric controls on soil moisture-boundary layer interactions: Three-dimensional wind effects. J. Geophys. Res.,108, 8385, doi:10.1029/2001JD001515.
Frich, P., Alexander L. V. , DellaMarta P. , Gleason B. , Haylock M. , Tank A. M. G. K. , and Peterson T. , 2002: Observed coherent changes in climatic extremes during the second half of the twentieth century. Climate Res., 19, 193–212, doi:10.3354/cr019193.
Geng, Y., and Coauthors, 2012: Soil respiration in Tibetan Alpine grasslands: Belowground biomass and soil moisture, but not soil temperature, best explain the large-scale patterns. PLoS One, 7, e34968, doi:10.1371/journal.pone.0034968.
He, Z., Zhao W. , Liu H. , and Chang X. , 2012: The response of soil moisture to rainfall event size in subalpine grassland and meadows in a semi-arid mountain range: A case study in northwestern China’s Qilian Mountains. J. Hydrol., 420–421, 183–190, doi:10.1016/j.jhydrol.2011.11.056.
Heisler-White, J. L., Knapp A. K. , and Kelly E. F. , 2008: Increasing precipitation event size increases aboveground net primary productivity in a semi-arid grassland. Oecologia, 158, 129–140, doi:10.1007/s00442-008-1116-9.
Heisler-White, J. L., Blair J. M. , Kelly E. F. , Harmoney K. , and Knapp A. K. , 2009: Contingent productivity responses to more extreme rainfall regimes across a grassland biome. Global Change Biol., 15, 2894–2904, doi:10.1111/j.1365-2486.2009.01961.x.
Huete, A., Didan K. , Miura T. , Rodriguez E. P. , Gao X. , and Ferreira L. G. , 2002: Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ., 83, 195–213, doi:10.1016/S0034-4257(02)00096-2.
Jankju, M., 2008: Individual performances and the interaction between arid land plants affected by the growth season water pulses. Arid Land Res. Manage., 22, 123–133, doi:10.1080/15324980801957986.
Jönsson, P., and Eklundh L. , 2004: TIMESAT—A program for analyzing time-series of satellite sensor data. Comput. Geosci., 30, 833–845, doi:10.1016/j.cageo.2004.05.006.
Kerr, Y. H., Waldteufel P. , Wigneron J.-P. , Martinuzzi J. , Font J. , and Berger M. , 2001: Soil moisture retrieval from space: The Soil Moisture and Ocean Salinity (SMOS) mission. IEEE Trans. Geosci. Remote Sens., 39, 1729–1735, doi:10.1109/36.942551.
Knapp, A. K., and Smith M. D. , 2001: Variation among biomes in temporal dynamics of aboveground primary production. Science, 291, 481–484, doi:10.1126/science.291.5503.481.
Knapp, A. K., and Coauthors, 2002: Rainfall variability, carbon cycling, and plant species diversity in a mesic grassland. Science, 298, 2202–2205, doi:10.1126/science.1076347.
Knapp, A. K., and Coauthors, 2008: Consequences of more extreme precipitation regimes for terrestrial ecosystems. Bioscience, 58, 811–821, doi:10.1641/B580908.
Koster, R. D., and Coauthors, 2004: Regions of strong coupling between soil moisture and precipitation. Science, 305, 1138–1140, doi:10.1126/science.1100217.
Legates, D. R., Mahmood R. , Levia D. F. , DeLiberty T. L. , Quiring S. M. , Houser C. , and Nelson F. E. , 2011: Soil moisture: A central and unifying theme in physical geography. Prog. Phys. Geogr., 35, 65–86, doi:10.1177/0309133310386514.
MacDonald, G. M., 2010: Water, climate change, and sustainability in the southwest. Proc. Natl. Acad. Sci. USA, 107, 21 256–21 262, doi:10.1073/pnas.0909651107.
Moran, M. S., and Coauthors, 2014: Change in the function of U.S. grasslands in response to the early 21st century drought. Ecology, 95, 2121–2133, doi:10.1890/13-1687.1.
NRCS, cited 2014: Web Soil Survey. [Available online at http://websoilsurvey.sc.egov.usda.gov/App/HomePage.htm.]
Ponce-Campos, G. E., and Coauthors, 2013: Ecosystem resilience despite large-scale altered hydroclimatic conditions. Nature, 494, 349–352, doi:10.1038/nature11836.
Schneider, J. M., Fisher D. K. , Elliott R. L. , Brown G. O. , and Bahrmann C. P. , 2003: Spatiotemporal variations in soil water: First results from the ARM SGP CART network. J. Hydrometeor., 4, 106–120, doi:10.1175/1525-7541(2003)004<0106:SVISWF>2.0.CO;2.
Seyfried, M. S., Grant L. E. , Du E. , and Humes K. , 2005: Dielectric loss and calibration of the hydra probe soil water sensor. Vadose Zone J., 4, 1070, doi:10.2136/vzj2004.0148.
Shaw, M. R., Zavaleta E. S. , Chiariello N. R. , Cleland E. E. , Mooney H. A. , and Field C. B. , 2002: Grassland responses to global environmental changes suppressed by elevated CO2. Science, 298, 1987–1990, doi:10.1126/science.1075312.
Suttle, K. B., Thomsen M. A. , and Power M. E. , 2007: Species interactions reverse grassland responses to changing climate. Science, 315, 640–642, doi:10.1126/science.1136401.
Wu, C., and Coauthors, 2012: An underestimated role of precipitation frequency in regulating summer soil moisture. Environ. Res. Lett., 7, 024011, doi:10.1088/1748-9326/7/2/024011.
Zhang, Y., and Coauthors, 2013: Extreme precipitation patterns and reductions of terrestrial ecosystem production across biomes. J. Geophys. Res. Biogeosci., 118, 148–157, doi:10.1029/2012JG002136.