National assessment of throughfall sensitivity to changes in storm magnitude for the forests of Iran

Aim of study: To understand throughfall (TF) sensitivity to variability in rainfall amount (Pg) for typical forest sites across the main climate types of Iran. Area of study: Nine forest stands of several common native and introduced tree species situated in all common Iranian climate types, but located primarily in northern Iran. Material and methods: A nondimensional relative sensitivity coefficient was employed to predict responses of TF to Pg changes. Projected Pg changes over the measurement sites for the period 2020-50 were estimated using one of the Coupled Model Intercomparison Project phase 5 (CMIP5) known as HadGEM2-ES under low and high emission scenarios (RCP 2.6 and 8.5). Main results: TF displayed strong positive linear relationships with Pg at all sites [TF=0.66 Pg -0.16; R 2=0.91]. The sensitivity coefficient ranged from 0.96-2.35 across the nine forest sites and large sensitivity coefficient differences were found between small (< mean annual Pg) and large (> mean annual Pg) storms for arid and Mediterranean plantations. Shifts in Pg and increased small storm frequency are predicted for these regions (2020-50) under low and high emission scenarios. Research highlights: TF sensitivity may be a useful variable when selecting tree species for afforestation to buffer expected shifts in Pg due to climate change. Additional keywords: climate change; forest ecosystems; precipitation projection; throughfall sensitivity. Authors ́ contributions: Pedram Attarod: Design and conception; Qiuhong Tang: Financial support; John T. Van Stan II: Language editor and scientific comments, and Xingcai Liu: Data preparation. Citation: Attarod, P., Tang, Q., Van Stan II, J-T., Liu, X. (2018). National assessment of throughfall sensitivity to changes in storm magnitude for the forests of Iran. Forest Systems, Volume 27, Issue 3, e019. https://doi.org/10.5424/fs/2018273-13857 Received: 27 Aug 2018. Accepted: 03 Dec 2018. Copyright © 2018 INIA. This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC-by 4.0) License. Funding: This research was supported by the University of Tehran (UT), Iran and by the International Partnership Program of Chinese Academy of Sciences (Grant No. 131A11KYSB20170113) and the National Natural Science Foundation of China (Grant No. 41790424). Competing interests: The authors have declared that no competing interests exist. Correspondence should be addressed to Pedram Attarod: attarod@ut.ac.ir Forest Systems 27 (3), e019, 12 pages (2018) eISSN: 2171-9845 https://doi.org/10.5424/fs/2018273-13857 Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria O. A., M. P. (INIA)


