iBet uBet web content aggregator. Adding the entire web to your favor.
iBet uBet web content aggregator. Adding the entire web to your favor.



Link to original content: https://doi.org/10.3390/rs16152753
Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland)
Next Article in Journal
Construction of High-Precision and Complete Images of a Subsidence Basin in Sand Dune Mining Areas by InSAR-UAV-LiDAR Heterogeneous Data Integration
Next Article in Special Issue
Climatic Influence on the Carotenoids Concentration in a Mediterranean Coastal Lagoon Through Remote Sensing
Previous Article in Journal
A New and Robust Index for Water Body Extraction from Sentinel-2 Imagery
Previous Article in Special Issue
Monitoring Harmful Algal Blooms and Water Quality Using Sentinel-3 OLCI Satellite Imagery with Machine Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland)

1
Department of Hydrology and Water Management, Adam Mickiewicz University, Krygowskiego 10, 61-680 Poznań, Poland
2
College of Hydraulic Science and Engineering, Yangzhou University, Yangzhou 225009, China
3
School of Energy and Environment, University of Phayao, Phayao 56000, Thailand
4
Atmospheric Pollution and Climate Research Unit, School of Energy and Environment, University of Phayao, Phayao 56000, Thailand
5
HUN-REN Balaton Limnological Research Institute, 8237 Tihany, Hungary
6
Department of Meteorology and Climatology, Adam Mickiewicz University, Krygowskiego 10, 61-680 Poznań, Poland
7
College of Fisheries, Huazhong Agricultural University, Wuhan 430070, China
8
Department of Land Improvement, Environmental Development and Spatial Management, Poznań University of Life Sciences, Piątkowska 94E, 60-649 Poznań, Poland
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(15), 2753; https://doi.org/10.3390/rs16152753
Submission received: 28 June 2024 / Revised: 22 July 2024 / Accepted: 23 July 2024 / Published: 27 July 2024
(This article belongs to the Special Issue Remote Sensing and GIS in Freshwater Environments)

Abstract

:
Water temperature is a fundamental parameter of aquatic ecosystems. It directly influences most processes occurring within them. Hence, knowledge of this parameter’s behavior, based on long-term (reliable) observations, is crucial. Gaps in these observations can be filled using contemporary methodological solutions. Difficulties in reconstructing water temperature arise from the selection of an appropriate methodology, and overcoming them involves the proper selection of input data and choosing the optimal modeling approach. This study employed the air2water model and Landsat satellite imagery to reconstruct the water temperature of Lake Miedwie (the fifth largest in Poland), for which field observations conducted by the Institute of Meteorology and Water Management—National Research Institute ended in the late 1980s. The approach based on satellite images in this case yielded less accurate results than model analyses. However, it is important to emphasize the advantage of satellite images over point measurements in the spatial interpretation of lake thermal conditions. In the studied case, due to the lake’s shape, the surface water layer showed no significant thermal contrasts. Based on the model data, long-term changes in water temperature were determined, which historically (1972–2023) amounted to 0.20 °C per decade. According to the adopted climate change scenarios by the end of the 21st century (SSP245 and SSP585), the average annual water temperature will be higher by 1.8 °C and 3.2 °C, respectively. It should be emphasized that the current and simulated changes are unfavorable, especially considering the impact of temperature on water quality. From an economic perspective, Lake Miedwie serves as a reservoir of drinking water, and changes in the thermal regime should be considered in the management of this ecosystem.

1. Introduction

Water temperature is one of the fundamental parameters of inland waters, as evidenced by extensive research on this topic. Water temperature has been analyzed in the context of the spawning process of Grass carp (Laurentian Great Lakes region). It was determined that spawning is most likely when the water temperature is below 25 °C and the discharge is above 10 m3/s [1]. Studies of several cascade reservoirs in China have shown that water temperature can directly or indirectly affect the diversity of microorganisms, where an increase in temperature caused an increase in the diversity of planktonic bacteria, while archaea showed the opposite situation. In the case of Lake Poyang, Yan et al. [2] demonstrated that water temperature regulates the dynamics of organophosphate esters in the food chain. According to Haddout et al. [3], water temperature in the lake is one of the key parameters in determining its ecological properties. Water temperature is closely related to qualitative parameters such as biochemical oxygen demand, chlorophyll-a, and transparency [4]. Due to the above, understanding the thermal characteristics of rivers or lakes is crucial for the proper interpretation of the processes and phenomena occurring in these ecosystems. This issue takes on particular significance in an era of rapidly changing climatic conditions, which have a clear impact on the components of the hydrosphere. Air temperature is commonly used in statistical models of water temperature as a surrogate variable for physical processes [5]. The currently observed increase in air temperature [6] has a clear translation into an increase in water temperature in lakes [7]. Changes in the thermal regime of lakes have further consequences for their functioning. Due to the complexity of the interactions between individual processes, their broader effects are often difficult to assess.
Originally, for many years, the primary source of information on water temperature in lakes was instrumental field measurements [8]. However, acquiring in situ data is not always satisfactory, i.e., it is limited both spatially and temporally. The first situation pertains to lakes in inaccessible areas, and the second relates to insufficiently long observation series or situations where measurements have been discontinued. The lack of knowledge about thermal conditions is particularly important in the case of large lakes, which have a real impact on the entire region in which they exist. An example of such a situation is Lake Miedwie—the fifth largest lake in Poland [9], for which water temperature measurements ended in 1989. Nowadays, in such cases, the environmental sciences are supported by a range of tools that allow for a detailed reconstruction and further prediction of various processes—including those concerning water temperature in lakes. Currently, an important role in studying the thermal characteristics of lakes is played by methodologies based on remote sensing [10], machine learning [11], or a combination of both approaches [12]. Difficulties in reconstructing water temperature arise from the selection of an appropriate methodology, and overcoming them involves the proper selection of input data and choosing the optimal modeling approach. According to research by Heddam et al. [11], the air2water model outperformed other water temperature modeling options. Therefore, this approach was used in the present article. The dataset obtained through this modeling method was enriched with Landsat satellite data, representing the first such combination compared to previous methods used in lake studies in Poland.
The goal of this study is the reconstruction and prediction of the water temperature of Lake Miedwie, followed by an assessment of changes in the thermal regime of the analyzed lake. The implementation of these objectives is based on the use of satellite imagery and modeling based on air temperature. Regarding other lakes in Poland where stationary water temperature measurements are lacking, this new approach, based on the proposed methodology, allows for a detailed understanding of one of the fundamental parameters of these ecosystems.

2. Materials and Methods

2.1. Study Area

In Poland, there are over 7000 natural lakes [13], primarily located in the northern part of the country [14]. However, the vast majority of these are small—nearly 3000 lakes fall within the 1–5 hectare range [13]. As mentioned previously, one of the largest is Lake Miedwie in the northwestern part of the country (Lon 14.88205116, Lat 53.26781949, 13.9 m a.s.l., Figure 1). Genetically, it is a post-glacial lake with a meridional trough shape [15]. Its surface area is 3500 ha; it has a maximum depth of 43.8 m and an average depth of 19.3 m, classifying it as a cryptodepression. The length of the lake’s shoreline is 39 km, and its surroundings are dominated by arable land. According to the climatic classification of Poland [16], Lake Miedwie is located in region 6—Szczecin Coastland. This region is characterized by an average annual temperature of 8.4 °C, with January being the coldest month (−0.7 °C) and July the warmest (17.6 °C). The average annual precipitation amounts to 558 mm. From an economic perspective, the lake serves as a water reservoir for the nearby city of Szczecin (with 400,000 residents), located approximately 20 km away. Additionally, there are recreational centers operating along the lake.

2.2. Materials

To achieve the study objective, a diverse dataset including in situ measurements and satellite images was utilized. Water temperature measurements at Lake Miedwie were carried out from 1972 to 1989. Temperature registration was conducted in the surface water layer (at a depth of 0.4 m) at a permanent hydrological station in the northern part of the lake (Figure 1). Additionally, to verify further model calculations, water temperature data from Lake Lubie during the years 1972–2022 were used, where measurements were performed in a similar way. Lake Lubie is located 63 km to the northeast. Meanwhile, air temperature, serving as an explanatory variable for water temperature, was recorded at the Szczecin station (1.0 m a.s.l., Figure 1). These included records obtained from a meteorological gauge at a height of 2 m above ground level. Both water and air temperature observations were conducted as part of the standard monitoring by the Institute of Meteorology and Water Management—National Research Institute.

