Changes in the Water Temperature of Rivers Impacted by the Urban Heat Island: Case Study of Suceava City

: Cities alter the thermal regime of urban rivers in very variable ways which are not yet deciphered for the territory of Romania. The urban heat island of Suceava city was measured in 2019 and its impact on Suceava River was assessed using hourly and daily values from a network of 12 water and air monitoring stations. In 2019, Suceava River water temperature was 11.54 ◦ C upstream of Suceava city (Mihoveni) and 11.97 ◦ C downstream (Ti¸s ă u¸ti)—a 3.7% increase in the water temperature downstream. After the stream water passes through the city, the diurnal thermal proﬁle of Suceava River water temperature shows steeper slopes and earlier moments of the maximum and minimum temperatures than upstream because of the urban heat island. In an average day, an increase of water temperature with a maximum of 0.99 ◦ C occurred downstream, partly explained by the 2.46 ◦ C corresponding di ﬀ erence between the urban ﬂoodplain and the surrounding area. The stream water diurnal cycle has been shifted towards a variation speciﬁc to that of the local air temperature. The heat exchange between Suceava River and Suceava city is bidirectional. The stream water diurnal thermal cycle is statistically more signiﬁcant downstream due to the heat transfer from the city into the river. This transfer occurs partly through urban tributaries which are 1.94 ◦ C warmer than Suceava River upstream of Suceava city. The wavelet coherence analyses and ANCOVA (analysis of covariance) prove that there are signiﬁcant (0.95 conﬁdence level) causal relationships between the changes in Suceava River water temperature downstream and the ﬂuctuations of the urban air temperature. The complex bidirectional heat transfer and the changes in the diurnal thermal proﬁles are important to be analysed in other urban systems in order to decipher in more detail the observed causal relationships. Minima of the air temperature, at Suceava Weather Station and in the ﬂoodplain, 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.


Introduction
An increase in water temperature of rivers all over the world due to global warming is depicted in numerous studies [1,2]. At the same time, the increasing urban population worldwide leads to a greater and more territorially focused impact of cities on the environment. This impact, enhanced by the current climate changes, is exerted on urban rivers through multiple anthropogenic stressors such as the increasing imperviousness of catchments or the diminishing of areas with riparian vegetation [3,4]. The land use changes associated with the urbanization of an area (e.g., widespread paving) [4,5] lead to alterations in the thermal regime of the air in cities, and urban heat islands occur [6]. The warmer

