Multiple Temporal Scales Assessment in the Hydrological Response of Small Mediterranean-Climate Catchments

Mediterranean-climate catchments are characterized by significant spatial and temporal hydrological variability caused by the interaction of natural as well human-induced abiotic and biotic factors. This study investigates the non-linearity of rainfall-runoff relationship at multiple temporal scales in representative small Mediterranean-climate catchments (i.e., <10 km2) to achieve a better understanding of their hydrological response. The rainfall-runoff relationship was evaluated in 43 catchments at annual and event—203 events in 12 of these 43 catchments—scales. A linear rainfall-runoff relationship was observed at an annual scale, with a higher scatter in pervious (R2: 0.47) than impervious catchments (R2: 0.82). Larger scattering was observed at the event scale, although pervious lithology and agricultural land use promoted significant rainfall-runoff linear relations in winter and spring. These relationships were particularly analysed during five hydrological years in the Es Fangar catchment (3.35 km2; Mallorca, Spain) as a temporal downscaling to assess the intra-annual variability, elucidating whether antecedent wetness conditions played a significant role in runoff generation. The assessment of rainfall-runoff relationships under contrasted lithology, land use and seasonality is a useful approach to improve the hydrological modelling of global change scenarios in small catchments where the linearity and non-linearity of the hydrological response—at multiple temporal scales—can inherently co-exist in Mediterranean-climate catchments.


Introduction
The complexity of Mediterranean fluvial systems is caused by the multiple temporal and spatial heterogeneities in the relationships between the natural and human-induced abiotic and biotic variables, especially in the Mediterranean Sea Region [1,2]. The Mediterranean climate lies between 32 • and 40 • N and S of the Equator and is characterized by a wet and mild winter, a warm and dry summer and a high inter-and intra-annual variability in rainfall patterns. Mediterranean climate regions comprise the Mediterranean Sea Region, the coast of California, Central Chile, the Cape region of South A temporal downscaling assessment of the inter-and intra-annual variability of the rainfall-runoff relationship was carried out in the Es Fangar Creek catchment, a headwater tributary of the Sant Miquel River catchment (151 km 2 ) located in the north-eastern part of Mallorca Island ( Figure 1b) and is representative of Mediterranean mid-mountainous catchments. The lithology is mainly composed of marl and marl-limestone formations from the medium-upper Jurassic and Cretaceous period in the valley bottoms. In the upper parts of the catchment, massive calcareous and dolomite materials from the lower Jurassic period and dolomite and marl formations from the Triassic period (Rhaetian) are dominant. The Es Fangar catchment has an area of 3.4 km 2 , with altitudes ranging from 72 m.a.s.l. to 404 m.a.s.l. (Figure 1c). The mean slope of the catchment is 26% and the length of the main channel is 3.1 km (average slope of 22%). The drainage network is natural in the headwater parts. In the bottom valley, flow lamination is applied, with transverse walls and also the straightening and diverting of the main stream, with the banks fixed with dry-stone walls for flood control and erosion prevention. In addition, subsurface tile drains are also installed to facilitate drainage due to the impervious materials which would impede agricultural activity during wet periods. As a result, 16% of the surface catchment is occupied by soil and water conservation structures. Since 1950, important socio-economic changes have caused a gradual abandonment of farmland in marginal areas, leading to afforestation. The land uses in 1956 were rainfed herbaceous crops (54%), forest (31%) and scrubland (15%). Nowadays, the main land uses (Figure 1c) are forest (63%), rainfed herbaceous crops (32%) and scrubland (5%). In addition, 54% of terraced land is currently covered by forests (Figure 1c), demonstrating the consolidation of the forest transition. The climate of the area is classified on the Emberger scale [49] as Mediterranean temperate sub-humid. The mean annual rainfall (1965-2016, Biniatró AEMET station) is 927 mm y −1 with a variation coefficient of 23%, and the mean annual temperature is 15.7 • C. A rainfall amount of 180 mm in 24 h is estimated to have a recurrence period of

Hydrological Response of Small Mediterranean-Climate Catchments
Bivariate statistical regressions were used to establish the correlations at the annual and event scales between rainfall and runoff in order to assess the hydrological response of small Mediterranean-climate catchments (i.e., <10 km 2 ; Figure 1c). These equations were used descriptively as least squares adjusted for the rainfall-runoff assessment. At the annual scale, data from the 43 representative catchments

Hydrological Response of Small Mediterranean-Climate Catchments
Bivariate statistical regressions were used to establish the correlations at the annual and event scales between rainfall and runoff in order to assess the hydrological response of small Mediterranean-climate catchments (i.e., <10 km 2 ; Figure 1c). These equations were used descriptively as least squares adjusted for the rainfall-runoff assessment. At the annual scale, data from the 43 representative catchments (Table A1 (a); hereinafter A1a) were collected to observe the influence of lithology on this response; i.e., the catchments were classified as pervious or impervious by using the information regarding the catchments' characteristics (e.g., soil type, soil texture or lithology materials) extracted from research papers. At the event scale, 203 events from 12 representative catchments were classified according to (a) seasonal occurrence (autumn, winter, spring or summer), (b) pervious or impervious lithology and (c) main land use (agricultural, agroforestry, forest or shrub) (Table A1 (b)). The main land uses of each catchment were divided into specific land uses when the information was available ( Figure A1). Most of these studies (79%) were located in the Mediterranean Sea Region and 47% of them studied catchments < 1 km 2 .

