Next Article in Journal
Optimal Allocation Model of Water Resources Based on the Prospect Theory
Next Article in Special Issue
Links between Teleconnection Patterns and Water Level Regime of Selected Polish Lakes
Previous Article in Journal
Evaluating the Hydrologic Benefits of a Bioswale in Brunswick County, North Carolina (NC), USA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Probabilistic Assessment of Correlations of Water Levels in Polish Coastal Lakes with Sea Water Level with the Application of Archimedean Copulas

Institute of Physical Geography and Environmental Planning, Adam Mickiewicz University, Krygowskiego 10, 61-680 Poznań, Poland
*
Author to whom correspondence should be addressed.
Water 2019, 11(6), 1292; https://doi.org/10.3390/w11061292
Submission received: 29 May 2019 / Revised: 14 June 2019 / Accepted: 18 June 2019 / Published: 20 June 2019
(This article belongs to the Special Issue Effects of Oceanic-Atmospheric Oscillations on Rivers)

Abstract

:
The hydrology of coastal lakes differs significantly from that of inland water bodies due to the influence of the neighboring sea. Observed climatic changes are expected to enhance the effect of the sea on coastal lake ecosystems, which makes research on sea–lake interactions even more significant. In this study, on the basis of maximum annual and monthly values of water level, dependencies among the water levels of six lakes located along the southern coast of the Baltic Sea in Poland, and the Baltic Sea water levels were analyzed. First, the Spearman rank correlation and the linear regression method were applied. Then, selected copulas were used to find joint distributions of the studied time series. In the next stage, the degrees of synchronous and asynchronous occurrences of maximum water levels in lakes and the sea were calculated. The study revealed that correlations between the maximum annual water levels in coastal lakes and in the Baltic Sea in the selected gauge stations were very strong and statistically significant. These results were confirmed by a synchronicity analysis carried out with the help of a copula function. The highest relationship was detected in the case of Lake Resko Przymorskie (correlation coefficient 0.86, synchronicity 75.18%), while the lowest were observed in Lakes Jamno (0.62 and 58.20%, respectively) and Bukowo (0.60 and 56.82%, respectively). The relation strength between maximum water levels of the sea and coastal lakes may increase in the future due to sea level rises caused by climate warming.

1. Introduction

A total of 7081 Polish lakes have a surface area of more than one hectare [1]. A large group of those are coastal lakes characterized by complex water circulation due to the intensive supply of river waters from inland areas and the impact of sea waters [2].
Unlike in the case of the remaining parts of Poland, lakes in the coastal zone in the north of the country are characterized by the occurrence of hydrological phenomena and processes that do not occur in inland areas, e.g., the phenomenon of backwater and sea water intrusion, wind-induced blocking of lake water outflows, and the phenomenon of water rise, which are particularly observed on the Baltic coast in the period of winter storms [3].
Lakes located at the interface of land and sea, i.e., areas different in terms of both physical and chemical properties, develop a specific geoecosystem in which the direction of physical, chemical, and biological processes depends on which of the two environments has a stronger impact at a given moment [4]. A study of the coastal zone carried out by Rotnicki [5] confirmed that changes in the sea level determine the existence of the coastal lakes. The maximum Baltic Sea levels in the years 1889–2006 grew both on the Swedish and Polish coasts. The trend of the maximum water levels has increased in the last 60 years [6].
Coastal lakes are characterized by complex water circulation related to their location at the interface of two “environments”, sea and land, and they supply marine waters and river waters from inland areas. Research on coastal lakes in Poland and in the world has covered issues concerning, among other things, the chemistry/quality of water (including its salinity) [7,8,9,10,11,12,13], the development of lakes in the Holocene [14], the trophic status of lakes in the Holocene [15], biodiversity [16,17,18,19,20], sea water intrusions [4,21,22], and ice phenomena [23]. All of these said studies have confirmed the effect of the sea on the aforementioned properties of lakes.
Water level fluctuations in lakes have various characteristics and causes. Water level changes, their scale, dynamics, and change tendencies depend on natural [24,25,26] and anthropogenic [27,28] factors [29]. In the case of coastal lakes, more dynamic changes in water levels are observed than in the remaining lakes. This is related to the effect of the Baltic Sea [30,31,32,33,34]. Research on the water level regime in lakes in North Poland has shown that coastal lakes constitute a separate type [35,36,37]. They are lakes with considerable mean annual amplitudes exceeding 100 cm [35], and according to Reference [36], they show a four-period rhythm of seasonal water level variability determined by atmospheric circulation and the hydrodynamic state of the Baltic Sea. In a classification of lakes through values of pentad water level coefficients in an average annual cycle, lakes are included in a separate group with the smallest range of the parameter. The lakes are the only ones to be characterized by a period of higher water levels in the winter (December, January) and in the summer–autumn period (July–September) and by a period of low water levels in spring (March–May) [37].
The observed strong effect of the sea on changes in the ecosystem of coastal lakes co-occurs with modern climatic changes. Changes in the thermal regime and water levels of the sea will undoubtedly affect lake ecosystems even more strongly, increasing the importance of this research. The possibility of forecasting and the accuracy of forecasts of one variable (water level in lakes) by means of another variable (sea water levels) depends on the degree of dependency between them. The dependency between variables is usually analyzed by means of the Pearson linear correlation coefficient, as well as Spearman and Kendall rank correlation coefficients. These correlation measures have their drawbacks, however, which in the case of analyses of correlations between variables may affect their credibility. According to literature on the subject, linear correlation is not a resistant measure [38]. This means that untypical observations (so-called outliers) can significantly affect the value of the correlation coefficient. The correlation also shows no invariance toward nonlinear transformations. Copulas and related measures of correlation strength have no such drawbacks. They are exceptionally useful and practical approaches in studies on correlations between random variables. The copula theory has already been applied in limnological research, e.g., in the determination of the probability of occurrence of hydrological draught in Lake Poyang and the Yangtse River system [39], in the determination of correlations between extreme water levels in Lake Ontario and the Saint Lawrence River [40], in the determination of correlations of water levels in the Great Lakes and the Saint Lawrence River [41], and in an analysis of the quantity of water supplied to dam reservoirs in the context of forecasting the occurrence of draught [42]. This paper points to the possibility of obtaining new information resulting from the applied new research approach. It concerns, for example, separate and combined probabilities of the occurrence of maximum water levels in lakes and the sea (probabilistic assessment of correlations of maximum water levels in lakes and the sea).
The objective of this paper was the determination of correlations of water levels in coastal lakes of the South Baltic with water levels in the sea (through the application of the correlation coefficient) and the consideration of the possibility of simultaneous modeling of correlations between the selected variables (with the application of copula theory), which provided more possibilities for interpretation of the detected relations than traditional statistical tools did.
The strength of the relationship between maximum water levels in lakes and the sea is presented in the form of the probability of occurrence of synchronous and asynchronous situations. This approach has so far been applied in research on, among other things, correlations between the volume of outflow and volume of supplied debris [43,44,45] and the volumes of outflow in rivers [46,47]. The results of this study provide a more precise view of connections between water levels in the Baltic Sea and the analyzed coastal lakes. This is particularly important in the context of projected climate warming, which through increased sea water levels may have an impact on coastal lakes, including fluctuations in their water levels, chemistry, and ecological status.

2. Materials and Methods

2.1. Materials and Study Area