Study Area
Suceava city is the centre of a metropolitan area of approximately 150,000 people and is located in the north-eastern part of Romania, in Suceava county. It is a city divided by Suceava River, which flows from NW towards SE and crosses the industrial area in the middle of the city. Terrain elevation in the city ranges from 265 m above sea level (a.s.l.), in the floodplain/south-eastern part of the city, to 435 m a.s.l., in the hills of the northern limits. Climate is temperate continental. Winters record small amounts of precipitations and frequent negative temperatures, while summers are warm and rainy, with frequent torrential rainfalls. Mean air temperature is 8.16 • C, whilst the average annual sum of precipitation is 621 mm.
Suceava River has an average flow rate of 16.87 m 3 /s at Iţcani station, located into the city [27]. It also has an annual flow regime imposed by the regional rainfalls and records frequent high waters (over 50 m 3 /s) in the late spring and early to mid-summertime interval. This river has well mixed waters due to its average maximum depth of under 1 m of water and due to numerous transversal hydraulic infrastructures, that generate hydraulic jumps (e.g., downstream of the Iţcani and Burdujeni bridges). The main stream water monitoring stations used in this study are positioned on Suceava River upstream and downstream of Suceava city ( Figure 1, Table 1). The other 3 monitoring stations are located on 3 small streamflows, typical for the tributaries received by Suceava River between Mihoveni (upstream) and Tişăuţi (downstream). The cumulative streamflow of Suceava River tributaries in the study area is roughly 0.5 m 3 /s. The Collector streamflow is a semi-artificial streamflow which represents a mixture of urban rainwater and groundwater drainage and, in a small proportion, wastewater. It was created in the communist era as part of the measures to drain water from Suceava floodplain. Untreated wastewaters are also present in the flow ofŞcheia and Cetăţii Creek, included in our monitoring network, and in the other streamflows, as a result of illicit domestic wastewater discharges. The wastewater treatment plant of Suceava city is located on the left bank of Suceava River and, in 2019, it discharged an average of 0.34 m 3 /s of water. This is also one of the 5 stations used to collect air temperature and relative humidity in Suceava floodplain ( Figure 1, Table 1). These latter stations are part of a wider monitoring network (19 stations) that we used prior to our hydrological analysis in order to obtain maps of the urban heat island and its meteorological effects that help us to better understand our study area (Figure 2, Supplementary Figure S1).
Major and static sources of heat are a few and are represented by the thermoelectric plant (near the wastewater treatment plant), the cellulose factory and 5 plants that belongs to the food industry. Many (tens of thousands) minor static sources are represented by domestic water and air heaters (such as boilers). Tens of thousands of cars are a daily source of mobile thermal pollution. In the official build-up area of Suceava city, the industrial areas represent 11.16%, other buildings (mostly domestic) represent 24.37%, roads and pavements represent 6.14% and the remaining space is constituted of semi-natural spaces, such as private yards, parks or bare lands [27].
The urban heat island of Suceava is easily identifiable in summer ( Figure 2b) and produces a distinct area of lower relative humidity in winter ( Figure 2c).
Water 2020, 4, x FOR PEER REVIEW 3 of 24 streamflow which represents a mixture of urban rainwater and groundwater drainage and, in a small proportion, wastewater. It was created in the communist era as part of the measures to drain water from Suceava floodplain. Untreated wastewaters are also present in the flow of Şcheia and Cetăţii Creek, included in our monitoring network, and in the other streamflows, as a result of illicit domestic wastewater discharges. The wastewater treatment plant of Suceava city is located on the left bank of Suceava River and, in 2019, it discharged an average of 0.34 m 3 /s of water. This is also one of the 5 stations used to collect air temperature and relative humidity in Suceava floodplain ( Figure 1, Table  1). These latter stations are part of a wider monitoring network (19 stations) that we used prior to our hydrological analysis in order to obtain maps of the urban heat island and its meteorological effects that help us to better understand our study area ( Figure 2, Supplementary Figure S1). Major and static sources of heat are a few and are represented by the thermoelectric plant (near the wastewater treatment plant), the cellulose factory and 5 plants that belongs to the food industry. Many (tens of thousands) minor static sources are represented by domestic water and air heaters (such as boilers). Tens of thousands of cars are a daily source of mobile thermal pollution. In the official build-up area of Suceava city, the industrial areas represent 11.16%, other buildings (mostly domestic) represent 24.37%, roads and pavements represent 6.14% and the remaining space is constituted of semi-natural spaces, such as private yards, parks or bare lands [27].

Data
We used hourly and daily time series with different lengths from 2018 and 2019 at UTC+2. The year 2019 was selected for in-depth analyses due to the broader data availability for all parameters. Data used in this study is divided into data obtained from our own instruments and data provided by various Romanian institutions. Data from our own instruments is represented by measurements of water temperature (5 sites) and, secondarily, water level (we used here, for a complementary analysis, only level data measured at Tişăuţi with AT500 instruments); also, from the air environment, we obtained hourly values of temperature and relative humidity in 18 sites, of which 4 (located in Suceava floodplain- Table 1) were selected for further integration in the meteorologicalhydrological analysis and comparison. Water temperature at Mihoveni (upstream) and Tişăuţi (downstream) was recorded from June 28, 2018 until the end of 2019. Water level at Tişăuţi used here was recorded for the entire year of 2019. Water temperature of the other streamflows was recorded from September 20, 2019 until the end of 2019. The sensors were placed in areas with persistent water flow. In the case of Suceava River, sensors were permanently under at least a 0.33 m column of water; for the tributaries, the minimum water column had a height of 0.1 m (maximum possible during low flow). Our meteorological data were recorded during the entire year of 2019 and above mostly grassy active surfaces. Water temperature data of Suceava River at Mihoveni (upstream) and Tişăuţi (downstream) and air temperature and relative humidity data of Suceava floodplain are available at http://water.usv.ro/data.php.
Data from Romanian institutions is represented by data provided by ANAR (Romanian Waters National Administration-discharge and temperature data of Suceava River and atmospheric precipitation, both at Iţcani-daily values for the entire year of 2019), ACET Suceava (the company responsible for water supply and sewerage in Suceava metropolitan area-discharge and temperature data of Suceava wastewater treatment plant effluent-daily data of 2019), ANM (National Meteorological Administration-air temperature, relative humidity and precipitation at Suceava Weather Station-hourly and daily data of 2019) and ANPM (National Environment Protection Agency-air temperature, relative humidity and solar radiation at SV2 station-hourly and daily data for 2019).  Urban heat island (UHI) above Suceava metropolitan area as revealed by hourly air temperature measurements in 2019 during winter-(a) (December, January, February)-and summer months-(b) (June-August)-UHI represents the area most impacted by the urban heat island.