Introduction
A significant portion of rainfall is intercepted on forest canopy surfaces, where it is either returned to the atmosphere as interception or redistributed to the ground as throughfall (TF) and stemflow.This rainfall partitioning by forest cover significantly alters hydrologic cycling along the soil-forest-atmosphere continuum (Miralles et al., 2010) and can impact feedbacks between the hydrologic cycle and global climate (Davies-Barnard et al., 2014).Stemflow is the proportion of rainfall that drains to the ground along the stem, usually accounting for <2% of annual rainfall in most forests (Van Stan & Gordon, 2018).The remaining majority (60-90%) of rainfall passes through canopy gaps or drips from the vegetation, reaching the surface as TF (Levia et al., 2011) where it can influence soil moisture (Raat et al., 2002), physicochemistry (Rosier et al., 2015), fine root distribution (Ford & Deans, 1978), and microbial processes (Moore et al., 2016).As these ecohydrological interactions play key roles in ecosystem functioning, understanding TF's response to natural and anthropogenic variability is critical to forest management.
Studies across forest ecosystems agree that rainfall amount (P g ) is the principal variable driving standscale TF amount (Levia & Frost 2006;Levia et al., 2011)-often explaining >90% of inter-storm TF variability-and that TF-P g relationships are shaped by canopy structures, like leaf area index, crown depth, etc. (Staelens et al., 2008;Toba & Ohta 2005).However, the authors are unaware of previous work quantifying TF sensitivity to changes in P g -a commonly anticipated climate response to hydrologic intensification (Huntington, 2006).This is surprising since hydrologic intensification has been linked to increased extreme P g in many regions (c.f., Tollefson, 2016).Moreover, afforestation, reforestation, urban forestry and other "greening" initiatives have grown in popularity and, therefore, increased the forest cover in most developed regions (i.e., McGovern & Pasher, 2016).In Iran, for example, the restoration of semi-arid and arid ecosystems through planting of low-demand and drought tolerant species has become a critical element of national ecosystem management plans (Attarod et al., 2015b).
Iran has invested in vast tree-plantings throughout its major cities for urban greening and air pollution mitigation, including the Chitgar and Lavizan Forest Parks (Sadeghi et al., 2016) and the ongoing "Jam Afforestation Project" which is tasked with largescale afforestation and reforestation in the Zagros region (FRWO, 2012).Restoration of the natural Caspian deciduous forests (that extend from the Alborz Mountains to the southern coast of the Caspian Sea) has also resulted in significant reforestation projects since the 1960s (Abbasian et al., 2015).Concerns have risen over the impact of these greening initiatives on the hydrological cycle (Sun et al., 2006;Wang et al., 2011), particularly for the 90% of Iran classified as arid or semi-arid (Ul Hassan et al., 2007).The first process in the rainfall-to-runoff pathway is the partitioning of P g by forest canopies (Savenije, 2004) and, accordingly, an improved understanding of TF sensitivity to climate change and P g is essential for addressing these concerns and quantifying the impacts of Iran's (and other nations') large-scale afforestation and reforestation efforts.
A sensitivity analysis is a technique used to determine how different values of an independent variable impact a particular dependent variable under a given set of assumptions.A method for estimating the "sensitivity coefficient" of a dependent variable (in this case, TF) on the relative changes of an independent variable (in this case, P g ) exists (McCuen, 1974) and has been rigorously applied to evapotranspiration and its principal meteorological drivers (Hupet & Vanclooster, 2001).However, the authors are unaware of its application to assess sensitivity of TF to its principal meteorological driver, P g , for multiple species of contrasting canopy structure.Thus, the objectives of this study are to (1) collect rainfall and TF across common forest species of differing canopy structures and climates in Iran, then (2) quantify and compare their TF sensitivity coefficients.Accomplishing these aims will provide novel information to complement existing data used by forest managers during species selection for restoration and afforestation activities.

Sites description
Data were collected in 9 forest stands of several common native and introduced tree species (Table 1) situated in all common Iranian climate types (Table 2), but located primarily in northern Iran (Fig. 1).The selected species represent a diversity of forest canopy architectures, ranging from the smooth-barked, broadleaved canopy of Fagus orientalis (FO) to the roughbarked, needle-leaved canopy of Pinus eldarica (PE).Leaf phenology also differs among the selected species as, for example FO is deciduous and PE is evergreen.The PE and Cupressus arizonica (CA) throughfall sites are in southern Alborz Mountain Range near the city of Tehran, while the remaining sites-FO, Quercus castaneifolia (QC 1 -QC 3 ), Acer velutinum (AV), Pinus brutia (PB), Cupressus sempervirens var.horizontalis (CS)-are in northern Alborz Mountain Range and Southern coasts of the Caspian Sea (Fig. 1).Stand structural (and TF) measurements for each forest stand were performed in 0.5 ha plots.Stand density ranged from 112 trees ha -1 in the FO forest to 1,600 trees ha -1 in the CS plantation (Table 1).Diameter at breast height (1.3 m, dbh) varied greatly among the measured forests, with the smallest values in CS (12 cm) and the largest values in QC 3 (65 cm) (Table 1).Mean canopy coverage also exhibited a wide range across sites, from 45% for CS and 95% for FO (Table 1).
Meteorological data were obtained from the nearest synoptic meteorological stations recording reliable long-term meteorological data (Fig. 1).The range in meteorological data records is from 1951 to 2015 (Table 2).There exists relatively long distances between the Mehr-Abad, Sari, and Gorgan meteorological stations (Fig. 1), but there is no significant topography between the measurement sites and the meteorological stations, and no closer meteorological stations exist in the region.Due to lack of meteorological station inside the Caspian forest, we used meteorological data recorded by Nou-Shahr meteorological station regardless of elevation difference.The ranges of annual precipitation (P) and temperature (T) in the weather station sites are 16.4-17.4ºC and 230-1291 mm, respectively (Table 2).The ''De Martonne'' climate classification, i.e., De Martonne aridity index (I DM ), as described by Baltas (2007), ranged roughly from 49 to 8.5, so that the study sites  were grouped into very humid (FO, QC 1 ), semi-humid (QC 2 , AV, PB), Mediterranean (QC 3 , CS), and arid (PE, CA) climates (Table 1).