2.3. Methods

2.3.1. Reconstruction and Prediction of Water Temperature with the Application of the Air2water Model

In the first stage, the water temperature of Lake Miedwie from 1990 to 2023 was reconstructed using the air2water model. Subsequently, predictions of water temperatures in Lake Miedwie for the period 2024 to 2100 were made, utilizing estimated daily air temperatures (Section 2.3.3).
The air2water model is a widely used tool to forecast lake surface water temperatures (LSWT) that has been well demonstrated in lakes worldwide [17,18,19,20]. Compared to other water temperature modeling options, it was shown that for the validation dataset, the air2water model had an average RMSE of 1.031, while the RMSE for ERT, MARS, M5Tree, MLPNN, and RF was 1.413, 1.520, 1.582, 1.726, and 1.499, respectively [11]. It is a lumped model coupling energy balance of the surface layer of the lake and stochastic calibration of model parameters. Among the available model versions, e.g., 3-para, 6-para, 8-para, and 9-para versions, the 6-para version has been widely used to model LSWT in different lakes worldwide [21,22], hence it is used in this study as well.
d T w d t = 1 δ a 1 + a 2 T a a 3 T w + a 5 c o s 2 π t t y a 6
δ     exp T w T h a 4                                 T w T h 1                                                     T w < T h                  
where Th is a reference value for the deep lake water temperature (4 °C for dimictic lakes and the minimum or maximum water temperature for the other lakes), δ is a dimensionless term representing the ratio between the volume of the surface lake layer and a reference volume, Tw is daily lake surface water temperature (°C), Ta is daily air temperature (°C), t is time step, and a1~a6 are model parameters that need to be determined during the calibration process.
For each lake, the available observed water temperatures are divided into two parts, around 75% for model calibration, and around 25% for model validation. To evaluate model performance, two widely used metrics are employed, namely root mean square error (RMSE) and the Nash–Sutcliffe efficiency coefficient (NSE). Once the models are sufficiently calibrated and the model performance is acceptable, the calibrated models will be used to project the daily water temperatures of lakes to 2100 considering two scenarios, namely SSP245 and SSP585.
R M S E = 1 n i = 1 n ( T o T s ) 2
N S E = 1 i = 1 n ( T o T s ) 2 i = 1 n ( T o T m ) 2
where To is observed water temperature (°C), Ts is the modeled water temperature (°C), Tm is the average value of To, and n is the number of samples.

2.3.2. Reconstruction of Water Temperature with the Application of Satellite Data

In addition to the air2water model, the reconstruction of the water temperature of Lake Miedwie also utilized a second source of data, specifically the Landsat Collection 2 Level-2 [23]. This dataset includes publicly available information on the Earth’s surface temperature in degrees Kelvin. Detailed information about the data used and the method of processing them to calculate the Earth’s surface temperature is presented in the instructions developed by the USGS for Landsat 4–7 and Landsat 8–9 [24,25]. The surface temperature data contained in Collection 2 are available for the period from 1982 to 2024 (Table 1). The Collection 2 data were released in 2020 by the USGS Earth Resources Observation and Science (EROS) Center. Data from Collection 2 have improved geometric accuracy and radiometric calibration compared with previous Collection 1 products. Due to the fact that the hydrological station performing in situ lake water temperature measurements is located within the shoreline range, it is not possible to obtain surface temperature data (Landsat 4–9) for this location specifically. That results in spatial resolution of the Landsat surface temperature science product. Pixels within the range of the hydrological station where the lake water temperature is measured also extend over land or shoreline vegetation. Therefore, surface temperature data were collected from a point with coordinates Lon 14.918857, Lat 53.347478, located 120 m towards the center of the lake in relation to the lake hydrological station. Similar procedures were used in an earlier study focused on lake water thermometry using Landsat 8 satellite data [12]. Water temperatures were determined based on all available data from Landsat 4, 5, 7, 8, and 9, when this point was not covered by clouds. This process resulted in a total of 900 surface temperature results, from which measurements from November to March were discarded due to possible ice cover on the lake. Consequently, 682 observations were retained for analysis. A script created in Google Earth Engine (GEE) [26] was used to gather the data, which were obtained in the form of a *.csv file arranged as follows: acquisition date; cloud cover; satellite name; and water temperature. Additionally, to determine water temperatures, GEEcode developed by Ermida et al. [27] was utilized, updated for use with Collection 2 and to determine the emissivity of snow and water surfaces. In this method, the surface water temperature is computed using the Statistical Mono-Window (SMW) algorithm developed by the Climate Monitoring Satellite Application Facility (CM-SAF). The results were obtained in the form of a *.csv file.
In the first stage of acquiring water temperature data for Lake Miedwie, quality statistics collected in the QA_PIXEL file were used, calculated based on image data and cloud mask information for the scene. In the second stage, the data underwent filtering based on water temperature results obtained from in situ measurements for Lake Lubie. For this purpose, minimum and maximum temperature values were calculated for each day from 1 April to 31 October based on long-term data from 1972 to 2023. These boundary values were used to filter the data in such a way that temperatures lower than the minimum values and higher than the maximum values were excluded from the analysis. The results obtained in this manner were compared with in situ measurement data from the period 1984–1989 and with results obtained from the air2water model for the period 1990–2023.
Additionally, satellite data were used to present the seasonal variability of surface water temperature in Lake Miedwie. For this purpose, visualizations were created using Landsat Collection 2 Level-2 data by selecting one image from each month of the year 2022 for the period from April to October.

2.3.3. Air Temperature Prediction