Instruments
From June 28 until September 28, 2018, water temperature data used in this study for the analysis of Mihoveni (upstream) and Tişăuţi (downstream) monitoring stations was measured with TruBlue 585 CTD instruments (accuracy: ±0.2 • C, resolution: 0.01 • C). From September 2018 until the end of 2019, water temperature and level measurements in our monitoring stations on Suceava River were obtained with 2 AquaTROLL500 instruments (accuracy: ±0.1 • C, resolution: 0.01 • C). The temperature measurements of the other streamflows were obtained with DS1922L-F5 iButton instruments (accuracy: ±0.5 • C, resolution: 0.0625 • C). Different stream waters, e.g., Suceava River or its tributaries, were measured with instruments with an accuracy identical for the entire streamflow type.
Air temperature and relative humidity (RH) in the stations used to create the maps in Figure 2 (with the exception of the official weather station of Suceava city) were measured with CEM DT-171 thermohygrometers (accuracy: ±1 • C and ±3.0% to ±5.0% RH, resolution: 0.1 • C and 0.1% RH) placed in miniature weather stations, made of wood, with good air circulation, preventing the direct solar radiation, with waterproof roof and fixed on wooden pillars at 2 m from the ground.

Data
We used hourly and daily time series with different lengths from 2018 and 2019 at UTC+2. The year 2019 was selected for in-depth analyses due to the broader data availability for all parameters. Data used in this study is divided into data obtained from our own instruments and data provided by various Romanian institutions. Data from our own instruments is represented by measurements of water temperature (5 sites) and, secondarily, water level (we used here, for a complementary analysis, only level data measured at Tişăuţi with AT500 instruments); also, from the air environment, we obtained hourly values of temperature and relative humidity in 18 sites, of which 4 (located in Suceava floodplain- Table 1) were selected for further integration in the meteorological-hydrological analysis and comparison. Water temperature at Mihoveni (upstream) and Tişăuţi (downstream) was recorded from 28 June 2018 until the end of 2019. Water level at Tişăuţi used here was recorded for the entire year of 2019. Water temperature of the other streamflows was recorded from 20 September 2019 until the end of 2019. The sensors were placed in areas with persistent water flow. In the case of Suceava River, sensors were permanently under at least a 0.33 m column of water; for the tributaries, the minimum water column had a height of 0.1 m (maximum possible during low flow). Our meteorological data were recorded during the entire year of 2019 and above mostly grassy active surfaces. Water temperature data of Suceava River at Mihoveni (upstream) and Tişăuţi (downstream) and air temperature and relative humidity data of Suceava floodplain are available at http://water.usv.ro/data.php.
Data from Romanian institutions is represented by data provided by ANAR (Romanian Waters National Administration-discharge and temperature data of Suceava River and atmospheric precipitation, both at Iţcani-daily values for the entire year of 2019), ACET Suceava (the company responsible for water supply and sewerage in Suceava metropolitan area-discharge and temperature data of Suceava wastewater treatment plant effluent-daily data of 2019), ANM (National Meteorological Administration-air temperature, relative humidity and precipitation at Suceava Weather Station-hourly and daily data of 2019) and ANPM (National Environment Protection Agency-air temperature, relative humidity and solar radiation at SV2 station-hourly and daily data for 2019).
In 2019, the average yearly flow rate of Suceava River was 19.5 m 3 /s at Iţcani station, which is only 15.6% higher than the multiannual discharge of 16.87 m3/s (6.9 m3/s standard deviation) [27] and, therefore, 2019 can be considered an average hydrological year. But, on the other side, the atmospheric environment above the city during the same year was hotter (+2.09 • C) and drier (−86 mm; −13.8%) than the last 60 years' average, the mark of a continuous climate change and increased urbanization of Suceava city metropolitan area. We suggest that 2019 can be taken into account as a case study year, which is valuable especially because studies using hourly data on thermal pollution of urban rivers are missing in Romania.