Field measurements
Field measurements were performed from July-2008 to March-2014 during growing and dormant seasons (Table 1).A discrete rain event was defined as a period with >0.1 mm of rainfall.The minimum inter-event dry period between discrete storms was 4-10 h, depending on the site.The effect of pre-storm canopy wetness was assumed negligible as the canopy is assumed to be dry after the minimum inter-event time, which is a common assumption among rainfall partitioning studies (Carlyle-Moses & Gash, 2011).Snowfall was ignored.P g at each forest site was measured by 3-10 funnel-type plastic collectors with 10 cm funnel diameter and 20-30 cm heights, placed in the nearest open area away from the forest stands.The average of water amount measured in all rain-gauges at a site was used to estimate P g for each site.Quantities of water in the collectors were measured manually using a graduated cylinder.Storms that reached only the weather stations but not the forest sites, or vice versa, were ignored.P g volumes were measured at the same time as TF volumes at each site, either immediately after a storm or at sunrise following a night storm.
TF was measured using 20-50 rain collectors of the same design as those used to quantify P g .TF collectors were randomly distributed beneath the canopy and fabric covered the neck of the collectors to avoid litter, needles, and debris from entering (Abbasian et al., 2015).TF data based on multiple field campaigns that used a different number of collectors is a common, often necessary, issue in national-to-international scale studies (i.e., Wallace et al., 2013;Návar, 2017); however, it is important to note that this methodological variability may introduce uncertainty to the stand-scale TF estimates (Vose et al., 2016).Variability about mean P g and TF will be expressed in standard error (SE) throughout.

Throughfall sensitivity coefficients
A practical method of presenting a sensitivity analysis is to plot relative changes of a dependent variable (in this case, TF) against relative change of an independent variable (i.e., P g ) as a curve (e.g., Singh & Xu, 1997;Goyal, 2004, Attarod et al., 2015a).This is different from a standard regression of TF versus P g (Fig. 2a) as the mean sensitivity coefficient is calculated as the slope of a correlation between the percent changes in P g against percent changes in TF (Fig. 2b).The sensitivity coefficient represents the fraction of change in P g transmitted to the change of TF, i.e. a sensitivity value of 0.1 would suggest that a 10% increase in P g may be predictable to increase TF by 1%.Negative coefficients would indicate that a decrease in TF would result from an increase of P g , which is not expected due to the universally reported positive relationship between P g and TF (Friesen et al., 2015).In the present study, the sensitivity coefficient of TF was determined in climate classifications only in response to changes in P g on an event-basis.

Projected P g changes over the measurement sites
General Circulation Models (GCMs) are numerical models representing physical processes in the atmosphere, ocean, cryosphere, and land surface.These models are the most advanced tools currently available for simulating the response of the global climate system to increasing greenhouse gas concentrations.GCM simulations for the fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC) have become available (Taylor et al., 2012;Miao et al., 2014).Comparing to the IPCC AR4, the GCMs in AR5 include a more varied set of model types (i.e., climate/ earth system models with more interactive components such as atmospheric chemistry, aerosols, dynamic vegetation, ice sheets and carbon cycle) (Liu et al., 2013).Several improvements in the physics, numerical algorithms and configurations are implemented in the IPCC AR5 models with a new set of scenarios called Representative Concentration Pathways (RCPs) for energy and industry CO 2 emissions (Moss et al., 2010).The RCPs span a large range of stabilization, mitigation and non-mitigation pathways.Phase 5 of the Coupled Model Intercomparison Project (CMIP5) is a standard experimental protocol for studying the output of GCMs which provides a community-based infrastructure in support of climate model diagnosis, validation, intercomparison, documentation and data access (Jones et al., 2011).Longer time-scale ("centennial") experiments have been performed at the Met Office Hadley Centre with the HadGEM2-ES Earth System model-one of the CMIP5 climate models (Collins et al., 2011;Jones et al., 2011;Miao et al., 2014).The HadGEM2-ES model used here is one of the state-ofthe-art GCMs and involves many typical and advanced representations of land and ocean processes (Jones et al., 2011).The HadGEM2-ES climate data has been widely used for climate studies (Huntingford et al., 2013).The r1i1p1 ensemble of HadGEM2-ES was used in this study.The r1i1p1 ensemble is the most accessible ensemble in the CMIP5 archive.
To understand P g variations under changing climate over the measurement sites, we focused on the precipitation projection for two scenarios: RCP 2.6 and RCP 8.5.RCP 2.6 represents a "low" emissions scenario featured by the radiative forcings of 2.6 Wm -2 and atmospheric CO 2 concentration of 421 ppm by 2100.RCP 8.5 represents a "high" emission scenario with the radiation forcing of 8.5 Wm -2 and CO 2 concentration of 936 ppm by 2100.The time resolution of projected storm magnitude (i.e., P g ) is daily.Future climate data in grids with 0.5°×0.5°horizontal resolution across measurement sites were obtained from the HadGEM2-ES model projections (Hempel et al., 2013).