This study employed multi-model ensembles (MMEs) utilizing the Bayesian model averaging (BMA) technique to obtain the surface air temperature data from multiple datasets from the Global Climatic Model (GCM) based on Coupled Model Intercomparison Project Phase 6 (CMIP6), as listed in Table 2. These datasets encompassed both historical and projected climate data. MMEs are a highly efficient method for reducing the uncertainty associated with individual models.
The NorESM2-MM dataset encompasses the time period from the pre-industrial era to the year 2100. It runs on a grid resolution of roughly 1 × 1 degree latitude/longitude, with a particular focus on marine biogeochemical processes. This model [28] highlights the significance of marine biogeochemical processes within its framework. The MPI-ESM1-2-HR dataset provides climate forecasts with a high level of detail on a 0.5-degree grid spanning the years 1850 to 2100. This dataset demonstrates exceptional proficiency in establishing connections between atmospheric observations and ocean models [29]. The EC-Earth3 dataset encompasses the time period from 1850 to future forecasts beyond 2100. It has a grid resolution of approximately 1 × 1 degree and is specifically designed to simulate the climate and variability of Europe [30]. The AWI-CM-1-1-MR dataset, which spans from the late 20th century to 2100, offers valuable insights into polar processes with a resolution of around 1 degree [31]. The BCC-CSM2-MR model covers a wide range of time periods, starting from 1850 and extending up to 2100 with a projected change of around 1 degree. It provides comprehensive information on the specific climate characteristics of East Asia and the West Pacific region [32]. The MRI-ESM2-0 dataset provides accurate simulations of atmospheric and oceanic dynamics, with a grid resolution ranging from 0.5 to 1.5 degrees [33]. The GFDL-ESM4 dataset encompasses the time period from the pre-industrial era to 2100. This dataset effectively combines the carbon cycle with physical climate aspects, hence offering valuable insights into the mechanisms of carbon feedback [34]. The CESM2-WACCM dataset encompasses time periods ranging from 1850 to 2100, featuring atmospheric resolution of up to 0.25 degrees. It includes a comprehensive model of the entire atmosphere [35]. The CMCC-CM2-SR5 encompasses a range of historical periods spanning from the 20th century to future possibilities in the 21st century, with a specific emphasis on the climate of the Mediterranean region [36].
To ensure consistency and facilitate direct comparisons, we employed a bilinear interpolation method to regrid various climate model datasets onto a standardized 1 × 1 degree latitude/longitude grid. Our study involved multiple datasets, including NorESM2-MM, MPI-ESM1-2-HR, EC-Earth3, AWI-CM-1-1-MR, BCC-CSM2-MR, MRI-ESM2-0, GFDL-ESM4, CESM2-WACCM, and CMCC-CM2-SR5, which span different time periods and geographic focuses. Bilinear interpolation calculates the interpolated value at a target grid point using values from the four nearest neighboring grid points on the source grid. Specifically, the interpolated value V′ at the target grid (xt, yt) is computed as follows:
V x t , y t = i = 1 4 w i · V ( x s i · y s i )
w i = 1 x · y · ( x d i s t i · y d i s t i )
V ( x x i · y x i ) represents the values of the variable at the four nearest neighboring grid points to the source grid that surrounds the target point (xt, yt), while the calculation of the weights, denoted as wi, is based on the distance that exists between each of the four locations and the target point. The variables denoted as Δx and Δy correspond to the intervals between adjacent points on the source grid. The variables x d i s t i and y d i s t i denote the horizontal and vertical distances, correspondingly, between the ith adjacent point and the target grid point. For instance, the NorESM2-MM dataset, with a grid resolution of roughly 1 × 1 degree, aligns directly with our standardized grid, while the MPI-ESM1-2-HR dataset, which has a 0.5-degree grid, is upscaled to match the 1 × 1 degree resolution using bilinear interpolation. Similarly, datasets with coarser resolutions, such as MRI-ESM2-0 (0.5 to 1.5 degrees), are regridded to achieve consistency. This interpolation process ensures that all datasets, regardless of their original resolution, are uniformly regridded, enabling a direct comparison of climate variables across different models and enhancing the robustness of our analyses.
The study employs Bayesian Model Averaging (BMA) as proposed by Hoeting et al. [37] to integrate data derived from several global climate models (GCMs). The Bayesian Model Analysis (BMA) is a statistical framework that facilitates the integration of predictions generated by several models, employing principles from probability theory. The underlying premise of this approach is the theorem of Bayes, which can potentially effectively be expressed mathematically within the framework of model averaging in the form that follows:
p M i y = p y M i · p ( M i ) p ( y )
p M i y denotes the posterior probability of model Mi given the data y, p y M i denotes the likelihood of the data y under model Mi, p ( M i ) denotes the prior probability of model Mi, and p(y) denotes the marginal likelihood of the data y and serves as a normalizing constant.
To determine the BMA prediction distribution for a new data y ~ , a weighted average is taken from the predictive distributions of each model. The equation incorporates the predictive distribution, denoted as p y ~ M i , which signifies the probability distribution of the new data y ~ to the model M i . Furthermore, the weight for model M i , denoted as p M i y , is computed under the previously elucidated methodology. It can be expressed mathematically in the form that follows:
p ( y ~ y i = 1 M p y ~ M i · p M i y )
The climate model ensemble provides a consolidated BMA prediction for the surface air temperature (T) at a particular grid point (x, y) and time (t) as follows:
T B M A x , y , t = i = 1 M w i x , y , t · T i ( x , y , t )
The prediction provided by model i is denoted as Ti (x, y, t), whereas the BMA weight assigned to model i at grid point (x, y) and time t is represented as wi (x, y, t). The weight in question is defined by the posterior probability p(Miy), which is specific to the location and time in question. The determination of uncertainty in the combined prediction of the BMA model is based on the standard deviation of the predictive distribution of the BMA model as follows:
σ B M A x , y , t = i = 1 M w i x , y , t · ( T i x , y , t T i x , y , t ) 2
σBMA (x, y, t) represents the degree of uncertainty at the given coordinates (x, y) and time t on the grid. The degree of uncertainty in the aggregated forecast and the level of agreement across the models can be inferred from the variability of the posterior distribution. The weighting factors from the BMA method for each model are listed in Table 3.
We utilized the nearest neighbor methodology [38] to choose surface air temperature data for a particular geographic location after applying the MME method. This approach involves calculating the distance between the target point and each location within the grid, and then choosing the closest grid point’s data as representative of the target location to ensure the selection of the most representative grid point. The notion of the nearest neighbor was articulated through the calculation of the distance between the target point and every location inside the grid.
The Euclidean distance d between a target point with geographic coordinates (xtarget, ytarget) and each grid point (xi, yi) in the dataset can be calculated as follows:
d i = ( x i x t a r g e t ) 2 + ( y i y t a r g e t ) 2
In this study, we applied a bias correction method to adjust the outputs of air temperature data from the MME method to better match observed historical climate data. This process corrects both the mean and the temporal variability in the GCM outputs, ensuring that the corrected data reflects the observed climate patterns more accurately. The corrected outputs are then used for future climate projections.
For the historical period up to 2014, we used observed temperature data and air temperature data from the MME method to calculate the necessary bias correction factors. Specifically, we computed the mean and standard deviation of the observed temperature data ( O R E F , σ O , R E F ) as well as the mean and standard deviation of the air temperature data from the MME method ( T R E F and σ T , R E F ). The bias correction was applied as follows [39,40]:
T B C t = O R E F + σ O , R E F σ T , R E F ( T R A W t T R E F )
where T B C t is the bias-corrected GCM output, O R E F is the mean of the observed data for the reference period, σ O , R E F is the standard deviation of the observed data for the reference period, T R A W t is the raw GCM output, T R E F is the mean of the GCM data for the reference period, and σ T , R E F is the standard deviation of the GCM data for the reference period.
This correction formula was applied to the GCM outputs for the future periods 2015–2023 and 2024–2100. By adjusting the raw GCM outputs based on the observed data’s mean and variability, the corrected data more closely align with historical observations. To evaluate the performance of the bias-corrected data, we calculated several statistical metrics for the period 2015–2023, including the correlation coefficient, root mean squared error (RMSE), and mean absolute error (MAE). The correlation coefficient measures the strength and direction of the linear relationship between observed and corrected temperatures, with values closer to 1 indicating a stronger positive correlation. RMSE and MAE measure the average magnitude of the errors and the average absolute difference between observed and corrected temperatures, respectively, with lower values indicating better performance.

2.3.4. Analysis of Trends of Changes in Water Temperature in Lake Miedwie

The study of long-term changes in air and water temperatures was conducted for different time periods. Initially, the analysis covered the entire period for which data were obtained from the air2water model, i.e., from 1990 to 2023. Subsequent analyses were performed for four 30-year periods. Firstly, the analysis was conducted for the past period from 1994 to 2023, and secondly for the near, mid, and far future periods, respectively from 2024 to 2053, 2047 to 2076, and 2071 to 2100. Future temperature changes were analyzed separately for the SSP245 and SSP585 scenarios. In order to determine the directions of long-term changes in average annual air and water temperatures in Lake Miedwie, non-parametric Mann–Kendall [41] and Sen [42] tests were used. In order to test the existence of a monotonic trend using the Mann–Kenadal test, two hypotheses were formulated. The null hypothesis (H0) states that the data show no monotonic trend, while the alternative hypothesis (HA) states that there is a monotonic trend in the data over time. To determine the magnitude of changes in air and water temperatures, the nonparametric Sen’s slope test was used. Sen’s test is insensitive to the presence of outliers in the dataset that significantly differ from the general population. Moreover, a trend-free prewhitening procedure was used to remove autocorrelation from the data series [43]. Meanwhile, to detect discontinuities in the data series, the nonparametric Pettitt test was applied. The Pettitt test allows for the detection of individual change points in time series data [44]. For the analysis of the directions and magnitudes of changes using the Mann–Kendall and Sen’s tests, the modifiedmk package developed by Patakamuri and O’Brien [45] was utilized. For change point detection, the trend package developed by Pohlert [46] was utilized.

3. Results

3.1. Historical Reconstruction of Water Temperature