The analysis of correlations of water levels in coastal lakes with water levels in the sea was based on the maximum annual and monthly values of water levels in lakes and the sea. They were calculated based on daily values of water levels from the years 1976–2015, which were obtained from the collections of the Institute of Meteorology and Water Management, National Research Institute. Due to the availability of measurement data, in the case of Lake Resko Przymorskie, the analysis covered the multiannual period 1981–2015. The analysis concerned 6 selected coastal lakes where observations of water levels are conducted, namely Resko Przymorskie (station Dźwirzyno), Jamno (station Unieście), Bukowo (station Bukowo Morskie), Gardno (station Gardna Wielka), Łebsko (station Izbica), Druzno (station Żukowo), and 4 water gauge stations in the Baltic Sea (Kołobrzeg, Ustka, Łeba, and Tolkmicko). The water gauge stations in the sea were assigned to lakes based on location. Each lake was ascribed to the nearest water gauge station on the sea (Figure 1).
Lakes Resko Przymorskie, Jamno, Bukowo, Gardno, and Łebsko are connected to the sea by means of canals, simultaneously constituting the mouth sections of rivers. Lake Jamno is connected to the sea by means of the Jamno Canal, Bukowo through the Szczucza Canal, Gardno through the Łupawa Canal, Łebsko through the Łeba Canal [32], and Lake Resko Przymorskie through the Resko Canal [48]. Lake Druzno, located in the east of the study area, has a connection to the Baltic Sea through the Vistula Lagoon. The analyzed lakes are shallow, with a mean depth of 1.2 m to 1.8 m (Table 1), but belong to the group of the largest lakes in Poland in terms of surface area. They are characterized by high values of mean, maximum, and extreme annual water level amplitudes. Extreme annual amplitudes in the lakes in the analyzed multiannual period exceed 100 cm. They are the highest in Lake Resko Przymorskie (200 cm) and the lowest in Lake Łebsko (127 cm) (Table 2). The highest water level amplitudes are observed in Kołobrzeg (Table 2). The lowest mean and maximum annual amplitudes occur in Łeba (132 cm and 151 cm), and the most extreme amplitudes are in Ustka (188 cm).
The variability coefficient of daily water levels varies from 0.036 in Lake Łebsko to 0.044 in Lake Bukowo. In the case of the Baltic Sea, the variability coefficient in all stations is 0.041 (Table 2).

2.2. Methods

2.2.1. Mann–Kendall Test

The first stage of the analysis involved the analysis of tendencies of fluctuations of mean annual and maximum annual water levels both in lakes and in the sea in particular stations. This involved the application of the nonparametric Mann–Kendall test, which detects trends in temporal sequences. The Mann–Kendall test was employed in the MAKESENS stencil in Microsoft Excel software. This is a Mann–Kendall test that was extended by scientists from the Finnish Meteorological Institute [50].

2.2.2. Correlation

Next, the Spearman rank correlation and the linear regression method were applied in the analysis of correlations of water levels in the lakes with water levels in the sea (similarly to the case of Pearson’s correlation coefficient). Spearman’s rank correlation test is a nonparametric test applied in analyses of correlations between quantitative parameters in the case of a low number of observations or a statistic distribution deviating from normal. Then, statistical significance was analyzed at a level of p = 0.05, p = 0.01, and p = 0.001. A statistical assessment of significance of correlation coefficients was performed by means of statistic t:
t = r n 2 1 r 2
Statistic t has a Student distribution with n − 2 degrees of freedom.

2.2.3. Application of the Copula Theory

The concept of copula was introduced by Reference [51], which defined copula as a joint distribution function of standard uniform random variables. Modeling joint distribution using copula relaxes the restriction of traditional flood frequency analysis by selecting marginals from different families of probability distribution functions for flood characteristics [52].
First, the best matching statistical distributions were selected for the analyzed data series. The following distributions reflect the maximum values in series of hydrological data: Weibull, Gamm, Gumbel, and log-normal [53]. Parameters of the distributions were estimated by means of the highest probability method [54]. For the purpose of an assessment of the goodness of fit of a given distribution in the data series, the Akaike information criterion (AIC) was applied [55], which was calculated from the following formulas:
AIC = Nlog (MSE) + 2(no. of fitted parameters),
where MSE is the mean square error, and N is the sample size, or
AIC = −2log(maximum likelihood for model) + 2(no. of fitted parameters).
The best model is the one that has the minimum AIC value [55].
In the next step, the copula method was used to construct the joint distribution of lake and sea stages. In general, a bivariate Archimedean copula can be defined as [56]
C θ ( u , v ) = ϕ 1 { ϕ ( u ) + ϕ ( v ) } ,
where the θ subscript of copula C is the parameter hidden in the generating function ϕ , and ϕ is a continuous function called a generator that strictly decreases and is convex from I = [0,1] to [0, ϕ (0)].
The Archimedean copula family is often applied in hydrological studies, for example in flood frequency analyses. It was found in References [57] and [58] that a copula-based flood frequency analysis performs better than a conventional flood frequency analysis, as joint distribution based on a copula fits the empirical joint distribution (i.e., from observed data using a plotting position formula) better than the established standard joint parametric distribution. Numerous successful applications of copula modeling have been achieved, most notably in survival analysis, actuarial science, and finance [52].
A large variety of copulas belong to the Archimedean copula family and can be applied when the correlation between hydrologic variables is positive or negative. The proofs of these properties have been reported by References [59] and [60]. For this reason, one-parameter Archimedean copulas, including the Clayton family, the Gumbel–Hougaard family, and the Frank family, were used in this study. The Gumbel–Hougaard and Clayton copula families are appropriate only for positively correlated variables, while the Frank family is appropriate for both negatively and positively correlated variables (Table 3).
where D k ( x ) is Debye function, for any positive integer k,
D k ( x ) = k k x 0 x t k e t 1 d t .
The best fitted joint distribution was selected through comparison to the empirical joint distribution using the Akaike information criterion (AIC), as mentioned earlier.
For each compared pairs of series, based on previously calculated parameters of statistical distribution, 5000 hypothetical points were generated. They were used for the selection of the best fitted family of copulas for a given pair of series and then for the development of an appropriate copula. Based on empirical pairs of values for particular years and generated hypothetical points, graphs with probability curves (expressed in return periods) were developed (Figure 2).
The next stage involved a calculation of the degree of synchronicity (synchronous occurrence) and asynchronicity (asynchronous occurrence) of maximum water levels in lakes and the sea. For each pair of stations, probability curves at a level of 62.5% (once in 1.6 years), 37.5% (once in approximately 2.7 years), 20% (once in 5 years), 10% (once in 10 years), 2% (once in 50 years), 1% (once in 100 years), 0.5% (once in 200 years), and 0.2% (once in 500 years) are presented (Figure 2).
The obtained data were then analyzed based on probabilities of 62.5% and 37.5% [44]. Nine sectors were designated, representing different relations between probable maximum water levels. Based on generated points with a distribution imitating the shared distribution of values from comparable water gauge stations and their participation in particular sectors (Figure 2), 3 sectors with synchronous occurrences of maximum water levels were designated:
  • Sector 1: LHWLS–LHWLL (XS62.5%, YL62.5%);
  • Sector 5: MHWLS–MHWLL (S62.5%< XS37.5%, L62.5% < YL37.5%);
  • Sector 9: HHWLS–HHWLL (X > S37.5%, Y > L37.5%);
as well as 6 sectors with asynchronous occurrence:
  • Sector 2: LHWLS–MHWLL (XS62.5%, L62.5% < YL37.5%);
  • Sector 3: LHWLS–HHWLL (XS62.5%, Y > L37.5%);
  • Sector 4: MHWLS–LHWLL (S62.5%< XS37.5%, YL62.5%);
  • Sector 6: MHWLS–HHWLL (S62.5%< XS37.5%, Y > L37.5%);
  • Sector 7: HHWLS–LHWLL (X > S37.5%, YL62.5%);
  • Sector 8: HHWLS–MHWLL (X > S37.5%, L62.5% < YL37.5%);