Delta change method
Because of the limitations of coarse-resolution GCM climate data (used in these predictions), the delta change (DC) method is generally used to derive scenarios of future climate (e.g., Fischer et al., 2007;Shahid, 2011;Chung & Nkomozepi, 2012).The method consists of simply scaling the observed climate data using monthly change factors calculated from the differences in climatology predicted by GCMs for the current and future periods.In this way, P g for future time periods was derived by scaling the historically observed climate data (OBS) by the GCM-computed change.This results in a new P g time series scaled according to the GCMs, but based on historical observations.In this study, the OBS from 1974 to 2004 (reference period) were scaled to derive the 30-yr future climate scenario for the period 2020-2050.Relative change factors (ΔVar) were then applied to the observed flux variable to calculate future P g .Following Leng & Tang (2014), the DC method are formulated as Eq., (1): (1) where Var Δ is the scaled flux variable using the DC method and Var OBS is the observed flux variable in the historic period.The suffixes i and j stand for the day and the month, respectively, and ΔVar is the monthly DC factor, which is calculated as follows: (2) where Var future (j) and Var current (j) are the mean values of the time series for month j for the future and current time periods by GCMs.

Overview of rainfall and throughfall events
There were 290 P g events recorded from 2008-2014 across all measurement sites and P g ranged from 0.5-54.7 mm.Mean P g was 10.8±0.7 mm across all forest sites, but site-specific mean P g ranged from 17-20 mm for forests in the very humid, semi-humid, and Mediterranean climate types (Table 3).Mean P g from the arid forest sites was 4.4 mm (Table 3).At arid sites, 67% of collected storms were < mean P g ; however, for other climate types, this fraction was generally around 50% (Table 3)."Large" storms (P g > 15 mm) were frequent at the very humid forest sites (85%), yet only 5% of storms exceeded 15 mm at arid sites (Table 3)."Small" storms (P g < 5 mm) were the norm for arid forests in this study (70%: Table 3).The respective percentages of P g < 5 mm or P g > 15 mm were similar for the semi-humid and Mediterranean forests (Table 3).
Strong positive linear relationships between event mean TF and P g were observed across all sites [TF=0.66(P g ) -0.16; R 2 = 0.91] and at each individual site (Table 4).Slopes of the regression lines between TF and P g (often assumed to be linked to varying canopy structures) were wide-ranging: 0.39 for CS to 0.81 for FO (Table 4).Relative TF (TF:P g ) varied roughly from 40% (Mediterranean CS) to 75% (very humid QC 1 ) (Table 4).In arid and Mediterranean climates, minimum TF:P g were measured at zero, however, in other climates the minimum was approximately 47% (Table 4).TF:P g was, surprisingly, nearly the same when averaged for all the needle-leaved forests (55.2%) and broad-leaved forests (50%) across the climate zones (Table 4).Marked differences in TF:P g on average, however, were observed for plantations versus natural forests, where plantations were found to produce 53.4% of TF compared to 71.9% for natural forests.For all sites, QC 1 produced the highest TF:P g in very humid conditions (QC 1 : 74.5%) and the lowest TF:P g under Mediterranean conditions (CS: 39.3%) (Table 4).Generally higher standard errors were observed for TF:P g in the Mediterranean (3.25%) and arid sites (3.35%) compared to forests in other climates (1.34%) (Table 4).