Monitoring and Data Acquisition in Es Fangar
The rainfall data since 2012 were obtained from the B696 Biniatró AEMET station (Figure 1d; 1 km away from the catchment). In October 2014, a rainfall gauge station (Míner Gran) was installed less than 2.5 km away from the catchment. The Míner Gran rainfall gauge is located 1 m above the ground and connected to a HOBO Pendant ® G Data Logger-UA-004-64 that records precipitation at 0.2 mm resolution. A linear regression was established (n = 978; R 2 : 0.88) for daily rainfall (2014-2017) between the Biniatró and Míner Gran stations to reconstruct rainfall data series from 2012 to 2014 for the Míner Gran station. Due to the lack of temperature data available in the studied catchment, the data of neighbouring AEMET weather stations (i.e., less than 8 km away from the catchment) were used to estimate the catchment's temperature by using the block kriging technique. With this information, the monthly evapotranspiration (i.e., ET o ) was estimated using the equation from Hargreaves and Samani [51].
The gauging station of the Es Fangar catchment was built in July 2012. Its cross section is formed by a rectangular broad-crested weir for low water stages to better measure low discharges. The water level was continuously (1 min time step) measured using a pressure sensor Campbell CS451 connected to a CR200X datalogger and average readings were kept every 15 min. The flow velocity was measured during baseflow conditions and flood events using an OTT MF Pro electromagnetic water flow meter. These flow velocity measurements (n = 17) were subsequently used to calibrate the stage-discharge relationship.

Rainfall-Runoff Relationship Assessment in Es Fangar
This study is based on data from 5 hydrological years (2012-2013 to 2016-2017). For each runoff event (i.e., when the water stage exceeds the low-flows channel; i.e., 0.036 m 3 s −1 ), a simple hydrograph separation between quickflow and baseflow components was performed through a visual technique based on the breakpoints detected on the logarithmic falling limb of the hydrograph [52]. This hydrograph separation method was used only to characterize the response of the catchment to a rainfall, with no aim of deriving any interpretation in terms of runoff processes. Over a 5-year period, hydrograph separation was conducted on 49 events. Based on the hydrograph separations, the baseflow and quickflow contributions for each hydrological year (October to September) were calculated.
Several variables were derived from the hyetograph and hydrograph from each rainfall-runoff event (Table 1). These variables aimed to characterize the pre-event conditions (antecedent precipitation in 24 h, AP1d; baseflow at the start of the event, Q 0 ) as well as the event characteristics (rainfall depth, P tot ; mean and maximum rainfall intensity, IP mean 30 and IP max 30; runoff depth and runoff coefficient, R and R c ; and peak-flow discharge, Q max ). The relationships between these variables were assessed Water 2020, 12, 299 6 of 27 through the Pearson correlation matrix. A more detailed analysis of rainfall-runoff relationships in the Es Fangar catchment was carried out in a second step to investigate the variability of the hydrological response to similar size rainfall events to assess the rainfall-runoff non-linearity. The widest range of events with a similar P tot (40-50 mm) generating the highest differences in R response was selected. Accordingly, 7 events with a P tot range from 41.8 to 49.8 mm but with different antecedent conditions or rainfall dynamics were selected and also classified in relation to wet, dry and transition periods, in accordance with Gallart et al. [53]. The annual rainfall ranged from 376 to 1794 mm yr −1 , whilst R ranged from 3.7 to 1200 mm within these representative catchments under Mediterranean-climate conditions. The relationship between annual rainfall and R ( Figure 2) showed a significant positive linear correlation (R 2 = 0.68; p < 0.01). However, some scattering was also apparent in the relationship because when catchments with pervious lithology were not included, the regression increased (R 2 = 0.82; p < 0.01; see Figure 2). The lowest annual R values (<10 mm) are related to catchments with annual rainfall < 500 mm yr −1 and a pervious lithology (Figure 2), located in the Eastern Mediterranean Sea (Lower Jordan River). Besides, catchments with pervious lithology with ca. 1500 mm yr −1 of annual rainfall showed large differences in their R amount, ranging from 70 to 974 mm.