where X = the values of x coordinates of generated points, Y = the values of y coordinates of generated points, S62.5% = the value of the maximum sea water level with a probability of exceedance of 62.5%, S37.5% = the value of the maximum sea water level with a probability of exceedance of 37.5%, L62.5% = the value of the maximum sea water level with a probability of exceedance of 62.5%, L37.5% is the value of the maximum sea water level with a probability of exceedance of 37.5%, WL = water level, LH = “low high”, MH = “mean high”, and HH = “high high”.
The percent contribution of points included in sectors 1, 5, and 9 permitted a determination of the degree of synchronicity of maximum water levels between the two analyzed water bodies in a given time unit.
Synchronous and asynchronous occurrences of maximum water levels were determined through a determination of threshold values of probability ranges:
  • Probable maximum water levels with a probability of occurrence of <62.5% were designated as LHWL;
  • Probable maximum water levels with a probability of occurrence in a range >62.5% and <37.5% were designated as MHWL; and
  • Probable maximum water levels with a probability of occurrence >37.5% were designated as HHWL.
For example, the occurrence of LHWL in a given lake is a synchronous event if LHWL also occurs in the Baltic Sea in a given time unit, and it is asynchronous if MHWL or HHWL is recorded there.
The total contribution of synchronous and asynchronous events is always 100%.
The mathematical and statistical processing of analysis results employed statistical procedures included in the following software programs: Excel (Microsoft), Statistica (TIBCO Software Inc.), and RStudio. The implementation of the graphic form employed QGIS (3.6.2. Noosa) and Publisher (Microsoft) software.

3. Results

3.1. Mann–Kendall Test

The study revealed that multiannual fluctuations of maximum annual water levels in three lakes (Bukowo, Gardno, Łebsko) showed statistically nonsignificant decreasing trends (Table 4). In the remaining lakes, increasing trends were observed, whereas a statistically significant trend was only determined for maximum water levels in Lake Jamno (p < 0.01). Maximum annual sea water levels in all stations in the analyzed multiannual period showed an increase, although it was not statistically significant.
In the case of multiannual fluctuations of maximum monthly water levels in lakes, an increasing tendency was observed for a major part of the year. Only Lake Resko Przymorskie showed no statistically significant trends. Maximum water levels in Lake Jamno showed positive, statistically significant trends throughout the year (p < 0.05). In the case of Lakes Bukowo and Druzno, significant trends were observed for maximum water levels from February to June, in the case of Lake Łebsko in May and June, and for Lake Gardno in June.
Fluctuations of maximum monthly water levels in the Baltic Sea showed positive, statistically nonsignificant trends for a major part of the year. In June, however, in the case of all stations with the exception of Łeba, an increasing statistically significant trend was observed (Table 4). Statistically significant increasing trends were also recorded in Kołobrzeg (1976–2015) and Ustka in May, and in Tolkmicko in August (p < 0.05).
In the majority of cases, significant trends of maximum sea water were accompanied by significant trends of maximum water levels in lakes in June (with the exception of Resko) and in May (with the exception of Resko, Bukowo, and Gardno).

3.2. Correlation

3.2.1. Maximum Annual Water Levels

This paper determined correlation coefficients and degrees of synchronicity of the occurrence of maximum annual water levels in the analyzed lakes and the sea. The highest coefficient of correlation of water levels was observed for the pair Kołobrzeg–Resko Przymorskie (0.86) and the lowest for the pair Ustka–Bukowo (0.60).
According to the analyses, correlations between maximum annual water levels in coastal lakes and maximum annual water levels in the Baltic Sea in selected water gauge stations were very strong and statistically significant (p < 0.001) (Table 5, Figure 3). In the case of Lake Resko Przymorskie, as much as 68% of the variation in maximum water levels in the lake was accounted for by the variability in sea water levels. The weakest correlations were observed in the case of Lakes Jamno and Bukowo, where only (respectively) 30% and 39% of the variation in maximum water levels in the lakes was accounted for by the variance of water levels in the Baltic Sea.

3.2.2. Maximum Monthly Water Levels

Coefficients of correlation between maximum monthly water levels in lakes and the sea were high and in the majority of cases statistically significant (p < 0.001) (Figure 4). The strongest and most statistically significant (p < 0.001) correlations with sea water levels were usually observed in the cool half-year, in the months from November to March. The highest coefficients of correlation between the analyzed lakes and the Baltic Sea (>0.90) were recorded in November, January (Gardno, Łebsko), and March (Resko Przymorskie, Łebsko). Water levels in three lakes, namely Gardno, Łebsko, and Resko Przymorskie, strongly correlated with sea water levels also in summer–autumn months. The weakest correlations of water levels in lakes with sea water levels were observed in April. Nonetheless, only the water levels in Lake Jamno correlated statistically nonsignificantly with sea water levels.

3.3. Synchronous–Asynchronous Encounter Probability

3.3.1. Maximum Annual Water Levels

Like the coefficient of correlation, modeling correlations by means of a copula function permitted the determination of the strength of correlations between water levels in lakes and the Baltic Sea. An additional advantage of the copula function was the possibility of quantifiable determination of the probability of occurrence of each analyzed variable separately and the probability of co-occurrence of maximum water levels in lakes and the sea, as well as the degree of their synchronicity.
The calculated synchronicity of maximum water levels in the lakes and the Baltic Sea oscillated between 56.82% and 75.18%, and asynchronicity between 24.82% and 43.18% (Table 5 and Table 6). The patterns of relations of water levels in lakes with sea water levels resulting from the synchronicity analysis were in accordance with results obtained from the correlation analysis (Table 5). The highest degree of synchronicity concerned the occurrence of maximum annual water levels in the sea and Lake Resko Przymorskie (75.18%), and the lowest was in Lake Bukowo (56.82%). This means that the maximum water level in Lake Resko Przymorskie in a hydrological year is in the same range of probability of occurrence as the maximum sea water level (i.e., one of three possible synchronous situations occurs: LHWLS–LHWLL, MHWLS–MHWLL, or HHWLS–HHWLL), with a probability of approximately 75%. On the contrary, i.e., in asynchronous situations, the maximum water level in Lake Resko Przymorskie, being in a different range of probability than the maximum sea water level in the same hydrological year, occurs with a probability of once in four years (p = 25%).
Next to Lake Resko Przymorskie, a relatively strong correlation (>65%) of maximum water levels was determined for Lakes Łebsko and Gardno, an average one (approximately 60%) for Lakes Druzno and Jamno, and a relatively low one (<60%) for Lakes Jamno and Bukowo. None of the lakes showed more than 50% asynchronicity. This means that the probability of occurrence of similar maximum water levels in a hydrological year in the sea and lakes is higher than the probability of occurrence of asynchronous events.

3.3.2. Maximum Monthly Water Levels