The average annual water temperature of Lake Miedwie during the period covered by field measurements was 9.3 °C, ranging from 8.5 °C (in 1980 and 1984) to 10.4 °C (in 1989). On a monthly basis, February was the coldest month (2.1 °C), and July was the warmest (18.1 °C). Due to the cessation of stationary measurements in 1989, further thermal characterization of Lake Miedwie was conducted by reconstructing water temperatures based on air temperatures up to 2023. For this purpose, the Air2Water model was used, which achieved calibration period results with an RMSE of 1.34 and an NSE of 0.95, and validation period results of 1.54 and 0.94, respectively. A similar analysis was also conducted for Lake Lubie (located 63 km away), which, due to ongoing continuous observations to the present day, serves as a reference base. The results of these statistics for Lake Lubie were similar, confirming the model’s assumptions as correct. The reconstructed water temperature of Lake Miedwie during the period with available air temperatures (1990–2023) was close to that for the period covered by field monitoring (1972–1989), with an average temperature of 9.6 °C, and the lowest and highest temperatures being 8.2 °C (in 1996) and 10.4 °C (in 2019 and 2020), respectively. Similarly, during this period, February remained the coldest month (2.1 °C), and July the warmest (18.5 °C). Long-term changes from 1972 to 2023 showed an increase in the average annual temperature by 0.20 °C per decade (statistically significant, p 0.05) (Table 4). Over the last 30 years, the thermal changes in the waters of Lake Miedwie were even greater, at 0.31 °C per decade (statistically significant, p 0.05). In the same period, the recorded increase in air temperature was 0.41 °C per decade (p 0.05).

3.2. Future Water Temperature

The analysis of the bias-corrected outputs compared to the observed temperatures for the period 2015–2023 demonstrates the effectiveness of the bias correction method applied in this study. The plot compares observed temperatures (black line, Figure 2) with bias-corrected GCM outputs for both scenarios SSP245 and SSP585. Several statistical metrics, including the correlation coefficient, Root Mean Squared Error (RMSE), and Mean Absolute Error (MAE), were used to evaluate the performance of the corrected data.
The spatial resolution of the GCMs used in this study, with a grid cell size of approximately 100 km, poses a limitation when studying Miedwie Lake, which has a size of about 10 to 20 km. This coarse resolution leads to the loss of fine-scale spatial details and the inability to accurately represent the lake’s dynamics and surrounding microclimate, potentially resulting in errors in temperature, precipitation, and other climate-related variables. The resolution of GCM outputs can be improved through dynamical downscaling based on Regional Climate Model (RCM) methods; however, these methods are computationally demanding and complex. Generally, RCMs can provide more detailed data and offer higher resolution. However, they are computationally intensive and require significant resources, making them less feasible for extensive temporal and spatial studies [47]. In order to address these constraints, we implemented a bias correction methodology that was inspired by the methodologies outlined by Hawkins et al. [39] and Ho et al. [40]. Bias correction methods are advantageous in this context because they can effectively adjust systematic biases in GCM outputs without the need for high computational resources [48,49]. This method aligns the GCM outputs more closely with observed data by adjusting both the mean and temporal variability. The efficacy of this methodology is illustrated by statistical metrics: the correlation coefficients for SSP245 and SSP585 are 0.79 and 0.77, respectively, indicating a robust positive correlation between the observed and corrected temperatures. The RMSE values for SSP245 and SSP585 are 6.28 and 6.62, respectively, indicating that the corrected outputs closely align with the observed data. Additionally, the MAE values for SSP245 and SSP585 are 4.95 and 5.17, respectively, further corroborating the accuracy of the bias-corrected data. Future research could explore advanced downscaling techniques, high-resolution models, and integrated approaches that combine multiple methods to maximize data accuracy and resolution.
Further analysis is based on the prediction of future climate changes (until the end of the 21st century) related to air temperature. Two commonly considered scenarios were adopted: SSP245 and SSP585. The study of the directions and pace of changes in water temperature of the analyzed lake is presented in Table 4. The data analysis was divided into three 30-year periods: near, mid, and far future (Figure 3). The results indicate that under the SSP245 scenario, an increase in water temperatures in Lake Miedwie is still projected, but the changes are expected to decrease over successive 30-year periods, amounting to 0.17, 0.09, and 0.08 °C per decade respectively. Considering the SSP585 scenario, an acceleration in the warming of the lake’s waters is expected in subsequent 30-year periods. For the period 2024–2053, the increase in water temperature at Lake Miedwie is forecasted to be 0.30 °C per decade. In the following 30-year periods, 2047–2076 and 2071–2100, the rate of change will be higher than in the last three decades (1994–2023), amounting to 0.35 and 0.42 °C per decade, respectively. Similar trends are observed for air temperatures. Assuming the SSP245 scenario, it is predicted that the pace of increase in air temperatures will be at a lower level than what occurred in the years 1994–2023 (Table 4). However, under the SSP585 scenario, the rate of increase in air temperatures in the near, mid, and far future will be higher than that observed in the years 1994–2023. Furthermore, the rate of change in air temperatures in subsequent 30-year periods is expected to increase.

3.3. Application of Satellite Images in Water Temperature Analysis

In the next phase, an attempt was made to use satellite images to reconstruct water temperatures in Lake Miedwie using data recorded by Landsat satellites 4, 5, 7, 8, and 9. For the period from April to October (ice-free) from 1984 to 2023, 682 observations were collected. After filtering using historical daily minimum and maximum temperatures, 511 observations remained for analysis. Comparing these results with in situ measurements (from the period 1984–1989) and results from the air2water model (from the period 1990–2023), deviations ranged from −4.9 to 5.2 °C. The RMSE value was 2.41 °C, MAE was 2.04 °C, and the coefficient of determination was 0.84 (Figure 3a). Using the SMW algorithm and code developed by Ermida et al. [27] and performing filtering based on long-term data for the reference Lake Lubie (as above), 516 observations remained for analysis. Comparing these results with in situ measurements (from the period 1984–1989) and results from the air2water model (from the period 1990–2023), deviations ranged from −5.5 to 4.9 °C. The RMSE value was 2.28 °C, MAE was 1.88 °C, and the coefficient of determination was 0.83 (Figure 4b).
The spatial variability in surface water temperatures in 2022, from April to October, is presented in Figure 5.
These results show that the surface water temperature in Lake Miedwie exhibits relatively low spatial variability in subsequent months. In the year analyzed, the water temperature within the reservoir ranged from approximately 3.5 °C to 27.8 °C. The relative homogeneity of the surface water temperature is due to the shape and catchment conditions. In the first case, there are no bays where water could heat up more quickly. Additionally, the absence of obstacles such as peninsulas means that water circulation evens out the thermal characteristics of the surface layer of water. Regarding the catchment conditions, it should be noted that the lake is not fed by a river with a high flow rate, which would alter the thermal system in the area of its mouth. Nonetheless, slight differences are visible (April, June), where the shallowest zones are characterized by higher water temperatures than the central sectors of the lake.

4. Discussion

