3. Results and Discussion
In the monitored period of 2018 and 2019, water temperature of Suceava River ranged from 0 to 28 °C. It became supercooled during the 2018/2019 winter (minimum: −0.19 °C at Mihoveni (upstream) and −0.13 °C at Tişăuţi (downstream)), when frazil ice was forming, and warm during the summers (maximum: 28.57 °C at Mihoveni and 28.65 °C at Tişăuţi), due to heat waves in the atmosphere (
Supplementary Figure S2). As the extreme values already revealed, Suceava River was generally warmer downstream of Suceava city (average: 11.97 °C upstream and 12.4 °C downstream). One can observe especially in the November–February months, that when water temperature upstream is near the freezing point, the same parameter downstream still exhibits diurnal cycles. These cycles in streamflows are created by similar cycles in air temperature and other mechanisms [
33] and are still present downstream of Suceava city during cold days due to the heat input that the river received inside the city.
In order to find out if this increase in temperature is caused by the warmer air above Suceava city, we compared water and air temperature during an entire calendar year, 2019 (
Figure 3). We selected two atmospheric datasets for this comparison. The first dataset is represented by the air temperature time series from Suceava Weather Station. Data from this monitoring station are less influenced by Suceava city due to its location outside the built-up area and are representative for the semi-natural environment that surrounds Suceava city and which represents the Suceava River catchment in Suceava Plateau.
The second dataset is represented by time series obtained by averaging hourly values from the selected monitoring stations located in the urban floodplain of Suceava River (SV2, ITC, TRN, AMB, ACT—
Table 1). This dataset is representative of the urban air of Suceava, as it was recorded in a complex residential, industrial and commercial area.
As expected, the thermal amplitude is higher in the air and the water temperature follows the air temperature (exceptions are the days with persistent negative air temperature, when water temperature stays around 0 °C due to the high energy loss necessary for freezing). In 2019, the average air temperature at Suceava Weather Station was 10.35 °C, while in the floodplain it was 11.05 °C. From the thermal difference of these stations (0.7 °C), only 0.5 °C might be explained by the vertical temperature gradient in the air due to the fact that the weather station is placed on a plateau at approximately 83 m above the mean elevation of the monitoring stations in the floodplain. Therefore, a surplus of 0.2 °C exists in the middle of the city as a result of the urban surfaces and activities. Water temperature was 11.54 °C upstream of Suceava city (Mihoveni) and 11.97 °C downstream (Tişăuţi). Studies on urban river temperature increase often find increases above 1 °C downstream of urban areas (e.g., 4–5 °C increase at some stations in North Carolina [
34]). The amplitude of the temperature increase of our studied river is more similar to those found by Zeiger and Hubbart [
35] in Missouri (0.2–0.7 °C). The difference between the mean values of the monitoring stations upstream and downstream in 2019 was 0.43 °C. The accuracy of the temperature sensors of the AT500 instruments used for measuring Suceava River water temperature is ±0.1 °C. The standard error of the mean is often used to assess the error of continuous measurements [
36,
37], and the standard error of the mean of Suceava River water temperature measurements is 0.085 °C for the upstream time series and 0.084 °C for the downstream time series. Therefore, we conclude that there is a relevant water temperature increase downstream.
We can observe that Suceava River is warmer than the air above the city at both moments, of entering and of exiting its metropolitan area. A calculus done for the warm semester (April–September, average value from hourly data), in order to exclude winter months with water persisting around 0 °C when the air temperature is frequently negative, has revealed a similar context: air at weather station = 16.67 °C, air in the floodplain = 17.7 °C, water upstream = 17.84 °C and water downstream = 18.19 °C. Suceava River is heated above the air temperature as a result of its warmer banks, which absorbed the solar radiation, and this reason is available for its tributaries too, whose shallower water is more prone to overheating. This heating is stronger in constructed areas (Suceava city in the floodplain is 0.7 °C warmer than at Suceava Weather Station) due to the numerous man-made structures (buildings, roads, channels) that transfer their excess heat to the water nearby. Small urban streams have a low inertia and are more susceptible to urban influence [
38].
During the communist era, a network of pipelines was built in order to transport hot water from the thermal power station towards the city inhabitants. Those pipelines are still used today not only for heating (during winter), but also for hot showers (in some homes). The pipelines transport warm water from the thermal power station (located in the floodplain, near the wastewater treatment plant) to various urban users. These pipelines are on land surface or underground and transfer some of their heat to the nearby soil and groundwater (hot water leaves the plant with 70–80 °C and arrives at the consumers, sometimes after more than 10 km of pipe, with 48 °C or less). The warm water is then discharged into the sewerage system in order to be collected at the wastewater treatment plant. The sewerage network often contains wastewater that is warmer than the groundwater. Warmer groundwaters diffusely recharge Suceava River and its urban tributaries in the floodplain (this recharge direction is proven by the small springs that appear along Suceava Riverbanks). No studies exist about the temperature of the urban groundwaters in our study area. Groundwater tends to be cooler than air in the summertime and warmer in wintertime, but even during summer, some groundwaters may become warmer because hot urban surfaces may transfer their heat through soils into phreatic waters.
It appears that the average heat transfer between Suceava River and Suceava city is from water to air after the water has been heated by the urban soil and constructions, but the analysis of the average diurnal profiles of water and air in our study area indicates that the water and air relationship is bidirectional. One can observe in
Supplementary Figure S3a,b, that the diurnal cycle of the air temperature in the floodplain has steeper slopes than at Suceava Weather Station and the moments of the maximum and minimum values occur earlier. During the late evening and the night, air is colder in the floodplain. Air temperature in the urban floodplain was warmest, compared to Suceava Weather Station air temperature, at 2 p.m. (2.46 °C difference), and air was colder by a maximum of 0.63 °C at 9 p.m. All these differences indicate that the air in the floodplain is strongly influenced by the urban land use (constructed surfaces heat up faster and cool more easily that a vegetated area). Similarly, after Suceava River flows through the city, the diurnal profile of Suceava River water temperature has steeper slopes and earlier moments of the maximum and minimum temperatures than upstream. Water downstream is warmer than upstream, especially in the afternoon when land and air become hotter. Water downstream was warmer than upstream with a maximum of 0.99 °C (4 p.m.) and a minimum of 0.05 °C (3 a.m.). Minima of the air temperature, at Suceava Weather Station and in the floodplain, and of water temperature, upstream and downstream, occur at 6 a.m., 5 a.m., 8 a.m. and 7 a.m., respectively. Moments of the average daily maxima for the same sites are 3 p.m., 2 p.m., 7 p.m. and 4 p.m., respectively. These values show that the peak position in the diurnal cycle of water has been strongly shifted towards the moment of the maximum air temperature because the heat gain during the day is an active process, in opposition to the passive process of heat loss during the night.
As previously observed, the heat transfer between air and water is bidirectional, both during an average day and during a year. This led to the question of whether this relationship alters the warming of Suceava River downstream. A simple statistical analysis revealed that 26.4% of the hourly values of Suceava River at the downstream station are actually lower than at the upstream station, and 11.2% of the days in 2019 had an average value that was lower downstream than upstream. We analysed the hourly differences in the diurnal regime of temperature during opposite seasons (summer and winter) and found that a constantly higher water temperature downstream is representative for winter, while, during summer night-time, Suceava River is colder downstream of the city (
Figure 4a,b).
During winter, the downstream–upstream thermal differences are lower than in summer. During the daytime, the heating of the active surface and of the air in the floodplain is much stronger than at Suceava Weather Station because of the built environment, but, during night-time, the floodplain becomes a space with frequent air temperature inversions of medium–low intensity (
Figure 4c,d). These inversions are able, during most days in summer, to decrease the river water temperature downstream. During both seasons, the water of Suceava River is a heating source for the atmosphere of Suceava city, especially in winter (during 3/4 of an average day—
Figure 4e,f).
Monthly average values in
Figure 5a shows that, during March–June, air temperature in the urban floodplain of Suceava River is warmer than the stream water in both stations (or similar to it, in April). This might be attributed to the melting snows in the mountain area of Suceava River catchment and to the important rainfalls specific to this period of the year—all these create the colder water of Suceava River in the first half of the year. It is to be taken into consideration that 2019 was a dry year, with an annual sum of precipitation of 535.03 mm at Suceava station and 486.5 mm at Iţcani station (in the floodplain). We can assume that, during years with average or excess precipitation, the impact of the urban heat island on Suceava River temperature is easier to detect.
Studies on other city–river pairs indicate that a river can be warmer than the urban air even during winters with frozen surface waters [
17]. Even small, linear water bodies such as a river may have a thermal effect on the surrounding environment, as proven by the 22 m wide UK river in the study of Hathway and Sharples [
39]. Suceava River is usually 50–60 m wide inside the city, during baseflow. A warm river enhances the urban heat island [
17].
On average, Suceava River receives heat from the urban air from ~9 a.m. until ~7 p.m. (
Figure 5b). During the daylight, Suceava River water temperature inside the city increases quickly, but, during the night, it decreases slowly due to the increasingly higher heat loss necessary to lose 1 °C at lower and lower temperatures (in this case, temperatures under 30 °C; see isobaric specific heat [
38]). In this way, during mid to late summer and during the autumn, the stream water gains energy from one day to another because the heat loss during the night-time is smaller than the energy gain during the daytime.
The most intense heat absorption from the atmosphere to the aquatic environment is done during the long summer daytimes by active evaporation processes from the water surface and the floodplain in an urban atmosphere drier than the surroundings. The atmosphere loses heat at the level of its basal layer—the heat is transferred to water through evaporation (
Supplementary Figure S1). The most intense heat transfer from the aquatic environment and floodplain to the urban atmosphere takes place during the long winter nights, foggy mornings and evenings, when the condensation or sublimation of excess vapor in the urban atmosphere of Suceava creates a caloric surplus given to the atmosphere. This contributes to mitigating night-time and morning radiative losses (
Supplementary Figure S4).
In both situations (both in the summer—through the remarkable intensity of the evaporation process—and in winter—through the long duration of the foggy/humid air interval in Suceava valley), the presence of water vapor in the atmosphere plays a very important and, at the same time, an ambivalent role in the thermo–caloric relationship between the urban atmosphere of Suceava city and the water of Suceava River. Similar strong thermal connections between water and air were also noted in previous studies [
40].
Boxplots of hourly water and air time series show the higher, urban variability of the air temperature in the constructed floodplain, especially in August (
Figure 6). In January, the interquartile interval of the water time series downstream of the city is bigger due to the less extreme negative temperatures in the floodplain air and the contribution of various heat sources. The observed distribution of air temperatures in 2019 is one of a warm year—the average temperature was 10.35 °C at Suceava Weather Station (in comparison, at the same station, the average temperature during 1950–2012 was 7.86 °C [
27]).
The average diurnal profiles of water and air vary in their shape from month to month and from one monitoring station to another (
Figure 7). The annual variation of the monthly average diurnal profiles is similar for the air time series (
Figure 7a,b), but the variation in water time series is different (
Figure 7c,d). Suceava River downstream recorded a high variability of the moment when the diurnal maximum value occurs. That moment varied from 3 p.m. to 8 p.m., and the urban air recorded only a fluctuation of ±1 h of the moment of the diurnal maximum. One can observe that the annual variation of Suceava River water has shifted from an original variation towards a variation that is specific to that of the local air temperature. In water time series upstream and downstream, the biggest differences in the position of maxima values occur during the warm semester. During winter, probably due to the atmospheric forcing that moves water temperatures to 0 °C for a long time, the differences in the hourly position of maxima values is minimal (the same is true for minima values).
A shape change of the stream water average diurnal profile from upstream to downstream of Suceava city similar to that observed in our study was also revealed by a previous study [
27] that analysed data from a shorter time interval (October 2012–January 2013): the hour of the maximum temperature shifted from the evening (upstream) to the afternoon inside the city and also downstream of Suceava city.
The scalograms of the continuous wavelet transform analysis indicate that the diurnal cycle occurs in more consecutive days in the monitoring station from Tişăuţi—downstream (
Figure 8). This persistence is sustained by heat transfer from the city into the river—this supplementary heat appears to attenuate the impact of external factors that tend to disrupt the repetition of a uniform diurnal cycle. The time interval with frequent significant (0.95 confidence level) diurnal cycles is longer downstream. During the cold semester, the diurnal cycles are less significant, especially in winter when, upstream, a diurnal periodicity may lack for many consecutive days. Also, the higher power of the diurnal periodicity downstream indicates that the diurnal cycle is easily shaped by the high thermal amplitude (day–night difference).
The wavelet coherence analysis of Suceava River water temperature at Mihoveni and Tişăuţi (
Figure 9a) reveals very frequent similar diurnal oscillations in both monitoring stations (that co-vary in a non-linear way). The diurnal periodicities are generally in phase, as indicated by the phase arrows pointing right. The scalograms of the wavelet coherence analysis of one station in air and the other station in water (
Figure 9b,c) show that there is a better covariance of the diurnal periodicities in the urban floodplain air and Suceava River water downstream, at Tişăuţi, than between the air relevant to the upstream catchment and Suceava River upstream, at Mihoveni. It is interesting to observe that the coherence between air and water downstream is stronger at the diurnal band of frequencies than the coherence between water upstream and water downstream. This indicates that the shape and timing of the diurnal water profile at Tişăuţi is imposed by the diurnal distribution of the air temperature, rather than the characteristics of water at Mihoveni. However, at lower frequencies, the relationship between water time series in both monitoring stations is stronger, as indicated by the power spectrum and the fact that the continuous significant coherent area (dark red area topped by the thin black line of the statistical significance area limit) starts from ~5 days (128 h), instead of ~21 days (512 h). The wavelet analysis is useful for the analysis of the non-linear covariation of the diurnal cycle. The higher variability of the diurnal cycle does not have an impact only on a scale of a few consecutive days, but also, often, at a monthly scale. Thus, for example, the hourly position of minima values in the average diurnal profiles of the air temperature at Suceava Weather Station and in the floodplain and of the water temperature at Mihoveni (upstream) and Tişăuţi (downstream) in July/August was 6/7 a.m., 5/6 a.m., 8/8 a.m. and 6/7 a.m., respectively.
In order to find punctual sources of heat input for Suceava River, we analysed the temperature of some of its urban tributaries and the discharge characteristics of the wastewater treatment plant of Suceava city (the wastewater treatment plants are well-known sources of thermal pollution in rivers, worldwide [
41] and in Suceava city [
27]). With the exception of Şcheia River, the urban tributaries of Suceava River recorded water temperatures significantly greater than that of Suceava River (
Figure 10) in the case study time interval.
During September 20–December 31, 2019, the average water temperature of (in upstream–downstream order; from hourly measurements) Suceava upstream, Şcheia, Cetăţii Creek, Collector and Suceava downstream was 8.54, 8.36, 11.64, 11.45 and 9.19 °C, respectively. Unlike the other urban tributaries, Şcheia flows mostly through a less urbanized landscape, before entering the city. This is why it has a temperature similar to that of Suceava River before entering the city. Şcheia lowers the temperature of Suceava River, as shown by measurements done ~0.4 km downstream, at Iţcani station, where the annual water temperature is 12.4 °C and the average of the case study time interval is 7.9 °C (calculated from daily data, used with precautions because the daily values are based on measurements done only twice per day).
The temperature of the other minor streamflows that discharge into Suceava River indicate how sensitive these waters are to the urban alteration of their streambed and catchment, to the urban heat island and to the disposal of wastewaters. Concrete embankments are frequent inside the city—these are an aquatic grey infrastructure that removes the natural processes from the rivers, including the interaction with the subsurface flow through the hyporheic exchange [
6].
In 2019, the wastewater treatment plant of Suceava city discharged approximately 0.34 m3/s of wastewaters with an average of 16.31 °C. During September 20–December 31, 2019, it discharged 0.32 m3/s with an average of 15.6 °C. The effluent of the wastewater treatment plant is colder than Suceava River only during summer months, when the tap water that becomes wastewater is colder than Suceava River (the main water supply for Suceava city is from Berchişeşti groundwater, located at the plateau border with the mountain area). The average temperature of the discharged wastewater was 20.71 °C during summer (with ~1.16 °C lower than Suceava River at Tişăuţi in the same time interval) and 12.24 °C during winter (with ~10.2 °C higher than Suceava River).
During September 20–December 31, 2019, we estimated the flow rate of Suceava River at Iţcani as 3.94 m3/s. The average temperature of the tributaries that Suceava River receives between Mihoveni and Tişăuţi is computed from the average of the three measured minor streamflows (10.48 °C). This average is combined in a weighted average with the wastewater temperature (cumulative discharge is 0.57 m3/s from 0.25 m3/s streamflow + 0.32 m3/s wastewater), resulting in a combined flow with a temperature of 13.35 °C (the discharge of tributaries is estimated as half of the 0.5 m3/s estimate by using the discharge of Şcheia during September 20–December 31, 2019, which is approximately half of the average specific to this time interval of the year). The combined flow temperature is added in a weighted average (discharge used for weighting) to Suceava River temperature upstream (8.54 °C) and results in a new temperature of Suceava River downstream of 9.14 °C.
The measured temperature in the field was 9.19 °C, resulting in this calculus/the surface water input into Suceava River explaining 0.6 °C (92.3%) from the total 0.65 °C difference between Tişăuţi and Mihoveni. The remaining heat gain has to be attributed to the urban heat island (directly) and the groundwater input. It is to be remembered that the final heat gain is the result of complex processes of gaining or losing heat (e.g., through water input or exchange with atmosphere) and the estimated values are specific to the studied period of the end of 2019 only. Elements not taken into consideration (e.g., temperature of the raindrops) may have some role in the variability of Suceava River temperature. The indirect impact of the urban heat island on Suceava River temperature is exerted through urban tributaries. The thermal impact of the tributaries represents 31.75% (part of the above-mentioned 92.3%). It remains to be discovered how much of the temperature increase of tributaries over Suceava River temperature is caused only by the urban heat island, but a rough estimate is at least 50%. We can also estimate that the 0.05 °C left unexplained in the calculus of water temperature increase from upstream to downstream is also caused at least 50% by the urban heat island (direct influence). The measured impact of the wastewaters on the water temperature of Suceava River is not a surprise because the municipal treatment plants are considered as being probably the most important heat source for urban rivers [
41]. Moreover, there was a constant increase in the temperature of wastewaters in cities reported, e.g., Tokyo [
42].
The discharge of the wastewater treatment plant varies considerably, and the highest values are created by stormwaters collected during important rainfalls (
Figure 11a). This is due to the fact that an exclusively pluvial drainage exists only in some small areas of the city. Most stormwaters are collected in the main (combined) sewerage and flow through the wastewater treatment plant. Therefore, the calculated impact of the wastewater is also representative of the impact of the stormwater runoff on the water quality of Suceava River. Temperature surges caused by stormwaters that washed warm city surfaces are difficult to detect, as most pluvial waters enter the combined sewers. Any detectable temperature surges caused by urban stormwater runoff must alter the average diurnal profile of Suceava River water temperature difference between the downstream station (Tişăuţi) and the upstream station (Mihoveni) (
Figure 11b). An analysis of the hourly temperature difference between these stations shows that the standard deviation is 0.76 °C and the maximum occurs at 2 p.m. In order to identify possible water temperature surges, we analysed the days with a temperature difference greater than 1 °C and with atmospheric precipitation greater than 10 mm (this amount generates consistent runoff, especially during high-intensity rains). From the 15 times when the rainfalls exceeded 10 mm/day, only 2 times generated a runoff that strongly impacted the stream water temperature downstream of Suceava city (
Figure 11c,d). During the surge event in August 2019, the temperature difference was 1.64 °C and the temperature increase induced by urban stormwater merged the diurnal peaks of 2 days into a larger peak with a duration of 24 h. The surge in July 2019 lasted 13 h and recorded the maximum temperature difference of 2019 (4.14 °C, which is 1.1 °C higher than the second highest value). The temperature surges occurred during warm summer days, similar to many cases reported in other studies [
3,
35]. Our surge values are similar to those of Rice et al. [
34] (1.9–3.27 °C), Zeigler and Hubbart [
35] (4 °C) or Nelson and Palmer [
3] (2–7 °C). The surge durations in our study are greater than those reported previously (up to 10–10.5 h) [
34,
35].
A correlation matrix (with Pearson coefficients) that took variables as daily averages during the entire year of 2019 into consideration (
Supplementary Table S1) revealed very good correlations (>0.95) between Suceava River water temperature and air temperature or wastewater temperature. A good correlation in this linear estimate is between water temperature and solar radiation (>0.6). The good correlation of the solar radiation indicates the potentially high impact of other, unknown/unmeasured parameters (such as water turbidity) on the variability of water temperature and highlights the limits of any mathematical model.
ANCOVA tests, applied to some relevant parameters in order to find out which natural and anthropogenic factors influence the water temperature of Suceava River at Mihoveni (upstream) and Tişăuţi (downstream), revealed the results displayed in
Table 2. The goodness-of-fit (R
2) of the linear model in explaining the upstream water temperature is 0.917, and it is 0.998 for the downstream case (in both cases, the risk in assuming that the selected explanatory variables explain the stream water temperature is <0.0001).
For the upstream case, the parameters that contribute more to the water temperature are air temperature, relative humidity and water level, but the contribution of the other parameters is not negligible. For the downstream case, the most important factors are water temperature from the upstream station and wastewater temperature, seconded by solar radiation and water level at the downstream monitoring station.
In order to detect non-linear dependencies between multiple parameters, we applied a multiple wavelet coherence analysis. This uses the same characteristics of the previously discussed (simple) wavelet coherence and is a reliable method for the analysis of the non-linearly dependent processes [
30]. The scalograms indicate the wavelet power of the first parameter explained by the other parameters (
Figure 12). As in the case of ANCOVA, at the upstream station, the strongest links were detected between water temperature upstream and air temperature at Suceava Weather Station + water level upstream. However, at the downstream station, the most important factors linked to the evolution of the water temperature are air temperature in the floodplain and water temperature upstream (the latter has a stronger role than the temperature of the wastewaters, which was detected as relevant by ANCOVA). One can observe that the explanatory variables better explain the explained variable at the downstream station (the dark red areas are dominant and indicate stronger coherence, especially at the weekly and monthly scales).
An analysis of the scaled correlation between various time series was also conducted (
Figure 13). The Pearson correlation coefficient was computed for every day (24 h/values generate 1 correlation coefficient and the result is a list of 365 correlation coefficients, 1 for each calendar day). We observed that, during 2019, the scaled correlation between water temperature at Tişăuţi (downstream) and air temperature in the floodplain is better (the average of the 365 correlation coefficients is 0.6188) than the scaled correlation between water temperature at Mihoveni (upstream) and air temperature at Suceava Weather Station (0.3407). The classical correlation between the entire time series gave an irrelevant difference (0.8819 versus 0.8898); therefore, the scaled correlation is more useful in the process of finding differences between upstream and downstream monitoring stations.
Figure 13b shows that the scaled correlation between the downstream/floodplain pair and air temperature in the floodplain is better when the air temperature is higher, for almost every month.
Figure 13c indicates that the same scaled correlation is lower when water level is high, especially during the warm season, when high waters occur frequently.
All data indicate that Suceava River temperature increases downstream of Suceava city because of numerous factors whose importance vary in time and in space. This increase is acquired effectively in the lower 2/3 of the distance between Mihoveni and Tişăuţi (11.6 km straight line) and is caused by Suceava city, including its urban heat island.
The increase in water temperature has an impact on other parameters of the river. For example, lower dissolved oxygen concentrations were reported downstream of Suceava city, at Tişăuţi (compared to Mihoveni), partly caused by the increase in water temperature [
31]. The lowered concentrations of dissolved oxygen are certainly diminishing the self-purification capacity of Suceava River that was attributed partly to this element [
43,
44].