Trends and correlations in annual extreme precipitation indices for

Theor Appl Climatol
DOI 10.1007/s00704-013-1079-6
ORIGINAL PAPER
Trends and correlations in annual extreme precipitation
indices for mainland Portugal, 1941–2007
M. Isabel P. de Lima & Fátima Espírito Santo &
Alexandre M. Ramos & Ricardo M. Trigo
Received: 7 March 2013 / Accepted: 19 December 2013
# Springer-Verlag Wien 2014
Abstract Precipitation extremes in mainland Portugal (southwestern Europe) using daily precipitation data recorded in the
period 1941–2007 (67 years) at 57 meteorological stations
scattered across the area are studied at an annual scale. Trends
in selected precipitation annual indices that are calculated
from these data are investigated, in particular trends in the
intensity, frequency and duration of extreme precipitation
events. Special attention is dedicated to local and regional
variability. The spatial correlations between the annual trends
in mean precipitation and in the extremes are analysed. Moreover, the relationships between the variability of the North
Atlantic Oscillation (NAO) index and several indices related
to the frequency and intensity of the precipitation at the 57
stations were also investigated. Results show that several
stations have predominantly negative tendencies in the precipitation indices, although the majority of stations did not
show statistically significant change over time in the 1941–
2007 period. At the regional level, the decreasing trend in the
M. I. P. de Lima (*)
Department of Civil Engineering, University of Coimbra,
Coimbra, Portugal
e-mail: [email protected]
M. I. P. de Lima
Institute of Marine Research, Marine and Environmental Research
Centre, Department of Civil Engineering, University of Coimbra,
Coimbra, Portugal
F. E. Santo
The Portuguese Sea and Atmosphere Institute—IPMA, I. P.,
Lisbon, Portugal
A. M. Ramos : R. M. Trigo
Instituto Dom Luiz, Universidade de Lisboa, Lisbon, Portugal
R. M. Trigo
Departamento de Engenharias, Universidade Lusófona,
1749 Lisbon, Portugal
simple daily precipitation intensity index is the only one
statistically significant at the 5 % level and appears to be
related to the predominance of the positive phase of the
NAO. For the period 1976–2007, the proportion of the total
precipitation attributed to heavy and very heavy precipitation
events increased and, consequently, daily precipitation events
show a tendency to become more intense. Moreover, correlation analysis show that the most extreme events could be
changing at a faster absolute rate in relation to the mean than
more moderate events.
1 Introduction
Over the twentieth century, changes in global and land precipitation have been observed across different time scales,
which are expected to result from the variability and change
in the climate. Particularly, changes in extreme precipitation
are of general concern because of the expected impact on
society and ecosystems; an extreme (weather or climate) event
is generally defined as the occurrence of a weather or climate
variable above (or below) a given threshold near the upper (or
lower) endpoints of the range of observed values of the
variable (e.g. IPCC 2012). Moreover, global climate variability and change are expected to be accompanied also by adjustments in other climate variables, which increase the complexity of the weather and climate systems. Water, soil and
energy are at the core of the discussions on this topic; they are
key factors for environmental and societal sustainability,
which depends much on the local conditions including resilience and adaptation capacities.
Insight into the properties of local and regional land precipitation in the recent past can be obtained by analysing
available ground-based point data; other alternative methods
(e.g. climate models) have often limited usefulness at small
time and space scales. However, precipitation measurement
M.I.P. de Lima et al.
limitations (e.g. data resolution, limited length of records) and
the sparse network of in situ ground-based precipitation monitoring contribute to the difficulties in detecting precipitation
low frequency fluctuations and change (both in space and
time), which is associated with large uncertainty. Findings
may also differ due to the different spatial and temporal scales
investigated and even the methodology used for each study.
Thus, the exploration of different scales and precipitation
parameters give opportunity to many studies to complement
previous analyses, towards a better understanding of those
systems and relevant processes.
Studies conducted at the global scale have usually been
consistently reporting increasing trends in extreme precipitation indices, which has also been found in some studies for
Europe, but at the regional level the findings may differ.
Despite general trends, the local precipitation structure, relevant generating mechanisms and other specificities can dictate
sometimes unexpected differences in precipitation characteristics and change over short distances that cannot be
disregarded. Examples of some of these studies are the following: for global scale, Frich et al. (2002), Alexander et al.
(2006) and IPCC (2007); for the countries of the western
Indian Ocean, Vincent et al. (2011); for China, Wang et al.
(2013) and Jiang et al. (2013); for Australia, Alexander et al.
(2007); for Central America and northern South America,
Aguilar et al. (2005); for southern and western Africa, New
et al. (2006); for central and southern Asia, Klein Tank et al.
(2006); at the European level, Klein Tank et al. (2002), Klein
Tank and Können (2003), Moberg and Jones (2005), Moberg
et al. (2006) and Karagiannidis et al. (2012); for the Mediterranean region and Iberian Peninsula, Xoplaki et al. (2004),
Norrant and Douguédroit (2006), Gallego et al. (2006),
Rodrigo and Trigo (2007), Rodrigo (2010), Gallego et al.
(2011), Acero et al. (2012), van den Besselaar et al. (2012)
and García-Barrón (2013); for mainland Portugal, Miranda
et al. (2002, 2006), de Lima et al. (2007, 2010a, b), Costa et al.
(2012), de Lima et al. (2013) and Espírito Santo et al. (2013b);
and particularly for the south of mainland Portugal, Costa and
Soares (2009), Durão et al. (2009) and Mourato et al. (2010)
and for the north, Santos and Fragoso (2013).
In general, no significant annual precipitation trends have
been reported for mainland Portugal, although some studies
found a decreasing tendency in extreme precipitation indices;
similar findings have been reported for the Mediterranean
region and Iberian Peninsula (see references, above). However, the majority of the studies that support these results have
relied on very few stations over western Iberia, where mainland Portugal is located; exceptions are the high density
station studies by Costa and Soares (2009) and Mourato
et al. (2010), but they are restricted to the southernmost region
of mainland Portugal. The strong precipitation spatial variability and gradients that are observed in the whole territory
requires that a higher number of locations is explored (e.g.
Trigo and DaCamara 2000; Miranda et al. 2002, 2006; BeloPereira et al. 2011). We note that some studies have examined
a restricted record period of only a few decades, which can
easily bias the estimation of trends (e.g. de Lima et al. 2010a, b).
Our study aims mainly to contribute to the increased understanding of precipitation in mainland Portugal by investigating daily precipitation data recorded at 57 meteorological
stations scattered across mainland Portugal. Trends in precipitation extremes in the period 1941–2007 (67 years) were
studied by testing time series of selected precipitation annual
indices over different multi-decadal periods: the indices were
calculated from daily data both at the local and regional levels.
These indices provide a view into the frequency of extreme
precipitation events and their intensity (see e.g. Frich et al.
2002; Moberg et al. 2006). Additionally, the correlations
between annual mean conditions and annual extreme precipitation indices, as well as correlations between trends in the
means and in the extremes were investigated. Moreover, the
correlation between precipitation indices and the most important pattern of large scale circulation in the northern hemisphere, i.e. the North Atlantic Oscillation (NAO), will be
analysed. With these analyses, we extended for Portugal the
work carried out by Klein Tank and Können (2003) for
Europe, thus providing a stronger basis for comparison of
trends in precipitation obtained for different locations and
European regions. Moreover, the analysis involves now more
indices, and longer and denser data sets than in previous
studies. These are the main contributions of our work, which
focuses on identifying trends in precipitation indices and
exploring whether such trends can be explained, in whole or
in part, by atmospheric circulation changes.
2 Study area and precipitation data
2.1 Brief description of the study area
Located in south-western Europe, mainland Portugal is confined between parallels 37°N and 42°N and within the relatively narrow meridional band that develops between 6.5 and
9.5°W. It lays in the transitional region between the subtropical anticyclone and the sub-polar depression zones. In
this territory, the latitude, orography and effect of the Atlantic
Ocean are the factors that dominate more the climate. Figure 1
shows a relief map of mainland Portugal; the highest altitude
is roughly 2,000 m.
Regional precipitation exhibits large inter-annual variability and very marked north–south and east–west gradients.
Mean annual precipitation varies from about 3,000 mm in
the north to roughly 500 mm in the southern part of the
country. The precipitation climate is also characterized by
strong seasonality: on average, about 40 % of the annual
precipitation falls in winter (December to February); summer
Annual extreme precipitation indices for mainland Portugal
1995/1996 (Santos et al. 2009), October–November 1997
(Carvalho 1997; Espírito Santo et al. 1998), autumn–winter
2000/2001 (Espírito Santo et al. 2001; Ventura et al. 2001;
WMO 2001) and winter 2009/2010 (Vicente-Serrano et al.
2011). In addition, on 18 February 2008, an active lowpressure system produced very intense precipitation that affected particularly Lisbon and the Beira Interior and Alentejo
regions (IM 2008). The station Lisbon/Geofísico recorded that
day a total of 118.4 mm, which surpassed the previous maxima of daily rainfall recorded at this station since 1864: 110.7
mm, on 5 December 1876, and 101.2 mm, on 30 January 2004
(Espírito Santo et al. 2004; IM 2008).
Long drought spells, which are more common in southern
Portugal, bring usually damages to agriculture and affect
water resources and availability for different uses. The decades of the 1940s, 1980s and 1990s are identified as being
particularly dry throughout the Iberian Peninsula territory (e.g.
Vicente-Serrano 2006a, b). Up to the present, the 2004/2005
episode is the worst drought event in the last 140 years and has
produced major socioeconomic impact not only in Portugal
but also in Spain (e.g. Espírito Santo et al. 2004; GarciaHerrera et al. 2007). Attention is drawn to the drought situation that was recorded in the last hydrological year 2011/2012,
which affected the whole territory and with greater severity in
the months of February and March.
2.3 Precipitation data set
Fig. 1 Relief map of mainland Portugal and location of the 57 climatological weather stations used in this study
(June to August) contributes to the total amounts with only 6 %,
approximately; in autumn and spring the amount of precipitation is highly variable (e.g. Miranda et al. 2002, 2006). The
small-scale precipitation variability has been less explored so
far and thus this process is still not fully characterized at small
spatial and temporal scales, in particular the extremes.
2.2 Records of precipitation extremes in the recent past
In mainland Portugal, the large temporal variability in precipitation leads to maximum precipitation extreme events (potentially leading to floods) and long drought spells, which are
relatively frequent. These types of events both affect different
regions, although the incidence and severity of the episodes
vary with location.
Major floods that occurred since the beginning of the
1980s, roughly in the last 30 years, were in December 1981
(Zêzere et al. 2005), November 1983 (catastrophic flooding in
the neighbouring Lisbon and Setúbal regions; e.g. Godinho
1984; Liberato et al. 2013), December 1983 (WMO 1984),
November and December 1989 (Zêzere et al. 2005), winter
The precipitation data cover the period 1941–2007 (67 years)
and were recorded at 57 climatological weather stations and
rain gauges from the networks of the Portuguese Sea and
Atmosphere Institute (Instituto Português do Mar e da
Atmosfera, IPMA, I.P.) and Portuguese Environment Agency
(Agência Portuguesa do Ambiente, APA). Data from 11 stations were provided by IPMA, while data from the other 46
stations were provided by the National Water Resources Information System (Sistema Nacional de Informação de
Recursos Hídricos), managed by APA.
The location of the stations selected in mainland Portugal is
shown in Fig. 1. These stations are representative, on the whole,
of the distribution of altitude of the territory; i.e. the percentage of
stations by classes of altitude is comparable to the area–altitude
distribution for mainland Portugal. Yet, there are differences in
the classes of lower altitude (i.e. the higher frequencies): more
stations in classes 0–100 m and 200–300 m and less in the class
100–200 m. We also note the lack of stations placed above
1,000 m: only one station is included in our data set.
2.4 Data processing and selection
The data were chosen from 200 precipitation series, which
were initially available; the selection was based on an exigent
combination of criteria related to the spatial distribution of the
M.I.P. de Lima et al.
series over mainland Portugal and data length, completeness,
quality and homogeneity:
(a) Completeness: Only stations with less than 2 % of missing values were used; any given month was considered
complete if no more than 3 days were missing from the
records; a year was considered complete if no more than
15 days were missing;
(b) Quality: Basic quality controls have been applied to all
the series with the purpose to identify errors in the data;
the daily data were searched for anomalous values (“outliers”), which were defined as having values that deviated more than 3 standard deviations from the climatologic
daily average, but no such extremes were identified;
(c) Homogeneity: Standard homogeneity tests (e.g. Wang
et al. 2007; Wang 2008a; Wang and Feng 2010) based
on the penalized t test (Wang et al. 2007), the penalized
maximal F test (Wang 2008b) and a Quantile Matching
algorithm (Wang 2009) were applied to the precipitation
data. The procedure used detected possible change points
in the monthly series that were checked against station
metadata records (when available). Series exhibiting evidence of discontinuities of non-climatic origin in the
study period were not used.
After the application of these tests, the number of stations’ data that were considered adequate to our study has
dropped to 57. The selection of the stations aimed also to
maximize the study period. The major limitations are not
only the data scarcity before 1941 (many series have started
being recorded at this time) but also and more importantly
the effective reduction of climate observations after 2007, as
a result of the reduction of stations in the monitoring network and the problems faced in the maintenance of the
existing ones; this resulted in gaps and a high percentage
of missing values. That justifies the difficulty of extending
this study beyond that date, while maintaining such a dense
network. More details about the data selection, quality control and homogenization can be found in de Lima et al.
(2013) and Espírito Santo et al. (2013b); these studies have
also used this data set.
3 Methods
Team Climate Change Detection and Indices (Peterson et al.
2001) and reviewed by Zhang et al. (2011). Here, the indices
were calculated at the annual scale for each individual station
and also for the study region (i.e. mainland Portugal) as a
whole. The regional indices were assessed as simple averages
over all 57 stations’ individual indices (i.e. giving all the
stations equal weight) and were used to obtain insight into
the inter-annual variability of precipitation in mainland Portugal. Anomalies are for the 1961–1990 reference period. A wet
day is a day with an accumulated precipitation of at least
1.0 mm.
The more relevant selected indices are the following:
–
–
–
–
–
Spell indices, consecutive wet days (CWD) and consecutive dry days (CDD), defined as the maximum number
of consecutive days with precipitation above/below
1 mm, respectively; thus, the CWD is the maximum
length of wet spells, which could intensity flooding,
whereas the CDD index can assess the region’s vulnerability to drought;
Absolute (fixed) thresholds, defined as the number of
days on which a precipitation value falls above a fixed
threshold: number of heavy precipitation days ≥10 mm
(R10), number of very heavy precipitation days ≥20 mm
(R20) and number of extremely heavy precipitation days
≥25 mm (R25);
The maxima of multi-day rainfall events’ indices, such as
the maximum precipitation on 1 and 5 days, RX1D and
RX5D; the RX5D index can be used as an indicator of
flood-producing events because, on the space scales considered here, severe floods are generally not caused by a
single heavy thunderstorm event but, more likely, by
long-lasting heavy precipitation events that extend over
a region;
Percentile (non-fixed) thresholds: the 90th, 95th and
99th percentiles of precipitation on wet, very wet and
extremely wet days, respectively, R90p, R95p and
R99p indices;
The R95pTot index, which is the ratio R95p/PrecTot
(precipitation fraction due to very wet days), represents
the percentage of annual total wet-day precipitation due
to events with precipitation above the 95th percentile; it
can be used to analyse the possibility of having a relatively larger variation in extreme precipitation events than
in total amount (Groisman et al. 1999).
3.1 Climate indices
A total of 13 indices of precipitation extremes calculated from
daily data were selected in this work for exploring changes in
the intensity, frequency and proportion of extremes in total
precipitation, in mainland Portugal. The description of the
selected indices is given in Table 1; these indices were firstly
defined by the joint CCl/WCRP-CLIVAR/JCOMM Expert
The indices selected here have also been used in other studies
for different regions around the world: e.g. Frich et al. (2002),
Klein Tank and Können (2003), Aguilar et al. (2005), Alexander
et al. (2006), Herrera et al. (2010) and Jiang et al. (2013); for
regions in mainland Portugal, see e.g. Costa et al. (2008), Costa
and Soares (2009), Durão et al. (2009), Costa et al. (2012), de
Lima et al. (2013) and Espírito Santo et al. (2013b).
Annual extreme precipitation indices for mainland Portugal
Table 1 Definition of the precipitation indices used in this study
Index
Description
Definition
PrecTot (mm)
SDII (mm)
CDD (days)
CWD (days)
R10 (days)
R20 (days)
R25 (days)
RX1D (mm)
RX5D (mm)
R90p (mm)
R95p (mm)
R99p (mm)
Annual total wet-day precipitation
Simple daily intensity index
Consecutive dry days
Consecutive wet days
Heavy precipitation days
Very heavy precipitation days
Extremely heavy precipitation days
Highest precipitation amount in one-day period
Highest precipitation amount in 5 consecutive days
Precipitation on wet days
Precipitation on very wet days
Precipitation on extremely wet days
R95pTot (%)
Precipitation fraction due to R95p
Annual total precipitation from days ≥1 mm
Annual total precipitation divided by number of wet days (≥1 mm)
Maximum length of dry spell (RR <1 mm)
Maximum length of wet spell (RR ≥1 mm)
Number of days per year with RR ≥10 mm
Number of days per year with RR ≥20 mm
Number of days per year with RR ≥25 mm
Annual maximum precipitation on 1-day intervals
Annual maximum precipitation sums on 5-day intervals
Precipitation amount per year above a site-specific threshold value for wet,
very wet and extremely wet days, calculated as the 90th, 95th and 99th
percentile of the distribution of daily precipitation amounts on days with
1 mm or more precipitation in the 1961–1990 reference period
Percentage of annual total wet-day precipitation due to events with
precipitation >95th percentile
3.2 North Atlantic Oscillation index
The NAO is the most important mode of atmospheric variability over the North Atlantic Ocean that affects significantly
Atlantic weather patterns, particularly in Europe and the Mediterranean Basin (e.g. Trigo et al. 2002, 2004; Vicente-Serrano
and López-Moreno 2008). The NAO index measures the
strength of the zonal flow across the North Atlantic and
consists in the pressure difference between Iceland and the
Azores (e.g. Hurrell 1995; Serreze et al. 1997; Hurrell and
Deser 2010) or Lisbon or Gibraltar (Jones et al. 1997).
The NAO index can be interpreted in terms of a large-scale
meridional exchange of atmospheric mass (van Loon and
Rogers 1978) or as the oscillation of a large-scale anomalous
pressure (or geopotential) pattern (Wallace and Gutzler 1981);
it has been found to correlate with surface climate in most of
the European region (e.g. Hurrell 1995; Hurrell and van Loon
1997; Trigo et al. 2002; Vicente-Serrano and López-Moreno
2008), in particular Portugal (e.g. Ulbrich et al. 1999; Trigo
et al. 2004; Santos et al. 2009).
In this work, we analyse the relationship between the NAO
and precipitation in mainland Portugal, using the precipitation
indices (Table 1) calculated from the daily data described in
Section 2.3. The annual NAO index used here was computed
as the difference between the normalized sea level pressure in
Ponta Delgada (Azores archipelago) and Stykkisholmur/
Reykjavik (Iceland) and the data were provided by the Climate Analysis Section, NCAR, Boulder, USA (Hurrell 1995).
The relationship between precipitation and the NAO was
already reported for the Iberian Peninsula in other studies: for
monthly or seasonal precipitation series, e.g. Trigo et al.
(2004) and Lopez-Bustins et al. (2008); for daily precipitation
series, e.g. Goodess and Jones (2002), Gallego et al. (2005)
and Rodrigo and Trigo (2007). Albeit in a very limited number
of studies, this relationship was also analysed for Portugal: at
the monthly or seasonal scales, by e.g. Zhang et al. (1997),
Ulbrich et al. (1999), Trigo and DaCamara (2000) and Miranda et al. (2002, 2006); and using daily data, by Santos et al.
(2009).
3.3 Trends and correlations
Linear trends in the time series of indices of precipitation
extremes were estimated by ordinary least squares fitting
(OLS). The t test was used to assess the statistical significance
of the trends with different levels of probability (see e.g.
Peterson et al. 2001; Klein Tank and Können 2003; Klein
Tank et al. 2009; van den Besselaar et al. 2012). In using this
methodology, we anchor on other studies that have reported
that the magnitude of the trends estimated by using the OLS
method and a more robust non-parametric method are very
similar (see e.g. Önöz and Bayazit 2003; Moberg and Jones
2005). Önöz and Bayazit (2003) showed that the parametric t
test is slightly more powerful than the non-parametric Mann–
Kendall test when the probability distribution is normal; they
also showed that such relative power decreases with the
increase in skewness but that for moderately skewed distributions the t test is almost as powerful as the Mann–Kendall test.
This implies that the two tests can be used interchangeably in
practical applications, with identical results in most cases.
Nevertheless, Nicholls (2001) recommended that in addition
to the test results at the 5 % significance level, also the results
at the 25 % level should be presented, which we took into
consideration here. The definition of the precipitation indices
listed in Table 1 shows that we are not examining extremely
rare precipitation events, for which the computation of
M.I.P. de Lima et al.
significant trends could be a priori hampered by the small
sample sizes (e.g. Tebaldi et al. 2006).
In this work, the analyses were carried out for the full 67year record period (1941–2007) and for two consecutive subperiods of approximately the same length: 1945–1975 and
1976–2007. In addition to the precipitation trends calculated
for individual stations, trends were also computed for the
region as a whole using the regional indices (see Section 3.1),
which were obtained by averaging over all stations. It is
expected that this averaging procedure will not decrease the
trend but will reduce the effect of natural variability (e.g. Klein
Tank et al. 2009), which anticipates more robust results.
Analysis of correlations between precipitation indices was
also used here as an additional tool to inspect the data closely;
the analysis was undertaken for mean conditions’ indices
(wet-day total precipitation (PrecTot) and simple daily precipitation intensity (SDII)) and all the other indices listed in
Table 1, which were calculated site-by-site, for all 57 precipitation stations, and at the regional level. The correlations were
calculated using the linear Pearson (product moment) correlation coefficient r and their statistical significance was assessed
at the 5 % level. This type of analysis might help to identify
stations’ data with anomalous behaviour (e.g. Moberg et al.
2006), which could result, for example, from a few individual
highly erroneous daily values in the data that could corrupt the
statistics for the entire series. The linear trends in the PrecTot
index and the precipitation extremes’ indices, for all the 57
stations, were also analysed: we aimed at investigating how
the trends in the mean and extremes vary at each station and
how well they are correlated. Because not all the precipitation
indices have the same units (see Table 1), the linear trends for
each index were calculated as a percentage of the average.
The Pearson’s correlation coefficient was also used to
explore the relationship between the precipitation indices
and the NAO index. The statistical significance of the correlation coefficients was calculated using the Student’s t test
(two-tailed test of the Student t distribution).
corresponding 95 % confidence interval is highlighted by
the width of the shaded area at each dot.
In the next sections, we dedicate special attention to some
of the results obtained. If not stated otherwise, the statistical
significance of the results given below is at the 5 % level.
4.1 Indices of wet-day total precipitation (PrecTot) and simple
daily precipitation intensity (SDII)
The mapping of the trends in the SDII and PrecTot indices,
which are indicators of mean conditions, is shown in Fig. 3 for
the three periods investigated.
The regional average anomalies for the SDII and PrecTot
indices are plotted in Fig. 4 (top), which shows the considerable inter-annual variability in these indices. The lowest SDII
regional indices (<8.0 mm) occurred in the years 2004, 2005
and 2007, which were the driest on record, and these occurrences coincide temporally with the lowest PrecTot regional
indices (<500 mm). The highest values of SDII (>11 mm)
were in 1989, 1997 and 1947, and of PrecTot (>1,200 mm)
happened in 1963 and 1960.
Overall, these two regional indices show decreasing trends
but there is only one result statistically significant in the 1941–
2007 period: SDII, which decreases on average
−0.13 mm day−1 decade−1 (Table 3); more than 70 % of the
stations’ data show small negative trends in this index that are
statistically significant in 40 % of the cases (see also Table 2).
The regional average series of PrecTot indicates a decreasing
trend of −13.21 mm decade−1 in the full 67-year period and a
more marked decrease of −44.60 mm decade−1 in the 1976–
2007 sub-period (Table 3). In this sub-period, more than 90 %
of the stations show decreasing trends, which range from −54
to −136 mm decade−1; however, only data from 11 % of the
stations show a statistically significant decrease. An increasing tendency is found for the 1945–1975 sub-period, but none
of the results is significant.
4.2 Indices of precipitation extremes
4 Observed trends and correlations in precipitation
indices
The results of trend analyses of 13 precipitation indices’ time
series are summarized in Tables 2 and 3, for the 1941–2007,
1945–1975 and 1976–2007 periods. Table 2 gives the percentage of precipitation stations that had positive/negative
trends in annual precipitation indices calculated individually
for the 57 stations’ data and the corresponding percentage of
statistically significant results at the 5 % level. And Table 3
shows the trends in annual regional precipitation indices and
the corresponding 95 % confidence intervals.
Figure 2 provides an overview of the trends in the regional
indices in the 1941–2007 period; for each index, the
The marked inter-annual variability, and overall change, found
in the regional extreme precipitation indices’s time series is
shown in Figs. 4 and 5. The highest values of the regional
indices were in 1997 for RX1D, RX5D, R90p, R95p and
R99p, and in 1963 for R10 and R20. In general, the lowest
values of the regional extreme precipitation indices were in the
drought years of 1980, 2004, 2005 and 2007, and in 1971 for
RX1D and R99p.
Overall, the results show statistically non-significant decreasing trends in regional extreme precipitation indices for
the 1941–2007 and 1976–2007 periods; the exceptions are the
RX1D and R99p regional indices, which show a nonsignificant increase in this last 32-year period.
Annual extreme precipitation indices for mainland Portugal
Table 2 Percentages of the 57 precipitation stations with linear positive (+)
and negative (−) trends in annual precipitation indices, and the respective
number of statistically significant trends (Sig+/Sig−) at the 5 % level.
Period
1941–2007
1945–1975
1976–2007
Sig+
Sig−
+
−
Sig+
Sig−
+
−
Sig+
Sig−
+
−
Positive (+) indicates significant wetting trends and negative (−) indicates
significant drying trends; it is the opposite for the CDD index. The data are
for the 1941–2007, 1945–1975 and 1976–2007 periods
PrecTot
SDII
CDD
CWD
R10
R20
R25
RX1D
RX5D
R90p
R95p
R99p
R95pTot
0
16
21
79
0
2
77
23
0
11
7
93
2
40
26
74
2
21
42
58
0
16
25
75
0
18
19
81
4
4
28
72
2
4
44
56
2
11
23
77
4
0
60
40
0
5
9
91
0
16
16
84
0
4
72
28
0
11
7
93
0
16
28
72
0
7
63
37
2
12
18
82
2
11
33
67
0
7
53
47
2
7
18
82
5
5
53
47
0
12
46
54
4
0
47
53
0
14
33
67
0
14
49
51
0
5
42
58
2
11
37
63
0
7
51
49
0
7
21
79
7
7
37
63
0
11
54
46
2
5
32
68
4
0
51
49
0
9
44
56
4
2
53
47
5
7
40
60
0
14
39
61
0
4
49
51
4.2.1 Indices of maximum length of dry and wet periods
The inter-annual variability of the regional CDD and CWD
indices is also shown in Fig. 4 (bottom) and trends are given in
Table 3.
The regional CDD index shows a decreasing trend of −1.3
day decade−1 in the 1941–2007 period. Negative trends were
found in 80 % of the stations; yet only roughly 18 % of them
exhibited a statistically significant decrease. In the 1976–2007
period, the decreasing trend is −1.2 day decade−1, but the trend
is not statistically significant.
The regional CWD annual indices take the lowest values
(<6 days) in 1982, 1944 and 2005, which also correspond to
the largest absolute values of the negative anomalies (Fig. 4).
In the 1941–2007 and 1976–2007 periods, CWD shows a
decreasing tendency for the majority of stations (77 and
91 %, respectively), but the statistical significance is low; only
a very small percentage of stations have statistically significant
trends, respectively, 10 and 5 %. This percentage of statistically
significant results is about the limit (assuming a 5 % rule of
thumb) of what could be considered random exceedance of the
significance level (see e.g. Livezey and Chen 1983).
4.2.2 Daily precipitation absolute threshold indices
Regional indices’ time series indicate fewer days with precipitation above 10, 20 and 25 mm (R10, R20 and R25, respectively, heavy, very heavy and extremely heavy precipitation
Table 3 Trends per decade (with the 95 % confidence intervals in brackets) in the annual regional indices of precipitation for the 1941–2007, 1945–1975
and 1976–2007 periods
Index
Unit
1941–2007
1945–1975
1976–2007
PrecTot
SDII
CDD
CWD
mm
mm/day
days
days
−13.21 (−37.42 to 11.01)
−0.13 (−0.25 to −0.01)
−1.32 (−2.92 to 0.28)
−0.18 (−0.46 to 0.10)
26.19 (−53.83 to 106.21)
−0.10 (−0.45 to 0.25)
−2.34 (−7.78 to 3.11)
0.20 (−0.71 to 1.12)
−44.60 (−118.83 to 29.64)
−0.19 (−0.63 to 0.26)
−1.16 (−5.86 to 3.54)
−0.50 (−1.37 to 0.36)
R10
R20
R25
RX1D
RX5D
R90p
R95p
R99p
R95pTot
days
days
days
mm
mm
mm
mm
mm
%
−0.60 (−1.48 to 0.27)
−0.30 (−0.72 to 0.12)
−0.17 (−0.47 to 0.13)
−0.14 (−1.04 to 0.76)
−1.56 (−4.24 to 1.13)
−6.44 (−18.84 to 5.96)
−3.00 (−11.72 to 5.71)
−0.86 (−2.61 to −4.33)
−0.14 (−0.73 to 0.46)
0.70 (−2.22 to 3.63)
0.18 (−1.21 to 1.56)
0.14 (−0.83 to 1.12)
−0.69 (−3.35 to 1.97)
−0.75 (−8.72 to 7.22)
2.14 (−37.02 to 41.30)
0.33 (−26.38 to 27.05)
−2.96 (−12.75 to 6.82)
−0.71 (−2.54 to 1.13)
−1.69 (−4.36 to 0.98)
−0.71 (−2.04 to 0.62)
−0.43 (−1.40 to 0.53)
0.25 (−2.90 to 3.40)
−2.29 (−11.23 to 6.66)
−13.54 (−54.77 to 27.69)
−5.90 (−35.65 to 23.84)
1.17 (−11.09 to 13.44)
0.30 (−1.67 to 2.28)
Significance level at 5 % is presented in bold, at 25 % in italics
M.I.P. de Lima et al.
Fig. 2 Overview of the trends in annual regional precipitation indices in
the 1941–2007 period: averages over the 57 stations (black dots); for each
index, the 95 % confidence intervals are highlighted by the width of the
shaded area at each dot. The indices on the horizontal axes are defined in
Table 1; those that belong to the same class (e.g. R10, R20, R25) are
connected with a line
daily events) in 1941–2007 and 1976–2007, and an increasing
trend in these indices in 1945–1975; however, trends are not
statistically significant (Table 3). Overall, this trend pattern is
similar to the one in the regional PrecTot index. But for the
1941–2007 period, the decrease in R10 and R20 is significant
only at the 25 % level.
The frequency of heavy, very heavy and extremely heavy
precipitation daily events in the station’s data generally
followed the tendency of their respective regional average.
More than 65 % of the stations show a decreasing trend in the
R10, R20 and R25 indices, but significant decreasing trends
are found in only 10 to 16 % of them, in the 67-year record
period. This behaviour is more pronounced in the 1976–2007
period: between roughly 80 and 90 % of the stations show
decreasing trends in these threshold indices (up to
5 days decade−1), but only between 7 and 12 % of the stations
have statistically significant results.
The percentage change in precipitation amount and number
of wet days have opposite signs at some stations (i.e. negative/
positive; when one decreases the other increases).
4.2.4 Daily precipitation percentile threshold indices
4.2.3 Indices of extreme precipitation events of 1- and 5-day
durations
The regional RX1D and RX5D indices exhibit decreasing
trends in the 1941–2007 and 1945–1975 periods but the trends
are not statistically significant (Table 3). Contrastingly, the
trend in the RX1D regional index is positive in 1976–2007. In
this period, the trends in the RX1D index range between −5
and +7 mm decade−1 across mainland Portugal; only 3.5 % of
the stations show significant positive trends and none of the
stations show significant negative trends.
Figure 6 maps the stations’ trends in the RX5D index in the
1941–2007, 1945–1975 and 1976–2007 periods: between 50
and 65 % of the stations show decreasing trends for the three
periods; 14 % of the stations report a significant decreasing
trend for the 1941–2007 and 1945–1975 periods; for the period
1976–2007, the trends range from −16 to +10 mm decade−1 and
only 5 % of the stations show significant negative trends.
The stations’ percentile-based R90p, R95p and R99p indices
do not show consistent trends. Only a few stations’ trends
(about 10 %) are statistically significant.
In general, the regional average series of these indices show
negative trends in the 1941–2007 period that accompanies the
PrecToc trend. But the trends in the regional R99p and PrecTot
indices have different sign (i.e. positive, negative) in the
1945–1975 and 1976–2007 sub-periods. The trend in the
R99p index is positive (1.2 mm decade−1) in the 1976–2007
period; the contribution of the extremely wet days to the total
precipitation amount is between 11 and 127 mm. The R95p
index shows a negative trend (−3.0 mm decade−1) in the full
67-year record period, while the contribution of the precipitation due to very wet days (above the 95th percentile) to the
total precipitation varies between 73 and 350 mm (approximately 10 and 30 %); this index exhibits a slight increase
(0.33 mm decade−1) in the 1945–1975 period and a decrease
(−5.9 mm decade−1) in the 1976–2007 period. The R90p
index has a trend pattern similar to R95p, although more
pronounced with an increase (2.1 mm decade−1) in 1945–
1975 and a large decrease (−13.5 mm decade−1) in the subsequent period (1976–2007); on average, the contribution of wet
days (above the 90th percentile) to the total annual precipitation varies between 120 and 500 mm in the 1941–2007
period.
On average, the trend in the annual precipitation fraction
due to very wet days (R95pTot) is −0.14 % decade−1 in 1941–
2007, which is not statistically significant (see Table 3). However, in the 1976–2007 period, the R95pTot trend is increasing
(0.30 % decade−1), which is accompanied by a decreasing
trend in the annual precipitation; these opposite trends may
indicate that the very wet days are less affected than the other
wet days.
Figure 7 shows the station’s trends in the R95p and
R95pTot indices in 1941–2007, which are similar. Positive
trends in the R95pTot index are found for 20 % of the stations
and are accompanied by increases in the total annual precipitation amount, but the result is not statistically significant.
Annual extreme precipitation indices for mainland Portugal
Fig. 3 Trends in 1941–2007 (left), 1945–1975 (centre) and 1976–2007 (right) for SDII (top) and PrecTot (bottom). The dots are scaled according to the
magnitude of the trend: blue for increasing (wetting) trends and yellow for decreasing (drying) trends
Only two stations that exhibit a statistically significant decrease in the annual precipitation show a significant change in
R95pTot. This may indicate a disproportionate large change in
the extremes relative to the total amount, but this is only apparent in areas associated with wetting trends, not drying trends.
For the 1976–2007 period, around 40 % of the stations (the
majority of them located in the Alentejo region, south of Portugal) show statistically non-significant positive trends in the
R95pTot index, and this coincides with the observation of a
decrease in the PrecTot index.
4.3 Relationship between annual mean precipitation
and extreme precipitation indices
4.3.1 Correlations between extreme and mean precipitation
indices
The correlation between eight extreme precipitation indices
and indices of wet-day precipitation (PrecTot) and intensity
(SDII) were investigated for the 1941–2007 period; the eight
extreme precipitation indices include the three threshold
M.I.P. de Lima et al.
Fig. 4 Regional annual anomalies for the SDII, PrecTot, CDD and CWD indices (Table 1) for 1941–2007. Superimposed are the piecewise trends
calculated for the sub-periods 1945–1975 and 1976–2007
indices (R10, R20 and R25), the two absolute indices (RX1D
and RX5D) and the three percentile-based indices (R90p,
R95p and R99p) listed in Table 1. For all the stations’ data,
the indices are positively and statistically significantly correlated; the average correlations for all the stations are equal to
or above 0.5.
On average, the correlations with PrecTot are stronger
(above 0.8) for R10, R20 and R90p; only four stations have
correlations between PrecTot and R90p below 0.8. For the
R25 and R95p indices, the correlations are also, on average,
above 0.76, so fairly strong. But for the R99p index, about half
of the stations have correlation coefficients below 0.5, and the
Fig. 5 Regional annual anomalies for the extreme precipitation R10, RX5D, R90p and R99p indices in 1941–2007. Superimposed to the time series are
the piecewise trends calculated for the sub-periods 1945–1975 and 1976–2007
Annual extreme precipitation indices for mainland Portugal
Fig. 6 Trends in 1941–2007 (left), 1945–1975 (centre) and 1976–2007 (right) for the RX5D index, in millilitre per decade. The dots are scaled according
to the magnitude of the trend: blue for increasing (wetting) trends and yellow for decreasing (drying) trends
spread in correlation coefficients among all the stations is
larger than for the other two percentiles indices, particularly
R90p.
The correlations between SDII and the R25, RX1D, RX5D
and R99p indices are on average stronger than the correlations
between PrecTot and those same indices. The average correlations between SDII and R10, R20 and R25 are strong (r>0.7).
For some stations located in the Alentejo region (low altitude),
the correlations are very close to 1 (r>0.95), and increase when
going from R10 to R25; it is for these stations that these correlations are stronger. For SDII and RX1D, only eight stations have
correlations coefficients above 0.8, and these stations are also
located in the southern region of mainland Portugal. The correlations between SDII and the percentile indices decrease the
nearer we get to the tail of the precipitation distribution, particularly for the precipitation on extremely wet days (R99p); this
pattern is similarly followed by the correlations between the
PrecTot index and the percentile indices.
4.3.2 Correlations between trends in precipitation extremes
and trends in mean precipitation
Here, we investigate the linear relationship between percent
trends (i.e. percentage of change over time) in the PrecTot
index and extreme precipitation indices, for all stations
individually. Table 4 and Fig. 8 show that percent trends
in precipitation extremes are, in general, highly spatially
correlated with the percent trend in PrecTot; moreover, all
the correlations are statistically significant. This is particularly true for the number of heavy to extremely heavy
precipitation days (R10 to R25) and precipitation on wet
and very wet days (R90p and R95p). Except for CDD,
CWD and R10, the slope of the line of best fit for all other
indices is above 1, indicating enhanced variability in the
trends of precipitation extremes compared to the mean
trends. The largest change (i.e. slope) is found for RX1D
and this result is statistically significant.
Table 4 Spatial correlations between annual percent trends in extreme precipitation indices and percent trends in annual total wet-day precipitation for
the 57 stations across mainland Portugal in the 1941–2007 period
Index
CDD
CWD
R10
R20
R25
RX1D
RX5D
R90P
R95P
R99P
s
r
R2
−0.3
−0.31
0.10
0.5
0.48
0.23
0.7
0.85
0.72
1.0
0.87
0.75
1.4
0.87
0.81
4.9
0.57
0.32
1.1
0.78
0.60
1.4
0.90
0.81
2.0
0.88
0.77
2.8
0.58
0.34
Slopes of the line of best fit (s) and correlations (r) that are statistically significant at the 5 % level are in bold; the coefficients of determination (R2 ) are
also given
M.I.P. de Lima et al.
Fig. 7 Trends (top) and regional annual anomalies (bottom) for the R95p
(left) and R95pTot (right) indices, in 1941–2007. Top: the dots are scaled
according to the magnitude of the trend—blue for increasing (wetting)
trends and yellow for decreasing (drying) trends. Bottom: superimposed to
the time series are the piecewise trends calculated for the sub-periods
1945–1975 and 1976–2007
We should highlight here the statistically significant results
obtained for R95p and R99p (Fig. 8): even though the correlations between PrecTot and R99p are lower than between
PrecTot and R95p (see also Table 4), the largest slope of the
line of best fit estimated for R99p (Fig. 8f) compared to R95p
(Fig. 8e) might suggest that the most extreme events could be
changing at a faster absolute rate in relation to the mean
conditions than more moderate events. A multiple OLS regression was carried out to quantify the relationship between
trends in annual precipitation and extreme precipitation indices, and geographical coordinates across mainland Portugal.
For elevation, latitude and longitude, the correlations are
weak; so the OLS results are unsatisfactory, indicating that
these attributes represent less than 6 % of the variation in the
trends. Only for R95 and R90, they represent 10 and 13 % of
the variation in their respective trends.
4.4 Correlation between the NAO index and precipitation
indices
Most research for the Northern Hemisphere on the relationship between the NAO and precipitation has been focused on
winter season (e.g. Hurrel 1995; Trigo et al. 2002, 2004);
recent works have stressed the need to look into summer
Annual extreme precipitation indices for mainland Portugal
(a)
(b)
(c)
(d)
(e)
(f)
Fig. 8 Trends (in percent per decade) in extreme precipitation indices
a R20; b R25; c RX1D; d RX5D; e R95p and f R99p plotted against
trends (in percent per decade) in annual total wet-day precipitation
(PrecTot) for the 57 stations across mainland Portugal. Each dot
represents a station. The line of best fit was calculated using total
least-squares regression, with s being the slope of this line and R2 the
corresponding coefficient of determination (see also Table 4)
NAO impacts but they are not significant for Portuguese
precipitation (e.g. Bladé et al. 2012). Starting from the early
1940s to the mid-1970s, the winter NAO exhibited a slight
downward trend followed by a significant trend towards positively values during the following decades until the mid1990s (e.g. Hurrell and van Loon 1997); since 1995, the
winter NAO has been declining but with large inter-annual
variance. The summer NAO has been characterized by less
decadal variability, recording a continuous increasing tendency between the 1950s and late 1970s followed by a decrease
afterwards. The correlations between the NAO and the regional average of the precipitation indices were analysed over the
full period of the records and the two sub-periods 1945–1975
and 1976–2006; the temporal correlations for these two subperiods are shown in Fig. 9.
Only the CDD index shows positive correlations with the
NAO for the three periods, although the correlations are not
statistically significant. For the 1941–2007 period, there is a
statistically significant anti-correlation between the precipitation indices and the NAO, with correlation coefficients varying between −0.5 and −0.3 (not shown). However, overall the
correlation coefficients are stronger (except for the regional
CDD index) in 1945–1975, ranging between −0.6 and −0.36.
It was in this period that the NAO stayed in a negative phase
M.I.P. de Lima et al.
Fig. 9 Correlations between the
NAO and 12 precipitation
regional indices for the 1945–
1975 and 1976–2007 subperiods. The horizontal dashed
lines indicate the 5 % significance
level. The labels on the horizontal
axes are defined in the text and in
Table 1
and this condition yielded also the strongest anti-correlation
with the RX5D index (an indicator of flood prone conditions).
These results are consistent with the occurrence of consecutive low pressure systems and associated fronts, which are
characteristics of an intensified Atlantic storm track (e.g. Trigo
et al. 2004). But in the following sub-period (1976–2007),
which is associated with a positive phase of the NAO (after the
mid-1970s), the respective anti-correlation is weaker (and
non-significant) in comparison with the other period. In addition, the other correlations are overall weaker (−0.4 to −0.2) in
this period and only five indices (PrecTot, R90p and R10,
R20, R25) show significant correlations.
Overall, the major differences between the trends in the
annual regional indices of precipitation (Table 3) in the two
sub-periods are the observation of a wetter period until the
mid-1970s, which reflects the tendency triggered by the negative phase of the winter NAO in that time, and a drier period
since then that might reflect the dominance of a NAO positive
phase in the last decades. For the same periods, the percentages of the 57 precipitation stations with positive and negative
correlations between the NAO and the precipitation indices
are given in Table 5. With the exception of the CDD index, all
other indices are dominated by negative correlations for the
majority of stations.
In the three periods, the anti-correlation with the NAO is
observed for all the stations for the PrecTot, R10 and R20
indices. However, it is clearly during the 1941–2007 period
that the anti-correlation is more frequently observed, including
those correlations that are statistically significant; in this period, all the stations show negative correlations between the
NAO and the other precipitation indices except for the RX1D,
R99p and R95pTOT indices. When positive correlations between the NAO and the precipitation indices are observed, any
of the results is statistically significant, except for the CDD
index. In the entire 1941–2007 period, for all the other indices,
their correlation with the NAO is negative for 84 to 100 % of
the stations, depending on the index; the results are statistically significant in 12 to 95 % of the stations. The strongest
results are observed for the PrecTot and R10 indices in 1941–
2007: negative correlations are found for all the stations,
which are significant in 95 % of them. The positive correlations between the NAO and the CDD index are more abundant in the 1976–2007 period; they are found in 68 % of the
stations and are significant in 7 % of them.
The correlations between the NAO index and the indices
that describe the number of heavy (R10) to extremely heavy
(R25) rainy days are negative for all the stations in the entire
1941–2007 period and also during 1945–1975, but the
Table 5 Percentages of the 57 precipitation stations with positive (+) and negative (−) correlations between the NAO and the precipitation indices, and
the respective statistically significant correlations (Sig+/Sig−) at the 5 % level. The data are for the 1941–2007, 1945–1975 and 1976–2007 periods
Period
1941–2007
1945–1975
1976–2007
PrecTot
SDII
CDD
CWD
R10
R20
R25
RX1D
RX5D
R90p
R95p
R99p
R95pTot
Sig+
Sig−
+
−
Sig+
0
95
0
100
0
0
61
0
100
0
5
0
53
47
0
0
58
0
100
0
0
95
0
100
0
0
89
0
100
0
0
72
0
100
0
0
12
16
84
0
0
33
0
100
0
0
74
0
100
0
0
54
0
100
0
0
9
12
88
0
0
14
14
86
0
Sig−
+
−
Sig+
Sig–
+
−
84
0
100
0
32
0
100
47
4
96
0
18
7
93
0
54
46
7
0
68
32
44
0
100
0
5
21
79
81
0
100
0
39
0
100
65
0
100
0
37
0
100
54
0
100
0
30
4
96
9
25
75
0
5
25
75
42
7
93
0
4
19
81
60
0
100
0
23
4
96
46
5
95
0
9
5
95
12
16
84
0
4
25
75
12
12
88
0
7
19
81
Annual extreme precipitation indices for mainland Portugal
number of stations having statistically significant correlations
decrease as the precipitation threshold increases from 10 to
25 mm. Likewise, for all periods, the number of stations with a
statistically significant negative correlation for the percentile
indices decreases as the percentile increases from R90p to
R99p. The strongest result found for the RX5D index is in
1945–1975, a period when the NAO stayed in a negative
phase (already discussed above): 93 % of the stations show
negative correlations with NAO, which are statistically significant in 42 % of them. Overall, the results for the RX1D index
are less significant; the strongest results are observed for the
entire period 1941–2007: 84 % of the stations show negative
correlations with NAO, which are statistically significant in
only 12 % of them. So, overall, there is anti-correlation
between NAO and the indicators of more intense rainy events,
but this is less significant for the extreme events. In fact, the
movement of the Atlantic storm track expressed by the NAO
is clearly more related with phenomena of frontal origin (see
e.g. Miranda et al. 2002), whereas the more intense events are
usually associated with local convective phenomena.
Figure 10 shows maps of the spatial distribution of the
Pearson’s correlation between the NAO and selected daily
precipitation indices (R10, R90p, SDII, PrecTot) in the
1941–2007 period, which confirm the existence of coherent
regions that are revealed by correlations of the same order of
magnitudes; the maps were constructed by interpolating the
correlations obtained for all the 57 stations over mainland
Portugal. They illustrate that the northeast and southeast have
precipitation regimes and patterns that are less influenced by
the NAO mode, when compared to the majority of areas in
mainland Portugal. It is likely that in these regions precipitation is less modulated by mid-latitude low pressure systems.
5 Discussion
The annual regional indices of precipitation extremes that
were used to characterize extreme wet and intense precipitation events in mainland Portugal do not suggest any marked
pattern of change over the 1941–2007 period, but the trends in
all the indices are negative, although the majority of trends are
not statistically significant and only the SDII index is statistically significant at the 5 % level.
The site-by-site assessment of the data from the 57 climatological stations shows that, in general, there is a mix of
increasing and decreasing trends in precipitation extremes
across mainland Portugal, but there are more stations
exhibiting decreasing trends than increasing trends in the
1941–2007 period, and this decreasing tendency is even more
notorious in the 1976–2007 sub-period.
We found a somewhat noticeable reduction in the regional
annual total wet-day precipitation index in the 1941–2007
period, which indicates a drying trend over mainland Portugal;
however, this result is not statistically significant. Nevertheless, at the regional scale, average annual precipitation (about
900 mm) presents a reduction over the full 67-year period of
about 10 %. The site-by-site analysis shows that the total wetday precipitation reduction is statistically significant for 16 %
of the stations. Klein Tank and Können (2003) analysed six
precipitation series across mainland Portugal (1946–1999)
and did not found significant trends in the annual total wetday precipitation (PrecTot) and the precipitation fraction due
to very wet days (R95pTot); only the data from one station,
Porto, showed significant decreasing trends in these indices.
However, the different smaller period and the limited data
examined by the authors in that study hamper the comparison
of results.
The regional CDD index showed also a decreasing trend
(statistically significant at the 25 % level) in this period; this
reduction is observed in 81 % of the stations, which is statistically significant in 17 % of them. Similar tendency towards
less extended dry periods has been reported by, e.g. Kiktev
et al. (2003), Alexander et al. (2006) and Tebaldi et al. (2006).
Conversely, the CWD regional index shows a decreasing
trend of −0.18 day decade−1 (statistically significant at the
25 % level) in the 1941–2007 period, which is more negative
(−0.50 day decade−1) in the 1976–2007 period (over 90 % of
the stations exhibit decreasing trends, but only 5 % of them are
statistically significant).
Decreasing trends in the 1941–2007 period were furthermore observed in RX1D and RX5D, and R10, R20 and R25
regional indices; however, only the trend in R10 (number of
days above 10 mm) is statistically significant at the 25 % level.
In addition, the results obtained for all the selected annual
extreme precipitation indices based on absolute precipitation
and relative (i.e. percentile) thresholds show that data from
only less than 16 % of the stations have statistically significant
decreasing trends; for RX1D, R95p and R99p less than 7 % of
the stations show significant increasing trends.
Despite all this, at the regional scale, the consistency of the
results for the different stations suggests less annual wet-day
precipitation along with fewer days with heavy precipitation.
Overall, these results hint a drying trend over mainland Portugal for the period 1941–2007. However, in this period, the
number of heavy and very heavy precipitation events (R10
and R20) decreased more significantly than the mean total
wet-day precipitation (PrecTot); so the proportion of the total
precipitation attributed to these events increased and, consequently, daily precipitation events tend to become more
intense.
The results also show that, on average, there is no statistically significant change in annual precipitation indices for
both sub-periods 1945–1975 and 1976–2007. But in general,
the precipitation indices detect a weak tendency towards more
extreme precipitation events in 1976–2007, particularly in the
southern region; this was shown by indicators of intensifying
M.I.P. de Lima et al.
Fig. 10 Correlations between the NAO and a R10, b R90p, c SDII and
d PrecTot indices across mainland Portugal in 1941–2007. Isolines are
for the Pearson correlation coefficients found for the 57 precipitation
stations; triangles indicate negative and statistically significant (5 %
level) correlations and crosses indicate non-significant correlations
Annual extreme precipitation indices for mainland Portugal
precipitation, such as RX1D and R99p. The percentage contribution of extremely wet days (upper 1 %) to the annual wetday precipitation is higher in 1976–2007 than in 1941–1975.
Correlation analysis suggests that precipitation extremes
are highly correlated with the mean precipitation and trends
in the precipitation extremes are highly correlated with
trends in mean precipitation, particularly for the R10, R20,
R25, R90 and R95 indices. We found evidence that the most
extreme events could be changing at a faster absolute rate in
relation to the mean conditions than more moderate events.
This result is consistent with Groisman et al. (1999) who
found a disproportionate change in precipitation intensities
whenever the mean precipitation changed; this was also shown
by Katz (1999). Furthermore, Tebaldi et al. (2006) point out that
there is a trend for intensified precipitation, with a greater frequency of heavy-precipitation and high-quantile events.
It is worth noting that, for the study area, other studies
identified a cooling tendency in the 1945–1975 sub-period
and a warming tendency in the 1976–2007 sub-period (see
e.g. Miranda et al. 2002, 2006; Ramos et al. 2011; de Lima
et al. 2013; Espírito Santo et al. 2013a). A similar overall
pattern was also found for Europe’s temperature (with the
warming period starting in the mid-1970s, see e.g. Klein Tank
et al. 2002; Klein Tank and Können 2003) and the global
average temperature (e.g. Rozelot and Lefebvre 2006). For
Portugal, the link between trends in air temperature and precipitation at the local and regional scales is outside the scope of this
work but should be investigated in future research. Also, the
correlation between the precipitation indices and geographic
attributes should be explored; Costa et al. (2010) and Durão
et al. (2010) have already analysed these relationships, but just
for southern Portugal and a limited number of indices and
record length. Nevertheless, we noticed that a larger percentage
of stations’ data from the central and southern regions of
Portugal showed positive trends in the extreme intensity, duration and percentile precipitation indices compared to the northern stations, for the entire 1941–2007 period and the last subperiod 1976–2007; but these trends are not statistically significant. This behaviour is consistent with the results reported for
the south by Costa and Soares (2009) and for the north by
Santos and Fragoso (2013), despite the different data and the
limited range of indices that were examined in those studies.
The correlation between the NAO index and precipitation
in western Iberian Peninsula highlights that precipitation in
this area, particularly mainland Portugal, is dominantly associated to mid-latitude low pressure systems and their fronts,
depending strongly on the exact location of the Atlantic storm
track (see e.g. Miranda et al. 2002). Strong positive phases of
the NAO tend to be associated with below-average precipitation over southern and central Europe whereas opposite
above-average precipitation anomalies are typically observed
during strong negative phases of the NAO (e.g. Trigo et al.
2002, 2004).
While, in general, spatially coherent regions of both
increasing and decreasing extreme precipitation in mainland Portugal did not emerge from the analysis of the
selected annual extreme indices, the result is different when
the precipitation is examined at sub-annual scales. We do
not focus on that behaviour here but other studies addressed
this issue: trends in monthly precipitation are discussed by,
e.g. de Lima et al. (2007, 2010a, b) based on the study of
long time series, some of them dating back to the nineteenth
century; Espírito Santo et al. (2013b) report a detailed
analysis of seasonal trends using a large number of indicators of daily precipitation extremes. According to these
studies, there are indications of the (re)distribution of precipitation during the year, with positive (increasing) trends
in some sub-annual periods (partially) offsetting negative
(decreasing) trends in other periods, leading in general to
non-significant trends at the annual scale.
6 Concluding remarks
The main purpose of the present study was to provide a more
comprehensive discussion of the precipitation structure and
changes in mainland Portugal, at the annual scale, using a
broader set of precipitation indices than in some previous studies.
Additionally, we aimed to use a higher number of stations, longer
records and to investigate relationships between the most important atmospheric circulation pattern and precipitation trends.
Analysis of a set of 13 annual precipitation indices, which
includes 11 extreme precipitation indices, derived from daily data
recorded at 57 stations distributed across mainland Portugal
identified several noticeable characteristics of precipitation in
the area in the 1941–2007 period. In general, the results found
in this study are in agreement with other studies that inspected
change in precipitation in western Iberia, where mainland Portugal lays (several references are given in Section 1); but we should
note that some of those studies were supported by monthly or
seasonal data whereas here we focussed on daily precipitation
data, which have been less explored.
The results of the site-by-site analysis highlight the
following:
–
–
–
Despite the presence of several stations with predominantly negative tendencies in the precipitation indices, the
majority of stations did not show statistically significant
change over time in the 1941–2007 period;
In the 1976–2007 period, there was a tendency towards
more extreme precipitation events, which was confirmed
by the highest daily precipitation amount (RX1D) and the
precipitation on extremely wet days (R99p), particularly
in the southern region of mainland Portugal;
Precipitation extremes are highly correlated with mean
precipitation and trends in the precipitation extremes are
M.I.P. de Lima et al.
–
highly correlated with trends in mean precipitation, particularly for the R10, R20, R25, R90 and R95 indices;
The most extreme precipitation events seem to be changing at a faster rate than the more moderate extreme events
and their intensity is increasing.
The results of the analysis of the regionally averaged indices show that:
–
–
–
–
–
There is an important but not statistically significant
decrease in regional average total precipitation;
The majority of precipitation regional indices, except
SDII, show small (non-significant) negative trends;
The regional averaged total wet-day precipitation has
declined, although this result is not statistically significant
at the 5 % level; however, since 1976 there is a large
reduction (−44 mm decade−1) in total wet-day precipitation, but this result is only statistically significant at the
25 % level; these drier conditions might reflect the predominance of warmer conditions in the last decades (e.g.
Klein Tank et al. 2005) but the relationship between
changes in temperature and precipitation needs to be
better understood;
Except for the CDD index, there is an overall significant
anti-correlation between the precipitation indices and the
NAO;
The decreasing trend found for SDII index may be related
to the predominance of the positive phase of the NAO
since the mid-1970s;
Although the average monthly and seasonal precipitation in
the Iberian Peninsula is highly influenced by the NAO (e.g.
Goodess and Jones 2002; Trigo et al. 2004), it depends also on
other modes of variability such as the Scandinavian and Eastern
Atlantic patterns; this was shown previously by, e.g. Trigo et al.
(2008) and Espírito Santo et al. (2013b). Moreover, at the daily
scale, smaller atmospheric circulation patterns play an important role favouring or damping the advection of moisture (Trigo
and DaCamara 2000) and may be better related (than simple
NAO indices) to the precipitation regime across the entire
Iberian Peninsula (e.g. Cortesi et al. 2013).
Acknowledgments The authors wish to thank Álvaro Silva and Sofia
Cunha (The Portuguese Sea and Atmosphere Institute, I. P.), for their help
in processing the maps in Figs. 1, 4, 6, 7 (top) and 10. Alexandre M.
Ramos was supported by the Portuguese Foundation for Science and
Technology (FCT) through grant FCT/DFRH/SFRH/BPD/84328/2012.
References
Acero FJ, Gallego MC, García JA (2012) Multi-day rainfall trends over
the Iberian Peninsula. Theor Appl Climatol 108:411–423. doi:10.
1007/s00704-011-0534-5
Aguilar E, Peterson TC, Ramírez Obando P, Frutos R, Retana JA, Solera
M, Soley J, González García I, Araujo RM, Rosa Santos A, Valle
VE, Brunet M, Aguilar L, Álvarez L, Bautista M, Castañón C,
Herrera L, Ruano E, Sinay JJ, Sánchez E, Hernández Oviedo GI,
Obed F, Salgado JE, Vázquez JL, Baca M, Gutiérrez M, Centella C,
Espinosa J, Martínez D, Olmedo B, Ojeda Espinoza CE, Núñez R,
Haylock M, Benavides H, Mayorga R (2005) Changes in precipitation and temperature extremes in Central America and northern
South America, 1961–2003. J Geophys Res 110, D23107. doi:10.
1029/2005JD006119
Alexander LV, Zhang X, Peterson TC, Caesar J, Gleason B, Klein Tank
AMG, Haylock M, Collins D, Trewin B, Rahimzadeh F, Tagipour
A, Rupa Kumar K, Revadekar J, Griffiths G, Vincent L, Stephenson
DB, Burn J, Aguilar E, Brunet M, Taylor M, New M, Zhai P,
Rusticucci M, Vazquez-Aguirre JL (2006) Global observed changes
in daily climate extremes of temperature and precipitation. J
Geophys Res 111. doi:10.1029/2005JD006290
Alexander LV, Hope P, Collins D, Trewin B, Lynch A, Nicholls N (2007)
Trends in Australia’s climate means and extremes: a global context.
Aust Meteorol Mag 56:1–18
Belo-Pereira M, Dutra E, Viterbo P (2011) Evaluation of global precipitation data sets over the Iberian Peninsula. J Geophys Res 116,
D20101. doi:10.1029/2010JD015481
Bladé I, Liebmann B, Fortuny D, van Oldenborgh GJ (2012) Observed
and simulated impacts of the summer NAO in Europe: implications
for projected drying in the Mediterranean region. Clim Dyn 39:709–
727. doi:10.1007/s00382-011-1195-x
Carvalho RC (1997) Nota sobre o temporal de 5/6 de Novembro de 1997.
IM - Instituto de Meteorologia, internal report. Lisboa, Portugal
Cortesi N, Trigo RM, Gonzalez-Hidalgo JC, Ramos AM (2013) High
resolution reconstruction of monthly precipitation of Iberian
Peninsula using circulation weather types. Hydrol Earth Syst Sci
17:665–678. doi:10.5194/hess-17-665-2013
Costa AC, Soares A (2009) Trends in extreme precipitation indices
derived from a daily rainfall database for the South of Portugal. Int
J Climatol 29(13):1956–1975. doi:10.1002/joc.1834
Costa AC, Durão R, Pereira MJ, Soares A (2008) Using stochastic
space-time models to map extreme precipitation in southern
Portugal. Nat Hazards Earth Syst Sci 8:763–773. doi:10.5194/
nhess-8-763-2008
Costa AC, Soares A, Pereira MJ, Durão R (2010) Evaluation of spacetime dynamics in extreme precipitation frequency using
geostatistical cosimulation with elevation. WSEAS Transactions
on Environment and Development 6:44–53
Costa AC, Santos JA, Pinto JG (2012) Climate change scenarios for
precipitation extremes in Portugal. Theor Appl Climatol 108:217–
234. doi:10.1007/s00704-011-0528-3
de Lima MIP, Marques ACP, de Lima JLMP, Coelho MFES (2007)
Precipitation trends in Mainland Portugal in the period 1941–
2000. IAHS publ no 310:94–102
de Lima MIP, Carvalho SCP, de Lima JLMP (2010a) Investigating annual
and monthly trends in precipitation structure: an overview across
Portugal. Nat Hazards Earth Syst Sci 10:2429–2440. doi:10.5194/
nhess-10-2429-2010
de Lima MIP, Carvalho SCP, de Lima JLMP, Coelho MFES (2010b)
Trends in precipitation: analysis of long annual and monthly time
series from mainland Portugal. Adv Geosci 25:155–160. doi:10.
5194/adgeo-25-155-2010
de Lima MIP, Espírito Santo F, Ramos AM, de Lima JLMP (2013)
Recent changes in daily precipitation and surface air temperature
extremes in mainland Portugal, in the period 1941–2007. Atmos Res
27:195–209. doi:10.1016/j.atmosres.2012.10.001
Durão R, Pereira MJ, Costa AC, Corte-Real JM, Soares A (2009) Indices
of precipitation extremes in Southern Portugal—a geostatistical
approach. Nat Hazards Earth Syst Sci 9:241–250. doi:10.5194/
nhess-9-241-2009
Annual extreme precipitation indices for mainland Portugal
Durão RM, Pereira MJ, Costa AC, Delgado J, del Barrio G, Soares A
(2010) Spatial–temporal dynamics of precipitation extremes in
southern Portugal: a geostatistical assessment study. Int J Climatol
30:1526–1537. doi:10.1002/joc.1999
Espírito Santo F et al. (1998) Caracterização climática do ano de 1997. IM
- Instituto de Meteorologia, internal report. Lisboa, Portugal
Espírito Santo F et al. (2001) Outono/Inverno de 2000/2001. Análise das
condições meteorológicas. IM - Instituto de Meteorologia, internal
report. Lisboa, Portugal
Espírito Santo F et al. (2004) Caracterização climática 2004. IM Instituto de Meteorologia, internal report. Ministério da Ciência,
Tecnologia e Ensino Superior, Lisboa, Portugal
Espírito Santo F, de Lima MIP, Ramos AM, Trigo RM (2013a) Trends in
seasonal surface air temperature in mainland Portugal, since 1941.
Int J Climatol. doi:10.1002/joc.3803
Espírito Santo F, Ramos AM, de Lima MIP, Trigo RM (2013b)
Seasonal changes in daily precipitation extremes in mainland
Portugal from 1941 to 2007. Reg Environ Chang. doi:10.
1007/s10113-013-0515-6
Frich P, Alexander LV, Della-Marta P, Gleason B, Haylock M, Klein Tank
AMG, Peterson T (2002) Observed coherent changes in climatic
extremes during the second half of the twentieth century. Clim Res
19:193–212
Gallego MC, García JA, Vaquero JM (2005) The NAO signal in daily
rainfall series over the Iberian Peninsula. Clim Res 29:103–109. doi:
10.3354/cr029103
Gallego MC, García JA, Vaquero JM, Mateos VL (2006) Changes in
frequency and intensity of daily precipitation over the Iberian
Peninsula. J Geophys Res 111:D24105. doi:10.1029/
2006JD007280
Gallego MC, Trigo RM, Vaquero JM, Brunet M, García JA, Sigró J,
Valente MA (2011) Trends in frequency indices of daily precipitation over the Iberian Peninsula during the last century. J Geophys
Res 116:D02109. doi:10.1029/2010JD014255
García-Barrón L, Morales J, Sousa A (2013) Characterisation of the intraannual rainfall and its evolution (1837–2010) in the southwest of the
Iberian Peninsula. Theor Appl Climatol. doi:10.1007/s00704-0130855-7
Garcia-Herrera R, Paredes D, Trigo RM, Trigo IF, Hernández H,
Barriopedro D, Mendes MT (2007) The outstanding 2004–
2005 drought in the Iberian Peninsula: associated atmospheric
circulation. J Hydrometeorol 8:483–498
Godinho SF (1984) Aspectos meteorológicos das inundações de 1983.
INMG – Instituto Nacional de Meteorologia e Geofísica, internal
report. Lisboa, Portugal
Goodess CM, Jones PD (2002) Links between circulation and changes in
the characteristics of Iberian rainfall. Int J Climatol 22(13):1593–
1615
Groisman PY, Karl TR, Easterling DR, Knight RW, Jamason PF,
Hennessy KJ, Suppiah R, Page CM, Wibig J, Fortuniak K,
Razuvaev VN, Douglas A, Forland E, Zhai P-M (1999) Changes
in the probability of heavy precipitation: important indicators of
climatic change. Clim Chang 43:243–283
Herrera S, Fita L, Fernández J, Gutiérrez JM (2010) Evaluation of the
mean and extreme precipitation regimes from the ENSEMBLES
regional climate multimodel simulations over Spain. J Geophys Res
115, D21117. doi:10.1029/2010JD013936
Hurrell JW (1995) Decadal trends in the North Atlantic Oscillation:
regional temperatures and precipitation. Science 269:676–679
Hurrell JW, Deser C (2010) North Atlantic climate variability: the role of
the North Atlantic Oscillation. J Mar Syst 79:231–244. doi:10.1016/
j.jmarsys.2009.11.002
Hurrell JW, van Loon H (1997) Decadal variations associated with the
North Atlantic Oscillation. Clim Chang 36:301–326
IM (Instituto de Meteorologia) (2008) Informação climática-relatório
mensal: Fevereiro 2008. IM, internal report, Março de 2008.
Ministério da Ciência, Tecnologia e Ensino Superior, Lisboa,
Portugal
IPCC (2007) Climate change 2007: the physical science basis.
Contribution of Working Group I to the Fourth Assessment Report
of the Intergovernmental Panel on Climate Change, Solomon S, Qin
D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller
HL (eds), Cambridge University Press, Cambridge, UK and New
York, NY, USA, 996 pp
IPCC (2012) Managing the risks of extreme events and disasters to
advance climate change adaptation. A Special Report of Working
Groups I and II of the Intergovernmental Panel on Climate Change
[Field CB, Barros V, Stocker TF, Qin D, Dokken DJ, Ebi KL,
Mastrandrea MD, Mach KJ, Plattner G-K, Allen SK, Tignor M,
Midgley PM (eds.)]. Cambridge University Press, Cambridge, UK,
and New York, NY, USA, 582 pp
Jiang F-Q, Hu R-J, Wang S-P, Zhang Y-W, Tong L (2013) Trends of
precipitation extremes during 1960–2008 in Xinjiang, the
Northwest China. Theor Appl Climatol 111(1–2):133–148. doi:10.
1007/s00704-012-0657-3
Jones PD, Jonsson T, Wheeler D (1997) Extension to the North Atlantic
Oscillation using early instrumental pressure observations from
Gibraltar and south-west Iceland. Int J Climatol 17:1433–1450
Karagiannidis AF, Karacostas T, Maheras P, Makrogiannis T (2012)
Climatological aspects of extreme precipitation in Europe, related
to mid-latitude cyclonic systems. Theor Appl Climatol 107:165–
174. doi:10.1007/s00704-011-0474-0
Katz RW (1999) Extreme value theory for precipitation: sensitivity analysis for climate change. Adv Water Resour 23:133–139
Kiktev D, Sexton DMH, Alexander L, Folland CK (2003) Comparison of
modeled and observed trends in indices of daily climate extremes. J
Clim 16:3560–3571
Klein Tank AMG, Können, GP (2003) Trends in indices of daily temperature and precipitation extremes in Europe, 1946–99. J Climate
16:3665–3680. doi:http://dx.doi.org/10.1175/1520-0442(2003)
016<3665:TIIODT>2.0.CO;2
Klein Tank AMG, Wijngaard JB, Können GP, Böhm R, Demarée G,
Gocheva A, Mileta M, Pashiardis S, Hejkrlik L, Kern-Hansen C,
Heino R, Bessemoulin P, Müller-Westermeier G, Tzanakou M, Szalai
S, Pálsdóttir T, Fitzgerald D, Rubin S, Capaldo M, Maugeri M, Leitass
A, Bukantis A, Aberfeld R, van Engelen AFV, Forland E, Mietus M,
Coelho F, Mares C, Razuvaev V, Nieplova E, Cegnar T, López JA,
Dahlström B, Moberg A, Kirchhofer W, Ceylan A, Pachaliuk O,
Alexander LV, Petrovic P (2002) Daily dataset of 20th-century surface
air temperature and precipitation series for the European Climate
Assessment. Int J Climatol 22:1441–1453. doi:10.1002/joc.773
Klein Tank AMG, Peterson TC, Quadir DA, Dorji S, Zou X, Tang H,
Santhosh K, Joshi UR, Jaswal AK, Kolli RK, Sikder AB,
Deshpande NR, Revadekar JV, Yeleuova K, Vandasheva S,
Faleyeva M, Gomboluudev P, Budhathoki P, Hussain A, Afzaal
M, Chandrapala L, Anvar H, Amanmurad D, Asanova VS, Jones
PD, New MG, Specktorman T (2006) Changes in daily temperature
and precipitation extremes in central and south Asia. J Geophys Res
111(D16105). doi:10.1029/2005JD006316
Klein Tank AMG, Zwiers FW, Zhang X (2009) Guideline on analysis of
extremes in a changing climate in support of informed decisions for
adaptation. Climate Data and Monitoring. WCDMP-No. 72, WMOTD No. 1500, Geneve, Switzerland
Liberato MRL, Ramos AM, Trigo RM, Trigo IF, Durán-Quesada AM,
Nieto R, Gimeno L (2013) Moisture sources and large-scale dynamics associated with a flash flood event. In: Lin J, Brunner D, Gerbig
C, Stohl A, Luhar A, Webley P (eds) Lagrangian modeling of the
atmosphere. American Geophysical Union, Washington, D. C. doi:
10.1029/2012GM001244
Livezey RE, Chen WY (1983) Statistical field significance and it’s
determination by Monte Carlo Techniques. Mon Weather Rev 111:
46–59
M.I.P. de Lima et al.
López-Bustins J-A, Martin-Vide J, Sanchez-Lorenzo A (2008) Iberia
winter rainfall trends based upon changes in teleconnection and
circulation patterns. Global Planet Chang 63:171–176
Miranda PMA, Coelho F, Tomé AR, Valente MA, Carvalho A, Pires
C, Pires HO, Cabrinha VP, Ramalho C (2002) 20th Century
Portuguese climate and climate scenarios. In: Santos FD,
Forbes K, Moita R (eds) Climate change in Portugal: scenarios,
impacts and adaptation measures. Gradiva, Lisbon, Portugal, pp
27–83
Miranda PMA, Valente MA, Tomé AR, Trigo R, Coelho MFES,
Aguiar A, Azevedo EB (2006) O Clima de Portugal nos
séculos XX e XXI. In: Santos FD, Miranda P (eds)
Alterações climáticas em Portugal. Cenários Impactos e
Medidas de Adaptação. Projecto SIAM II. Gradiva, Lisbon,
Portugal, pp 47–113
Moberg A, Jones PD (2005) Trends in indices for extremes in daily
temperature and precipitation in central and western Europe 1901–
1999. Int J Climatol 25:1149–1171
Moberg A, Jones PD, Lister D, Walther A et al. (2006) Indices for daily
temperature and precipitation extremes in Europe analyzed for the
period 1901–2000. J Geophys Res 111(D22106). doi:10.1029/
2006JD007103
Mourato S, Moreira M, Corte-Real J (2010) Interannual variability of
precipitation distribution patterns in Southern Portugal. Int J
Climatol 30:1784–1794. doi:10.1002/joc.2021
New M, Hewitson B, Stephenson DB, Tsiga A, Kruger A, Manhique A,
Gomez B, Coelho B, Caio AS, Masisi DN, Kululanga E,
Mbambalala E, Adesina F, Saleh H, Kanyanga J, Adosi J, Bulane
L, Fortunata L, Mdoka ML, Lajoie R (2006) Evidence of trends in
daily climate extremes over Southern and West Africa. J Geophys
Res 111:D14102. doi:10.1029/2005JD006289
Nicholls N (2001) The insignificance of significance testing. Bull Am
Meteorol Soc 82:981–986
Norrant C, Douguédroit A (2006) Monthly and daily precipitation trends
in the Mediterranean (1950–2000). Theor Appl Climatol 83(1–4):
89–106. doi:10.1007/s00704-005-0163-y
Önöz B, Bayazit M (2003) The power of statistical tests for trend
detection. Turk J Eng Environ Sci 27:247–251
Peterson TC, Folland C, Gruza G, Hogg W, Mokssit A, Plummer N
(2001) Report on the Activities of the Working Group on Climate
Change Detection and Related Rapporteurs 1998–2001. World
Meteorological Organization, WCDMP–No. 47/WMO–TD
No.1071: Geneva, Switzerland
Ramos A, Trigo RM, Santo FE (2011) Evolution of extreme temperatures
in Portugal: reporting on recent changes and future scenarios. Clim
Res 48:177–192. doi:10.3354/cr00934
Rodrigo FS (2010) Changes in the probability of extreme daily precipitation observed from 1951 to 2002 in the Iberian Peninsula. Int J
Climatol 30:1512–1525. doi:10.1002/joc.1987
Rodrigo FS, Trigo RM (2007) Trends in daily rainfall in the Iberian
Peninsula from 1951 to 2002. Int J Climatol 27:513–529. doi:10.
1002/joc.1409
Rozelot JP, Lefebvre S (2006) Is it possible to find a solar signature in the
current climatic warming? Phys Chem Earth 31(1–3):41–45. doi:10.
1016/j.pce.2005.01.005
Santos M, Fragoso M (2013) Precipitation variability in Northern
Portugal: data homogeneity assessment and trends in extreme precipitation indices. Atmos Res 131:34–45
Santos JA, Andrade C, Corte-Real J, Leite S (2009) The role of
large-scale eddies in the occurrence of winter precipitation deficits in Portugal. Int J Climatol 29(10):1493–1507. doi:10.1002/
joc.1818
Serreze MC, Carse F, Barry RG, Rogers JC (1997) Icelandic low cyclone
activity: climatological features, linkages with the NAO, and relationships with recent changes in the Northern Hemisphere circulation. J Clim 10:453–464
Tebaldi C, Hayhoe K, Arblaster JM, Meehl GA (2006) Going to the
extremes: an intercomparison of model-simulated historical and
future changes in extreme events. Clim Chang 79:185–211
Trigo RM, DaCamara CC (2000) Circulation weather types and their
influence on the precipitation regime in Portugal. Int J Climatol 20:
1559–1581
Trigo RM, Osborn TJ, Corte-Real J (2002) The North Atlantic Oscillation
influence on Europe: climate impacts and associated physical mechanisms. Clim Res 20:9–17
Trigo RM, Pozo-Vazquez D, Osborn TJ, Castro-Diez Y, Gámis-Fortis S,
Esteban-Parra MJ (2004) North Atlantic Oscillation influence on
precipitation, river flow and water resources in the Iberian
Peninsula. Int J Climatol 24:925–944. doi:10.1002/joc.1048
Trigo RM, Valente MA, Trigo IF, Miranda PM, Ramos AM, Paredes D,
García-Herrera R (2008) North Atlantic wind and cyclone trends
and their impact in the European precipitation and Atlantic significant wave height. Ann N Y Acad Sci 1146:212–234. doi:10.1196/
annals.1446.014
Ulbrich U, Christoph M, Pinto JG, Corte-Real J (1999) Dependance of
winter precipitation over Portugal on NAO and Baroclinic Wave
Activity. Int J Climatol 19:379–390
van den Besselaar EJM, Klein Tank AMG, Buishand TA (2012) Trends
in European precipitation extremes over 1951–2010. Int J Climatol.
doi:10.1002/joc.3619
van Loon H, Rogers JC (1978) The seesaw in winter temperatures
between Greenland and Northern Europe. Part I: general description.
Mon Weather Rev 106:296–310
Ventura JE, Espírito Santo F, Moita S, Silva A (2001) Le caractère
exceptionnel des précipitations d’hiver de 2000/2001 au Portugal
Continental: Caractérisation, Causes et Conséquences. XIV
Colloque International de Climatologie, Sevilla, Spain
Vicente-Serrano SM (2006a) Differences in spatial patterns of drought on
different time scales: an analysis of the Iberian Peninsula. Water
Resour Manag 20:37–60
Vicente-Serrano SM (2006b) Spatial and temporal analysis of
droughts in the Iberian Peninsula (1910–2000). Hydrol Sci J
51:83–97
Vicente-Serrano SM, López-Moreno JI (2008) Nonstationary influence
of the North Atlantic Oscillation on European precipitation. J
Geophys Res 113, D20120. doi:10.1029/2008JD010382
Vicente-Serrano SM, Trigo RM, López-Moreno JI, Liberato MLR,
Lorenzo-Lacruz J, Beguería S, Morán-Tejeda E, El Kenawy A
(2011) Extreme winter precipitation in the Iberian Peninsula in
2010: anomalies, driving mechanisms and future projections. Clim
Res 46:51–65
Vincent LA, Aguilar E, Saindou M, Hassane AF, Jumaux G, Roy D,
Booneeady P, Virasami R, Randriamarolaza LYA, Faniriantsoa FR,
Amelie V, Seeward H, Montfraix B (2011) Observed trends in
indices of daily and extreme temperature and precipitation for the
countries of the western Indian Ocean, 1961–2008. J Geophys Res
116, D10108. doi:10.1029/2010JD015303
Wallace JM, Gutzler DS (1981) Teleconnections in the geopotential
height field during the Northern Hemisphere winter. Mon Weather
R e v 1 0 9 : 7 8 4 – 8 1 2 . d o i :1 0 . 11 7 5 / 1 5 2 0 - 0 4 9 3 ( 1 9 8 1 )
109<0784:TITGHF>2.0.CO;2
Wang XL (2008a) Accounting for autocorrelation in detecting meanshifts in climate data series using the penalized maximal t or F test.
J Appl Meteorol Climatol 47:2423–2444. doi:10.1175/
2008JAMC1741.1
Wang XL (2008b) Penalized maximal F-test for detecting undocumented
mean-shifts without trend-change. J Atmos Ocean Technol 25(3):
368–384. doi:10.1175/2007/JTECHA982.1
Wang, XL (2009) A quantile matching adjustment algorithm for Gaussian
data series. Climate Research Division, Science and Technology
Branch, Environment Canada. 5 pp. http://cccma.seos.uvic.ca/
ETCCDMI/software.shttm
Annual extreme precipitation indices for mainland Portugal
Wang XL, Feng Y (2010) RHtestsV3 User Manual. Climate Research
Division, Science and Technology Branch, Environment Canada. 27
pp. http://cccma.seos.uvic.ca/ETCCDMI/RHtest/RHtestsV3_
UserManual.doc
Wang XL, Wen QH, Wu Y (2007) Penalized maximal t-test for detecting
undocumented mean change in climate data series. J Appl Meteorol
Climatol 46(6):916–931. doi:10.1175/JAM2504.1
Wang H, Chen Y, Chen Z (2013) Spatial distribution and temporal trends
of mean precipitation and extremes in the arid region, northwest of
China, during 1960-2010. Hydrol Process 27:1807–1818. doi:10.
1002/hyp.9339
WMO (World Meteorological Organization) Bulletin (1984) Unusual
weather in 1983—Part I. July 1984, 33(3):216–220
WMO RA–VI Bulletin (2001) Annual Bulletin on the Climate in WMO
Region VI–2001–Europe and Middle East, 58 pp
Xoplaki E, González-Rouco JF, Luterbacher J, Wanner H (2004) Wet
season Mediterranean precipitation variability: influence of largescale dynamics and trends. Clim Dyn 23:63–78. doi:10.1007/
s00382-004-0422-0
Zêzere JL, Trigo RM, Trigo IF (2005) Shallow and deep landslides
induced by rainfall in the Lisbon region (Portugal): assessment of
relationships with the North Atlantic Oscillation. Nat Hazards Earth
Syst Sci 5:331–344
Zhang X, Wang XL, Corte-Real J (1997) On the relationships between
daily circulation patterns and precipitation in Portugal. J Geophys
Res 12:2474–2489
Zhang X, Alexander L, Hegerl C, Jones P, Klein Tank A, Peterson C,
Trewin B, Zwiers Francis W (2011) Indices for monitoring changes
in extremes based on daily temperature and precipitation data.
WIREs Clim Change 2011. doi:10.1002/wcc.147