The ability to interpret long-term, and therefore reliable, changes that various components of the natural environment undergo requires the most accurate data possible. This becomes particularly significant in the era of global warming for assessing the scale of current and future transformations. According to Obrębska-Starkel and Starkel [50], direct data provide the most accurate information about climate changes. However, it is not always possible to rely on such compilations, which can be supplemented by various methodological approaches. Rapid progress in the field of modeling and the continuous development of remote sensing techniques create new possibilities for interpreting hydrological processes. This approach was used for one of the largest lakes in Poland, where field measurements of water temperature were completed in the late 1980s. The results obtained in the study enrich the widely developed contemporary research direction concerning changes in the thermal regime of lakes. Referring to specific situations, it was noted, for instance, that in the case of Lake Nam Co (China), the water temperature during the summer increased by 0.52 °C per decade, and at the same time, the onset of stratification accelerated and its duration extended [51]. Between 1984 and 2019, two morphometrically different lakes in Italy showed an increase in water temperature by 0.05 °C/year (Lake Bracciano, larger and deeper) and 0.06 °C/year (Lake Martignano, smaller and shallower) [52]. The largest of the lakes in Poland (Śniardwy Lake), over the last five decades, experienced a temperature rise of 0.44 °C per decade [53]. In this context, Lake Miedwie is similar, i.e., from 1972 to 2023 (during which period there were stationary measurements and taking into account model data based on air temperature), there was a statistically significant increase in average annual water temperature, at a rate of 0.20 °C per decade. At the same time, it can be observed that despite maintaining a uniform direction of changes similar to that of Lake Śniardwy, the scale is lower. This difference may result from the morphometry of both lakes, with Miedwie being significantly deeper (the average depth is 5.8 m and 19.3 m, respectively). The inflow of cooler water masses from the extensive hypolimnion zone may cause the observed difference. Such a situation has been presented, among others, in the case of Lake Hańcza (the deepest in Poland), where the rate of water temperature increase was lower compared to other cases [54].
The approach based on satellite imagery in the considered case yielded worse results than modeling using the air2water model. It can be therefore said that with access to appropriate input data, i.e., continuous recordings of air temperature as a variable explaining water temperature trends, the use of such methodology is more justified. Nonetheless, it is important to emphasize that in the absence of reliable data, namely data related to areas not covered by any field observations (including mountainous areas), satellite data are crucial for understanding the thermal conditions of lakes [55,56,57]. Another issue is the impact of local conditions on the thermodynamics of lakes, where air–water relationships can be obscured by the morphometric parameters of the lakes or their watershed characteristics [58]. Furthermore, issues of spatial variability in water temperature are significant, where point field measurements may not reflect the dynamics of changes within a single lake. For example, such an illustration for Lake Drawsko (located 85 km from Lake Miedwie) showed significant spatial variability, with the warmest water in narrow, sheltered bays. In the case of Lake Miedwie, its relatively compact shape (without bays and peninsulas) influenced the absence of major thermal contrasts in different months (Figure 5).
The rapid progress of civilization necessitates a priority understanding of future changes, which also pertains to aquatic ecosystems. So far, studies on this matter do not look optimistic, and the recorded increase in water temperature is evident. By the end of the 21st century, the water temperature in Southeast Asian lakes will rise by 0.38 °C in the best-case scenario, while the most extreme simulation indicates an increase of 2.29 °C [59]. Lake temperature projections are based on one-dimensional, vertically resolved hydrodynamic simulations for 29 lakes. A prediction of water temperature for 29 lakes in Switzerland by the end of the 21st century showed an average increase of 3.3 °C (RCP8.5) [60]. Referring to predictions of water temperature changes in Lake Miedwie, two commonly used scenarios were applied (SSP245 and SSP585). According to them, the analyzed lake will be warmer by the end of the 21st century (compared to the initial period of study) by an average of 1.8 °C and 3.2 °C, respectively. Piccolroaz et al. [61] conducted studies for 25 lakes and found that the average increase in water temperature in the discussed region will be 0.15 °C/decade in the case of the RCP4.5 pathway and 0.34 °C/decade for the RCP8.5 pathway. Therefore, in the first case, Lake Miedwie achieved the same result, and in the second, it is lower, which may be due to the presence of shallower lakes in the set analyzed by Piccolroaz et al. [61]. In the case of polymictic lakes, due to multiple circulations, the water exhibits similar thermal characteristics from the surface to the bottom. The SSP585 pathway shows a continuous increase in temperature (Figure 3b), consequently leading to greater differences between shallow and deep lakes.
These are therefore significant changes in the thermal regime that will be reflected in the course of most processes occurring in the lake. As noted by Chen et al. [62], water temperature is a fundamental parameter of water quality in lakes and reservoirs. Water quality issues are particularly important for Lake Miedwie given that it serves as a source of drinking water for Szczecin. As Stepanowska et al. [63] note, Lake Miedwie was classified as eutrophic at the turn of the 20th and 21st centuries and has been classified as mesotrophic in recent years. Wilmański [64] addresses plankton issues, noting that current blooms are less intense. Maintaining the lake’s waters in a mesotrophic state depends on maintaining the appropriate level of various components that determine water quality. According to the data obtained, meeting such conditions will be more difficult in the future due to changes in one of the fundamental characteristics of water, namely its temperature. For example, based on the RCP4.5 scenario, an analysis of water temperature changes in the case of Fairy Lake (USA) showed that there will be an extension of the duration of cyanobacterial blooms by more than 39% [65]. Simulations of changes in water temperature and dissolved oxygen in the case of Lake Mississippi show that an increase in the former by 9.2% and 18.4% (RCP4.5 and RCP8.5) will be noticeable alongside a decrease in oxygen by 6.2% and 14.3%, respectively [66]. Reduced oxygen levels will impact the limited possibilities for the oxidation of organic matter [67]. In this light, the institutions managing Lake Miedwie will face new challenges, where alongside actions aimed at reducing pollution at the watershed level, it will be necessary to consider threats to water quality on a macroscale level—associated with progressing global warming.

5. Conclusions

The article reconstructs and subsequently predicts changes in the water temperature of one of the largest lakes in Poland, Lake Miedwie, for which stationary measurements had ceased over 30 years ago. Utilizing contemporary models that define the interaction between water and the atmosphere, as well as thermal satellite imagery, it was possible to obtain results characterizing the thermal regime of Lake Miedwie, including an analysis of long-term changes. Using field measurements of water temperature and subsequently modeling results based on air temperature and Landsat images (1972–2023), it was determined that over the past five decades, the analyzed lake has been warming at an average rate of 0.19 °C per decade. Further analysis concerning future changes indicated that the increase in water temperature by 2100 will average 0.15 °C per decade (SSP245) and as much as 0.25 °C per decade (SSP585). Both contemporary and projected future changes are significant, which, given the fundamental importance of water temperature, will impact the ecosystem discussed in the article. Given that Lake Miedwie serves as a source of drinking water, issues concerning the maintenance or improvement of water quality in the era of global warming will be crucial for water management authorities.

Author Contributions

Conceptualization, M.P.; methodology, M.P., M.S., S.Z. and T.A.; software, M.S., S.Z., T.A. and H.L.; validation, M.S., S.Z. and T.A.; formal analysis, M.S., S.Z. and T.A.; investigation, M.P., M.S., S.Z. and T.A.; resources, M.P., M.S. and T.A.; data curation, M.P., M.S. and T.A.; writing—original draft preparation, M.P., M.S., S.Z., T.A., H.L., K.S.-P., S.J. and L.W.; writing—review and editing, M.P., M.S., S.Z., T.A., H.L., K.S.-P., S.J. and L.W.; visualization, M.P., M.S., S.Z. and T.A.; supervision, M.P.; project administration, M.P.; funding acquisition, M.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

www.dane.imgw.pl (accessed on 27 June 2024).

Acknowledgments