Historical annual storm characteristics
Historical records at the meteorological stations show that mean annual precipitation decreased from very humid (1291 mm y -1 ) to arid climate (230 mm y -1 ) (Table 5).Compared with other climates, the arid climate had relatively constant annual precipitation, varying only by 9 mm y -1 between 1951-2015 (Table 5).The ratio of maximum P g to mean annual precipitation was respectively 11%, 16%, 18%, and 22% for the semi-humid, very humid, Mediterranean, and arid climates (Table 5).Although mean P g had a descending trend from very humid (10.5 mm) to arid (3.9 mm), the proportion of mean P g to maximum P g in different climates ranged from 8.5% in the semi humid climate to 5.0% in the very humid climate (Table 5).

Throughfall sensitivity to changes in storm magnitude
TF exhibited varying degrees of sensitivity to P g , showing large fluctuations (i.e., doubling in sensitivity) Table 4. Relationships between throughfall (TF) and rain event magnitude (P g ), percent of event based average relative throughfall (TF: P g ), and related statistics.See Table 1 for the tree species represented by the location codes.n refers to the number of recorded events.for species planted in different climates (Fig. 3).TF from the arid PE forest was most sensitive to fluctuations in P g (= 2.35) and TF from the Mediterranean QC 3 forest was least sensitive (= 0.96) (Fig. 3).Needle-leaved plantations (i.e., PE, CA) generally had higher sensitivity (> 1.5) to changing P g compared to broadleaved forests (~ 1.0) (Fig. 3).
For broadleaved forests, QC 1 was less sensitive than FO in a very humid climate (Fig. 3).QC 1 , QC 2 , and QC 3 had similar sensitivity coefficients regardless of climate: very humid (=1.09) to semi-humid (=1.07) and Mediterranean climates (= 0.96) (Fig. 3).TF sensitivity was tested for larger and smaller storm sizes than the mean P g in each climate (Table 6).
Excluding AV, PB, and QC 1 , TF was found to be more sensitive to small storms (Table 6).Plantations in arid or Mediterranean climates had large differences in sensitivity coefficient between small and large events, e.g.2.95 vs. 1.14 for PE and 1.78 vs. 0.74 for CS plantation (Table 6).Interestingly, the FO forest had nearly identical sensitivity values for storms smaller (1.06) and larger (1.05) than mean P g (Table 6).Sensitivity of the very humid QC 1 natural forest was greater for large events (1.07 against 1.24) in comparison with the semi-humid QC 2 and Mediterranean QC 3 forests (Table 6).

Projected changes in annual precipitation and storm magnitude
During the 2020-2050 period, the low emissions (RCP 2.6) scenario predicted increased annual precipitation on average in all climate types, excluding the semi-humid climate, yet the high emissions (RCP 8.5) scenario predicted decreased annual precipitation on average in all climates excluding the arid climate (Table 7).Projected changes to mean P g from the low emission scenario was +16.0%, +3.2%, 0%, and -1.8% for Mediterranean, very humid, arid, and semihumid forests, respectively (Tables 5 and 7).The high emissions scenario predicts that mean P g will change by -5% in very humid forests and -10% in arid forests, but mean P g will not change for forests in the Mediterranean and semi-humid climates (Tables 5  and 7).Table 1 for the tree species represented by the location codes.

Storm magnitude and throughfall
We collected a wide range of storms across sites hosting common forests in all of Iran's climate types (Tables 1-2; Fig. 1).Stands for which comparable TF data exist (e.g., FO and QC 1 ) are comparable to those reported elsewhere, for example: FO versus Fagus sylvatica (e.g., Staelens et al., 2008) and QC versus Quercus serrata and Quercus acutissima (e.g., Toba & Ohta, 2005).Relationships between P g and TF for individual forest often indicate differences in the canopy-meteorological interactions that control TF production by plants (Zhang et at., 2016).In this study, as in others, P g accounted for most (nearly all) of TF variability regardless the species, yet the slope and intercept of the P g -TF correlation varied across the forests of differing species-climate combination (Table 4).Relative TF:P g in our study, 39.3-74.5%,appears to be strongly tied to forest structure (e.g., density, seasonal change, vegetation area index, gap fraction, and canopy storage capacity), and climate conditions (P g and, perhaps storm intensity and wind conditions) as shown by many others (Crockford & Richardson, 2000;2011;Staelens et al., 2008;Muzylo et al., 2012;Motahari et al., 2013).Results indicated that, when the same species experiences different climate conditions, they can differentially partition P g into TF-something particularly evident QC in three different climates (Table 4).