Analysis Methods
The maps that used interpolation of stations/points with measurements of air temperature and relative humidity were created in ArcGIS with a hybrid method that combines the multiple regression Water 2020, 12, 1343 6 of 24 method with the surface interpolation of stations/points using an inverse distance weighted (IDW) method [28]. The spatial tendency of a parameter is calculated by applying the multiple regression. The difference between the real values and those estimated by regression (residues) are then interpolated using the IDW. Finally, the raster with the interpolated residues is added to the spatial trend in order to correct the spatial representation of the analysed parameter.
In this study, boxplots, seasonal adjustment, smoothing, normalising, surface plots and scalograms are done in MATLAB. The surfaces plot was created in MATLAB by using the "surf" function. The surfaces plot in this study used normalised matrices composed of vectors extracted from time series that had the seasonal oscillations removed by subtracting from the raw time series the smoothed variant of the selected time series. The subtraction acted as a high-pass filter. The smoothing used a moving average filter with a span of 49 values (in order to include 2 consecutive days). The vectors/columns of the matrices represent average hourly values and the normalisation process computed the z-score values of each vector (centre 0, standard deviation 1). The normalisation was necessary for representing on the same scale (for comparison purposes) the diurnal thermal amplitudes from both semesters of 2019 and from both environments, even if the real amplitude is higher in the warm semester and in the air.
Scalograms represent the graphic result of continuous wavelet analysis or wavelet coherence analysis. The wavelet transform involved the Morlet mother wavelet, and the continuous wavelet transform analysis used a rectified power spectrum [29]. We have chosen Morlet as the mother wavelet because it is frequently used in the analysis of time series in geosciences due to the good time-frequency localisation of events. The wavelet coherence analysis used the methodology of Grinsted et al. [30], which continued to implement the Morlet mother wavelet in the study of hydrological time series. For a detailed description of the wavelet analysis and equations, see our complementary studies [31,32].
XLSTAT (Statistical Software for Excel) was used for ANCOVA (analysis of covariance), which involved one quantitative dependent variable (parameter) and multiple quantitative explanatory variables (other parameters) for calculating the goodness-of-fit coefficients of the linear regression model, obtaining the Fisher's F test result (estimates the risk in assuming that the null hypothesis is wrong/that there is no effect of the explanatory variables) and producing the Type III SS Table (Sum of Squares analysis). This table estimates the contribution of each explanatory variable to the model by removing one variable at a time and evaluating its effect on the quality of the computed model (this removes the undesirable effect of the order in which the variables are selected in the process of model calculation). ANCOVA assumptions are: the dependent variable and covariate variables are measured on a continuous scale, there are no significant outliers, residuals are approximately normally distributed for the independent variable, there is a homogeneity of variances and there is a homogeneity of regression slopes.
An energy balance of the water system in this study was calculated according to the following formula: temperature SVT = (mass SVM * temperature SVM + mass TR * temperature TR + mass WW * temperature WW )/ (mass SVM + mass TR + mass WW ) where mass is the mass of the water flows calculated from discharge, SVT is Suceava River downstream (Tişăuţi), SVM is Suceava River upstream (Mihoveni), TR is the sum of Suceava River tributaries between the upstream and downstream monitoring stations on Suceava River and WW is the wastewater from the wastewater treatment plant of Suceava city.

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.
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 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 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. 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].  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 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. Water 2020, 4, x FOR PEER REVIEW 13 of 24 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. 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. Water 2020, 4, x FOR PEER REVIEW 14 of 24 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. 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).  During September 20-December 31, 2019, the average water temperature of (in upstreamdownstream 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 m 3 /s of wastewaters with an average of 16.31 °C. During September 20-December 31, 2019, it discharged 0.32 m 3 /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 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 m 3 /s of wastewaters with an average of 16.31 • C. During September 20-December 31, 2019, it discharged 0.32 m 3 /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 m 3 /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 m 3 /s from 0.25 m 3 /s streamflow + 0.32 m 3 /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 m 3 /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]. 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 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. 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. Figure 12. Multiple wavelet coherence analysis scalograms between Suceava river water temperature and air parameters: (a) water temperature upstream versus air temperature at Suceava Weather Station and water level upstream, (b) water temperature downstream versus air temperature in the floodplain and water temperature upstream. Figure 12. Multiple wavelet coherence analysis scalograms between Suceava river water temperature and air parameters: (a) water temperature upstream versus air temperature at Suceava Weather Station and water level upstream, (b) water temperature downstream versus air temperature in the floodplain and water temperature upstream.
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].