Synchronicity for maximum monthly water levels was calculated for Lake Jamno and stations in Kołobrzeg and for Lake Gardno and a station in Ustka (Table 7). The lakes were selected because the results of correlation and synchronicity pointed to their different relations with maximum annual sea water levels (Table 6).
The averaged synchronicity of maximum monthly water levels in Lake Jamno and the sea (55.56%) was inconsiderably lower than the synchronicity of maximum annual water levels in these water bodies (58.20%). Values above the average were obtained for the period from November to March. In the remaining months, they were below the average. In April, August, and October, the occurrence of asynchronous situations was more probable than synchronous ones. The highest synchronicity was characteristic of February and was slightly lower in December, and the lowest was in April. Concerning the correlation coefficient, the same pattern was observed (Figure 5).
The analysis of the results of synchronicity and asynchronicity for the pair Lake Jamno–station Kołobrzeg in particular sectors (Table 7; compare Figure 2) showed that the highest probability (approximately 10%) of occurrence of the situation MHWLS–MHWLL occurred in February. The probability (approximately 28%) of occurrence of situation HHWLS–HHWLL was in January, and LHWLS–LHWLL was in December (approximately 27%).
Averaged synchronicity of maximum monthly water levels in Lake Gardno and the Baltic Sea (72.39%) was, unlike in the case of the correlations described earlier for Jamno, higher than the synchronicity of maximum annual water levels (65.92%, Table 7). The strongest correlations (>75% synchronicity) were determined for maximum water levels in November, January, February, March, and October (Table 7). Weaker correlations were particularly observed in the period from April to September, and the weakest, as in the case of Lake Jamno, was in April (62.00%). Unlike in the case of Lake Jamno, however, no month showed a prevalence of total shares in asynchronous sectors. In other months, the highest synchronicity values also occurred (MHWLS–MHWLL and HHWLS–HHWLL) in November (17.28% and 33.92%, respectively), and LHWLS–LHWLL in February (32.82%).

4. Discussion

According to earlier research concerning the correlations of water levels in Polish coastal lakes with sea water levels, the strongest correlations occur in the cool season of the year, from November to April (Lakes Gardno and Łebsko). Weaker correlations are observed in the warm season, from May to October, particularly in the case of Lakes Jamno and Bukowo [32].
The highest correlations of monthly water levels in the lakes and the sea were usually observed in autumn, winter, and spring (Figure 3). This is generally related to lake alimentation in the autumn–winter period, where marine waters enter lakes through canals as a result of storms [33,61]. Storms on the Baltic Sea are usually observed from October to March and last from four to seven days [62].
In spring, in the melt season, alimentation with river waters is dominant. Months in which strong correlations were observed also corresponded to the course of water levels in lakes in the annual cycle. This was exemplified by the presented graph on the course of maximum pentad water levels in Lake Jamno and the Baltic Sea at station Kołobrzeg in the analyzed multiannual period (Figure 6). The maximum pentad water level was calculated as the averaged value of maximum water levels in the lake for a consecutive five days. The highest maximum water levels in the average annual cycle in Lake Jamno and Kołobrzeg were observed from December to January and from September to October.
The strongest correlation concerned water levels in Lakes Resko Przymorskie, Łebsko, and Gardno. In the case of the last two lakes, similar results were obtained by References [32,33] (the studies did not involve an analysis of the effect of the sea on water levels in Lake Resko Przymorskie). According to the author, strong correlations of water levels in the lakes and the sea can be justified by almost undisturbed water flow, which contributes to relatively fast evening of the water levels between the two water bodies. Fishing ports are located in the mouths of the Łupawa and Łeba canals (Rowy and Łeba), where dredging is performed on a current basis. The mouth of the Łupawa Canal is very rarely obstructed, and in the Łeba Canal, flow is practically undisturbed. In the mouth section of the Resko Canal (Błotnica), a marine port is also located in the municipality of Dźwirzyno. This may explain high correlations of water levels in Lake Resko Przymorskie with sea water levels. The weakest correlation with sea water levels concerned water levels in Lakes Jamno and Bukowo. According to Reference [32], this is related to inhibited water exchange with the sea, which is caused by the obstruction of the mouth sections of the canals to the sea by sand and frazil-ice banks. Blocking the outflow causes relatively high water levels in the lakes in comparison to the sea water level.
Moreover, as a result of climate change, positive trends of the maximum water levels in the Baltic Sea were detected. It was found that in the years 1889–2006, the maximum water levels increased on average from 0.39 cm/year in Świnoujście and Kołobrzeg to 0.54 cm/year in Gdańsk [6]. It is worth noting that correlations between the studied coastal lakes and the Baltic Sea refer not only to their water levels: the coast of the Baltic Sea is a place where high qualitative and quantitative variability of biotic [16,17,18,19,20] and abiotic elements [7,8,9,10,11,12,13,23] is observed, too.
Processes and phenomena occurring in coastal lakes are also characterized by enormous strength and dynamics. This was very evident, among other things, in the analysis of water level fluctuations in an average annual cycle. The highest mean monthly amplitudes in Polish lakes were recorded in coastal lakes [63]. The highest mean amplitude values were particularly observed in Lakes Druzno and Resko Przymorskie. In Lake Druzno, they were observed from October to April, varying from 45 cm in April to 53 cm in January. The highest mean amplitudes of water levels among the analyzed lakes were recorded in Lake Resko Przymorskie in the months from May to September. They reached values from 45 cm in May to 58 cm in August.
Another dynamically occurring process in the coastal zone of the South Baltic is the phenomenon of sea water intrusions [3,4]. The periodic inflow of marine waters into lakes results in changes in the water quality of the lake, particularly its salinity. A hydrodynamic comparative analysis of the waters of Lakes Jamno and Bukowo showed that they differ in terms of concentrations of indicators related to the inflow of sea waters. Considerably higher values were recorded in the case of Lake Bukowo. The waters of Lake Jamno were characterized by higher concentrations of indicators related to the inflow of waters from the catchment [64]. Similar patterns were observed in the analysis of the correlation of mean monthly water levels in the lakes with those in the Baltic Sea [33]. Stronger correlations of water level fluctuations were observed in the case of Lake Bukowo, and the weakest in Lake Jamno.
Over the last 20 years, copula functions have been found to have broad applications in hydrological, meteorological, and climatological data analyses. One of the methods of their application is the analysis of synchronicity and asynchronicity of occurrences of phenomena. According to other authors, copula functions have proven efficient in comprehensive analyses of, e.g., dependencies between the discharge volume and volume of supplied debris [44,45]. They also permit the quantitative determination of the frequency of occurrence of phenomena or disasters and can be used to prevent undesirable events (i.e., drought and floods) and limit the risk of their occurrence, i.e., their use in shaping the rational management of water resources in a given area [43,45,46,47]. This paper presents the possibility of applications of synchronicity analyses in the analysis of the dependency of maximum water levels in coastal lakes on sea water levels, which can also contribute to an increase in the level of water management in coastal areas.
In a broader approach, copula functions have also been used in limnological analyses, as shown in the introduction. The research was primarily conducted in the context of the possibility of the occurrence of forecasting drought. As in this paper, other publications [40,41] have presented analyses of maximum water levels in lakes. The analyses, however, were performed particularly for the purpose of determining changes in the value of average water levels and their variability in a given multiannual period. The authors did not focus on the direct analysis of correlations of water levels between water bodies. Coastal lakes are an exceptional group of lakes. Their analysis permits a determination of the strength of impact of the sea on the water conditions of the coast: Similar research has not been sufficiently presented in global literature.

5. Conclusions