Annual Scale: Lithology Influence
The annual rainfall ranged from 376 to 1794 mm yr −1 , whilst R ranged from 3.7 to 1200 mm within these representative catchments under Mediterranean-climate conditions. The relationship between annual rainfall and R ( Figure 2) showed a significant positive linear correlation (R 2 = 0.68; p < 0.01). However, some scattering was also apparent in the relationship because when catchments with pervious lithology were not included, the regression increased (R 2 = 0.82; p < 0.01; see Figure 2). The lowest annual R values (<10 mm) are related to catchments with annual rainfall < 500 mm yr −1 and a pervious lithology (Figure 2), located in the Eastern Mediterranean Sea (Lower Jordan River). Besides, catchments with pervious lithology with ca. 1500 mm yr −1 of annual rainfall showed large differences in their R amount, ranging from 70 to 974 mm. To investigate the variability of the hydrological response at the event scale, the rainfall-runoff relationships of 203 events from 12 small Mediterranean-climate catchments were analysed (Table A1b and Figure 3). Within Figure 3, the outliers are marked with an ellipsoid and excluded from the  To investigate the variability of the hydrological response at the event scale, the rainfall-runoff relationships of 203 events from 12 small Mediterranean-climate catchments were analysed (Table A1 (b) and Figure 3). Within Figure 3, the outliers are marked with an ellipsoid and excluded from the correlations because outliers enlarged the rating curves of these correlations, obtaining significant relationships in all cases.  The seasonal distribution of the events was winter (34%), autumn (32%), spring (22%) and summer (12%)-most of them occurred between November and February (47%). Similar seasonal median rainfall was observed, ranging from 26.2 (winter) to 29.0 mm (autumn). The highest seasonal median of event R occurred in winter (4.2 mm). However, those events occurred during the transition periods depicted a similar median R (autumn = 2.3 mm; spring = 2.2 mm) being the lowest value for summer events (0.6 mm). In the case of rainfall-runoff relationships, the highest correlation was obtained in spring (R 2 = 0.63), followed by winter, autumn and summer (Figure 3a).
Small differences between the median of R events in catchments with pervious lithology (2.2 mm) and catchments without pervious lithology (3.1 mm) were observed. Nevertheless, significant differences in rainfall-runoff relationships were detected as events in catchments with pervious lithology showed the highest correlation and the lowest scattering ( Figure 3b). Rainfall events > 55 mm generated a R > 4 mm, excepting events occurred in late spring (Figure 3a) in catchments with predominance of pervious lithology (Figure 3b).
Event rainfall-runoff relationship was carried out considering the differences in the main land uses of studied catchments (Figure 3c). The median event R in forest (37.6 mm) was higher than agricultural (5.0 mm), shrub (2.0 mm) and agroforest (1.5 mm) catchments. However, the highest correlations between rainfall and R were obtained in catchments with a predominance of forest (R 2 = 0.96), agricultural (R 2 = 0.61) and agroforest (R 2 = 0.58) land uses.

Annual Scale
The mean annual rainfall (840.8 mm yr −1 ± 213 mm y −1 ) calculated over the five hydrological years (2012-2017) was broadly representative (−9%) of the long term mean annual rainfall (927 mm y −1 ± 215 mm y −1 : 1965-2016), but showed a high range of variation (coefficient of variation 25%, being 553 mm y −1 to 1090 mm y −1 ), characteristic of Mediterranean conditions. Figure 4 shows the annual variability of rainfall, baseflow and quickflow contributions as well the annual R c and the number of days with observed flow at the gauging station. Linear positive relationships were observed between annual rainfall, R, R c , baseflow and quickflow (R 2 ≥ 0.84, data not shown). The annual R c ranged from 2.9% to 14.2% (mean value = 10.4%) and quickflow contribution from 9.9% to 45.0% (mean value = 33.0%). During the study period, flow was observed 42.8% of the time. From year to year, the cumulated number of days with flow ranged from 37.8% (138 days) to 65.2% (238 days) of the total days ( Figure 4). An inverse relation between annual R and the number of days with flow was established. On the one hand, hydrological years with the largest R (2013-2014, 2014-2015 and 2016-2017) showed fewer days with flow and a lower baseflow contribution (<60%). In these years, 50% to 62% of the annual R was reached in 5 or fewer days (Table A2). In addition, days with more R were always during autumn and winter. On the other hand, the hydrological years with the lowest R (2012-2013 and 2015-2016) showed 208 and 238 days with flow. The baseflow contributions represented 70% and 90%, respectively, and 50% of the annual R was reached in 10 and 17 days, also respectively. The contribution of the 5 days with more R did not exceed 28% and 37%, respectively, and these days with more R were distributed among the four seasons.   Figure 5 shows the monthly mean values of rainfall, reference evapotranspiration and minimum, median and maximum R values. The mean monthly rainfall values were broadly comparable to those observed for the long-term period, except for October, when rainfall during the study period was much lower (56% lower than long-term rainfall).

Seasonal Scale
The seasonal dynamics of rainfall and evapotranspiration controlled the R response. Characteristic wet (winter) and dry (summer) periods alternated throughout the year, separated by transition periods (last autumn and early spring). Autumn was the rainiest season, with low evapotranspiration and flow observed at the outlet for 52.7% of the time, on average. The autumn mean Rc was 9.1% and ranged from 4.1% to 11.4%. During November and December, a continuous baseflow was generated in response to a large rainfall amount after the pre-filling of the initial water storage in October. From October to December, the mean Rc increased from 1.3% to 13.1%. Finally, the minimum, median, average and maximum monthly R amounts in December were higher than in November, even with lower rainfall amounts ( Figure 5). Although differences between the hydrological years were observed, autumn was the season with less inter-variability in terms of R response.
During winter, flow occurred during 90.6% of the time, although the mean rainfall amount (294 mm) was lower than in autumn (326 mm). This more frequent presence of flow was caused by the low evapotranspiration demand, thus keeping the hydrological pathways active, as already established in autumn. Consequently, in winter, a R amount similar to that of autumn can be generated from a lower rainfall amount. The winter mean R coefficient was 16.9%, ranging from 1.5% to 27%. During winter, the mean monthly Rc decreased from 18.9% to 14.6% between January and March as a consequence of a decreasing rainfall amount and an increasing evapotranspiration demand.  Figure 5 shows the monthly mean values of rainfall, reference evapotranspiration and minimum, median and maximum R values. The mean monthly rainfall values were broadly comparable to those observed for the long-term period, except for October, when rainfall during the study period was much lower (56% lower than long-term rainfall).