Throughfall sensitivity to storm amount
Large fluctuations in TF sensitivity were observed between different species, climate types, and storm sizes (greater or less than mean annual P g ) common to Iran (Table 6).For forest managers, it is apparent that different species selected for planting in different climates will exhibit differing TF sensitivities to shifting P g (Fig. 3).Practically, managers can use these sensitivity coefficients (Fig. 3 or Table 6) and projected changes to P g (Tables 5 versus 7) and estimate the potential shift in TF supply to the forest surface.For example, PE's sensitivity coefficient (2.35) indicates that a 10% decrease in P g could approximately reduce TF by 23%.Although we identified little change in the total yearly precipitation in the arid climate for the most recent decade (where these PE plantations are generally situated), an approximate estimation of 8% decrease in the storm size observed in the recent decade can induce an 18% decrease in TF in these plantations.Greater sensitivity coefficients for PE (and CA) plantations may be a result of (1) these species having larger storage capacities than other tree species in this study (Sadeghi et al., 2015), (2) their arid climate allowing their canopies to dry more efficiently between storms, and (3) most-to-nearly all the storms being small (Table 6).The higher sensitivity of these species to P g has implications for forest managers dealing with climate change, since both are the most widely-used species for afforestation in the arid and semiarid regions of Asia (e.g., Iran, Lebanon, Syria, Pakistan, Iraq, and Afghanistan) and are frequently selected due to their greater tolerance to drought, high/low temperature extremes, and being faster-growing than native broador needle-leaved tree species (Jazirei, 2009).Under a changing climate, water resource management in arid regions is therefore complicated by afforestation initiatives and forest managers may be able to use TF sensitivity as a variable in their decision to choose a Table 6.Throughfall (TF) sensitivity coefficients classified for storm magnitudes (P g ) smaller and larger than mean P g recorded in the measurement sites.See Table 1 for the tree species represented by the location codes.species with the least sensitivity to expected shifts in P g .
The remaining species located in the Caspian forests of northern Iran with very humid, semi-humid and Mediterranean climates exhibited roughly the same sensitivity (mean = 1.13;Fig. 3).However, according to our results, the replacement of QC 3 natural forest with CS man-made in the Mediterranean areas of the eastern Caspian region will likely increase TF sensitivity by 0.31 (Fig. 3).This increased TF sensitivity with converting Mediterranean QC forests to CS may be a product of greater tree density of young CS (Table 1) or, more theoretically, greater roughness lengths and zero-plane displacement heights (Rutter et al., 1975;Valente et al., 1997).Moreover, CS's scale-like needle-leaves have different leaf shedding habits compared to QC 3 's broadleaves, which will differentially interact with P g to alter TF and its sensitivity (Pypker et al., 2011).
Natural broad-leaved forests in the Caspian region (FO and QC 1 ) showed roughly the same sensitivity (Fig. 3).The primary function of the Caspian forests, other than wood production, is conservation of soil and water resources (Sagheb Talebi et al., 2014).Thus, it is expected that the low TF sensitivity of FO and QC 1 to shifts in P g will help buffer the region against climate changes.But, restoration of the Caspian deciduous forest of northern Iran by planting species, like CS (which had a greater TF sensitivity coefficient), were extensively performed by Iran's Forest, Rangeland, and Watershed organization, and may have considerable effects on ecosystem ecohydrology through altered TF supply.In contrast, introduction of PB in the semi humid climate of the central Caspian region for restoration of degraded forests showed roughly the same, even lower, sensitivity compared to native species (AV and QC 2 ) (Fig. 3).Thus, understanding the relationship between hydrologic cycling and the impact of afforestation projects on these variables can be useful for forest management and selection of suitable species for reforestation of degraded forest ecosystems.
Although the consistently-observed strong correlation between stand-scale TF and storm amount across forest types and climates confirms that storm amount is the primary factor affecting TF generation after canopy structures are saturated (Levia et al., 2011), the degree to which stand-scale TF responds to storm amount has been linked to vegetation structure (leaf area index, crown depth, etc.) (Staelens et al., 2008;Toba & Ohta 2005).Consequently, it is advisable to incorporate vegetation characteristics with TF sensitivity in response to changes in P g responses in future investigations.

