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 and temporal 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 [
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 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 climatic elements is a phenomenon in which certain climatic elements or features tend to separate into two extreme or opposite categories [
31,
32]. This means that over time, certain climatic features become more clearly concentrated in the two opposite extremes, which can lead to increased variability and diversity in atmospheric conditions. Polarisation can manifest itself as a clear separation between extremes of temperature, precipitation or other climatic parameters, with implications for long-term patterns and trends in the climate of an area. The property of bimodality emphasises the concentration of features in two opposite extremes [
8,
33,
34,
35]. The polarisation of precipitation events can be defined as a process in which the frequency and intensity of extreme precipitation events, such as rainstorms, 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,
36,
37] is also due to the consequences of floods and periods of hydrological drought [
38]. Studies of the variability of precipitation and temperature over different periods [
5,
9,
14,
39,
40,
41] have made it possible to analyse statistical trends and assess the signal of these changes [
37,
42,
43,
44]. The results of these analyses provide information about the moments of trend change and the form of them. In the context of the importance of the variability of climatic factors, it is important to apply statistical techniques [
4,
45] that allow the detection of changes and trends [
2,
43], 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,
46,
47]. According to the Intergovernmental Panel on Climate Change [
48], 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,
36,
39]. It is clearly emphasised that human influence significantly outweighs the action of natural forces. The warming of the Earth is evident in data that show a 0.6°C increase in global temperature from 1901 to 2001 [
48]. 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 [
40], who found increasing trends mainly in flood magnitude 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) [
49,
50]. Statistical analyses must be interpreted in the context of observed physical [
3,
7,
51], social and economic phenomena [
3,
14,
19,
52]. Therefore, predicting temporal trends in precipitation and temperature is very useful in social and urban planning [
53].
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 [
54,
55]. The shaping of river flow depends on many factors, such as precipitation, temperature, evaporation, soil conditions, landforms and land use changes [
56]. 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 [
57]. 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 [
57].
When assessing climate polarisation, a process in which certain climatic features become concentrated in two opposing categories, there are strong links between precipitation and temperature [
39,
57,
58]. An increase in global temperatures can lead to an increase in the amount of water vapour in the atmosphere. More water vapour in the atmosphere increases the potential for cloud formation and precipitation, especially in areas with higher temperatures where there is access to more thermal energy As a result, more water vapour may in turn increase precipitation, meaning that areas with higher temperatures may experience more 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,
50,
59].
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 link between important atmospheric phenomena, such as El Niño and La Niña cycles [
60,
61], 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 [
50]. 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 [
62]. 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.
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 [
53,
74,
75]. Other, simpler methods are based on statistical characteristics [
4,
10,
33,
45]. 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,
76,
77]. 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 [
33].
Further methods of building indicators are based on the idea of unevenness [
78,
79], which can also be a measure of polarisation. One such indicator is kurtosis [
80], 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 [
77,
81,
82,
83,
84], built using the Lorentz curve, which describes the degree of concentration of a one-dimensional distribution of a random variable with non-negative values [
82]. 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 [
62]. 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 [
85] 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 [
86] 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 [
84,
87] 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 [
78,
82] 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 [
77] 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 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 describe 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.
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. The main reason for these changes is the economic and demographic development in the area. In 1947, the town of Egilsstadir was founded in eastern Iceland on the banks of the Lagarfljót River and became the largest town in eastern Iceland and its main service, transport and administrative centre. It is now the largest settlement in the eastern region. In the 1950s, Iceland saw a surge in the use of oil as a heating medium. This resulted in increased CO
2 emissions, rising temperatures and may have caused changes in rainfall patterns [
99].
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 URUGUAY (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).
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 km
2, 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. Changes in trends can be linked to processes of decolonisation, independence and economic development in the area. The apogee of this process occurred between 1960 and 1968, and 1960 is referred to as the 'Year of African Independence' [
100]. 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. The 1950s in Russia saw intensive economic development, industrialisation and intensive agricultural development [
101]. 1949 was the beginning in China of land reform and nationalisation of the main economic sectors. The 1970s was the Cultural Revolution where economic development collapsed [
102]. 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. The 1960s and 1970s were a watershed period in the economic development of Latin American countries [
103]. 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. This is a period of industrial and demographic boom [
102]. 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. This is a period of intense industrial development and demographic growth, as well as the decolonisation of Oceania and its economic development [
102].
The reasons for the change in precipitation trends in the analysed catchments may be related to a variety of climatic, geographical and anthropogenic factors. Potential causes of changes in rainfall trends by region are outlined below:
Africa: changes in atmospheric circulation, including fluctuations in the rainfall belt and belt patterns, El Niño and La Niña phenomena, deforestation, overgrazing of animals and changes in land use, land degradation, [
32,
36]
Asia: diversity of topography can affect air mass movements and precipitation formation, changes in atmospheric circulation such as monsoons, increasing urbanisation and infrastructure development can create the so-called 'heat island effect', which affects the microclimate, air pollution including particulate matter which affects condensation and cloud formation [
36,
104],
South America: atmospheric circulation fluctuations, including those associated with El Niño and La Niña, topography of the region, including the Andes, deforestation, conversion of land for agriculture and urbanisation [
32],
North America: changes in atmospheric circulation such as the North Atlantic Oscillation, deforestation, urbanisation and land use change, greenhouse gas emissions [
36,
105],
Australia and Oceania: El Niño - La Niña cycle , changes in ocean circulation - variation in ocean surface temperature and atmospheric pressure between western and eastern areas of the Indian Ocean, deforestation, change in water use and air pollution, intensive agriculture and overexploitation of water resources [
32,
36].
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.
Changes in temperature trends in continental catchments are the result of a complex interaction of both natural climatic and anthropogenic factors. The causes of changing trends are outlined below:
Africa: changes in atmospheric circulation: fluctuations in the belts and monsoons, changes in the oceanic system: phenomena such as El Niño and La Niña. Greenhouse gas emissions, land use change: urbanisation, deforestation and agriculture [
8,
32,
36,
106].
Asia: monsoon changes: changes in intensity and timing of monsoons, topography: hills, mountains and valleys, air pollution: dust and gas emissions can affect cloud formation and solar absorption, urbanisation: creation of heat island effect in urbanised areas [
32,
35,
36].
South America: changes in ocean circulation: fluctuations in ocean surface temperature, changes in land use: excessive deforestation and conversion of land to agriculture, [
8,
32,
36].
North America: changes in atmospheric circulation: influence of circulations such as AO and NAO, greenhouse gas emissions: increase in greenhouse gas concentrations, land use change: urbanisation and deforestation, [
8,
32,
36].
Australia and Oceania: changes in ocean cycles: El Niño and La Niña fluctuations, GHG emissions: increase in GHG concentrations, changes in water management: over-intensive use of water resources [
8,
32,
36,
106].
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 P2 measure factors are positive then there is a positive sign (+), if the trends of P2 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 [
107,
108].
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 [
107,
108].
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.