Seasonal Scale
The seasonal dynamics of rainfall and evapotranspiration controlled the R response. Characteristic wet (winter) and dry (summer) periods alternated throughout the year, separated by transition periods (last autumn and early spring). Autumn was the rainiest season, with low evapotranspiration and flow observed at the outlet for 52.7% of the time, on average. The autumn mean R c was 9.1% and ranged from 4.1% to 11.4%. During November and December, a continuous baseflow was generated in response to a large rainfall amount after the pre-filling of the initial water storage in October. From October to December, the mean R c increased from 1.3% to 13.1%. Finally, the minimum, median, average and maximum monthly R amounts in December were higher than in November, even with lower rainfall amounts ( Figure 5). Although differences between the hydrological years were observed, autumn was the season with less inter-variability in terms of R response.
During winter, flow occurred during 90.6% of the time, although the mean rainfall amount (294 mm) was lower than in autumn (326 mm). This more frequent presence of flow was caused by the low evapotranspiration demand, thus keeping the hydrological pathways active, as already established in autumn. Consequently, in winter, a R amount similar to that of autumn can be generated from a lower rainfall amount. The winter mean R coefficient was 16.9%, ranging from 1.5% to 27%. During winter, the mean monthly R c decreased from 18.9% to 14.6% between January and March as a consequence of a decreasing rainfall amount and an increasing evapotranspiration demand.
In spring, flow was observed 41.3% of the time. The monthly rainfall decreased during the season and was much lower (127.4 mm) than in autumn and winter. Although monthly evapotranspiration losses increased and were higher than rainfall, the water reserves accumulated during the autumn and winter months sustained the flow contribution. The spring mean R c was 5.1%, ranging from 0.4% to 8.4% during the study period. The average monthly R decreased during the season from April (7.7%) to June (0.3%).
Finally, in summer, flow was observed only 0.9% of the time. The stream was most often dry in this season due to the high negative balance between rainfall and evapotranspiration. R was only ephemeral in response to the episodic rainfall events more frequent in September, when 80% of rainfall occurred on average ( Figure 5). The summer R c was 1.4%. In spring, flow was observed 41.3% of the time. The monthly rainfall decreased during the season and was much lower (127.4 mm) than in autumn and winter. Although monthly evapotranspiration losses increased and were higher than rainfall, the water reserves accumulated during the autumn and winter months sustained the flow contribution. The spring mean Rc was 5.1%, ranging from 0.4% to 8.4% during the study period. The average monthly R decreased during the season from April (7.7%) to June (0.3%).
Finally, in summer, flow was observed only 0.9% of the time. The stream was most often dry in this season due to the high negative balance between rainfall and evapotranspiration. R was only ephemeral in response to the episodic rainfall events more frequent in September, when 80% of rainfall occurred on average ( Figure 5). The summer Rc was 1.4%.

Non-linearity Assessment at Event Scale
During the study period, the number of flood events per hydrological year was between 2 (2015-2016) and 14 (2014-2015). Their seasonal distribution was 13 events in autumn, 25 in winter, 7 in spring and 4 in summer. The events' characteristics are detailed in Table A3. An investigation of the variability and the non-linearity of the hydrological response at the event scale was carried out by means of a comparative assessment of seven rainfall-runoff events with similar Ptot, ranging from 41.8 to 49.8 mm, but different antecedent conditions and rainfall dynamics ( Figure 6).
The event on 21 January 2017 ( Figure 6a) occurred under wet conditions following a runoff event on the day before. AP1d was 35.4 mm and Q0 was high (0.229 m 3 s −1 ). The rainfall characteristics