Combining throughfall sensitivity and climate projections
Changes in P g to an area are expected to be compounded by the TF sensitivity of each forest type.An example "rough" estimation using the 16% increase in P g predicted by RCP 2.6 scenario in the Mediterranean climate results in a 20% and 15% increase in TF generation in CS and QC 3 forests, respectively.The projected 3.2% increase (by RCP 2.6) and 5% decrease (by RCP 8.5) in mean P g , however, in the very humid climate of Iran where the Caspian forests are located (Table 7) is not anticipated to significantly alter TF supply to the forest floor due to the tree species' low sensitivity (Fig. 3).Despite the highest sensitivity of plantations in the arid climate (~2), very slight changes in mean P g are predicted by the RCP 2.6 scenario which may not influence TF receipt at the surface.However, under the RCP 8.5 scenario, a 10% decrease in mean P g in the arid climate could decrease TF by 23% and 15% per event beneath PE and CA plantations, respectively.Fewer impacts may be realized in the semi-humid climate where the dominant tree species showed low TF sensitivity and not significant changes in mean event size are predicted by both RCPs.Clearly, these climate change scenarios (and others) can be a reference for setting minimum and maximum configurations of forest cover in water management and planning associated with adaptation.There will be large uncertainty among GCM models, most GCM models may have different performances across regions, and, as a consequence of lacking large scale observations, downscaling from GCM grid data to local areas is difficult.Moreover, the delta change factor strategy used in this study does not adjust climate projections, but assumes that the signal or changes are reasonably projected by climate models, even though the models are biased (Leng & Tang, 2014).A key feature of this approach is that, because the method uses historical precipitation as its basis and GCM data only to change the magnitude of the historical precipitation, it fails to account for changes in P variability predicted by different GCM models (Leng & Tang, 2014).Still, the projections and downscaling of changes in P g under the most conservative (low emission) and the highest greenhouse gas emissions pathway (high emission) scenarios (Table 7) underscores the necessity of work to understand TF sensitivity of common forest types in a region.

Conclusion
This national assessment of throughfall beneath common forest types in Iran showed large fluctuations in throughfall sensitivity for species planted in different climates.Arid needle-leaved plantations exhibited generally higher sensitivity responses to changing storm amount, while very humid, natural broad-leaved forests had low sensitivity responses.Projections of mean storm magnitude under a low emission scenario indicated that, compared to historical data, modest changes would be expected for most climates (excluding arid) on average during 2020-2050.However, under a high emission scenario, larger alterations are expected in all climate types except the arid climate.Results indicate that any projected change in storm size will not simply be translated directly to a change in throughfall.Rather, the increase or decrease in throughfall amount may be better predicted as a product of the projected change in storm magnitude and the forest-specific throughfall sensitivity.These findings have implications for selection of the most appropriate and adapted species for afforestation under climate change.

Figure 2 .
Figure 2. Example regression of throughfall (TF) against storm magnitude (P g ) for data collected at the Fagus orientalis site.Through shifting P g by any percentage (ΔP g ), the regression equation (in panel a) can be used to calculate a corresponding percent change in TF (ΔTF).These shifts are (b) plotted and a regression calculated where the slope is a nondimensional coefficient of "sensitivity."

Figure 3 .
Figure3.Mean throughfall (TF) sensitivity coefficients for the measurement sites.See Table1for the tree species represented by the location codes.

Table 1 .
Forest sites where throughfall (TF) data were collected, including their climate, geographical coordinates, elevation, and stand characteristics."G" and "D" refer to growing and dormant seasons, respectively, and dbh refers to diameter at breast height.

Table 2 .
: Baltas, 2007)her stations included in this study, their climate types according to the De Martonne aridity index (I DM: Baltas, 2007), and the duration of their measurement record.

Table 3 .
Characteristics of recorded rainfall events (P g ) in the measurement sites with respect to the different climate types.n refers to the number of recorded events.

Table 5 .
Historical characteristics of annual precipitation, number of storms, and other related rain event magnitude (P g ) characteristics recorded in the meteorological stations nearest to the measurement sites.

Table 7 .
Projections of annual precipitation, and mean storm magnitude (P g ) during the period of 2020-2050 predicted by the GCM under RCP 2.6 (low emission), and RCP 8.5 (high emission) scenarios.