This paper confirmed the effect of water level fluctuations in the Baltic Sea on water levels in coastal lakes. The analysis of maximum correlations of maximum annual water levels in the lakes confirmed strong correlations with sea water levels. The highest correlation coefficients for maximum annual water levels were recorded on Lakes Resko Przymorskie, Gardno, and Łebsko. In the case of maximum monthly water levels, the highest correlations were observed for Lakes Łebsko, Gardno, and Resko Przymorskie in autumn–winter months, and the lowest for Lakes Jamno and Bukowo in spring and summer. The application of a copula function permitted the calculation of probable maximum water levels in the lakes and the sea and a determination of the degree of synchronicity of occurrence of maximum water levels in the lakes and the sea.
The research showed the usefulness of advanced statistical tools in the analysis of correlations between water levels in coastal lakes and the sea. Considering the fact that correlation coefficients (also rank correlation coefficients) provide incomplete information on correlation, copulas were applied as a tool for modelling the correlation structure. The application of copulas proved to be a flexible tool for modeling the correlation structure in both distribution tails. Dependencies in the tails determine higher or lower degrees of synchronicity for differently correlated variables.
The analysis of correlations of maximum annual water levels in the lakes and the sea revealed that the highest water levels in all of the analyzed cases showed higher synchronicity than the lowest water levels. In the annual cycle, maximum monthly water levels in the lakes and the sea showed that in the cool half-year from November to April, the correlation between the highest water levels, particularly in the case of Lake Jamno, was stronger than between the lowest water levels, as evidenced by higher synchronicity values (Table 7). In the summer half-year (May–October), the situation was the opposite. A stronger correlation and higher synchronicity occurred in the lowest water levels in lakes and the sea than in the highest water levels. Despite doubts and limitations [65,66], copulas and related measures of correlation strength are an exceptionally useful approach in the issue of the analysis of correlations between hydrological and hydroclimatic variables.
Given the relatively small number of studies of that kind in the world scientific literature, this research is of great importance both on a regional and global scale. Due to climate warming, the sea level has been increasing, which may have a significant impact on the functioning of coastal lakes, including fluctuations of their water levels, chemistry, and ecological status. Thus, more detailed studies on these issues are necessary, and this research can be the starting point for such studies in the future.

Author Contributions