We would like to thank Copernicus Climate Change Service, Climate Data Store, (2021): CMIP6 climate projections. Copernicus Climate Change Service (C3S) Climate Data Store (CDS) for providing the CMIP6 dataset.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jaffe, S.; Qian, S.S.; Mayer, C.M.; Kocovsky, P.M.; Gouveia, A. Assessing the probability of grass carp (Ctenopharyngodon idella) spawning in the Sandusky River using discharge and water temperature. J. Great Lakes Res. 2024, 50, 102303. [Google Scholar] [CrossRef]
  2. Yan, Z.; Feng, C.; Xu, Y.; Wang, J.; Huang, N.; Jin, X.; Wu, F.; Bai, Y. Water temperature governs organophosphate ester dynamics in the aquatic food chain of Poyang Lake. Environ. Sci. Ecotechnol. 2024, 21, 100401. [Google Scholar] [CrossRef] [PubMed]
  3. Haddout, S.; Priya, K.L.; Casila, J.C.C.; Hoguane, A.M.; Ljubenkov, I. Modelling of depth profiles of the water temperature in the Lake Sidi Ali (Morocco). Int. J. River Basin Manag. 2024, 22, 23–30. [Google Scholar] [CrossRef]
  4. Yang, K.; Yu, Z.; Luo, Y.; Yang, Y.; Zhao, L.; Zhou, X. Spatial and temporal variations in the relationship between lake water surface temperatures and water quality—A case study of Dianchi Lake. Sci. Total Environ. 2018, 624, 859–871. [Google Scholar] [CrossRef] [PubMed]
  5. Issak, D.J.; Horan, D.L.; Wollrab, S.P. Air temperature data source affects inference from statistical stream temperature models in mountainous terrain. J. Hydrol. X 2024, 22, 100172. [Google Scholar] [CrossRef]
  6. Twardosz, R.; Walanus, A.; Guzik, I. Warming in Europe: Recent Trends in Annual and Seasonal temperatures. Pure Appl. Geophys. 2021, 178, 4021–4032. [Google Scholar] [CrossRef]
  7. Ptak, M.; Sojka, M.; Kozłowski, M. The increasing of maximum lake water temperature in lowland lakes of Central Europe: Case study of the Polish Lakeland. Ann. De Limnol. -Int. J. Limnol. 2019, 55, 11. [Google Scholar] [CrossRef]
  8. Piccolroaz, S.; Zhu, S.; Ladwig, R.; Carrea, L.; Oliver, S.; Piotrowski, A.P.; Ptak, M.; Shinohara, R.; Sojka, M.; Woolway, R.I.; et al. Lake water temperature modeling in an Era of climate change: Data sources, models, and future prospects. Rev. Geophys. 2024, 62, e2023RG000816. [Google Scholar] [CrossRef]
  9. Duda, L.; Duklas, K. Retencja Jeziorna Jeziora Miedwie w Latach 1902/03-1939/40; Szczecińskie Towarzystwo Naukowe; Wydział Nauk Matematycznych i Technicznych; Państwowe Wydawnictwo Naukowe Oddział w Poznaniu: Szczecin, Poland, 1970; p. 24. [Google Scholar]
  10. Wloczyk, C.; Richter, R.; Borg, E.; Neubert, W. Sea and lake surface temperature retrieval from Landsat thermal data in Northern Germany. Int. J. Remote Sens. 2006, 27, 2489–2502. [Google Scholar] [CrossRef]
  11. Heddam, S.; Ptak, M.; Zhu, S. Modelling of daily lake surface water temperature from air temperature: Extremely randomized trees (ERT) versus Air2Water, MARS, M5Tree, RF and MLPNN. J. Hydrol. 2020, 588, 125130. [Google Scholar] [CrossRef]
  12. Dyba, K.; Ermida, S.; Ptak, M.; Piekarczyk, J.; Sojka, M. Evaluation of Methods for Estimating Lake Surface Water Temperature Using Landsat 8. Remote Sens. 2022, 14, 3839. [Google Scholar] [CrossRef]
  13. Choiński, A. Katalog Jezior Polski; Wydawnictwo Naukowe UAM: Poznań, Poland, 2006. [Google Scholar]
  14. Ptak, M.; Choiński, A.; Strzelczak, A.; Targosz, A. Disappearance of Lake Jelenino since the end of the XVIII century as an effect of anthropogenic transformations of the natural environment. Pol. J. Environ. Stud. 2013, 22, 191–196. [Google Scholar]
  15. Pasławski, Z. Zarys limnologii fizycznej jeziora Miedwie. Prace PIHM 1969, 96, 57–72. [Google Scholar]
  16. Woś, A. Klimat Polski w Drugiej Połowie XX Wieku; Wydawnictwo Naukowe UAM: Poznań, Poland, 2010. [Google Scholar]
  17. Piccolroaz, S.; Toffolon, M.; Majone, B. A simple lumped model to convert air temperature into surface water temperature in lakes. Hydrol. Earth Syst. Sci. 2013, 17, 3323–3338. [Google Scholar] [CrossRef]
  18. Zhu, S.; Ptak, M.; Yaseen, Z.M.; Dai, J.; Sivakumar, B. Forecasting surface water temperature in lakes: A comparison of approaches. J. Hydrol. 2020, 585, 124809. [Google Scholar] [CrossRef]
  19. Zhu, S.; Piotrowski, A.P.; Ptak, M.; Napiorkowski, J.J.; Dai, J.; Ji, Q. How does the calibration method impact the performance of the air2water model for the forecasting of lake surface water temperatures? J. Hydrol. 2021, 597, 126219. [Google Scholar] [CrossRef]
  20. Wang, W.; Shi, K.; Wang, X.; Zhang, Y.; Qin, B.; Zhang, Y.; Woolway, R.I. The impact of extreme heat on lake warming in China. Nat. Commun. 2024, 15, 70. [Google Scholar] [CrossRef] [PubMed]
  21. Toffolon, M.; Piccolroaz, S.; Majone, B.; Soja, A.M.; Peeters, F.; Schmid, M.; Wüest, A. Prediction of surface temperature in lakes with different morphology using air temperature. Limnol. Oceanogr. 2014, 59, 2185–2202. [Google Scholar] [CrossRef]
  22. Zhu, S.; Ptak, M.; Sojka, M.; Piotrowski, A.P.; Luo, W. A simple approach to estimate lake surface water temperatures in Polish lowland lakes. J. Hydrol. Reg. Stud. 2023, 48, 101468. [Google Scholar] [CrossRef]
  23. Landsat Collection 2 Level-1 Product Courtesy of the U.S. Geological Survey. Available online: https://www.usgs.gov/centers/eros/science/usgs-eros-archive-landsat-archives-landsat-4-5-thematic-mapper-collection-2 (accessed on 27 June 2024).
  24. Landsat 8-9 Collection 2 (C2) Level 2 Science Product (L2SP) Guide Version 5.0; United States Geological Survey: Asheville, NC, USA, 2023; pp. 1–43.
  25. Landsat 4-7 Collection 2 (C2) Level 2 Science Product (L2SP) Guide Version 4.0; United States Geological Survey: Asheville, NC, USA, 2021; pp. 1–44.
  26. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  27. Ermida, S.L.; Soares, P.; Mantas, V.; Göttsche, F.-M.; Trigo, I.F. Google Earth Engine open-source code for Land Surface Temperature estimation from the Landsat series. Remote Sens. 2020, 12, 1471. [Google Scholar] [CrossRef]
  28. Seland, Ø.; Bentsen, M.; Olivié, D.; Toniazzo, T.; Gjermundsen, A.; Graff, L.S.; Debernard, J.B.; Gupta, A.K.; He, Y.C.; Kirkevåg, A.; et al. Overview of the Norwegian Earth System Model (NorESM2) and key climate response of CMIP6 DECK, historical, and scenario simulations. Geosci. Model Dev. 2020, 13, 6165–6200. [Google Scholar] [CrossRef]
  29. Müller, W.A.; Jungclaus, J.H.; Mauritsen, T.; Baehr, J.; Bittner, M.; Budich, R.; Bunzel, F.; Esch, M.; Ghosh, R.; Haak, H.; et al. A higher-resolution version of the Max Planck Institute Earth System Model (MPI-ESM1.2-HR). J. Adv. Model. Earth Syst. 2018, 10, 1383–1413. [Google Scholar] [CrossRef]
  30. Döscher, R.; Acosta, M.; Alessandri, A.; Anthoni, P.; Arneth, A.; Arsouze, T.; Bergmann, T.; Bernadello, R.; Bousetta, S.; Caron, L.P.; et al. The EC-Earth3 Earth system model for the Coupled Model Intercomparison Project 6. Geosci. Model Dev. 2022, 15, 2973–3020. [Google Scholar] [CrossRef]
  31. Semmler, T.; Danilov, S.; Gierz, P.; Goessling, H.F.; Hegewald, J.; Hinrichs, C.; Koldunov, N.; Khosravi, N.; Mu, L.; Rackow, T.; et al. Simulations for CMIP6 with the AWI climate model AWI-CM-1-1. J. Adv. Model. Earth Syst. 2020, 12, e2019MS002009. [Google Scholar] [CrossRef]
  32. Xin, X.-G.; Wu, T.-W.; Zhang, J.; Zhang, F.; Li, W.-P.; Zhang, Y.-W.; Lu, Y.-X.; Fang, Y.-J.; Jie, W.-H.; Zhang, L.; et al. Introduction of BCC models and its participation in CMIP6. Adv. Clim. Change Res. 2019, 15, 533. [Google Scholar]
  33. Yukimoto, S.; Kawai, H.; Koshiro, T.; Oshima, N.; Yoshida, K.; Urakawa, S.; Tsujino, H.; Deushi, M.; Tanaka, T.; Hosaka, M.; et al. The Meteorological Research Institute Earth System Model version 2.0, MRI-ESM2.0: Description and basic evaluation of the physical component. J. Meteorol. Soc. Jpn. Ser. II 2019, 97, 931–965. [Google Scholar] [CrossRef]
  34. Dunne, J.P.; Horowitz, L.W.; Adcroft, A.J.; Ginoux, P.; Held, I.M.; John, J.G.; Krasting, J.P.; Malyshev, S.; Naik, V.; Paulot, F.; et al. The GFDL Earth System Model version 4.1 (GFDL-ESM 4.1): Overall coupled model description and simulation characteristics. J. Adv. Model. Earth Syst. 2020, 12, e2019MS002015. [Google Scholar] [CrossRef]
  35. Lauritzen, P.H.; Nair, R.D.; Herrington, A.; Callaghan, P.; Goldhaber, S.; Dennis, J.; Bacmeister, J.; Eaton, B.; Zarzycki, C.; Taylor, M.A.; et al. NCAR release of CAM-SE in CESM2.0: A reformulation of the spectral element dynamical core in dry-mass vertical coordinates with comprehensive treatment of condensates and energy. J. Adv. Model. Earth Syst. 2018, 10, 1537–1570. [Google Scholar] [CrossRef]
  36. Cherchi, A.; Fogli, P.G.; Lovato, T.; Peano, D.; Iovino, D.; Gualdi, S.; Masina, S.; Scoccimarro, E.; Materia, S.; Bellucci, A.; et al. Global mean climate and main patterns of variability in the CMCC-CM2 coupled model. J. Adv. Model. Earth Syst. 2019, 11, 185–209. [Google Scholar] [CrossRef]
  37. Hoeting, J.A.; Madigan, D.; Raftery, A.E.; Volinsky, C.T. Bayesian model averaging: A tutorial. Stat. Sci. 1999, 14, 382–417. [Google Scholar]
  38. Cover, T.M.; Hart, P.E. Nearest Neighbor Pattern Classification. IEEE Trans. Inf. Theory 1967, 13, 21–27. [Google Scholar] [CrossRef]
  39. Hawkins, E.; Osborne, T.M.; Ho, C.K.; Challinor, A.J. Calibration and bias correction of climate projections for crop modelling: An idealised case study over Europe. Agric. For. Meteorol. 2013, 15, 19–31. [Google Scholar] [CrossRef]
  40. Ho, C.K.; Stephenson, D.B.; Collins, M.; Ferro, C.A.; Brown, S.J. Calibration strategies: A source of additional uncertainty in cli-mate change projections. Bull. Am. Meteorol. Soc. 2012, 93, 21–26. [Google Scholar] [CrossRef]
  41. Kendall, M.G.; Stuart, A. The Advanced Theory of Statistics, 3rd ed.; Charles Griffin Ltd.: Cheshire, UK, 1968. [Google Scholar]
  42. Gilbert, R.O. Statistical Methods for Environmental Pollution Monitorin; Van Nos-trand Reinhold Co.: New York, NY, USA, 1987. [Google Scholar]
  43. Yue, S.; Pilon, P.; Phinney, B.; Cavadias, G. The influence of autocorrelation on the ability to detect trend in hydrological series. Hydrol Process 2002, 16, 1807–1829. [Google Scholar] [CrossRef]
  44. Pettitt, A.N. A non-parametric approach to the changepoint problem. Appl. Stat. 1979, 28, 126–135. [Google Scholar] [CrossRef]
  45. Patakamuri, S.K.; O’Brien, N. 31 October 2022. Modified Versions of Mann Kendall and Spearman’s Rho Trend Tests, Version 1.6. Available online: https://cran.r-project.org/web/packages/modifiedmk/modifiedmk.pdf (accessed on 27 June 2024).
  46. Pohlert, T. 10 October 2023. Non-Parametric Trend Tests and Change-Point Detection, Version 1.1.6. Available online: https://cran.r-project.org/web/packages/trend/trend.pdf (accessed on 27 June 2024).
  47. Christensen, J.H.; Christensen, O.B. A summary of the PRUDENCE model projections of changes in European climate by the end of this century. Clim. Chang. 2007, 81, 7–30. [Google Scholar] [CrossRef]
  48. Ehret, U.; Zehe, E.; Wulfmeyer, V.; Warrach-Sagi, K.; Liebert, J. HESS Opinions “Should we apply bias correction to global and regional climate model data?”. Hydrol. Earth Syst. Sci. 2012, 16, 3391–3404. [Google Scholar] [CrossRef]
  49. Teutschbein, C.; Seibert, J. Bias correction of regional climate model simulations for hydrological climate-change impact studies: Review and evaluation of different methods. J. Hydrol. 2012, 456, 12–29. [Google Scholar] [CrossRef]
  50. Obrębska-Starkel, B.; Starkel, L. The greenhouse effect and global environmental changes. Zeszyty IGiPZ PAN. Inst. Geogr. I Przestrz. Zagospod. (IGiPZ) PAN 1991, 4, 1–72. [Google Scholar]
  51. Huang, L.; Wang, J.; Zhu, L.; Ju, J.; Daut, G. The Warming of Large Lakes on the Tibetan Plateau: Evidence from a Lake Model Simulation of Nam Co, China, During 1979–2012. J. Geophys. Res. Atmos. 2017, 122, 13–095. [Google Scholar] [CrossRef]
  52. De Santis, D.; Del Frate, F.; Schiavon, G. Analysis of Climate Change Effects on Surface Temperature in Central-Italy Lakes Using Satellite Data Time-Series. Remote Sens. 2022, 14, 117. [Google Scholar] [CrossRef]
  53. Ptak, M.; Sojka, M.; Nowak, B. Effect of climate warming on a change in thermal and ice conditions in the largest lake in Poland–Lake Śniardwy. J. Hydrol. Hydrodyn. 2020, 68, 260–270. [Google Scholar] [CrossRef]
  54. Ptak, M.; Sojka, M. Different reactions to climate change as observed in selected elements of hydrological regime of the deepest lake on the Central European Plain (Lake Hańcza). Hydrol. Sci. J. 2021, 66, 1083–1095. [Google Scholar] [CrossRef]
  55. Zhang, G.; Yao, T.; Xie, H.; Qin, J.; Ye, Q.; Dai, Y.; Guo, R. Estimating surface temperature changes of lakes in the Tibetan Plateau using MODIS LST data. J. Geophys. Res. Atmos. 2014, 119, 8552–8567. [Google Scholar] [CrossRef]
  56. Wan, W.; Zhao, L.; Xie, H.; Liu, B.; Li, H.; Cui, Y.; Ma, Y.; Hong, Y. Lake Surface Water Temperature Change Over the Tibetan Plateau From 2001 to 2015: A Sensitive Indicator of the Warming Climate. Geophys. Res. Lett. 2018, 45, 11177–11186. [Google Scholar] [CrossRef]
  57. Liu, B.; Wan, W.; Xie, H.; Li, H.; Zhu, S.; Zhang, G.; Wen, L.; Hong, Y. A long-term dataset of lake surface water temperature over the Tibetan Plateau derived from AVHRR 1981–2015. Sci. Data 2019, 6, 48. [Google Scholar] [CrossRef]
  58. Ptak, M.; Sojka, M.; Choiński, A.; Nowak, B. Effect of environmental conditions and morphometric parameters on surface water temperature in Polish lakes. Water 2018, 10, 580. [Google Scholar] [CrossRef]
  59. Virdis, S.G.P.; Kongwarakom, S.; Juneng, L.; Padedda, B.M.; Shrestha, S. Historical and projected response of Southeast Asian lakes surface water temperature to warming climate. Environ. Res. 2024, 247, 118412. [Google Scholar] [CrossRef]
  60. Michel, A.; Råman Vinnå, C.L.M.; Bouffard, D.; Epting, J.; Huwald, H.; Schaefli, B.; Schmid, M.; Wüest, A.J. Evolution of Stream and Lake Water Temperature under Climate Change; Hydro-CH2018 Project; FOEN: Bern, Switzerland, 2021. [Google Scholar] [CrossRef]
  61. Piccolroaz, S.; Zhu, S.; Ptak, M.; Sojka, M.; Du, X. Warming of lowland Polish lakes under future climate change scenarios and consequences for ice cover and mixing dynamics. J. Hydrol. Reg. Stud. 2021, 34, 100780. [Google Scholar] [CrossRef]
  62. Chen, R.; Wang, B.; Qian, L.; Luo, Y.; Ma, M. Seasonal Variation of Water Quality of Taiping Lake in Eastern China. J. Risk Anal. Crisis Response 2024, 14, 70–79. [Google Scholar] [CrossRef]
  63. Stepanowska, K.; Biernaczyk, M.; Machula, S.; Kubiak, J. Struktura po3owów rybackich na tle czynników abiotycznych wód w jeziorze Miedwie. Komunikaty Rybackie 2014, 5, 1–4. [Google Scholar]
  64. Wilmański, K. Rozwój technologii oczyszczania wody z użyciem pylistego węgla aktywnego w Polsce. Inżynieria I Ochr. Sr. 2016, 19, 265–275. [Google Scholar]
  65. Bhatti, M.; Singh, A.; McBean, E.; Vijayakumar, S.; Fitzgerald, A.; Siwierski, J.; Murison, L. Climate Change Impacts on Water Temperatures in Urban Lakes: Implications for the Growth of Blue Green Algae in Fairy Lake. Water 2024, 16, 587. [Google Scholar] [CrossRef]
  66. Yaghouti, M.; Heidarzadeh, N.; Nakhaei, N.; Ulloa, H.N. The impacts of climate change on thermal stratification and dissolved oxygen in the temperate, dimictic Mississippi Lake, Ontario. Ecol. Inform. 2023, 75, 10208. [Google Scholar] [CrossRef]
  67. Ptak, M.; Nowak, B. Variability of oxygen-thermal conditions in selected lakes in Poland. Ecol. Chem. Eng. 2016, 23, 639–650. [Google Scholar] [CrossRef]
