1. Introduction
Analyses of historical observations of climate variables clearly indicate the anthropogenic causes of climate change [
1,
2,
3]. Numerous studies at various spatial scales have thoroughly explored the impact of projected changes in the climate system on hydrology and water resources [
1,
3,
4,
5,
6]. After air temperature, fresh water is proving to be the second most reactive environmental variable to climate change [
1,
7,
8,
9,
10]. Despite the numerous existing studies and analyses of the temporal and spatial characteristics of precipitation and temperature, the timeliness of these studies is becoming increasingly appreciated appreciated [
11,
12,
13,
14,
15,
16,
17,
18], especially in the context of documenting the impact of human activities on regional climate factors [
19,
20,
21]. Increasing computational capabilities and comprehensive analyses based on long measurement series [
8,
22,
23] are attracting our attention, especially in the context of the search for the justification of theses on the polarisation of extreme climatic events [
1,
24,
25,
26,
27,
28,
29,
30].
The polarisation of precipitation events can be defined as a process in which the frequency and intensity of extreme precipitation events, such as downpours, hail, storms or floods, increases. Simultaneously, the "normal" precipitation condition is becoming less frequent. In the case of temperature phenomena, polarisation can be defined as a process in which the frequency and intensity of extreme temperatures, such as heat waves or extreme cold, increases. Again, the "normal" state of temperatures is becoming increasingly rare. Both of these phenomena are part of a broader process of climate change that affects atmospheric circulation and processes, leading to increased weather variability and a greater risk of extreme weather. Polarisation can be interpreted as a change in the nature of the occurrence of extreme events, where extreme events, specifically extreme temperatures and precipitation, are becoming more frequent, while "normal" conditions are becoming less frequent.
The increased interest in climatic factors such as precipitation and temperature [
13,
31,
32] is also due to the consequences of floods and periods of hydrological drought [
33]. Studies of the variability of precipitation and temperature over different periods [
5,
9,
14,
34,
35,
36] have made it possible to analyse statistical trends and assess the direction of these changes [
32,
37,
38,
39]. The results of these analyses provide information about the moments of trend change and the form of these changes. In the context of the importance of the variability of climatic factors for water resources, it is important to apply statistical techniques [
4,
40] that allow the detection of changes and trends [
2,
38], in among other sources, data relating to monthly precipitation totals and monthly average temperatures.
Climate change predicted by the IPCC (Intergovernmental Panel on Climate Change) occurs as a result of increases in greenhouse gas emissions, which are mainly caused by the industrialisation of the world, the burning of fossil fuels and changes in land use [
3,
10,
41,
42]. According to the Intergovernmental Panel on Climate Change [
43], both natural forces and human activities contribute to changes in climate patterns, including increases in land and ocean surface temperatures, changes in spatial and temporal precipitation patterns, increases in the frequency of extreme events, rising sea levels and increased El Niño [
3,
4,
31,
34]. The warming of the Earth is evident in data that show a 0.6°C increase in global temperature from 1901 to 2001 [
43]. To see changes in precipitation and temperature trends, anomalous values and statistical models are used to measure the timing of the trend change and its magnitude. The change in trends of various climate variables can be quantified using global atmospheric circulation models (GCMs) or, for example, using statistical models that identify the moment of change of the trend and its magnitude [
8].
Many studies have shown that climatic variability has a significant impact on river flow values, which are sensitive to both precipitation and temperature. The polarisation of climatic phenomena has been demonstrated by researchers [
35], who found increasing trends mainly in flood climaxes in German river catchments, which was attributed to climate change. In addition, seasonal analyses have shown that these changes were greater in winter than in summer. Rapid changes in both mean and variance can be attributed to both climate (e.g., changes in climate regimes) and anthropogenic effects (e.g., construction of dams and reservoir systems, changes in land use and vegetation cover, changes in the location of measurement points) [
44,
45]. Statistical analyses must be interpreted in the context of observed physical [
3,
7,
46], social and economic phenomena [
3,
14,
19,
47]. Therefore, predicting temporal trends in precipitation and temperature is very useful in social and urban planning [
48].
The world's rivers play a key role in linking the atmosphere, hydrosphere and biosphere by transporting 40% of precipitation runoff and 95% of sediment to the oceans [
49,
50]. The shaping of river flow depends on many factors, such as precipitation, temperature, evaporation, soil conditions, landforms and land use changes [
51]. With the available long-term precipitation and temperature data, it is very important to conduct long-term analyses of this data in order to understand the polarisation of climate phenomena and its potential consequences [
52]. In addition, such analyses provide valuable information that can be used to develop decision-making tools for the purposes of effectively managing extreme events, such as droughts and floods, and to mitigate ecological damage in regions particularly vulnerable to these hazards [
52].
When assessing climate polarisation, precipitation and temperature show strong relationships [
34,
52,
53]. Increased global temperatures can lead to increased water vapour in the atmosphere, which in turn can result in higher precipitation. Conversely, lower temperatures can result in reduced precipitation. In extreme cases, such as prolonged droughts or floods, fluctuations in temperature and precipitation can have catastrophic consequences, including the failure of crops and damage to property [
1,
45,
54]. Therefore, the accurate monitoring of precipitation and temperature is crucial for understanding climate change and assessing the polarisation of climate phenomena [
55].
Long-term data series make it possible to identify temperature and precipitation trends and help to determine whether a region is experiencing changes in polarity. Furthermore, the analysis of precipitation-temperature interactions provides a better understanding of important atmospheric phenomena, such as El Niño and La Niña cycles [
56], which affect global climate patterns. An example of this is the fact that rising temperatures can affect the occurrence of intense precipitation, while reduced precipitation can contribute to extreme heat waves [
45]. It is import to consider the amount of precipitation and temperature along with other climate phenomena as their interaction can affect the overall effect of climate polarisation [
57]. Understanding magnitudes and the relationship between precipitation and temperature is crucial for assessing the polarisation of climate phenomena and its consequences for the environment, economy and society. This knowledge is essential with regard to taking preventive measures in order to mitigate the negative effects of climate change.
The main topic of this study is the analysis of long-term sequences of precipitation and temperature data in terms of the possibility of introducing the concept of polarisation of climate extremes. Section one reviews the literature and presents the background of the variability of extreme climate events that provide the basis for analysing the polarisation phenomenon. Section two presents the temporal and spatial extent of the data that was used in the analysis. Section three characterizes the phenomenon of the polarisation of extreme climate events in the areas of precipitation and temperature, and discusses the impact of polarisation on important consequences for ecosystems, economies and infrastructure. The fourth section presents the approaches used in constructing a measure of polarisation and presents a proposal for the introduction of a static measure based on a stationary time series and a dynamic time series based on trends. The fifth and sixth sections present known statistical techniques in trend analysis and trend change point recognition. The seventh section discusses and evaluates the results of the performed analyses. The eighth section presents conclusions.
2. Preparation of Data for Analysis
The great interest in analysing long precipitation sequences stems from the need to assess climate change and its impacts on all spatial scales. This demand has led to the initiation and support of many research and monitoring programs conducted by international organisations. In this context, the Global Precipitation Climatology Centre (GPCC) was established on behalf of the World Meteorological Organization (WMO) in 1989. This institution is supported and operated by the Deutscher Wetterdienst (DWD, Germany's National Meteorological Service) as a German contribution to the World Climate Research Program (WCRP). The main objective of the GPCC is the global analysis of monthly precipitation on the Earth's surface based on "in situ" precipitation station data [
31,
58]. An equally important institution is the NOAA (National Oceanic and Atmospheric Administration) the US government agency for the study of the atmosphere, oceans and climate. It conducts research and monitoring of climate change, including collecting, processing, analysing, updating and making long-term precipitation and temperature sequences available [
59,
60].
The goal of the aforementioned institutions is to meet user requirements for the accuracy of the data and analysis results made available, as well as the timeliness and availability of the product. Among the various types of GPCC and NOAA products, gridded sets of long-term precipitation and temperature data are made available [
58,
59]. These items of data are not made available in real time. This paper relies on grid data of monthly precipitation totals from the GPCC products made available and grid data of monthly mean temperatures from NOAA products. This data corresponds to a spatial resolution of 0.5°x 0.5°, are consistent in spatial and temporal extent and cover the years 1901 to 2010. Products from both the GPCCC and NOAA are made available via the Internet.
The NOAA (National Oceanic and Atmospheric Administration) and the GPCC (Global Precipitation Climatology Centre) are two institutions that provide the opportunity, based on the available long-term data, to perform analysis and evaluation of the polarisation of climatic phenomena, precipitation and temperatures. The NOAA is engaged in the research and monitoring of the atmosphere, oceans and climate around the world. The NOAA collects data from various sources such as satellites, ocean buoys, aircraft and weather stations to assess the state of the climate and predict changes. The GPCC is involved in assessing precipitation around the world. The GPCC uses a variety of data sources, such as weather stations, radar and satellites, to develop global precipitation maps. The two institutions use their knowledge and data to assess changes in climate phenomena, including precipitation and temperatures around the world, and predict how they will change in the future. They are also working with other climate institutions to increase understanding of the climate and to help make decisions about climate change. This paper examines global trends in monthly precipitation totals and monthly mean temperatures from the area of 377 river basins distributed over all of the continents. Assuming 509.9 10
6 square kilometres of land area on the earth, the analysis covers 12.76% of the total land area of the globe.
Table 1 shows the areas covered by the analysis. The regions in the text are identified by code
where
(
Table 1,
Figure 1).
GPCC data relating to rainfall totals by month from 1901 to 2010, calculated from grid data with a spatial resolution of 0.5°x 0.5° longitude and latitude, were converted to catchment areas. In this way, a sequence of monthly precipitation was obtained, which became the subject of the analyses presented in this article. GIS interpolation mechanisms were used in the spatial analysis of data preparation. The calculated values of the sequences were subjected to a simple statistical analysis, determining the basic statistics: the minimum and maximum value, the mean, the standard deviation of the sample and the value of the coefficient of variation. The data was analysed in monthly cross sections as well as calendar years. The analyses of monthly precipitation covered the years 1901 to 2010.
NOAA monthly mean temperature data, from 1901 to 2010, calculated from grid data with a spatial resolution of 0.5°x 0.5° longitude and latitude, were converted to catchment areas. In this way, sequences of monthly average temperatures were obtained, which became the subject of the analyses presented in this article. GIS interpolation mechanisms were used in the spatial analysis of data preparation. The area of one of the smallest catchments in this study: SKJERN A (Europe, Denmark) is, according to calculations, 1090 km2. An area of 0.5°x 0.5° at the latitude of Denmark covers an area of 3100 km2. The location of this catchment area, for calculations of precipitation or temperature, requires the consideration of five neighbouring areas of 0.5°x 0.5°. With even such a small catchment area, precipitation and temperature characteristics are averaged. The calculated values of the sequences were subjected to a simple statistical analysis, determining the basic statistics: the minimum and maximum value, the mean and, the standard deviation of the sample. The data was analysed in monthly as well as calendar year cross sections. Analyses of temperature data covered the years 1901 to 2010.
4. The polarisation of precipitation and temperature phenomena
The occurrence of temporal variability and extreme events in air temperature and precipitation is usually assessed by analysing a set of indicators that define variability and extreme conditions [
39]. The common understanding of an extreme event is based on the assumption that a "normal" condition exists. In the context of the common understanding of an extreme event, polarisation can be interpreted as a change in the nature of the occurrence of extremes, with extreme events – extreme temperatures and precipitation – becoming more frequent and the "normal" state becoming less frequent. The polarisation of extreme events can result from a variety of factors, including human activity and climate change, which can affect precipitation cycles, the intensity of extreme events and temperature conditions. In this way, polarisation can be seen as a process in which climate characteristics change and extreme events become more frequent, while the "normal" state becomes increasingly unstable. Polarisation can have serious consequences for the functioning of ecosystems, the economy and the quality of human life.
Analysis of the polarisation of precipitation and temperature is essential because of their key role in global energy and water cycles and their impact on climate change. Accurate knowledge of precipitation and temperature is particularly important for assessing the amount of available freshwater and managing water resources, which is essential for reducing the risk of floods and droughts. In addition, there is growing scientific evidence that human activities are influencing climate change, contributing to shorter periods of intense precipitation and longer periods of high temperature and low precipitation. The polarisation of extreme events, such as floods and droughts, is increasingly evident and can be attributed to the unevenness and intensity of human activity. Therefore, studying the polarisation factors associated with monthly precipitation and monthly average temperatures is crucial to the understanding of climate change at the regional level and developing strategies to manage water resources and reduce the risk of extreme climate events.
This article analyses long-term sequences of precipitation and temperature to assess the polarisation of climate phenomena. The term "polarisation" in terms of precipitation or temperature is used in climate science to describe the process by which some areas of the planet become more extreme in terms of temperature and/or precipitation. This can lead to significant changes in the local environment, including changes in the distribution of plant and animal species, changes in weather patterns and changes in sea level [
61]. Long-term sequences of precipitation and temperature are an essential tool for studying climate polarisation. These sequences provide a detailed record of an area's climate over many years, allowing trends and patterns in the data to be identified. By analysing these sequences, it is possible to identify areas where the climate is becoming increasingly polarized and track changes over time. One of the key advantages of using long-term precipitation and temperature sequences is that they provide a high degree of accuracy and precision. The time series are created using advanced measurement techniques, such as data from satellites and weather stations, and are subject to rigorous quality control measures [
27,
58,
62]. As a result, they provide a reliable record of climate data in a given area, enabling precise observations and measurements.
The impact of polarisation can be observed in various areas of life on earth: in ecosystems, the economy and infrastructure, in changes in the size and location of water resources. The impact of polarisation on ecosystems in areas of extreme temperature and precipitation changes affects biodiversity and ecosystem functioning [
61]. The polarisation of climatic factors can lead to changes in the distribution of plant and animal species [
63,
64], which can have negative consequences for entire ecosystems. The consequences of polarisation for agriculture can affect crops and food production [
17,
36]. In some regions, climate polarisation may lead to a reduction in the amount of available water resources, which in turn affects agricultural production and the health of the industry. The consequences of polarisation for infrastructure can affect roads, bridges, buildings, water and sewage networks and water management systems. In extreme cases, severe damage or destruction can occur. Due to the complex and dynamic nature of the climate [
65,
66], polarisation is variable and difficult to predict. It is also difficult to predict exactly what the consequences of polarisation will be in a given region. In addition, climate variability means that some regions may experience polarisation in different years, complicating the process of predicting and monitoring changes. The interaction between polarisation and other climate phenomena can trigger certain climatic events, such as hurricanes, tornadoes or droughts. In this way, one phenomenon can amplify or weaken another, complicating the process of understanding the scale of climate change. The influence of anthropogenic factors on the polarisation of climate factors can occur through greenhouse gas emissions [
17,
42,
52,
57], among other factors.
Polarisation of climate events refers to changes in the intensity and frequency of influential weather events. NOAA and the GPCC monitor these changes [
58,
67] to determine the factors influencing the polarisation of climate phenomena, their effects, and ways to reduce the negative consequences associated with these events. Assessing temperature variability is crucial to understanding climate change and its impact on ecosystems and people [
68]. All of this information is essential for developing strategies and actions to manage change and to adapt and minimise its negative impacts. NOAA and GPCC are key sources of information for scientists, policymakers and other stakeholders who make decisions about climate change.
5. The concept of polarisation measure
The concept of measuring climate polarisation is based on measuring the degree of diversity and variability of climate elements in a given area. This measure is used to determine the level of polarisation of climate elements in the region of interest. Polarisation refers to the extremes of parameters such as temperature, precipitation, humidity, etc. occurring in a specific area. A high level of polarisation indicates significant differences between maximum and minimum values, which can lead to an increased risk of extreme weather events. The measure of polarisation can be determined from various climatic parameters, such as temperature, precipitation, humidity and atmospheric pressure. This can be calculated on different spatial and temporal scales, for example, for the entire country over the course of a year, for a specific region over the course of a month, or for individual weather stations over the course of a day. Depending on the purpose of the study and the availability of data, the measure of polarisation can be determined for different combinations of parameters. It is worth noting that a high level of polarisation does not necessarily mean unfavourable climatic conditions. For example, some regions have significant temperature differences between summer and winter, which can have a positive impact on seasonal tourism. However, in other cases, high polarisation can lead to serious problems, such as droughts, floods or storms. Therefore, a measure of polarisation can be useful in identifying areas that require special attention when planning climate change adaptation activities.
Different approaches are possible for constructing measures of polarisation in climate phenomena, depending on the type and scale of the data under analysis. One of the basic methods is to use probabilistic techniques based on classical theories of one- or multi-dimensional random variables with distributions built using copula functions [
48,
69,
70]. Other, simpler methods are based on statistical characteristics [
4,
10,
40,
71]. A simple concept is to construct indicators measuring the degree of extremity that a phenomenon can take, for example, by measuring the dispersion relative to the long-term average. An example of such an indicator is the coefficient of variation. Another way to create a polarisation indicator is to define it as the degree of extremity compared to the mean value relative to the variability measured by the dispersion measure [
12,
72,
73]. Such indicators can be based on classical measures, such as variance, standard deviation, mean deviation and the coefficient of variation, or on positional measures, such as range, quartile deviation and the coefficient of variation [
71].
Further methods of building indicators are based on the idea of unevenness [
74,
75], which can also be a measure of polarisation. One such indicator is kurtosis, which is influenced by the intensity of extreme values, so it measures what is happening in the "tails" of the distribution. Other types of indicators are measures that have their genesis in economics and econometrics. One such example is the Gini index, also called the Gini coefficient [
73,
76,
77,
78,
79], built using the Lorentz curve, which describes the degree of concentration of a one-dimensional distribution of a random variable with non-negative values [
77]. The Gini index ranges from a minimum value of zero, when all values are equal, to a theoretical maximum of one in an infinite population in which every element except one has a magnitude of zero. In the context of climate change, the Gini index can be used to measure inequality in exposure to the effects of climate change, such as droughts, floods or sea level rises [
57]. Higher values of the Gini index indicate greater inequality in the occurrence of these phenomena. Several other examples of measures can be cited, including the following:
The concentration ratio [
80] determines the degree of concentration of values at one end of the distribution and is similar to the Gini coefficient. However, it should be noted that the Gini coefficient may be less useful in analysing asymmetric distributions, which means that other indicators such as the concentration ratio or Lorenz curve should be considered in such cases.
The GMD (Gini Mean Difference) index [
81] is an inequality measure used in statistical and econometric analysis to measure polarisation or inequality in a sample distribution. Unlike the Gini coefficient, which measures unevenness, the GMD index enables the analysis of unevenness in the distribution of any variable, such as income, age, weight, height, precipitation and temperature. The GMD index ranges from zero to one, where zero indicates complete evenness in the distribution and one indicates the concentration of all values in one class. The higher the GMD index value, the greater the unevenness in the variable distribution.
The Theil index [
79,
82] is a measure of inequality in the distribution of quantitative variables, based on the idea of information entropy, taking into account differences between groups of values in the distribution, similar to the Gini coefficient, but with a greater emphasis on extreme values.
The Lorenz indicator [
74,
77] is an inequality indicator in a distribution, which is based on the Lorenz curve. It is often used to measure income inequality but can also be used to measure inequality in other quantitative variables, including climate change studies. Higher values of the Lorenz curve indicate greater inequality in the occurrence of climate change effects such as droughts, floods, or sea level rises, meaning that some regions or social groups are more vulnerable to the effects of climate change than others.
The Atkinson index [
73] is a measure of inequality in the distribution of quantitative variables, which is based on the idea of absolute deviations. It takes into account the differences between groups of values in the distribution; it is similar to the Gini index, but focuses more on average values than extreme values.
Range relates to values calculated as max-min usually refering to the difference between the maximum and minimum values of a given variable in a specific period of time. In the case of assessing the polarisation of precipitation and temperature, max-min can be used as a measure of the amplitude of these variables in a given period.
In this study, two measures are adopted to evaluate the phenomenon of polarisation, the first P1 is built based on a stationary time series, the second is based on calculated trends.
The first measure is defined as follows:
Where:
- the amplitude of change, range, the evaluated changes or ranges are taken over a selected period of time (110 years), the difference between the maximum and minimum value of a characteristic - is a measure that characterizes the empirical area of variation of the studied characteristic;
- standard deviation.
The amplitude of change values can be useful in identifying extreme climate conditions, such as periods of drought or heat waves, as well as periods of intense precipitation or extremely low temperatures. However, max-min as a single measure may be limited in its usefulness because it does not take into account other factors such as the length and intensity of the period, as well as other variables that affect weather conditions. Due to the fact that it ignores all data except for two extreme values, it does not provide information about the diversity of individual feature values in the population. Therefore, along with max-min values, other measures such as mean values, standard deviations, or cumulative indices are typically used to obtain a more comprehensive and diverse view of climate variability. For example, max-min values can be calculated for individual months, seasons or years, and then compared with mean values, standard deviations, or other measures to better understand climate variability over time and space.
The adopted measure, which is the ratio of the max-min values to the standard deviation , can be used as a measure to evaluate the polarisation of precipitation and temperature. The values refer to the amplitude of changes of a given variable over a period of time, while the standard deviation σ refers to the degree of variability of these values around their mean. By applying this measure, we can see how large the amplitude of changes is in relation to the variability around the mean. Values of this measure greater than 1 suggest that the amplitude of changes is greater than the variability around the mean. If the variability is relatively low compared to the amplitude, it may indicate the occurrence of periods of extreme climatic conditions, such as periods of drought or heatwaves, as well as periods of intense rainfall or extremely low temperature. However, the use of a single measure may be limited because it does not take into account other factors affecting climate variability. Therefore, it is worth using it together with other measures that show the character of climate variability. Note that the values of this measure for precipitation and temperature will always take a non-negative value.
The second form of adopted measures could be dynamic indicators showing the variability of the indicator over a longer period of time. The simplest indicators are based on the idea of a trend coefficient. Trend is a long-term feature referring to the general direction of changes over time in the value of a variable. Trend can indicate an increase, decrease, or stabilization of the variable over time. In time series analysis, the trend coefficient is one of the basic elements that help to forecast future values of the variable. The trend can be analysed for different time scales, from short-term cycles to long-term trends [
2]. Trend analysis is applied in many fields, such as economics, finance, meteorology, earth sciences, medicine, sociology, and marketing. In scientific research, trend analysis is often used to study changes in long-term time series, such as climate change, demographic changes, or evolution of species.
The second measure is defined as:
where
is the trend of the amplitude of changes, range or difference between the maximum and minimum value of the variable - it is a measure characterizing the empirical range of variability of the studied feature;
is the standard deviation to be evaluated are taken over a selected period of time (110 years).
The above measure can be interpreted as the ratio of the trend in variability of the amplitude to the trend in variability of the standard deviation. The trend in variability of the amplitude refers to the direction and speed of changes in the maximum and minimum values of a given variable over time, while the trend in variability of the standard deviation refers to the direction and speed of changes in the variability of the values of a given variable around their mean over time. Values of this measure greater than 1 suggest that the trend in the increase of amplitude changes is stronger than the trend in the increase of variability around the mean. This means that changes in the maximum and minimum values of a given variable are increasing faster than the variability around their mean, which may indicate the occurrence of extreme weather conditions. In the proposed measure, the trend of the ratio of the indicators was not calculated, but the trend of the numerator and denominator was separately maintained. Both measures are used to assess trends in data, but they differ in how they relate to data variability.
The measure: was not used because max-min changes are normalised by the standard deviation, meaning that high data changes in one direction can be offset by changes in the other direction.
The measure , by contrast, uses two indicators of data variability: trend amplitude and trend standard deviation. The use of amplitude in this formula means that more emphasis is placed on the extreme values of the data, rather than its overall variability. Therefore, when one is interested in more extreme values in the data, the measure adopted may be a better choice.
Note that the values of this measure for precipitation and temperature can be negative as well as positive. In the case of negative values, this means that the trends of polarisation factors (amplitude and variability) will be opposite, in the case of a positive value, the trends of polarisation factors will be consistent.
If we wanted a more balanced evaluation of the trend, taking into account the overall variability of the data, the measure might be more appropriate. Therefore, the proposed measure here was not applicable.
6. Detecting a change point in the trend
Various methods can be applied to determine change points in a time series [
8,
83,
84,
85,
86]. In this analysis, the non-parametric Pettitt change point test (PCPT) [
87] was used to detect the occurrence of changes. The Pettitt test (PCPT) is a non-parametric test for detecting sudden changes in a time sequence. It is used to detect the turning point where a sudden change, known as a "jump", occurs in the time series. The Pettitt test (PCPT ) involves comparing the sum of ranks of two subsets of data, which are divided by a threshold value, to determine whether there is a statistically significant change in the time sequence.
This test can be applied to analyse data with any distribution, and the test result is not dependent on the data having a normal distribution. The result of the Pettitt test (PCPT ) is a test statistic value, which is compared with the critical value for the level of significance to determine whether the null hypothesis of no sudden changes in the time sequence can be rejected. The rank is given after sorting and is dependent on the variable that gives the order of the records in the set.
The PCPT has been widely used to detect changes in observed climatic and hydrological time series [
8,
16,
88,
89]. The Pettitt test is also applicable to investigate an unknown change point by considering a sequence of random variables
, which have a change point at
. As a result
has a common distribution
, but
have a different distribution
, where
. The null hypothesis
: no change (but
); is tested against the alternative hypothesis
: change (lub
); using the non-parametric statistic
where:
for the downward shift and
or the upward shift [
86]. The confidence level associated with
lub
is approximately determined by:
When
is smaller than the specified confidence level (for example, in this study, 0.95 was adopted), the null hypothesis is rejected. The approximate probability of significance for the change point is defined as:
It is obvious that in the case of a significant change point, the series is segmented at the change point into two sub-time series.
The main aim of this study is to investigate the existence of change points in the time series of monthly sum precipitation and monthly average temperature characteristics. For time series showing a significant change point, the trend test is applied to partial series, and if the change point is not significant, the trend test is applied to the entire time series [
8].
7. Trend test
To examine the trend in a given time series, the Mann-Kendall test (MKT) can be applied. Originally, this test was used by Mann [
90], and later Kendall in 1975 derived the distribution of the test statistic [
86]. This test is independent of the type of distribution and we do not need to adopt any specific form of the data distribution function [
91]. This test was widely recommended by the World Meteorological Organization for public applications, and has been used in many scientific studies to assess trends in water resource data [
2,
8,
86]. Therefore, the MKT has been recognised as an excellent tool for trend detection by other scholars in similar applications. It should also be noted that the MKT test considers only the relative values of all elements in the series
. for analysis. The test statistic for the MKT test is given by:
where
are the consecutive values of the data, and n is the number of elements in the data. Under the null hypothesis of no trend, and assuming that the data are independent and identically distributed with a zero mean and a variance
denoted by
, which is calculated as (
- the repeated values in the analyzed sequence):
To test a hypothesis, the standard normal distribution is used, denoted as the test statistic
, which is defined as follows for the trend test:
In the two-sided test for trend, the null hypothesis is represented as: : there is no trend in the dataset, which will be rejected if the calculated statistic is greater than the critical value of this statistic obtained from the standard normal distribution table corresponding to the previously established level of significance. A positive value of indicates an increasing trend, and a negative value indicates a decreasing trend.
The magnitude of the trend is estimated using a non-parametric slope estimator based on the median proposed by Sen [
92] and extended by Hirsch [
93]. The slope estimator is given by:
Where , and is treated as the median of all possible pairs of combinations for the entire dataset.
9. Results and Discussion
In the present study, the following properties of long-term sequences of monthly precipitation totals and monthly mean temperatures were investigated for 377 catchments from the area of 6 WMO regions.
The scope of research and calculations included:
the determination of the values of long-term monthly strings over a period of 110 years;
the calculation of statistics relating to average values for each calendar month, the minimum value, the maximum value, the mean value and the standard deviation;
the determination of the values of these trends with the Mann Kendall test (MKT) being used to evaluate the trend values;
the examination of whether the long-term series showed change points using the Pettitt test (PCPT);
the examination of whether the sub-series (after the point of change up to 2010) has a significant trend in cases in which the long-term series showed points of change and to what extent this trend has changed.
Of all the catchments analysed in terms of sequences of monthly precipitation totals, only for one catchment of the LAGARFLJOT River located in Iceland (EUROPE) was it shown that four tests (MKT for RANGE, MKT for STD, PCPT for trend change and MKT for new trend) were simultaneously satisfied. This means that at the 5% significance level for this catchment, the measure trend was recognised, the year of the trend change was determined, and additionally the new trend was recognised. The RANGE trend changed not only in value, but also in direction, from negative to positive, from a value of (-0.136) to a value of (0.310) in 1953, and the STD trend changed from a value of (-0.533) to a value of (0.889) also in 1953. In the case described for the LAGARFLJOT catchment, both the values and directions of the RANGE and STD trends changed. The change in the direction of the RANGE trend from negative to positive means that in 1953, there was an increase in the maximum value of monthly precipitation in the study area, while the STD trend indicates that the variability of this phenomenon has also changed.
Based on the analysis of monthly precipitation data in the analysed catchments, it should be noted that for 93/377 catchments, i.e. in 25% of the analysed catchments, the trends of polarisation coefficients (RANGE and STD) in monthly precipitation totals were confirmed at the significance level of 5%, while in the remaining 75%, the trends were not confirmed.
Analysing the values (
Figure 2) of the polarisation coefficient
sequentially in the ten largest catchments ranging from 244,000 – 2,900,000 km
2 showed (
): NILE (Africa, Egypt): 4.26/0.71, YENISEI (Asia, Russian Federation): 4.25/0.64, JANGCY RIVER: 5.28/ 0.66, GANGES: 3.81/1.24, HUANG HE (Asia, China): 6.11/1.01, BRAHMAPUTRA (Asia, Bangladesh): 5.42/0.94, AMU DARYA (Asia, Uzbekistan): 5.10/0.83, EUPHRATES (Asia, Iraq): 5.93/0.84, SENEGAL (Africa, Senegal): 4.95/1.34, and URUGWAY (South America, Uruguay): 7.15/1.34. The highest value of this assumed polarisation index among the ninety-three catchments was obtained at 12.78 for the COPPER RIVER catchment (North America, United States), and the lowest was at 3.81 for the GANGES River catchment (Asia, India).
for monthly precipitation totals of the period 1901 to 2010.
Analysing the values (Figure 3) of the trends of polarisation factors, (trend(RANGE) and trend(STD)) were shown at a significance level of 5%. The following values were obtained from the next ten largest catchments ranging from 244,000 – 2,900,000 km2, namely NILE (Africa, Egypt): 2.28, YENISEI (Asia, Russian Fed.): 3.66, YANGTZE RIVER (Asia, China): 4.69, GANGES (Asia, India): 3.20, HUANG HE (Asia, China): 4.02, BRAHMAPUTRA (Asia, Bangladesh): 5.17, AMU DARYA (Asia, Uzbekistan): 3.31, EUPHRATES (Asia, Iraq): 3.81, SENEGAL (Africa, Senegal): 3.00, and URUGUAY (South America, Uruguay): 3.55. The highest value of this assumed polarisation index, among the ninety-three catchments, was obtained at 5.17 for the BRAHMAPUTRA (Asia, Bangladesh) river catchment, the lowest was at 2.28 for the NILE (Africa, Egypt): river catchment.
Figure 4 shows the spatial location of the analysed catchments, where significant trend change points in the
polarisation coefficients of monthly precipitation were recognised at the 5% significance level. For the catchments of the rivers of Africa, the years of changes in rainfall trends mainly occur in the late nineteen-sixties to the nineteen-seventies. For the catchments of the rivers of Asia, the years of changes in rainfall trends mainly occur in the range of the nineteen-fifties to the nineteen-sixties in the area of Russia, the nineteen-forties to the nineteen-sixties in the area of China, and from the nineteen-thirties to the nineteen-seventies in the area of India. For the catchment areas of South American rivers, years of changes in precipitation trends mainly occur in the period from the nineteen-forties to the nineteen-seventies. In the area of North America, this period is from the nineteen-forties to the nineteen sixties on the west coast and through the nineteen-sixties for the east coast. For the catchments of the rivers of the Australia-Oceania area, trend changes have been shown in the years from the nineteen-fifties to the nineteen-seventies.
Based on the analysis of the data of monthly average temperatures in the analysed catchments, it should be noted that for 46/377 catchments, i.e. in 12.2% of the analysed catchments, the trends of polarisation coefficients in the range of temperatures were proven at the significance level of 5%, and in the remaining 88.8, they were not confirmed.
Analysing the values (
Figure 5) of the trends of polarisation factors (
) sequentially in the ten largest catchments with an area of 190,000 – 1,730,000 km
2, the following values apply: AMUR: 3.42/-10.07, ORINOCO (South America, Venezuela): 6.10/0.04, GANGES (Asia, India): 3.74/0.25, INDUS (Asia, India): 3.34/0.47, BRAHMAPUTRA (Asia, Bangladesh): 3.39/0.63, SAO FRANCISCO (South America, Brazil): 5.24/0.07, VOLTA (Africa, Ghana): 4.67/0.06, RIO PARNAIBA (South America, Brazil): 5.65/0.04, GODAVARI (Asia, India): 4.29/0.15, URAL (Asia, Kazakhstan): 3.87/3.87. The highest value of this assumed polarisation index among the ninety-three catchments was obtained at 6.97 for the ESMERALDAS River (South America, Ecuador) catchment, and the lowest was at 6.10 for the ORINOCO (South America, Venezuela) River catchment.
Analysing the values (
Figure 6) of the trends of the
polarisation factors for the next ten largest catchments ranging from 190,000 – 1,730,000 km
2, showed: AMUR (Asia, Russian Fed.): 3.42, ORINOCO (South America, Venezuela): 6.10, GANGES (Asia, India): 3.74, INDUS (Asia, India): 3.34, BRAHMAPUTRA (Asia, Bangladesh): 3.39, SAO FRANCISCO (South America, Brazil): 5.24, VOLTA (Africa, Ghana): 4.67, RIO PARNAIBA South America, Brazil): 5.65, GODAVARI (Asia, India): 4.29, URAL (Asia, Kazakhstan): 3.87. The highest value of this assumed polarisation index among the ninety-three catchments was obtained at 6.97 for the ESMERALDAS (South America, Ecuador ) River catchment, the lowest was at 6.10 for the ORINOCO (South America, Venezuela) River catchment.
Figure 7 shows the spatial location of the analysed catchments, where significant trend change points in the
polarisation coefficients of monthly mean temperatures were recognised at the 5% significance level. For the catchments of African rivers, the years of changes in precipitation trends mainly occur in the nineteen-sixties to nineteen-eighties. For the catchments of Asian rivers, the years of changes in precipitation trends mainly occur in the nineteen-forties to the nineteen-sixties. For the catchments of South American rivers, years of changes in rainfall trends mainly occur in the nineteen-thirties to the nineteen-sixties. For the catchment area of North America, this period for the west coast is in the years following the start of the nineteen-eighties. For the catchments of the rivers of the Australia-Oceania area, trend changes have been identified in the nineteen-fifties.
9. Conclusions
This paper presents an analysis of monthly precipitation totals based on the GPCC database and monthly mean temperatures (NOAA data) for 377 catchments distributed around the world, representing 12.76% of the globe's total area. Data sequences covering 110 years from 1901 to 2010 calculated from grid data with a spatial resolution of 0.5°x0.5° longitude and latitude were analysed. The data was analysed in monthly periods of the calendar year. The study created and analysed 377 catchments x 110 years x 12 months = 497,640 precipitation data strings and the same number of temperature data strings, for a total of about 1 million long-term strings characterising climatic phenomena in terms of precipitation and temperature variability. Statistical characteristics of calendar months were calculated and the values of min, max, and standard deviation were determined. The indices of the adopted measures of polarisation were calculated on the basis of the coefficients of amplitude and standard deviation, as well as on the basis of the trends of these characteristics. Non-parametric MKT and PCPT trends were used in the analysis.
The polarisation index taking into account both precipitation and temperature phenomena is shown in
Figure 8. On the horizontal axis is the polarisation index
of monthly precipitation totals, while on the vertical axis is the polarisation index
of monthly average temperatures. The points depicting the polarisation phenomenon farthest from the center of the system indicate a high intensity of polarisation.
The situations of sign compatibility of trend(RANGE) and trend(STD) were used here. It was assumed that if the trends of P
2 measure factors are positive then there is a positive sign (+), if the trends of P
2 measure factors are negative then there is a negative sign (-). The concordance of the sign of trend(RANGE) and trend(STD) for monthly precipitation and temperature is due to the fact that in most cases, the variability of precipitation is correlated with changes in its maximum and minimum values. Similarly, in the area of temperature, the concordance of the sign of the trend(RANGE) and trend(STD) for monthly temperatures is due to the fact that in most cases, the variability of temperatures is correlated with changes in their maximum and minimum values. The characteristics of the river basins in which polarisation was found in both precipitation and temperature areas also included locations relating to Köppen-Geiger classification (K-G) areas [
94,
95].
The highest values are found in the VINCES (South America) catchment. This is located on the border of areas: As (equatorial, dry summer), Am (equatorial, monsoonal) and Cfc (warm temperature, fully humid and cool summer). It is characterised by positive values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates an increase in precipitation factor anomalies with a decrease in temperature anomalies. This means that the amount of precipitation may be increasing or there may be more variability in precipitation in the area. In contrast, temperatures are less variable or stable, with decreasing temperature anomalies. An increase in precipitation factor anomalies while temperature anomalies are decreasing indicates a certain inverse relationship between the two factors in the catchment. This means that with higher precipitation, there may be a tendency towards lower temperatures or less extreme temperature changes.
The catchment area of the MANGOKY River (Africa) is located on the island of Madagascar in the regions: Bsh (arid, summer dry, hot arid), Aw (equatorial, winter dry). It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors. The decrease in the anomalies of precipitation factors while the decrease in temperature anomalies indicates a certain compatible relationship between the two factors in the catchment. This means that as precipitation anomalies decrease and stabilise, there may be a tendency toward lower temperatures or less extreme temperature changes.
The YONGDING HE river catchment is located in Asia in the Bsk (arid, summer dry) and Dwb (snow, winter dry and warm summer) regions. It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors.
The catchment area of the BRAHMAPUTRA River is located in Asia in the regions: ET (polar, polar tundra) and Cwa (warm temperature, winter dry, hot summer). It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors.
The TAPTI RIVER catchment is located in Asia in the Aw (equatorial, monsoonal) and As (equatorial, summer dry) regions. It is characterised by positive values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates an increase in precipitation factor anomalies with a decrease in temperature anomalies.
The ATRATO river catchment is located in South America in the Af (equatorial, fully humid) region. It is characterised by positive values of both trends of precipitation polarisation factors and positive values of both trends of temperature polarisation factors. This indicates an increase in precipitation factor anomalies and temperature anomalies. If the region has positive values of both trends of precipitation polarisation factors and positive values of both trends of temperature polarisation factors, this indicates an increase in precipitation factor anomalies and temperature anomalies. Positive values of the precipitation polarisation factor trends suggest an increase in precipitation and greater variability of precipitation in the region. This could mean that precipitation becomes more abundant or more frequent and intense rainfall occurs. Positive trend values of temperature polarisation factors indicate an increase in temperature and greater temperature variability in the region.
The MIRA river catchment is located in South America in the ET (polar, polar tundra ), Cfb (warm temperature, fully humid, warm summer) region. It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors.
The ELLICE RIVER catchment is located in North America, on the border of regions: ET (polar, polar tundra) and Dfc (snow, fully humid, cool summer). It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors.
The catchment area of the ST. MAURICE river is located in North America, on the border of the Dfb (snow, fully humid, cool summer) and Dfc (snow, fully humid, warm summer) regions. It is characterised by positive values of both precipitation polarisation factor trends and decreasing values of temperature polarisation factor trends. This indicates an increase in precipitation factor anomalies with a decrease in temperature anomalies.
The PURARI River catchment is located in Australasia and Oceania, in the Af (equatorial, fully humid) region. It is characterised by negative values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates a decrease in the anomalies of precipitation and temperature factors.
The DE GREY RIVER catchment is located in Australasia and Oceania, in the Bwh (arid, winter dry, hot arid) region. It is characterised by positive values of both trends of precipitation polarisation factors and negative values of both trends of temperature polarisation factors. This indicates an increase in precipitation factor anomalies with a decrease in temperature anomalies.
Calming of precipitation and temperature anomalies (negative trends in both precipitation and temperature factors) are expected in the BRAHMAPUTRA and YONGDING HE catchments. For the PURARI (Australia and Oceania) and MANGOKY - Africa catchments, calculations indicate a calming of anomalies in both the temperature and precipitation variability area. Temperature and precipitation anomalies are to be expected in the ATRATO: South America catchments.
The analysis of K-G areas and catchments where trends in polarisation coefficients in both the area of monthly precipitation and monthly average temperatures have been recognised at the 5% significance level does not show associations of these catchment areas with the zones proposed by Köppen-Geiger [
94,
95].
In the calculations performed to determine the polarisation index, in addition to cases of trends with compatible signs (i.e. both positive and negative polarisation coefficients), there were also cases in which catchments showed opposite signs of polarisation coefficient trends in both precipitation and temperature sequence analysis. However, adopting a 5% significance level for the MKT test resulted in the rejection of these cases. The concordance of the sign of the trend (max-min) and trend (STD) for precipitation is due to the fact that in most cases, the variability of precipitation is correlated with changes in its maximum and minimum values. In other words, when there are periods of increased precipitation, we usually also observe higher maximum and minimum values of precipitation, and thus an increase in variability relative to the average value. Similarly, when there are periods of increased drought, there are usually lower maximum and minimum precipitation values, and variability relative to the average is also lower. Note, however, that there are also periods in which maximum and minimum precipitation values may increase or decrease, but variability relative to the average remains constant or changes in the opposite direction.
Analysis of the polarity of temperature and precipitation is necessary because of their impact on many aspects of the natural and human environment. Based on long-term sequences of precipitation and temperature, the article shows that the polarity process is present and can lead to significant changes in the local environment. Such changes can affect the distribution of plant and animal species, weather patterns, agriculture and the economy, and more. In addition, precipitation and temperature play a key role in global energy and water cycles, and their variability can lead to floods, droughts and other natural disasters. For this reason, knowledge of the nature of change and time scales is used to reduce the risk of extreme climate conditions, such as periods of drought or heat waves, as well as periods of intense precipitation and flooding or extremely cold temperatures.
A proper understanding of the polarisation of extreme events is also key to developing strategies related to mitigating and minimizing the impact of anthropogenic factors. This paper has presented opportunities to assess polarisation in precipitation and temperature variability, which can help develop such strategies. The following of climate changes suggest that polarisation is becoming more entrenched, and the associated extremes are becoming more intense and unevenly distributed. Therefore, analysing the polarisation of temperature and precipitation is important for understanding climate change and its impact on our environment as well as for developing effective strategies to manage risks and minimise human impact on the environment.
Funding
This research received no external funding.
Data Availability Statement
Data available upon request due to necessary commentary.
Conflicts of Interest
The author declares no conflict of interest.
References
- Romanowicz, R.J.; Bogdanowicz, E.; Debele, S.E.; Doroszkiewicz, J.; Hisdal, H.; Lawrence, D.; Meresa, H.K.; Napiórkowski, J.J.; Osuch, M.; Strupczewski, W.G.; et al. Climate Change Impact on Hydrological Extremes: Preliminary Results from the Polish-Norwegian Project. Acta Geophys. 2016, 64, 477–509. [Google Scholar] [CrossRef]
- Palaniswami, S.; Muthiah, K. Change Point Detection and Trend Analysis of Rainfall and Temperature Series over the Vellar River Basin. Pol. J. Environ. Stud. 2018, 27, 1673–1681. [Google Scholar] [CrossRef] [PubMed]
- Groves, D.G.; Yates, D.; Tebaldi, C. Developing and applying uncertain global climate change projections for regional water management planning. Water Resour. Res. 2008, 44, 1-16. [Google Scholar] [CrossRef]
- R. Katz, “Statistics of Extremes in Climatology and Hydrology,” Adv. Water Resour., vol. 25, pp. 1287–1304, 2002.
- Herschy, R.W. The world’s maximum observed floods. Flow Meas. Instrum. 2002, 13, 231–235. [Google Scholar] [CrossRef]
- Blöschl, G.; Bierkens, M.F.; Chambel, A.; Cudennec, C.; Destouni, G.; Fiori, A.; Kirchner, J.W.; McDonnell, J.J.; Savenije, H.H.; Sivapalan, M.; et al. Twenty-three unsolved problems in hydrology (UPH) – a community perspective. Hydrol. Sci. J. 2019, 64, 1141–1158. [Google Scholar] [CrossRef]
- Lewis, S.C.; King, A.D. Evolution of mean, variance and extremes in 21st century temperatures. Weather Clim. Extrem. 2017, 15, 1–10. [Google Scholar] [CrossRef]
- Jaiswal, R.K.; Lohani, A.K.; Tiwari, H.L. Statistical Analysis for Change Detection and Trend Assessment in Climatological Parameters. Environ. Process. 2015, 2, 729–749. [Google Scholar] [CrossRef]
- Heim, R.R. An overview of weather and climate extremes – Products and trends. Weather. Clim. Extremes 2015, 10, 1–9. [Google Scholar] [CrossRef]
- Sillmann, J.; Thorarinsdottir, T.; Keenlyside, N.; Schaller, N.; Alexander, L.V.; Hegerl, G.; Seneviratne, S.I.; Vautard, R.; Zhang, X.; Zwiers, F.W. Understanding, modeling and predicting weather and climate extremes: Challenges and opportunities. Weather. Clim. Extremes 2017, 18, 65–74. [Google Scholar] [CrossRef]
- Młyński, D.; Cebulska, M.; Wałęga, A. Trends, Variability, and Seasonality of Maximum Annual Daily Precipitation in the Upper Vistula Basin, Poland. Atmosphere 2018, 9, 313. [Google Scholar] [CrossRef]
- Młyński, D.; Wałęga, A.; Petroselli, A.; Tauro, F.; Cebulska, M. Estimating Maximum Daily Precipitation in the Upper Vistula Basin, Poland. Atmosphere 2019, 10, 43. [Google Scholar] [CrossRef]
- C. M. Twardosz R., “Temporal variability of maximum monthly precipitation totals in the Polish Western Carpathian Mts during the period 1951 – 2005,” pp. 123–134, 2012. [CrossRef]
- Ziernicka-Wojtaszek, A.; Kopcińska, J. Variation in Atmospheric Precipitation in Poland in the Years 2001–2018. Atmosphere 2020, 11, 794. [Google Scholar] [CrossRef]
- Venegas-Cordero, N.; Kundzewicz, Z.W.; Jamro, S.; Piniewski, M. Detection of trends in observed river floods in Poland. J. Hydrol. Reg. Stud. 2022, 41. [Google Scholar] [CrossRef]
- Z. W. Kundzewicz and M. Radziejewski, “Methodologies for trend detection,” IAHS-AISH Publ., no. 308, pp. 538–549, 2006.
- Z. W. Kundzewicz and A. Robson, “Detecting Trend and Other Changes in Hydrological Data,” World Clim. Program. - Water, no. May, p. 158, 2000, [Online]. Available online: http://water.usgs.gov/osw/wcp-water/detecting-trend.pdf.
- Berezowski, T.; Szcześniak, M.; Kardel, I.; Michałowski, R.; Okruszko, T.; Mezghani, A.; Piniewski, M. CPLFD-GDPT5: High-resolution gridded daily precipitation and temperature data set for two largest Polish river basins. Earth Syst. Sci. Data 2016, 8, 127–139. [Google Scholar] [CrossRef]
- B. Twaróg, “Characteristics of multi-annual variation of precipitation in areas particularly exposed to extreme phenomena. Part 1. the upper Vistula river basin,” E3S Web Conf., vol. 49, 2018.
- Yu, Y.; He, F.; Vavrus, S.J.; Johnson, A.; Wu, H.; Zhang, W.; Yin, Q.; Ge, J.; Deng, C.; Petraglia, M.D.; et al. Climatic factors and human population changes in Eurasia between the Last Glacial Maximum and the early Holocene. Glob. Planet. Chang. 2023, 221. [Google Scholar] [CrossRef]
- Chen, T.; Zou, L.; Xia, J.; Liu, H.; Wang, F. Decomposing the impacts of climate change and human activities on runoff changes in the Yangtze River Basin: Insights from regional differences and spatial correlations of multiple factors. J. Hydrol. 2022, 615. [Google Scholar] [CrossRef]
- Szolgayova, E.; Parajka, J.; Blöschl, G.; Bucher, C. Long term variability of the Danube River flow and its relation to precipitation and air temperature. J. Hydrol. 2014, 519, 871–880. [Google Scholar] [CrossRef]
- Pechlivanidis, I.G.; Olsson, J.; Bosshard, T.; Sharma, D.; Sharma, K. Multi-Basin Modelling of Future Hydrological Fluxes in the Indian Subcontinent. Water 2016, 8, 177. [Google Scholar] [CrossRef]
- Mudelsee, M.; Börngen, M.; Tetzlaff, G.; Grünewald, U. Extreme floods in central Europe over the past 500 years: Role of cyclone pathway “Zugstrasse Vb”. J. Geophys. Res. Atmos. 2004, 109. [Google Scholar] [CrossRef]
- Vavrus, S.J.; Notaro, M.; Lorenz, D.J. Interpreting climate model projections of extreme weather events. Weather. Clim. Extremes 2015, 10, 10–28. [Google Scholar] [CrossRef]
- Angélil, O.; Perkins-Kirkpatrick, S.; Alexander, L.V.; Stone, D.; Donat, M.G.; Wehner, M.; Shiogama, H.; Ciavarella, A.; Christidis, N. Comparing regional precipitation and temperature extremes in climate model and reanalysis products. Weather. Clim. Extremes 2016, 13, 35–43. [Google Scholar] [CrossRef] [PubMed]
- Michaelides, S.; Levizzani, V.; Anagnostou, E.; Bauer, P.; Kasparis, T.; Lane, J. Precipitation: Measurement, remote sensing, climatology and modeling. Atmospheric Res. 2009, 94, 512–533. [Google Scholar] [CrossRef]
- Das, N.; Bhattacharjee, R.; Choubey, A.; Ohri, A.; Dwivedi, S.; Gaur, S. Time series analysis of automated surface water extraction and thermal pattern variation over the Betwa river, India. Adv. Space Res. 2021, 68, 1761–1788. [Google Scholar] [CrossRef]
- Reinking, R. An approach to remote sensing and numerical modeling of orographic clouds and precipitation for climatic water resources assessment. Atmospheric Res. 1995, 35, 349–367. [Google Scholar] [CrossRef]
- López-Bermeo, C.; Montoya, R.D.; Caro-Lopera, F.J.; Díaz-García, J.A. Validation of the accuracy of the CHIRPS precipitation dataset at representing climate variability in a tropical mountainous region of South America. Phys. Chem. Earth, Parts A/B/C 2022, 127. [CrossRef]
- Becker, A.; Finger, P.; Meyer-Christoffer, A.; Rudolf, B.; Schamm, K.; Schneider, U.; Ziese, M. A description of the global land-surface precipitation data products of the Global Precipitation Climatology Centre with sample applications including centennial (trend) analysis from 1901–present. Earth Syst. Sci. Data 2013, 5, 71–99. [Google Scholar] [CrossRef]
- J. D. Gómez, J. D. J. D. Gómez, J. D. Etchevers, A. I. Monterroso, C. Gay, J. Campo, and M. Martínez, “Spatial estimation of mean temperature and precipitation in areas of scarce meteorological information,” Atmosfera, vol. 21, no. 1, pp. 35–56, 2008.
- Easterling, D.R.; Kunkel, K.E.; Wehner, M.F.; Sun, L. Detection and attribution of climate extremes in the observed record. Weather. Clim. Extrem. 2016, 11, 17–27. [Google Scholar] [CrossRef]
- Nobre, G.G.; Jongman, B.; Aerts, J.; Ward, P.J. The role of climate variability in extreme floods in Europe. Environ. Res. Lett. 2017, 12, 084012. [Google Scholar] [CrossRef]
- Petrow, T.; Merz, B. Trends in flood magnitude, frequency and seasonality in Germany in the period 1951–2002. J. Hydrol. 2009, 371, 129–141. [Google Scholar] [CrossRef]
- Singh, P.; Gupta, A.; Singh, M. Hydrological inferences from watershed analysis for water resource management using remote sensing and GIS techniques. Egypt. J. Remote. Sens. Space Sci. 2014, 17, 111–121. [Google Scholar] [CrossRef]
- Zhang, X.; Vincent, L.A.; Hogg, W.D.; Niitsoo, A. Temperature and precipitation trends in Canada during the 20th century. Atmos.-Ocean 2000, 38, 395–429. [CrossRef]
- N. N. Karmeshu Supervisor Frederick Scatena, “Trend Detection in Annual Temperature & Precipitation using the Mann Kendall Test – A Case Study to Assess Climate Change on Select States in the Northeastern United States,” Mausam, vol. 66, no. 1, pp. 1–6, 2015, [Online]. Available online: http://repository.upenn.edu/mes_capstones/47.
- R. Dankers and R. Hiederer, “Extreme Temperatures and Precipitation in Europe: Analysis of a High-Resolution Climate Change Scenario,” JRC Sci. Tech. Reports, p. 82, 2008.
- Katz, R.W.; Parlange, M.B. Overdispersion Phenomenon in Stochastic Modeling of Precipitation. J. Clim. 1998, 11, 591–601. [Google Scholar] [CrossRef]
- Roy, S.S.; Balling, R.C., Jr. Trends in extreme daily precipitation indices in India. Int. J. Climatol. 2004, 24, 457–466. [Google Scholar] [CrossRef]
- Colmet-Daage, A.; Sanchez-Gomez, E.; Ricci, S.; Llovel, C.; Estupina, V.B.; Quintana-Seguí, P.; Llasat, M.C.; Servat, E. Evaluation of uncertainties in mean and extreme precipitation under climate change for northwestern Mediterranean watersheds from high-resolution Med and Euro-CORDEX ensembles. Hydrol. Earth Syst. Sci. 2018, 22, 673–687. [Google Scholar] [CrossRef]
- G. Y. Lenny Bernstein, Peter Bosch, Osvaldo Canziani, Zhenlin Chen, Renate Christ, Ogunlade Davidson, William Hare, Saleemul Huq, David Karoly, Vladimir Kattsov, Zbigniew Kundzewicz, Jian Liu, Ulrike Lohmann, Martin Manning, Taroh Matsuno, Bettina Menne, Bert M, “Climate Change 2007 : An Assessment of the Intergovernmental Panel on Climate Change,” Change, vol. 446, no. November, pp. 12–17, 2007, [Online]. Available online: http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf.
- Chapman, J.; Eckley, I.A.; Killick, R. A nonparametric approach to detecting changes in variance in locally stationary time series. Environmetrics 2019, 31, e2576. [Google Scholar] [CrossRef]
- Thornton, P.K.; Ericksen, P.J.; Herrero, M.; Challinor, A.J. Climate variability and vulnerability to climate change: a review. Glob. Chang. Biol. 2014, 20, 3313–3328. [Google Scholar] [CrossRef] [PubMed]
- Swanson, K.L.; Tsonis, A.A. Has the climate recently shifted? . Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef]
- Balhane, S.; Driouech, F.; Chafki, O.; Manzanas, R.; Chehbouni, A.; Moufouma-Okia, W. Changes in mean and extreme temperature and precipitation events from different weighted multi-model ensembles over the northern half of Morocco. Clim. Dyn. 2022, 58, 389–404. [Google Scholar] [CrossRef]
- Mesbahzadeh, T.; Miglietta, M.M.; Mirakbari, M.; Sardoo, F.S.; Abdolhoseini, M. Joint Modeling of Precipitation and Temperature Using Copula Theory for Current and Future Prediction under Climate Change Scenarios in Arid Lands (Case Study, Kerman Province, Iran). Adv. Meteorol. 2019, 2019, 1–15. [Google Scholar] [CrossRef]
- Gerten, D.; Rost, S.; von Bloh, W.; Lucht, W. Causes of change in 20th century global river discharge. Geophys. Res. Lett. 2008, 35. 1-5. [Google Scholar] [CrossRef]
- Walling, D. Human impact on land–ocean sediment transfer by the world's rivers. Geomorphology 2006, 79, 192–216. [Google Scholar] [CrossRef]
- Hunter, T.S.; Clites, A.H.; Campbell, K.B.; Gronewold, A.D. Development and application of a North American Great Lakes hydrometeorological database — Part I: Precipitation, evaporation, runoff, and air temperature. J. Great Lakes Res. 2015, 41, 65–77. [Google Scholar] [CrossRef]
- Chai, Y.; Yue, Y.; Zhang, L.; Miao, C.; Borthwick, A.G.; Zhu, B.; Li, Y.; Dolman, A. Homogenization and polarization of the seasonal water discharge of global rivers in response to climatic and anthropogenic effects. Sci. Total. Environ. 2020, 709, 136062. [Google Scholar] [CrossRef] [PubMed]
- Li, Z.; Shi, Y.; Argiriou, A.A.; Ioannidis, P.; Mamara, A.; Yan, Z. A Comparative Analysis of Changes in Temperature and Precipitation Extremes since 1960 between China and Greece. Atmosphere 2022, 13, 1824. [Google Scholar] [CrossRef]
- Jahn, M. Economics of extreme weather events: Terminology and regional impact models. Weather. Clim. Extremes 2015, 10, 29–39. [Google Scholar] [CrossRef]
- Zhang, X.; Alexander, L.; Hegerl, G.C.; Jones, P.; Tank, A.K.; Peterson, T.C.; Trewin, B.; Zwiers, F.W. Indices for monitoring changes in extremes based on daily temperature and precipitation data. Wires Clim. Chang. 2011, 2, 851–870. [Google Scholar] [CrossRef]
- Allan, R.P.; Soden, B.J. Atmospheric Warming and the Amplification of Precipitation Extremes. Science 2008, 321, 1481–1484. [Google Scholar] [CrossRef]
- Change, I.P.O.C. Climate Phenomena and their Relevance for Future Regional Climate Change. 2013, 1217–1308. [CrossRef]
- B. Rudolf, C. B. Rudolf, C. Beck, J. Grieser, and U. Schneider, “Global Precipitation Analysis Products of the GPCC,” Internet Pbulication, pp. 1–8, 2005. Available online: ftp://ftp-anon.dwd.de/pub/data/gpcc/PDF/GPCC_intro_products_2008.pdf.
- T. C. P. and R. S. Vose, An overview of the global historical climatology network-daily database. Bull. Am. Meteorol. Soc. 1997, 78, 897–910. [CrossRef]
- Donat, M.G.; Alexander, L.V.; Yang, H.; Durre, I.; Vose, R.; Dunn, R.J.H.; Willett, K.M.; Aguilar, E.; Brunet, M.; Caesar, J.; et al. Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: The HadEX2 dataset. J. Geophys. Res. Atmos. 2013, 118, 2098–2118. [Google Scholar] [CrossRef]
- Persson, J.; Blennow, K.; Gonçalves, L.; Borys, A.; Dutcă, I.; Hynynen, J.; Janeczko, E.; Lyubenova, M.; Martel, S.; Merganic, J.; et al. No polarization–Expected Values of Climate Change Impacts among European Forest Professionals and Scientists. Sustainability 2020, 12, 2659. [Google Scholar] [CrossRef]
- Arheimer, B.; Pimentel, R.; Isberg, K.; Crochemore, L.; Andersson, J.C.M.; Hasan, A.; Pineda, L. Global catchment modelling using World-Wide HYPE (WWH), open data, and stepwise parameter estimation. Hydrol. Earth Syst. Sci. 2020, 24, 535–559. [Google Scholar] [CrossRef]
- Iverson, L.R.; McKenzie, D. Tree-species range shifts in a changing climate: detecting, modeling, assisting. Landsc. Ecol. 2013, 28, 879–889. [Google Scholar] [CrossRef]
- J. Franklin. Mapping Species Distributions: Spatial Inference and Prediction. Oryx 2010, 44, 615–615. [CrossRef]
- Viner, D.; Ekstrom, M.; Hulbert, M.; Warner, N.K.; Wreford, A.; Zommers, Z. Understanding the dynamic nature of risk in climate change assessments—A new starting point for discussion. Atmospheric Sci. Lett. 2020, 21. [Google Scholar] [CrossRef]
- Rosenzweig, C.; Parry, M.L. Potential impact of climate change on world food supply. Nature 1992, 367, 133–138. [Google Scholar] [CrossRef]
- Smith, T.M.; Reynolds, R.W.; Peterson, T.C.; Lawrimore, J. Improvements to NOAA’s Historical Merged Land–Ocean Surface Temperature Analysis (1880–2006). J. Clim. 2008, 21, 2283–2296. [Google Scholar] [CrossRef]
- Kosanic, A.; Anderson, K.; Harrison, S.; Turkington, T.; Bennie, J. Changes in the geographical distribution of plant species and climatic variables on the West Cornwall peninsula (South West UK). PLOS ONE 2018, 13, e0191021. [Google Scholar] [CrossRef]
- Chen, L.; Guo, S. Copulas and Its Application in Hydrology and Water Resources. In Copulas and its application in hydrology and water resources; Springer: Singapore, 2019. [Google Scholar] [CrossRef]
- Hou, W.; Yan, P.; Feng, G.; Zuo, D. A 3D Copula Method for the Impact and Risk Assessment of Drought Disaster and an Example Application. Front. Phys. 2021, 9. [Google Scholar] [CrossRef]
- S. M. Ross, Introduction to Probability and Statistics, no. 5. Academic Press is an imprint ofElsevier, 2014.
- Santos-Gómez, J.D.; Fontalvo-García, J.S.; Osorio, J.D.G. Validating the University of Delaware’s precipitation and temperature database for northern South America. DYNA 2015, 82, 86–95. [Google Scholar] [CrossRef]
- Amiel, Y.; Cowell, F. Thinking about Inequality. 1999. [CrossRef]
- Lorenz, M.O. Methods of Measuring the Concentration of Wealth. Publ. Am. Stat. Assoc. 1905, 9, 209–219. [Google Scholar] [CrossRef]
- M. G. and F. V. F. Francesco Nicolli, “Inequality and Climate Change: Two Problems, One Solution,” no. 2022, 2023, [Online] Available: https://wwwjstororg/stable/resrep44832.
- Ogwang, T. Calculating a Standard Error for the Gini Coefficient: Some Further Results: Reply. Oxf. Bull. Econ. Stat. 2004, 66, 435–437. [Google Scholar] [CrossRef]
- Damgaard, C.; Weiner, J. Describing Inequality in Plant Size or Fecundity. Ecology 2000, 81, 1139–1142. [Google Scholar] [CrossRef]
- T. Panek, “Polaryzacja ekonomiczna w Polsce,” Wiadomości Stat., vol. 1, no. 1, pp. 41–62, 2017.
- Sitthiyot, T.; Holasut, K. A simple method for measuring inequality. Palgrave Commun. 2020, 6, 1–9. [Google Scholar] [CrossRef]
- Samson, P.-M. Concentration of Measure Principle and Entropy-Inequalities. 2017, 55–105. [CrossRef]
- S. Yitzhaki, “Gini ’ s Mean Difference : A Superior Measure of Variability for Non-Normal Gini ’ s Mean difference : a superior measure of variability for non-normal distributions,” no. November, 2016.
- P. F. Pedro Conceição, “Young Person’s Guide to the Theil Index: Suggesting Intuitive Interpretations and Exploring Analytical Applications,” World, pp. 1–54, 2000, [Online]. Available online: https://hal.science/hal-01365314.
- Buishand, T. Some methods for testing the homogeneity of rainfall records. J. Hydrol. 1982, 58, 11–27. [Google Scholar] [CrossRef]
- Chen, J.; Gupta, A.K. Parametric Statistical Change Point Analysis. 2021, 951–952, 2013–2015. [CrossRef]
- Radziejewski, M.; Bardossy, A.; Kundzewicz, Z.W. Detection of change in river flow using phase randomization. Hydrol. Sci. J. 2000, 45, 547–558. [Google Scholar] [CrossRef]
- Salarijazi, M. Trend and change-point detection for the annual stream-flow series of the Karun River at the Ahvaz hydrometric station. Afr. J. Agric. Res. 2012, 7. [Google Scholar] [CrossRef]
- Pettitt, A.N. A Non-Parametric Approach to the Change-Point Problem. J. R. Stat. Soc. Ser. C (Appl. Stat.) 1979, 28, 126–135. [CrossRef]
- Verstraeten, G.; Poesen, J.; Demarée, G.; Salles, C. Long-term (105 years) variability in rain erosivity as derived from 10-min rainfall depth data for Ukkel (Brussels, Belgium): Implications for assessing soil erosion rates. J. Geophys. Res. Atmos. 2006, 111. [Google Scholar] [CrossRef]
- Conte, L.C.; Débora, M.B.; Fábio, M.B. Bootstrap Pettitt test for detecting change points in hydroclimatological data: Case study of Itaipu Hydroelectric Plant, Brazil. Hydrolog. Sci. J. 2019, 64, 1312–1326. [Google Scholar] [CrossRef]
- Mann, H.B. Nonparametric tests against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
- Yue, S.; Pilon, P.; Phinney, B.; Cavadias, G. The influence of autocorrelation on the ability to detect trend in hydrological series. Hydrol. Process. 2002, 16, 1807–1829. [Google Scholar] [CrossRef]
- Sen, P.K. Estimates of the regression coefficient based on Kendall’s Tau. J. Am. Stat. Assoc. 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
- Hirsch, R.M.; Slack, J.R.; Smith, R.A. Techniques of trend analysis for monthly water quality data. Water Resour. Res. 1982, 18, 107–121. [Google Scholar] [CrossRef]
- Beck, H.E.; Zimmermann, N.E.; McVicar, T.R.; Vergopolan, N.; Berg, A.; Wood, E.F. Present and future Köppen-Geiger climate classification maps at 1-km resolution. Sci. Data 2018, 5, 180214. [Google Scholar] [CrossRef] [PubMed]
- Kourtzanidis, K.; Angelakoglou, K.; Giourka, P.; Tsarchopoulos, P.; Nikolopoulos, N.; Ioannidis, D.; Kantorovitch, J.; Formiga, J.; Verbeek, K.; de Vries, M.; et al. Technical and Innovation Management Plans. POCITYF 2020. https://pocityf.eu/ Kottek, M.; Grieser, J.; Beck, C.; Rudolf, B.; Rubel, F. World map of the Köppen-Geiger climate classification updated. Meteorol. Z. 2006, 15, 259–263. [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).