Conceptualization, D.W.; methodology, A.P. and L.S.; software, K.P., A.P., L.S.; validation, K.P., A.P.; formal analysis, K.P., A.P.; investigation, K.P., A.P.; resources, K.P.; writing—original draft preparation, K.P., A.P.; writing—review and editing, K.P., A.P., D.W., L.S.; visualization, K.P., A.P.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Choiński, A. Limnologia Fizyczna Polski; Wydawnictwo Naukowe UAM: Poznań, Poland, 2007. [Google Scholar]
  2. Mikulski, Z.; Bojanowicz, M. Bilans wodny jeziora znajdującego się pod wpływem morza (na przykładzie jeziora Druzno). Przegląd Geofiz. 1967, 12, 3–4. [Google Scholar]
  3. Cieśliński, R. Zróżnicowanie typologiczne i funkcjonalne jezior w polskiej strefie brzegowej południowego Bałtyku. Probl. Ekol. Kraj. 2010, 26, 135–144. [Google Scholar]
  4. Cieśliński, R.; Drwal, J.; Chlost, I. Sea water intrusions to the Lake Gardno. Balt. Coast. Zone 2009, 13 Pt I, 85–98. [Google Scholar]
  5. Rotnicki, K. (Ed.) Przemiany Środowiska Geograficznego Nizin Nadmorskich Południowego Bałtyku w Vistulianie i Holocenie; Bogucki Wydawnictwo Naukowe UAM: Poznań, Poland, 2001. [Google Scholar]
  6. Wolski, T.; Wiśniewski, B. Changes of maximum sea levels at selected gauge stations on the Polish and Swedish Baltic coast. Stud. Prac. WNEIZ 2012, 29, 209–227. [Google Scholar]
  7. Schallenberg, M.; Larned, S.T.; Hayward, S.; Arbuckle, C. Contrasting effects of managed opening regimes on water quality in two intermittently closed and open coastal lakes. Estuar. Coast. Shelf Sci. 2010, 86, 587–597. [Google Scholar] [CrossRef]
  8. Dye, A.; Barros, F. Spatial patterns of macrofaunal assemblages in intermittently closed/open coastal lakes in New South Wales, Australia. Estuar. Coast. Shelf Sci. 2005, 64, 357–371. [Google Scholar] [CrossRef]
  9. Trojanowski, J.; Antonowicz, J. Heavy metals in surface microlayer in water of Lake Gardno. Arch. Environ. Prot. 2011, 37, 75–88. [Google Scholar]
  10. Jarosiewicz, A. Seasonal changes of nutrients concentration in two shallow estuarine lakes Gardno and Lebsko; Comparison. Balt. Coast. Zone J. Ecol. Prot. Coastline 2009, 13, 121–133. [Google Scholar]
  11. Mudryk, Z.J. Antibiotic Resistance among Bacteria Inhabiting Surface and Subsurface Water Layers in the Estuarine Lake Gardno. Pol. J. Environ. Stud. 2002, 11, 401–406. [Google Scholar]
  12. Cieśliński, R. Zmiany zasolenia i poziomu wody jeziora Jamno wynikające z budowy wrót przeciwsztormowych. Inżynieria I Ochr. Środowiska 2016, 19, 517–539. [Google Scholar] [CrossRef]
  13. Sheela, A.M.; Letha, J.; Sabu, J.; Ramachandran, K.K.; Justus, J. Detection of extent of sea level rise in a coastal lake system using IRS Satellite Imagery. Water Resour. Manag. 2013, 27, 2657–2670. [Google Scholar] [CrossRef]
  14. Wunsam, S.; Schmidt, R.; Müller, J. Holocene lake development of two Dalmatian lagoons (Malo and Veliko Jezero, Isle of Mljet) in respect to changes in Adriatic sea level and climate. Palaeogeogr. Palaeoclimatol. Palaeoecol. 1999, 146, 251–281. [Google Scholar] [CrossRef]
  15. García-Rodríguez, F.; Sprechmann, P.; Metzeltin, D.; Scafati, L.; Melendi, D.L.; Volkheimer, W.; Mazzeo, N.; Hiller, A.; von Tümpling, W.; Scasso, F. Holocene trophic state changes in relation to sea level variation in Lake Blanca, SE Uruguay. J. Paleolimnol. 2004, 31, 99–115. [Google Scholar] [CrossRef] [Green Version]
  16. Timms, B.V. A limnological survey of the freshwater coastal lakes of east Gippsland, Victoria. Mar. Freshw. Res. 1973, 24, 1–20. [Google Scholar] [CrossRef]
  17. Paturej, E. Zooplankton Przymorskich Jezior Pobrzeża Bałtyckiego; Wydawnictwo Uniwersytetu Warmińsko-Mazurskiego: Olsztyn, Poland, 2005. [Google Scholar]
  18. Paturej, E.; Goździejewska, A. Zooplankton-based assessment of the trophic state of three coastal lakes–Łebsko, Gardno, and Jamno. Bull. Sea Fish. Inst. 2005, 3, 7–25. [Google Scholar]
  19. Paturej, E. Assessment of the trophic state of the coastal Lake Gardno based on community structure and zooplankton-related indices. Electron. J. Pol. Agric. Univ. 2006, 9, 17. [Google Scholar]
  20. Morozińska-Gogol, J.; Bachowska, A. Parasites of roach Rutilus rutilus (L.) from the Lake Lebsko. Balt. Coast. Zone 2002, 7, 75–80. [Google Scholar]
  21. Drwal, J.; Cieśliński, R. Coastal lakes and marine intrusions on the southern Baltic coast. Oceanol. Hydrobiol. Stud. 2007, 36, 61–75. [Google Scholar] [CrossRef]
  22. Grudzinska, I.; Vassiljev, J.; Saarse, L.; Reitalu, T.; Veski, S. Past environmental change and seawater intrusion into coastal Lake Lilaste, Latvia. J. Paleolimnol. 2017, 57, s10917–s10933. [Google Scholar] [CrossRef]
  23. Girjatowicz, J. Charakterystyki zlodzenia polskich jezior przybrzeżnych. Inżynieria Morska I Geotech. 2001, 2, 73–76. [Google Scholar]
  24. Li, X.-Y.; Xu, H.-Y.; Sun, Y.-L.; Zhang, D.-S.; Yang, Z.-P. Lake-level change and water balance analysis at lake Qinghai, West China during recent decades. Water Resour. Manag. 2007, 21, 1505–1516. [Google Scholar] [CrossRef]
  25. Duo, B.; Bianbaciren; Li, L.; Wang, W.; Zhaxiyangzong. The response of lake change to climate fluctuation in north Qinghai-Tibet Plateau in last 30 years. J. Geogr. Sci. 2009, 19, 131–142. [Google Scholar] [CrossRef]
  26. Singh, C.R.; Thompson, J.R.; French, J.R.; Kingston, D.G.; MacKay, A.W. Modelling the impact of prescribed global warming on runoff from headwater catchments of the Irrawaddy River and their implications for the water level regime of Loktak Lake, northeast India. Hydrol. Earth Syst. Sci. 2010, 14, 1745–1765. [Google Scholar] [CrossRef]
  27. Konatowska, M.; Rutkowski, P. The changes of the area and of the water-level of Kamińsko Lake (Zielonka Experimental Forestry Division) in the period of recent 150 years. Studia I Mater. Cent. Edukac. Przyr. Leśnej 2008, 10, 205–217. [Google Scholar]
  28. Zhuang, C.; Ouyang, Z.; Xu, W.; Bai, Y.; Zhou, W.; Zheng, H.; Wang, X. Impacts of human activities on the hydrology of Baiyangdian Lake, China. Environ. Earth Sci. 2011, 62, 1343–1350. [Google Scholar] [CrossRef]
  29. Wrzesiński, D.; Ptak, M. Water level changes in Polish lakes during 1976–2010. J. Geogr. Sci. 2016, 26, 83–101. [Google Scholar] [CrossRef]
  30. Chlost, I.; Cieśliński, R. Change of level of waters Lake Łebsko. Limnol. Rev. 2005, 5, 17–26. [Google Scholar]
  31. Girjatowicz, J.P. Związek między poziomem wody w jeziorach przybrzeżnych i wodami morskimi polskiego wybrzeża Morza Bałtyckiego. Przegląd Geofiz. 2008, 2, 141–153. [Google Scholar]
  32. Girjatowicz, J.P. Miesięczne i sezonowe charakterystyki poziomów wody wybranych polskich jezior przybrzeżnych. Inżynieria Morska Geotech. 2008, 1, 27–32. [Google Scholar]
  33. Girjatowicz, J.P. Wpływ Morza Bałtyckiego na poziomy wód polskich jezior przybrzeżnych. Inżynieria Morska Geotech. 2011, 1, 18–22. [Google Scholar]
  34. Fac-Beneda, J. Charakterystyka hydrologiczna jeziora Druzno. In Monografia Jeziora Druzno, Monografia Przyrodnicza; Nitecki, C., Ed.; Wyd. Mantis: Olsztyn, Poland, 2013. [Google Scholar]
  35. Choiński, A. Zarys Limnologii Fizycznej Polski; Wydawnictwo Naukowe UAM: Poznań, Poland, 1975. [Google Scholar]
  36. Borowiak, D. Reżimy Wodne i Funkcje Hydrologiczne Jezior Niżu Polskiego; Katedra Limnologii Uniwersytetu Gdańskiego: Gdańsk, Poland, 2000. [Google Scholar]
  37. Plewa, K. Typy przebiegu pentadowych współczynników stanu wody jezior Niżu Polskiego. Bad. Fizjogr. Ser. A – Geogr. Fiz. 2018, A69, 161–177. [Google Scholar]
  38. Gurgul, P.; Syrek, R. Zastosowanie mieszanki kopul do modelowania współzależności pomiędzy wybranymi sektorami gospodarki. Ekon. Menedżerska 2009, 6, 129–139. [Google Scholar]
  39. Zhang, D.; Chen, P.; Zhang, Q.; Li, X. Copula-based probability of concurrent hydrological drought in the Poyang lake-catchment-river system (China) from 1960 to 2013. J. Hydrol. 2017, 553, 773–784. [Google Scholar] [CrossRef]
  40. Assani, A.A.; Landry, R.; Biron, S.; Frenette, J.-J. Analysis of the interannual variability of annual daily extreme water levels in the St Lawrence River and Lake Ontario from 1918 to 2010. Hydrol. Process. 2014, 28, 4011–4022. [Google Scholar] [CrossRef]
  41. Assani, A.A. Analysis of the Impacts of Man-Made Features on the Stationarity and Dependence of Monthly Mean Maximum and Minimum Water Levels in the Great Lakes and St. Lawrence River of North America. Water 2016, 8, 485. [Google Scholar] [CrossRef]
  42. Kim, K.; Lee, S.; Jin, Y. Forecasting quarterly inflow to reservoirs combining a copula-based bayesian network method with drought forecasting. Water 2018, 10, 233. [Google Scholar] [CrossRef]
  43. Zhou, N.-Q.; Zhao, L.; Shen, X.-P. Copula-based Probability Evaluation of Rich-Poor Runoff and Sediment Encounter in Dongting Lake Basin. Sci. Geogr. Sin. 2014, 34, 242–248. [Google Scholar] [CrossRef]
  44. Zhang, J.; Ding, Z.; You, J. The joint probability distribution of runoff and sediment and its change characteristics with multi-time scales. J. Hydrol. Hydromech. 2014, 62, 218–225. [Google Scholar] [CrossRef] [Green Version]
  45. You, Q.; Jiang, H.; Liu, Y.; Liu, Z.; Guan, Z. Probability Analysis and Control of River Runoff–sediment Characteristics based on Pair-Copula Functions: The Case of the Weihe River and Jinghe River. Water 2019, 11, 510. [Google Scholar] [CrossRef]
  46. Gu, H.; Yu, Z.; Li, G.; Ju, Q. Nonstationary Multivariate Hydrological Frequency Analysis in the Upper Zhanghe River Basin, China. Water 2018, 10, 772. [Google Scholar] [CrossRef]
  47. Chen, J.; Shixiang, G.; Zhang, T. Synchronous-Asynchronous Encounter Probability Analysis of High-Low Runoff for Jinsha River, China, using Copulas. Matec Web Conf. 2018, 246, 01094. [Google Scholar] [CrossRef]
  48. Burandt, P.; Kobus, S.; Sidoruk, M.; Glińska-Lewczuk, K. Hydrographic and hydrological characteristics Part I: Liwia Łuża, Resko Przymorskie, Jamno, Kopań and Wicko. In Hydroecological Determinants of Fuctioning of Southern Baltic Coastal Lakes; Obolewski, K., Ed.; Polish Scientific Publishers PWN: Warszawa, Poland, 2017. [Google Scholar]
  49. Choiński, A. Katalog Jezior Polski; Wydawnictwo Naukowe Uniwersytetu im. Adama Mickiewicza: Poznań, Poland, 2006. [Google Scholar]
  50. Salmi, T.; Määttä, A.; Anttila, P.; Ruoho-Airola, T.; Amnell, T. Detecting Trends of Annual Values of Atmospheric Pollutants by the Mann-Kendall Test and Sen’s Slope Estimates—The Excel Template Application Makesens; Publications on Air Quality 31; Finnish Meteorological Institute: Helsinki, Finland, 2002; p. 35.
  51. Sklar, M. Fonctions de répartition à n dimensions et leurs marges. Publ. L’inst. Stat. L’université De Paris 1959, 8, 229–231. [Google Scholar]
  52. Karmakar, S.; Simonovic, S.P. Water resources research report. In Flood Frequency Analysis Using Copula with Mixed Marginal Distributions; The University of Western Ontario Department of Civil and Environmental Engineering: London, ON, Canada, 2007. [Google Scholar]
  53. Koutsoyiannis, D. Probability and Statistics for Geophysical Processes; National Technical University of Athens: Athens, Greece, 2008. [Google Scholar]
  54. Zhao, L.; Xia, J.; Sobkowiak, L.; Wang, Z.; Guo, F. Spatial Pattern Characterization and Multivariate Hydrological Frequency Analysis of Extreme Precipitation in the Pearl River Basin, China. Water Resour. Manag. 2012, 26, 3619–3637. [Google Scholar] [CrossRef]
  55. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–722. [Google Scholar] [CrossRef]
  56. Nelsen, R.B. An Introduction to Copulas; Springer: New York, NY, USA, 1999. [Google Scholar]
  57. Zhang, L.; Singh, V.P. Bivariate flood frequency analysis using the copula method. J. Hydrol. Eng. 2006, 11, 150–164. [Google Scholar] [CrossRef]
  58. Grimaldi, S.; Serinaldi, F. Asymmetric copula in multivariate flood frequency analysis. Adv. Water Resour. 2006, 29, 1115–1167. [Google Scholar] [CrossRef]
  59. Genest, C.; Rivest, L.-P. Statistical inference procedure for bivariate Archimedean copulas. J. Am. Stat. Assoc. 1993, 88, 1034–1043. [Google Scholar] [CrossRef]
  60. Genest, C.; Favre, A.-C. Everything you always wanted to know about copula Modeling but were afraid to ask. J. Hydrol. Eng. 2007, 12, 347–368. [Google Scholar] [CrossRef]
  61. Cyberski, J.; Jędrasik, J. Wymiana i cyrkulacja wód w jeziorze Gardno. In Zlewnia Przymorskiej Rzeki Łupawy i jej Jeziora; Korzeniewski, K., Ed.; WSP: Słupsk, Poland, 1992. [Google Scholar]
  62. Rokiciński, K. Geograficzna i hydrometeorologiczna charakterystyka Morza Bałtyckiego jako obszaru prowadzenia działań asymetrycznych. Zesz. Nauk. Akad. Mar. Wojennej 2007, 48, 65–82. [Google Scholar]
  63. Mikosch, T. Copulas: Tales and facts. Extremes 2006, 9, 3–20. [Google Scholar] [CrossRef]
  64. Mikosch, T. Copulas: Tales and facts—Rejoinder. Extremes 2006, 9, 55–62. [Google Scholar] [CrossRef]
  65. Plewa, K.; Wrzesiński, D.; Baczyńska, A. Przestrzenne i czasowe zróżnicowanie amplitud stanów wody jezior w Polsce w latach 1981-2015. Bad. Fizjogr. Ser. A Geogr. Fiz. 2017, A68, 115–126. [Google Scholar] [CrossRef]
  66. Cieśliński, R. Hydrochemiczna ocena porównawcza wód jeziora Jamno i Bukowo. Przegląd Geol. 2005, 53, 1066–1067. [Google Scholar]