Conclusions
In 2019, the average air temperature was 10.35 • C at Suceava Weather Station and 11.05 • C in Suceava River floodplain/centre of the city. Suceava River water temperature was 11.54 • C upstream of Suceava city (Mihoveni) and 11.97 • C downstream (Tişăuţi). The analysis of the average diurnal profiles of water and air in our study area indicates that the heat transfer between water and air is bidirectional. The peak position in the diurnal cycle of water has been strongly shifted towards the moment of the maximum air temperature of city. Monthly average values show that, in the first half of the year, air temperature in the urban floodplain of Suceava River was higher than the stream water temperature in both stations.
During a detailed spatial analysis (last quarter of 2019), the water temperature of some urban tributaries was~2 • C higher than the water temperature of Suceava River. These tributaries are strongly impacted by the urban heat island and, together with the effluent of the wastewater treatment plant, increase the temperature of the main river.
The wavelet coherence analysis shows that there is a better covariance of the diurnal periodicities in the temperatures of urban floodplain air and Suceava River water downstream, at Tişăuţi, than between the air temperature nearby the city and Suceava River upstream, at Mihoveni. Simple and multiple wavelet coherence analyses and ANCOVA indicate that the evolution of Suceava River water temperature downstream of Suceava city is controlled mainly by urban air temperature, water temperature upstream, wastewater temperature and solar radiation.
Last, but not least, 2019 was a dry and warm year (atmospheric environment) and Suceava River had warm waters. Further analyses using high temporal resolution data are needed during average or wetter/colder years in order to better understand the relationship between the urban heat island and Suceava River. For a better understanding of the distinct influence of the urban heat island and wastewaters on the changes of Suceava River, a monitoring station on Suceava River immediately upstream of the wastewater treatment plant is necessary for future studies. Data provided from this new station will help in assessing the gradual changes in the diurnal thermal profile of Suceava River from the uppermost monitoring station towards downstream of Suceava city. Data collected in this and future studies could be integrated into computer models that might help local authorities mitigate the impact of the urban heat island on urban stream waters. Best management practices should be promoted for Suceava floodplain. Restoration treatments may lower a river temperature and increase its dissolved oxygen concentration, as proven in other urban areas [6]. Because riparian and floodplain forests keep the river temperature low [45], we recommend a forestation with willows (Salix sp.) of Suceava River floodplain inside Suceava city in order to mitigate the urban stream syndrome.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4441/12/5/1343/s1, Figure S1: Urban heat island above Suceava city as revealed by hourly air relative humidity measurements in 2019 during winter -a -and summer -b -(values used in maps are averages per season) -UHI represents the area most impacted by the urban heat island. Figure S2: Comparative evolution of Suceava River water temperature upstream (Mihoveni) and downstream (Tişăuţi) of Suceava city (hourly data, June 28th, 2018 -December 31st, 2019). Figure S3: Comparisons of water temperature at Mihoveni (upstream) and Tişăuţi (downstream) and air temperature in Suceava city area in 2019: a. average diurnal profiles of water temperature, b. average diurnal profiles of air temperature. Figure S4: Relative humidity diurnal regime of the air at Suceava Weather Station and in the floodplain in 2019 during summer (a) and winter (b) (∆RH is the difference between the floodplain values and those of Suceava Weather Station). Table S1: Correlation matrix of selected water and air parameters in Suceava city in 2019.