Non-Linearity Assessment at Event Scale
During the study period, the number of flood events per hydrological year was between 2 (2015-2016) and 14 (2014-2015). Their seasonal distribution was 13 events in autumn, 25 in winter, 7 in spring and 4 in summer. The events' characteristics are detailed in Table A3. An investigation of the variability and the non-linearity of the hydrological response at the event scale was carried out by means of a comparative assessment of seven rainfall-runoff events with similar P tot , ranging from 41.8 to 49.8 mm, but different antecedent conditions and rainfall dynamics ( Figure 6).
(duration and intensity) were similar to those observed for other events of the same magnitude ( Figure  6b,e,g), but the hydrological response was an order of magnitude higher in terms of R (28.9 mm), Rc (65%) and Qmax (3.942 m 3 s −1 ) due to wet antecedent conditions (Table A4). This event had a similar duration to events that occurred under higher rainfall intensities (11 November 2012, 29 September 2014). However, wet antecedent conditions were much more favourable than rainfall intensity to reach a high Qmax and especially a higher R contribution. Figure 6. Selected events for non-linearity analysis at Es Fangar Creek. Events with a total precipitation between ca. 40 and 50 mm. The numbers located at the peak of each hydrograph indicate the ID of the events recorded in Es Fangar Creek during the study period 2012-2017, further exposed in Table A3. . Selected events for non-linearity analysis at Es Fangar Creek. Events with a total precipitation between ca. 40 and 50 mm. The numbers located at the peak of each hydrograph indicate the ID of the events recorded in Es Fangar Creek during the study period 2012-2017, further exposed in Table A3.
The event on 21 January 2017 ( Figure 6a) occurred under wet conditions following a runoff event on the day before. AP1d was 35.4 mm and Q 0 was high (0.229 m 3 s −1 ). The rainfall characteristics (duration and intensity) were similar to those observed for other events of the same magnitude (Figure 6b,e,g), but the hydrological response was an order of magnitude higher in terms of R (28.9 mm), R c (65%) and Q max (3.942 m 3 s −1 ) due to wet antecedent conditions (Table A4). This event had a similar duration to events that occurred under higher rainfall intensities (11 November 2012, 29 September 2014. However, wet antecedent conditions were much more favourable than rainfall intensity to reach a high Q max and especially a higher R contribution. The event on 4 April 2012 (Figure 6b) occurred in early spring, ergo at the beginning of the transition period from wet to drier conditions. The antecedent conditions were favourable to R generation, due to the presence of a baseflow maintained by water reserves accumulated in the autumn and winter months.  (Table A4). Nevertheless, R and R c were slightly larger for the event on 8 February 2017, likely due to the large rainfall amount accumulated during the two events (more than 400 mm of rainfall between 19 December 2016 and 8 February 2017).
The event on 7 May 2016 (Figure 6g) was similar to the event on 19 December 2016 (Figure 6e) in terms of P tot , duration and intensity. Nevertheless, R response started earlier (i.e., shorter response time) for the event on 19 December 2016. This difference was caused because the December 2016 event occurred during the wetting-up period when catchment water reserves were increasing. As a result, whilst the R response was relatively small (4.4 mm), R c (9%) and Q max (0.394 m 3 s −1 ) were significantly larger than for the event on 7 May 2016 (Table A4), which occurred during the transition period towards dry conditions. With an AP1d of 1.6 mm and Q 0 of 0.006 m 3 s −1 , the lowest R response was generated (see total runoff, R c and Q max in Table A4). This R response was also the most delayed, as shown on Figure 6g, starting when 95% of the rain (39.8 mm) had fallen on the catchment.
The event on 29 September 2014 (Figure 6f) occurred during the dry season (i.e., late summer, early autumn) but with an AP1d of 11.8 mm. The rainfall event was the shortest of these seven selected events (3.2 h), with the second highest IP max 30 (56.4 mm h −1 ) of the study period. This very high intensity triggered a R response during dry conditions. Similar rainfall but lower intensities during dry conditions always produced lower responses (R c ≤ 3%). Although the R response was relatively small (4.2 mm) with R c < 10%, Q max was relatively high (2.356 m 3 s −1 ), promoted by the rainfall intensity.

Small Mediterranean-Climate Catchments: Lithology Influences
Catchments with annual rainfall ranging from 500 to 900 mm yr −1 showed annual R ca. 20-350 mm, where the role played by lithology on the annual R response was important. Accordingly, the highest annual R values (>350 mm) were observed in catchments with badland areas [34] and a high subsoil clay content promoting lateral flow and a perennial regime [35]. Swarowsky et al. [54] also reported an annual R 17% lower than in an adjacent catchment to that investigated by Lewis et al. [35] and related this difference to greater deep infiltration. Similarly to our results, Merheb et al. [46] obtained a significant rainfall-runoff correlation (R 2 = 0.69) using 160 catchments from the Mediterranean Region (0.35 to 21,700 km 2 ). These authors highlighted that the scattering observed was due to karstic catchments or snowmelt contributions.
Also, scattering can be observed in the relationship for catchments with annual rainfall >1000 mm yr −1 . For this group of catchments, the annual R ranged from 135 to 1200 mm, except for the Santa Magdalena catchment, which is entirely on a limestone lithology and yielded an annual R value of 70 mm [47]. Large annual R values usually corresponded to watertight catchments with R values > 1000 mm yr −1 for two small Mediterranean-climate catchments on granite bedrock [36]. However, karstic environments can also promote very high annual R values, e.g., Ref. [55], as a result of incoming karstic spring sources. The conceptual model developed by Borg Galea et al. [2], which is based in exogenous and endogenous variables that link natural and human-derived variables and their influences in the system on a spatial-temporal scale, emphasized that climate is the main exogenous driver of the rainfall-runoff relationship, although their processes through precipitation, temperature and wind are mediated by catchment geology (endogenous variable).
Different seasonal rainfall-runoff relationships in relation to the influence of rainfall and evapotranspiration [56] are related to the alternation of rainfall and reference evapotranspiration throughout the year reproducing wet (winter), dry (summer) and transition periods (last autumn and early spring) [57]. Accordingly, winter is the season with the highest R c , from 17% to 56%, due to rainfall being accumulated during autumn-and maintained in winter-and low evapotranspiration demand [10,15,56]. The same pattern was also observed during the five hydrological years assessed in Es Fangar Creek. In addition, all these studies also emphasized that in summer the hydrological response is limited or null, with R c < 10%. Serrano-Muela [56] observed runoff coefficients < 5% in a range of rainfall amounts of 15-200 mm during summer in a small catchment characterized by a large forest cover. The driest environments in these studies obtained a summer R c < 2%. R c in autumn and spring were similar, ranging from 9% to 25% and 5% to 28% respectively. However, a variability throughout the seasons can be observed, especially in late autumn and early spring, as these are transition periods. Additionally, Lana-Renault [10] observed that the wetting-up period was steeper than the drying-down period, which was more progressive in a catchment located in the Pyrenees. Autumn was the season after the driest period, and it was when the evapotranspiration reached the maximum values. The beginning and the end of this season can be quite different in relation to R c . The catchments have a null or limited response until a succession of rainfall events that fill the water reserves and generate favourable conditions for R generation [56]. The findings of Lana-Renault [10] in spring concluded that the accumulated rainfall during autumn and winter maintain high water reserves, allowing high runoff coefficients, even though spring was not the rainiest season.

Es Fangar
In the Es Fangar catchment, a low mean annual R value (87 mm) was recorded during the 2012-2017 period, despite its mean annual rainfall value (black dot on Figure 2). This result is most likely related to the presence of massive karstic dolomites and breccias in the lithology of the catchment headwaters.
The reported annual R in Es Fangar is within the R range of small catchments located in Mallorca, with annual rainfall ranging from 500 to 900 mm yr −1 , which yielded annual R between 36 and 130 mm. From these catchments, the highest R value was observed in an agricultural small lowland catchment prone to receiving subsurface contributions [15]. Spatial differences in the annual R between the headwater (102 mm) and catchment outlet (36 mm) due to transmission losses along the main channel were detected in the Sa Font de la Vila River, a mid-mountainous terraced catchment affected by forest fires [58]. Within this catchment, temporal differences were observed as R c ranged from 1% to 22% from one year to year.

Seasonality, Lithology and Land Use Influences on Small Mediterranean-Climate Catchments Hydrological Response
The seasonal dynamics of rainfall and evapotranspiration during the year generate wet, dry and transition periods, which control the R response [10,53]. Therefore, seasonally, the highest correlation observed in spring could be related to the water reserves accumulated during the autumn and winter seasons, which is favourable for R generation [13]. According to the seasonal distribution of the 203 events, other studies in Mediterranean-climate catchments observed an event seasonality in the hydrological response [8,10,59] where autumn rainfall events produced a limited or null R response, especially at the beginning of the season. In addition, in winter and spring, the R c values were always larger than 3%, being the maximum in spring (70%).
The role of transmission losses due to lithology has been investigated [32,60], depicting how low runoff values were recorded due to infiltration. In this study (see Section 3.2), events assessed in catchments with a predominance of pervious lithology showed a higher scattering than events in impervious catchments. In addition, rainfall events >55 mm generated R > 4 mm, except in events occuring in late spring in catchments with pervious lithology. In this case, the transmission losses due to pervious lithology have more influence on the hydrological response than the seasonal role of the water catchment reserves. Ries et al. [61] identified a low and large R response in rainfall events >50 mm based on the physiographic catchment and rainfall characteristics. Events occurring in catchments, which even had different land uses, with less bedrock permeability and less soil water storage but with higher values of rainfall intensity generate a greater R response than catchments with more bedrock permeability and soil water storage under low rainfall intensities. Contrarily, other types of lithology (i.e., badland areas) are characterized by a low infiltration capacity [62]. Latron and Gallart [63] observed that frequently badlands are the only active area to R response. According to these findings, Nadal-Romero et al. [64] observed an event stormflow coefficient of up to 20% in a Mediterranean-climate catchment with badland areas, even when the headwater catchment was forested (30% of the catchment area).
Events occurring in agricultural catchments obtained the second highest values in the median R and in the rainfall-runoff correlation, as arable land and its spatial distribution within a catchment have been demonstrated to be a driver for R generation [65]. However, events in the forest catchment had the highest correlation and median R, as the water yield and rainfall-runoff relationship decrease in forests catchments when compared to agricultural catchments, due to forest cover [66]. However, this correlation (R 2 : 0.96; n: 3) and median value is related to the presence of karstic springs in the forest catchment, which was the only forest catchment found in the literature. The median R of forest catchments could be similar to the values obtained in afforested ones, as both land uses trigger a reduction in the R generation. This reduction was quantified in a loss of water yield ca. 10-40% [38,[66][67][68][69]. Nevertheless, larger flows and Q max were observed in afforested catchments than forest ones. The shrub median R (2.0 mm) was higher than afforested catchments, as shrub coverage is characterized by lower vegetation cover than forest and afforested land uses. In addition, the second main land use in shrub catchments is agricultural land use ( Figure 1A), which could explain the higher median R compared with the afforested catchments. The same results were observed in García-Ruiz et al. [66], where catchments under shrub land use had lower R contributions than catchments with forest land use due to a lower plant cover density in the shrub catchment.

Rainfall-Runoff Relationship at the Es Fangar Catchment
The identification of driving factors explaining the variability of R c can be explored through the relationships between R c and R and P tot , rainfall, maximum intensity and baseflow at the beginning of the event, as Latron et al. [8] carried out in a mountainous Mediterranean-climate catchment in North Eastern Spain.
Only the correlation between R and R c was significant in the Es Fangar catchment (Figure A2a), explaining 73% of the variance. As in Latron et al. [8] this relationship had the highest correlation. However, contrarily to Latron et al. [8], the relationship between rainfall and R c was highly non-lineal in Es Fangar, presenting a huge scattering, especially for P tot between 20 and 60 mm (Figure A2b), which yielded R c ranging from 1% to 65%. For larger rainfall events (i.e., 60 to 100 mm), the R c values were high (>20%), but low R responses (with R c ranging 1% to 3%) were also observed when large rainfall events were recorded after or during dry conditions. In line with the results from Latron et al. [8], R c was not related with rainfall intensity ( Figure A2c).
The relationship between Q 0 and R c was not linear in the Es Fangar catchment ( Figure A2d). A threshold was identified leading to larger R c (>10%) for events with Q 0 above 0.04 m 3 s −1 , which always occurred between November and April. The relationship between Q 0 and R c in Latron et al. [8] was stronger, because of the presence of a baseflow during most of the year. In general terms, the frequent absence of a baseflow in the Es Fangar catchment between May and October worsens most of the relationships presented in Figure A2 in comparison to a catchment with an almost permanent baseflow, as in Latron et al. [8] and Nadal-Romero et al. [40]. Then, baseflow-R c correlation increased as the environment got wetter, from R 2 = 0.19 in Es Fangar to R 2 = 0.70 in a humid catchment [70].
In order to help in identifying the main factors involved in the hydrological response in the Es Fangar catchment, multiples relationships were carried out to investigate the influence of rainfall, Q 0 and IP max 30 on R and Q max (Figure 7). The hydrological response by using R and Q max was highly non-linear, as observed in other Mediterranean-climate catchments [8,11] because the soil moisture antecedent conditions are crucial for R generation. Figure 7a shows how the R response could not be explained by Q 0, as the highest R values presented a large scattering along the Y axis. This pattern confirmed that large R amounts (i.e., >10 mm) may occur from large rainfall events occurring in dry or wet conditions (i.e., with low or high Q 0 values). Therefore, P tot was a key factor for R response, as shown by the significant positive correlation (p < 0.01) between both variables ( Table 2). The largest R values were always recorded during last autumn or winter, with contributions of >17 mm ( Figure A3 and Table A5), and are characterized by P tot , R, R c , Q max and AP1d values much higher than median values (Table A3). Contrarily, large rainfall events occurring in spring resulted in a small R contribution because of the dry antecedent conditions (i.e., low or null Q 0 values). Other authors also assessed relationships between seasonality and R generation. In a mountainous catchment, Lana-Renault [10] observed that the highest rainfall-runoff correlations at the event scale were depicted in winter and spring due to higher water reserves than in autumn and summer. In Es Fangar, a group of events during autumn and winter characterized by Q 0 > 0.080 m 3 s −1 always generated R c greater than 14% ( Figure A2d). The P tot of these events ranged from 5.2 to 61.4 mm but the AP1d was larger than 20 mm, except in one event (AP1d: 3.4 mm), which was characterised as the highest IP max 30 (i.e., 77.5 mm h −1 ) of the study period. These event characteristics suggest that P tot was the main factor controlling an effective hydrological response as antecedent conditions (i.e., soil moisture degree of the catchment) have been demonstrated to play an important role in R generation [18]. Similarly, Lana-Renault [10] identified a cluster of events characterized with the highest Q 0 that generated a larger R response, but the convective rainfall occurring in the Pyrenees during summer blurred the seasonal pattern described in the Es Fangar catchment. Figure 7b shows that high Q max values were observed in response to large rainfall amounts (above 50 mm) or for flood events with a Q 0 value higher than 0.080 m 3 s −1 . Both factors had a combined effect on the generation of high Q max values as shown by the significant correlation observed between P tot and Q 0 and Q max . The role of wet conditions over Q max was also observed in other Mediterranean-climate catchments where AP1d correlated significantly with Q max [15,17,71]. (c) IP max 30-rainfall-runoff and (d) IP max 30-rainfall-Q max . The points with an ID are depicting the selected events for the non-linearity assessment, as it can be further consulted in Figure 6 and Table A3.  (Table 2). Figure 7c,d show how high values of R and Q max resulted for events with low IP max 30 when rainfall was > 50 mm or high Q 0 . This absence of correlation between rainfall and the catchment response (R and Q max ) indicates that R generation is most likely the result of saturation excess processes [72], where P tot and antecedent conditions played a key role. For events between June and September, saturation excess processes were, however, probably less dominant and combined with infiltration excess processes in response to highest IP max 30 values ranging from 19 to 56.4 mm h −1 . As a result, the events observed between June and September were shorter (5.7 to 14 h duration compared to 8 to 65.5 h for wet conditions events). Therefore, R mechanisms co-exist within a Mediterranean catchment [13] such as Es Fangar, with a marked seasonality, with saturation processes being dominant during the wet period and Hortonian R during the dry period. Consequently, the alternation of the R mechanisms generated lower and larger R contributions in events occurring during the dry and wet seasons, respectively. These processes caused the non-linearity of the hydrological response because the same amount of P tot generated a non-lineal R response due to different soil moisture conditions [15,64].

Non-Linearity Assessment at Es Fangar Catchment
The events analysed in Figure 6 presented a huge variability in their hydrological response, despite recording a similar P tot , ranging from 41.8 to 49.8 mm. R c ranged from 1% to 65%, depending on the catchment moisture conditions, rainfall intensities and seasonality characteristics. The highest hydrological response occurred under marked wet soil moisture conditions in the winter period, even with low rainfall intensities (Figure 6a). A similar complexity in R responses has already been observed in other Mediterranean-climate catchments where seasonality and antecedent soil moisture conditions played a key role in R generation [11,57]. Gallart et al. [73] also showed that different R generation behaviours could be observed during the year due to varying catchment wetness conditions and changing rainfall events characteristics. The results obtained by these authors are clearly illustrated by comparing Events 2 and 20 with Events 18 and 37 in Figure 6. Indeed, R events that occurred during wet conditions with low rainfall intensities (18 and 37) showed similar R c than events that occurred during dry conditions with high rainfall intensities (2 and 20). However, during dry conditions, runoff events resulting from low intensity rainfall yielded the lowest R responses (e.g., Event 36). Our results are also in agreement with the observations made by Schnabel et al. [74], who observed higher R contributions for low rainfall intensities events occurring in periods with high soil water content than for events with high rainfall intensities. Besides, those authors pointed out that this condition was common during years with above average rainfall (i.e., during wet years) but was rarely observed during dry years.

Conclusions
The evaluation of multiple temporal scales in contrasting small Mediterranean-climate catchments has improved the understanding of the role played by lithology and land use in the hydrological response. The assessment of rainfall-runoff relationships at the annual scale in small Mediterranean-climate catchments showed a significant linearity in the hydrological response due to the importance of the annual rainfall amount. Nevertheless, lithology effects on R generation explained an increase in the scattering in the rainfall-runoff relationship because pervious and impervious materials triggered larger and lower R contribution, respectively. Despite the significant correlation between rainfall and R, the Es Fangar Creek dataset illustrated a huge intra-annual variability of the rainfall-runoff relationship during the five hydrological years analysed, as seasonal rainfall and evapotranspiration dynamics controlled the R response. These dynamics were clearly observed in the average seasonal R c , decreasing from winter to summer. As a result, these seasonal differences should be considered as a starting point of the non-linearity generation in the rainfall-runoff relationships at the event scale.
At the event scale, lineal and non-lineal relationships were observed in the rainfall-runoff relationships in small Mediterranean-climate catchments, suggesting that different factors conditioned the R response. The total rainfall was the most significant driving factor, although the interaction between seasonality and the spatial diversity of lithology and land uses at the catchment scale also played an important role in R generation. Thus, the highest correlations at the seasonal scale were observed in those events which occurred in winter and spring when the highest water reserves favoured the R response. Lithology caused a higher dispersion in the rainfall-runoff relationships at the event scale in the set of small Mediterranean catchments because pervious materials required higher antecedent wetness conditions. Different land uses promoted a decrease in R generation, comparing agricultural with scrubland uses, because agriculture promoted the highest correlation in the rainfall-runoff relationships due to lower vegetation cover. Therefore, human-induced alterations related to land uses in the R generation must be assessed in those countries were the abandonment of upland agricultural catchments have allowed natural vegetation to expand (i.e., southern Europe) and also in countries were agriculture and irrigated areas are increasingly replacing forest and shrub lands (i.e., North African Maghreb). However, events under agricultural land use or which occurred in the winter season independently of the land use and lithology catchments were the situations which were able to generate a higher linearity in the rainfall-runoff relationships.
This temporal downscaling from the annual to event scale elucidated how different R mechanisms can co-exist in small Mediterranean-climate catchments, considering the main temporal and spatial factors that govern the river catchment connectivity. Despite this, controls on R generation in Mediterranean-climate catchments require further attention to assess the role of lithology, land use and seasonality and their combined effects on the hydrological response for going beyond in the comprehension of highly sensitive areas to global change, such as the Mediterranean region. Considering that the performance gaps of hydrological processes in hydrological models due to the results of rainfall-runoff at daily scales are smoothed in small Mediterranean-climate catchments, the findings of this study improve the comprehension of hydrological processes as practical field knowledge needed to complete hydrological modelling at the event scale. Acknowledgments: Part of the meteorological data were provided by the Spanish Meteorological Agency (AEMET).

Conflicts of Interest:
The authors declare no conflict of interest.       Creek during the study period 2012-2017, further exposed in Table A3. Figure A3. Highest R event contribution during the study period  at Es Fangar Creek. The numbers located at the peak of each hydrograph indicate the ID of the events recorded in Es Fangar

Appendix A
Creek during the study period 2012-2017, further exposed in Table A3.