Figure 1. Location of the research object.
Figure 1. Location of the research object.
Remotesensing 16 02753 g001
Figure 2. The bias-corrected General Circulation Model (GCM) outputs for scenarios SSP245 (blue dashed line) and SSP585 (red dashed line) compared to the observed temperatures for the period 2015–2023.
Figure 2. The bias-corrected General Circulation Model (GCM) outputs for scenarios SSP245 (blue dashed line) and SSP585 (red dashed line) compared to the observed temperatures for the period 2015–2023.
Remotesensing 16 02753 g002
Figure 3. Prediction of annual water temperature in Lake Miedwie by the end of the 21st century based on the SSP245 (a) and SSP585 (b) scenarios.
Figure 3. Prediction of annual water temperature in Lake Miedwie by the end of the 21st century based on the SSP245 (a) and SSP585 (b) scenarios.
Remotesensing 16 02753 g003
Figure 4. Differences between in situ water temperature measurements and the results from the air2water modeling, as well as temperatures obtained from Landsat Collection 2 Level-2 using the SMW algorithm according to the GEE code developed by Ermida et al. [27] for Lake Miedwie (a) and Lake Lubie (b).
Figure 4. Differences between in situ water temperature measurements and the results from the air2water modeling, as well as temperatures obtained from Landsat Collection 2 Level-2 using the SMW algorithm according to the GEE code developed by Ermida et al. [27] for Lake Miedwie (a) and Lake Lubie (b).
Remotesensing 16 02753 g004
Figure 5. Spatial variability in thermal conditions of water of Lake Miedwie in the period from April to October 2022.
Figure 5. Spatial variability in thermal conditions of water of Lake Miedwie in the period from April to October 2022.
Remotesensing 16 02753 g005
Table 1. Access to Landsat data.
Table 1. Access to Landsat data.
SatelliteInstrumentOperation PeriodNumber of Acquired Data PointsNumber of Used Data Points (April–October)
Landsat 4Thematic Mapper (TM)August 1982 to December 199388
Landsat 5Thematic Mapper (TM)March 1984 to May 2012343280
Landsat 7Enhanced Thematic Mapper Plus (ETM+)July 1999 to April 2022323233
Landsat 8Operational Land Imager (OLI)April 2013 to present195142
Landsat 9Operational Land Imager 2 (OLI-2)February 2022 to present3119
Table 2. List of Global Climate Models (GCMs) and their corresponding institutions.
Table 2. List of Global Climate Models (GCMs) and their corresponding institutions.
GCMInstitutions
NorESM2-MMNorwegian Climate Centre (NCC)
MPI-ESM1-2-HRMax Planck Institute for Meteorology (MPI-M)
EC-Earth3EC-Earth-Consortium
AWI-CM-1-1-MRAlfred Wegener Institut (AWI)
BCC-CSM2-MRBeijing Climate Center (BCC)
MRI-ESM2-0Meteorological Research Institute (MRI)
GFDL-ESM4Geophysical Fluid Dynamics Laboratory (GFDL) of the
National Oceanic and Atmospheric Administration (NOAA)
CESM2-WACCMNational Center for Atmospheric Research (NCAR)
CMCC-CM2-SR5Euro-Mediterranean Centre on Climate Change (CMCC) Foundation
Table 3. BMA weighting factors for each model.
Table 3. BMA weighting factors for each model.
ModelSSP245SSP585
AWI-CM-1-1-MR0.110.12
BCC-CSM2-MR0.080.08
CESM2-WACCM0.120.11
CMCC-CM2-SR50.090.10
EC-Earth30.120.13
GFDL-ESM40.100.09
MPI-ESM1-2-HR0.120.12
MRI-ESM2-00.140.13
NorESM2-MM0.120.12
Table 4. Water temperature changes (Miedwie Lake) and air temperature changes (Szczecin).
Table 4. Water temperature changes (Miedwie Lake) and air temperature changes (Szczecin).
PeriodMann–Kendal and Sen TestsPettitt Test
Sz-Valuep-ValueSen Slope Value
(°C per Decade)
Change Pointp-Value
Water temperature
1972–20235454.300.0000.2019880.003
1994–20231803.210.0010.3120140.004
2024–2053 a2374.270.0000.1720360.006
2047–2076 a1332.410.0160.09 0.214
2072–2100 a1502.740.0060.0820790.035
2024–2053 b2785.020.0000.3020370.000
2047–2076 b3175.680.0000.3520610.000
2072–2100 b3446.150.0000.4220840.000
Air temperature
1972–20236805.210.0000.41280.000
1994–20232013.570.0000.49200.007
2024–2053 a2254.000.0000.28110.007
2047–2076 a1242.190.0280.14 0.195
2072–2100 a1552.750.0060.14 0.077
2024–2053 b2815.000.0000.49110.000
2047–2076 b2985.300.0000.55140.000
2072–2100 b3436.100.0000.67130.000
a SSP245; b SSP585.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ptak, M.; Zhu, S.; Amnuaylojaroen, T.; Li, H.; Szyga-Pluta, K.; Jiang, S.; Wang, L.; Sojka, M. Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland). Remote Sens. 2024, 16, 2753. https://doi.org/10.3390/rs16152753

AMA Style

Ptak M, Zhu S, Amnuaylojaroen T, Li H, Szyga-Pluta K, Jiang S, Wang L, Sojka M. Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland). Remote Sensing. 2024; 16(15):2753. https://doi.org/10.3390/rs16152753

Chicago/Turabian Style

Ptak, Mariusz, Senlin Zhu, Teerachai Amnuaylojaroen, Huan Li, Katarzyna Szyga-Pluta, Sun Jiang, Li Wang, and Mariusz Sojka. 2024. "Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland)" Remote Sensing 16, no. 15: 2753. https://doi.org/10.3390/rs16152753

APA Style

Ptak, M., Zhu, S., Amnuaylojaroen, T., Li, H., Szyga-Pluta, K., Jiang, S., Wang, L., & Sojka, M. (2024). Utilizing Multi-Source Datasets for the Reconstruction and Prediction of Water Temperature in Lake Miedwie (Poland). Remote Sensing, 16(15), 2753. https://doi.org/10.3390/rs16152753

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