Figure 1. Location of lakes and water gauge stations in the Baltic Sea.
Figure 1. Location of lakes and water gauge stations in the Baltic Sea.
Water 11 01292 g001
Figure 2. Example of combined accumulated curves of the probability of exceedance of maximum water levels and the determination of sectors with various degrees of synchronicity or asynchronicity of maximum water levels.
Figure 2. Example of combined accumulated curves of the probability of exceedance of maximum water levels and the determination of sectors with various degrees of synchronicity or asynchronicity of maximum water levels.
Water 11 01292 g002
Figure 3. Correlations of maximum annual water levels of the analyzed lakes with maximum annual water levels in the Baltic Sea.
Figure 3. Correlations of maximum annual water levels of the analyzed lakes with maximum annual water levels in the Baltic Sea.
Water 11 01292 g003
Figure 4. Coefficients of correlation of maximum water levels in the analyzed lakes and in the Baltic Sea and their levels of significance.
Figure 4. Coefficients of correlation of maximum water levels in the analyzed lakes and in the Baltic Sea and their levels of significance.
Water 11 01292 g004
Figure 5. Total accumulated curves of probability of exceedance of maximum annual water levels for the correlation between water levels in the lakes and water level in the Baltic Sea. Red points are empirical pairs of values of maximum annual water levels. Gray points are generated theoretical points based on calculated parameters of joint distribution.
Figure 5. Total accumulated curves of probability of exceedance of maximum annual water levels for the correlation between water levels in the lakes and water level in the Baltic Sea. Red points are empirical pairs of values of maximum annual water levels. Gray points are generated theoretical points based on calculated parameters of joint distribution.
Water 11 01292 g005
Figure 6. Course of maximum pentad water levels in Lake Jamno and the Baltic Sea at station Kołobrzeg.
Figure 6. Course of maximum pentad water levels in Lake Jamno and the Baltic Sea at station Kołobrzeg.
Water 11 01292 g006
Table 1. Morphometric data of the studied lakes.
Table 1. Morphometric data of the studied lakes.
No.LakeArea (ha)Volume
(thousand m3)
Average Depth (m)Max Depth (m)
1Resko Przymorskie (P.)559.07703.41.32.5
2Jamno2231.531,528.01.43.9
3Bukowo1644.032,071.71.82.8
4Gardno2337.530,950.51.32.6
5Łebsko7020.0117,521.01.66.3
6Druzno1147.517,352.01.22.5
Source: Reference [49].
Table 2. Water levels and their amplitudes in the analyzed lakes and in the Baltic Sea.
Table 2. Water levels and their amplitudes in the analyzed lakes and in the Baltic Sea.
Lake/Gauge StationWater Level (m)Amplitudes of Water Level (m)Coefficient of Variation (Cv)
Annual MeanMaxMinAnnual MaxAnnual MeanExtreme Multiyear
Lakes
Resko P.5025714601751112000.043
Jamno533578506104731310.039
Bukowo529577494126821580.044
Gardno520573485122881340.040
Łebsko509561475120861270.036
Druzno520576479138971760.041
Baltic Sea
Kołobrzeg5035794361811322200.041
Ustka5055734431611191880.041
Łeba5055704461511141960.041
Tolkmicko5095814561681242010.041
Table 3. Copula function, parameter space, generating function ϕ (t), and functional relationship of Kendall’s τθ with a copula parameter for selected single-parameter bivariate Archimedean copulas.
Table 3. Copula function, parameter space, generating function ϕ (t), and functional relationship of Kendall’s τθ with a copula parameter for selected single-parameter bivariate Archimedean copulas.
Copula Family C θ ( u , v ) Generator ϕ ( t ) Parameter θ Kendall’s τ θ (5)
Clayton m a x ( ( u θ + v θ 1 ) 1 θ , 0 ) 1 θ ( t θ 1 ) [ 1 , ) \ { 0 } τ = θ / ( 2 + θ )
Gumbel–Hougaard e x p { [ ( l n   u ) θ + ( l n   v ) θ ] 1 θ } ( ln t ) θ [ 1 , ) ( θ 1 ) / θ
Frank 1 θ l n [ 1 + ( e θ u 1 ) ( e θ v 1 ) e θ 1 ] l n e θ t 1 e θ 1 ( , ) \ { 0 } 1 + 4 [ D 1 ( θ ) 1 ] / θ
Table 4. Results of Mann–Kendall test for maximum annual and monthly water levels.
Table 4. Results of Mann–Kendall test for maximum annual and monthly water levels.
Lake/GaugesPeriodYearNovDecJanFebMarAprMayJunJulAugSepOct
Lakes
Resko P.1981–20150.612−0.724−0.2271.0530.6971.151−0.4421.0101.0960.355−0.2700.7120.555
Jamno1976–20153.0322.937 3.403 ¯ 2.0052.342 3.334 ¯ 3.742 ¯ 3.605 ¯ 4.340 ¯ 3.207 3.964 ¯ 3.1832.682
Bukowo1976–2015−0.944−0.723−1.5160.1752.0062.7292.0301.1442.3810.0701.2020.443−0.816
Gardno1976–2015−0.642−0.234−1.353−0.1750.3730.7590.6070.3392.662−0.233−0.467−1.295−0.455
Łebsko1976–2015−0.0820.233−0.4660.1751.0491.6210.6652.217 3.918 ¯ 0.9580.969−0.0580.793
Druzno1976–20141.0061.017−0.2181.0652.3252.4701.2732.967 4.515 ¯ 1.9772.3150.6540.921
Baltic Sea
Kołobrzeg1981–20150.483−0.6390.241−0.810−0.1421.364−0.1141.0942.0031.1220.0000.170−0.156
Kołobrzeg1976–20150.874−0.3500.606−0.0701.5031.4330.0002.4242.3071.2580.804−0.0230.466
Ustka1976–20150.513−0.1520.3960.1281.2131.2710.3382.4501.9951.4500.945−0.0580.023
Łeba1976–20150.000−0.2680.012−0.0350.8751.2700.3851.9021.8900.0580.187−0.6650.000
Tolkmicko1976–20141.393−0.0610.0480.4361.6341.659−0.1451.805 3.445 ¯ 1.1161.9620.1570.242
Statistically significant trends at levels of p < 0.001 ¯ , p < 0.01, p < 0.05.
Table 5. Results of analysis of correlation and synchronicity between the Baltic Sea and the analyzed lakes.
Table 5. Results of analysis of correlation and synchronicity between the Baltic Sea and the analyzed lakes.
No.Gauge stations
(Sea–Lake)
PeriodAnnual Max Water Level
Correlation CoefficientSynchronicity (%)Asynchronicity (%)
1Kołobrzeg–Resko P.1981–2015 0.86 ¯ 75.1824.82
2Kołobrzeg–Jamno1976–2015 0.62 ¯ 58.2041.80
3Ustka–Bukowo1976–2015 0.60 ¯ 56.8243.18
4Ustka–Gardno1976–2015 0.72 ¯ 65.9234.08
5Łeba–Łebsko1976–2015 0.76 ¯ 69.8230.18
6Tolkmicko–Druzno1976–2014 0.62 ¯ 60.4239.58
Statistically significant correlation at a level of p < 0.001 ¯ .
Table 6. Synchronous and asynchronous occurrence of maximum annual water levels in lakes and the sea (division into nine sectors) (% share).
Table 6. Synchronous and asynchronous occurrence of maximum annual water levels in lakes and the sea (division into nine sectors) (% share).
SectorKołobrzeg–JamnoUstka–BukowoUstka–GardnoŁeba–ŁebskoTolkmicko–DruznoKołobrzeg–Resko
124.5623.8026.7829.4825.8630.64
58.548.3010.1410.428.4613.12
925.1024.7229.0029.9226.1031.42
28.929.888.066.628.466.48
48.889.346.907.208.766.14
87.686.928.047.226.925.52
67.267.626.946.667.085.32
94.524.681.741.304.320.78
74.544.742.401.184.040.58
Syn.58.2056.8265.9269.8260.4275.18
Asyn.41.8043.1834.0830.1839.5824.82
Syn. = synchronicity, Asyn. = asynchronicity.
Table 7. Synchronicity and asynchronicity of maximum monthly water levels in Lake Jamno and the sea in Kołobrzeg and Lake Gardno and the sea in Ustka (division into nine sectors) (% share).
Table 7. Synchronicity and asynchronicity of maximum monthly water levels in Lake Jamno and the sea in Kołobrzeg and Lake Gardno and the sea in Ustka (division into nine sectors) (% share).
SectorNovDecJanFebMarAprMayJunJulAugSepOct
Jamno–Kołobrzeg
124.3026.6425.1624.8425.7017.2624.7625.4225.1023.1223.0622.72
58.649.069.0610.069.227.407.647.687.947.546.647.12
925.6027.0028.0427.8626.2818.1222.0022.0622.0619.2222.3620.00
28.327.788.328.368.349.789.527.186.307.968.507.80
48.947.728.768.709.009.288.566.907.587.668.667.08
87.167.427.146.966.808.868.9610.369.189.508.3810.22
67.088.027.506.767.668.868.589.5410.289.988.7811.12
95.503.003.083.023.5411.084.845.385.927.506.726.72
74.463.362.943.443.469.365.145.485.647.526.907.22
Syn.58.5462.7062.2662.7661.2042.7854.4055.1655.1049.8852.0649.84
Asyn.41.4637.3037.7437.2438.8057.2245.6044.8444.9050.1247.9450.16
Gardno–Ustka
132.6631.6431.4232.8231.3825.830.4228.6429.4630.2628.7232.94
517.2810.7415.4816.2213.049.6211.2011.0612.549.169.3014.34
933.9227.5232.9032.7430.8226.5826.0630.6830.4226.3025.7428.90
23.825.365.364.686.328.785.126.926.746.106.144.32
44.264.925.624.646.068.665.367.426.205.766.544.46
84.028.083.944.605.146.909.185.866.348.968.787.38
63.829.004.604.125.767.308.726.126.448.989.706.90
90.101.220.360.040.883.042.121.540.962.262.380.46
70.121.520.320.140.603.321.821.760.902.222.700.30
Syn.83.8669.9079.8081.7875.2462.0067.6870.3872.4265.7263.7676.18
Asyn.16.1430.1020.2018.2224.7638.0032.3229.6227.5834.2836.2423.82

Share and Cite

MDPI and ACS Style

Plewa, K.; Perz, A.; Wrzesiński, D.; Sobkowiak, L. Probabilistic Assessment of Correlations of Water Levels in Polish Coastal Lakes with Sea Water Level with the Application of Archimedean Copulas. Water 2019, 11, 1292. https://doi.org/10.3390/w11061292

AMA Style

Plewa K, Perz A, Wrzesiński D, Sobkowiak L. Probabilistic Assessment of Correlations of Water Levels in Polish Coastal Lakes with Sea Water Level with the Application of Archimedean Copulas. Water. 2019; 11(6):1292. https://doi.org/10.3390/w11061292

Chicago/Turabian Style

Plewa, Katarzyna, Adam Perz, Dariusz Wrzesiński, and Leszek Sobkowiak. 2019. "Probabilistic Assessment of Correlations of Water Levels in Polish Coastal Lakes with Sea Water Level with the Application of Archimedean Copulas" Water 11, no. 6: 1292. https://doi.org/10.3390/w11061292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop