Next Article in Journal
The Impact of Reforestation Induced Land Cover Change (1990–2017) on Flood Peak Discharge Using HEC-HMS Hydrological Model and Satellite Observations: A Study in Two Mountain Basins, China
Next Article in Special Issue
The Assessment of Hydrological Availability and the Payment for Ecosystem Services: A Pilot Study in a Brazilian Headwater Catchment
Previous Article in Journal
Effect of Plastic Film Residue on Vertical Infiltration Under Different Initial Soil Moisture Contents and Dry Bulk Densities
Previous Article in Special Issue
Tracking Lake and Reservoir Changes in the Nenjiang Watershed, Northeast China: Patterns, Trends, and Drivers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Geography, Ștefan cel Mare University of Suceava, 720229 Suceava, Romania
2
Computers, Electronics and Automation Department, Ștefan cel Mare University of Suceava, 720229 Suceava, Romania
3
Suceava Weather Station, Regional Meteorological Centre of Moldova, National Meteorological Administration of Romania, 720237 Suceava, Romania
*
Author to whom correspondence should be addressed.
Water 2020, 12(5), 1343; https://doi.org/10.3390/w12051343
Submission received: 26 March 2020 / Revised: 3 May 2020 / Accepted: 6 May 2020 / Published: 9 May 2020

Abstract

:
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şăuţi)—a 3.7% increase in the water temperature downstream. After the stream water passes through the city, the diurnal thermal profile 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 difference between the urban floodplain and the surrounding area. The stream water diurnal cycle has been shifted towards a variation specific 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 significant 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 significant (0.95 confidence level) causal relationships between the changes in Suceava River water temperature downstream and the fluctuations of the urban air temperature. The complex bidirectional heat transfer and the changes in the diurnal thermal profiles are important to be analysed in other urban systems in order to decipher in more detail the observed causal relationships.

1. 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 urban surfaces generate temperature surges in streams after localized rainstorms [3,7]. The impervious urban areas reduce rainwater infiltration in the catchment and, as a consequence, the high waters are more intense and the baseflow is weaker [8,9]. The shallower waters during baseflow are more prone to being heated by the urban air [3]. The higher temperature of urban stream waters reduce their self-purification capacity by affecting the aquatic biota and the amount of dissolved oxygen in the water [10]. Changes in fish populations’ structure and diversity are easily identifiable in rivers switched to an urban regime, and species that better tolerate warmer waters replace the previous ones [11,12].
Rivers passing through cities may have mean temperatures above or below the corresponding mean air temperature [13]. In both cases, discussions may arise concerning not only the influence of the city on the river, but also the inverse case. Colder rivers can have a cooling effect on a city, which may be enhanced by air circulation along the stream corridor [14,15,16]. At the opposite side, the warmer rivers have an effect of enhancing the urban heat island [17].
Due to the important role of water temperature in aquatic ecosystems, river water temperature is included in water quality indices [18]. Also, some software solutions were proposed to model the thermal impact of cities on rivers in order to find pollution causes and improve city management [19,20]. Overall, the increase in the urban stream water temperature is recognised as an essential aspect of the urban stream syndrome [21].
As an effect of climate change, river temperatures in north-eastern Romania are expected to increase by 1–2 °C at the end of this century, according to some models [13]. The predicted increases in air temperature in this part of Europe [22] will be higher in urban areas. The temporal evolution of stream water temperature of rivers in Romania was analysed only in a few studies [23,24,25,26,27]. Fewer are the studies about river water temperature in cities [26,27]. In this study, we aim to describe the thermal impact of Suceava city on its stream waters, especially Suceava River, by using high temporal resolution data and multiple analysis methods. The syntheses obtained from our data are relevant for a usual city–river relationship in Romania. Specific objectives of this study are (1) to find if and how the urban heat island alters the water temperature of Suceava River and (2) to identify the role played by urban tributaries and wastewaters in the observed variations of Suceava River temperature.

2. Materials and Methods

2.1. 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 m3/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 m3/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 m3/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 m3/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).

2.2. 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.

2.3. 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 m3/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.

2.4. 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 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:
temperatureSVT = (massSVMtemperatureSVM + massTRtemperatureTR + massWWtemperatureWW)/
(massSVM + massTR + massWW)
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.

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 (R2) 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].

4. 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 https://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.

Author Contributions

Conceptualization, A.-E.B.; Data curation, A.-E.B. and D.M.; Investigation, A.-E.B. and D.M.; Methodology, A.-E.B. and D.M.; Resources, A.-E.B., A.P., D.I.O. and P.I.B.; Supervision, A.G.; Visualization, D.I.O., A.P. and P.I.B.; Writing—original draft, A.-E.B. and D.M.; Writing—review and editing, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CNCS–UEFISCDI, project number PN-III-P1-1.1PD-2016-2106.

Acknowledgments

Due to important field and laboratory work, Dumitru Mihăilă is to be considered, together with Andrei-Emil Briciu, as one of the principal authors. Wavelet coherence software was provided by A. Grinsted. Thanks for the data provided is addressed to three Romanian institutions, ANAR (Siret Water Basin Administration), ANM, ANPM, and to ACET Suceava. This study was supported by data obtained within the research project SQRTDA (Streamwater Quality Real-Time Data Analysis). This work was supported by a grant of Ministry of Research and Innovation, CNCS–UEFISCDI, project number PN-III-P1-1.1PD-2016-2106, within PNCDI III.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hannah, D.M.; Garner, G. River water temperature in the United Kingdom: Changes over the 20th century and possible changes over the 21st century. Prog. Phys. Geogr. 2015, 39, 68–92. [Google Scholar] [CrossRef] [Green Version]
  2. Wagner, T.; Midway, S.R.; Whittier, J.B.; DeWeber, J.T.; Paukert, C.P. Annual Changes in Seasonal River Water Temperatures in the Eastern and Western United States. Water 2017, 9, 90. [Google Scholar] [CrossRef] [Green Version]
  3. Nelson, K.C.; Palmer, M.A. Stream temperature surges under urbanization and climate change: Data, models, and responses. J. Am. Water Resour. Assoc. 2007, 43, 2. [Google Scholar] [CrossRef]
  4. Van Sickle, J.; Baker, J.; Herlihy, A.; Bayley, P.; Gregory, S.; Haggerty, P.; Ashkenas, L.; Li, J. Projecting the Biological Condition of Streams under Alternative Scenarios of Human Land Use. Ecol. Appl. 2004, 14, 368–380. [Google Scholar] [CrossRef] [Green Version]
  5. Paul, M.J.; Meyer, J.L. Streams in the Urban Landscape. Annu. Rev. Ecol. Syst. 2001, 32, 333–365. [Google Scholar] [CrossRef]
  6. Abdi, R.; Endreny, R.; Nowak, D.A. A model to integrate analysis of urban river thermal cooling in river restoration. J. Environ. Manag. 2020, 258, 110023. [Google Scholar] [CrossRef]
  7. Galli, J.; Dubose, R. Water Temperature and Freshwater Stream Biota: An Overview; U.S. Environmental Protection Agency: Fort Meade, MD, USA, 1990. [Google Scholar]
  8. Pluhowski, E.J. Urbanization and Its Effect on the Temperature of the Streams on Long Island; Professional Paper 627-D; U.S. Government Printing Office: Washington, DC, USA, 1970. [Google Scholar]
  9. Brabec, E.; Schulte, S.; Richards, P.L. Impervious Surfaces and Water Quality: A Review of Current Literature and Its Implications for Watershed Planning. J. Plan. Lit. 2002, 16, 499–514. [Google Scholar] [CrossRef]
  10. Zaltman, R.; Reid, G.W. Theoretical study on thermal pollution. Geofís. Int. 1972, 12, 1. [Google Scholar]
  11. Daufresne, M.; Boët, P. Climate change impacts on structure and diversity of fish communities in rivers. Glob. Chang. Biol. 2007, 13, 2467–2478. [Google Scholar] [CrossRef]
  12. Demars, B.O.; Russell Manson, J.; Ólafsson, J.S.; Gíslason, G.M.; Gudmundsdóttir, R.; Woodward, G.; Reiss, J.; Pichler, D.E.; Rasmussen, J.J.; Friberg, N. Temperature and the metabolic balance of streams. Freshw. Biol. 2011, 56, 1106–1121. [Google Scholar] [CrossRef]
  13. Van Vliet, M.T.H.; Franssen, W.H.P.; Yearsley, J.R.; Ludwi, F.; Haddeland, I. Global river discharge and water temperature under climate change. Glob. Environ. Chang. 2013, 23, 450–464. [Google Scholar] [CrossRef]
  14. Sugawara, H.; Narita, K.; Kim, M.S. Cooling effect by urban river. In Proceedings of the seventh International Conference on Urban Climate, Yokohama, Japan, 29 June–3 July 2009. [Google Scholar]
  15. Qi, J.; Liu, J.; Song, X.; Guo, L. Field measurement of the influence of large urban river on urban thermal climate. J. Harbin Inst. Tech. 2011, 43, 56–59. [Google Scholar]
  16. Cheng, L.; Guan, D.; Zhou, L.; Zhao, Z.; Zhou, J. Urban cooling island effect of main river on a landscape scale in Chongqing, China. Sustain. Cities Soc. 2019, 47, 101501. [Google Scholar] [CrossRef]
  17. Moyer, A.N.; Hawkins, T.W. River effects on the heat island of a small urban area. Urban Clim. 2017, 21, 262–277. [Google Scholar] [CrossRef]
  18. Kannel, P.R.; Lee, S.; Lee, Y.S.; Kanel, S.R.; Khan, S.P. Application of water quality indices and dissolved oxygen as indicators for river water classification and urban impact assessment. Environ. Monit. Assess. 2007, 132, 93–110. [Google Scholar] [CrossRef]
  19. Abdi, R.; Endreny, T. A River Temperature Model to Assist Managers in Identifying Thermal Pollution Causes and Solutions. Water 2019, 11, 1060. [Google Scholar] [CrossRef] [Green Version]
  20. Abdi, R.; Endreny, T.; Nowakc, D. i-Tree Cool River: An open source, freeware tool to simulate river water temperature coupled with HEC-RAS. MethodsX 2020, 7, 100808. [Google Scholar] [CrossRef]
  21. Somers, K.A.; Bernhardt, E.S.; Grace, J.B.; Hassett, B.A.; Sudduth, E.B.; Wang, S.; Urban, D.L. Streams in the urban heat island: Spatial and temporal variability in temperature. Freshw. Sci. 2013, 32, 309–326. [Google Scholar] [CrossRef] [Green Version]
  22. Mihăilă, D.; Briciu, A.-E. Actual climate evolution in the NE Romania. Manifestations and consequences. In Proceedings of the SGEM2012 Conference Proceedings, Albena, Bulgaria, 1 Augest 2012; Volume 4, pp. 241–252. [Google Scholar]
  23. Ajeagah, G.; Cioroi, M.; Praisler, M.; Constantin, O.; Palela, M.; Bahrim, G. Bacteriological and environmental characterisation of the water quality in the Danube River Basin in the Galati area of Romania. Afr. J. Microbiol. Res. 2012, 6, 292–301. [Google Scholar]
  24. Cojoc, G.M.; Romanescu, G.; Tirnovan, A. The importance of water temperature fluctuations in relation to the hydrological factor. case study—Bistrita river basin (Romania). PESD 2014, 8, 2. [Google Scholar] [CrossRef] [Green Version]
  25. Briciu, A.-E.; Oprea-Gancevici, D.I. Diurnal thermal profiles of selected rivers in Romania. In Proceedings of the SGEM2015 Conference Proceedings, Albena, Bulgaria, 18 June 2015; Volume 1, pp. 221–228. [Google Scholar]
  26. Briciu, A.-E.; Mihaila, D.; Mihaila, D. Short, medium and long term stochastic analysis of the Suceava River pollution evolution in the homonymous city. In Proceedings of the SGEM 2012 Conference Proceedings, Albena, Bulgaria, 1 Augest 2012; Volume 3, pp. 809–816. [Google Scholar]
  27. Briciu, A.-E. Studiu de hidrologie urbană în arealul municipiului Suceava; (Urban hydrology study in Suceava municipality area); Ştefan cel Mare University Publishing House: Suceava, Romania, 2017. [Google Scholar]
  28. Lu, G.; Wong, D. An adaptive inverse-distance weighting spatial interpolation technique. Comput. Geosci. 2008, 34, 1044–1055. [Google Scholar] [CrossRef]
  29. Ng, E.K.W.; Chan, J.C.L. Geophysical applications of partial wavelet coherence and multiple wavelet coherence. J. Atmos. Oceanic Technol. 2012, 29, 1845–1853. [Google Scholar] [CrossRef]
  30. Grinsted, A.; Moore, J.C.; Jevrejeva, S. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Process. Geophys. 2004, 11, 5–6. [Google Scholar] [CrossRef]
  31. Briciu, A.-E.; Graur, A.; Oprea, D.I.; Filote, C. A Methodology for the Fast Comparison of Streamwater Diurnal Cycles at Two Monitoring Stations. Water 2019, 11, 2524. [Google Scholar] [CrossRef] [Green Version]
  32. Briciu, A.-E. Changes in Physical Properties of Inland Streamwaters Induced by Earth and Atmospheric Tides. Water 2019, 11, 2533. [Google Scholar] [CrossRef] [Green Version]
  33. Gribovszki, Z.; Szilágyi, J.; Kalicz, P. Diurnal fluctuations in shallow groundwater levels and streamflow rates and their interpretation—A review. J. Hydrol. 2010, 385, 371–383. [Google Scholar] [CrossRef] [Green Version]
  34. Rice, J.S.; Anderson, W.P., Jr.; Thaxton, C.S. Urbanization influences on stream temperature behavior within low-discharge, headwaters streams. Hydrol. Res. Lett. 2011, 5, 27–31. [Google Scholar] [CrossRef] [Green Version]
  35. Zeiger, S.J.; Hubbart, J.A. Urban Stormwater Temperature Surges: A Central US Watershed Study. Hydrology 2015, 2, 193–209. [Google Scholar] [CrossRef]
  36. Orr, H.G.; Simpson, G.L.; des Clers, S.; Watts, G.; Hughes, M.; Hannaford, J.; Dunbar, M.J.; Laizé, C.L.R.; Wilby, R.L.; Battarbee, R.W.; et al. Detecting changing river temperatures in England and Wales. Hydrol. Process. 2015, 29, 752–766. [Google Scholar] [CrossRef] [Green Version]
  37. Jonkers, A.R.; Sharkey, K.J. The Differential Warming Response of Britain’s Rivers (1982–2011). PLoS ONE 2016, 11, e0166247. [Google Scholar] [CrossRef] [Green Version]
  38. Cox, R.A.; Smith, N.D. The specific heat of sea water. Proc. Roy. Soc. Lond. 1959, 252, 51–62. [Google Scholar]
  39. Hathway, E.A.; Sharples, S. The interaction of rivers and urban fromin mitigating the urban heat island effect: A UK case study. Build. Environ. 2012, 58, 14–22. [Google Scholar] [CrossRef] [Green Version]
  40. Lolli, S.; D’Adderio, L.P.; Campbell, J.R.; Sicard, M.; Welton, E.J.; Binci, A.; Rea, A.; Tokay, A.; Comerón, A.; Barragan, R.; et al. Vertically Resolved Precipitation Intensity Retrieved through a Synergy between the Ground-Based NASA MPLNET Lidar Network Measurements, Surface Disdrometer Datasets and an Analytical Model Solution. Remote Sens. 2018, 10, 1102. [Google Scholar] [CrossRef] [Green Version]
  41. Kinouchi, T.; Yagi, H.; Miyamoto, M. Increase in stream temperature related to anthropogenic heat input from urban wastewater. J. Hydrol 2007, 335, 78–88. [Google Scholar] [CrossRef]
  42. Kinouchi, T. Impact of long-term water and energy consumption in Tokyo on wastewater effluent: Implications for the thermal degradation of urban streams. Hydrol. Process. 2007, 21, 1207–1216. [Google Scholar] [CrossRef]
  43. Briciu, A.-E.; Toader, E.; Romanescu, G.; Sandu, I. Urban Streamwater Contamination and Self-purification in a Central-Eastern European City. Part I. Rev. Chim. 2016, 67, 1294–1300. [Google Scholar]
  44. Briciu, A.-E.; Toader, E.; Romanescu, G.; Sandu, I. Urban Streamwater Contamination and Self-purification in a Central-Eastern European City. Part, B. Rev. Chim. 2016, 67, 1583–1586. [Google Scholar]
  45. Sun, N.; Yearsley, J.; Voisin, N.; Lettenmaier, D.P. A spatially distributed model for the assessment of land use impacts on stream temperature in small urban watersheds. Hydrol. Process. 2015, 29, 2331–2345. [Google Scholar] [CrossRef]
Figure 1. Map of the study area (centred on Suceava city) and location of the monitoring stations.
Figure 1. Map of the study area (centred on Suceava city) and location of the monitoring stations.
Water 12 01343 g001
Figure 2. 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.
Figure 2. 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.
Water 12 01343 g002
Figure 3. Comparison of water temperature at Mihoveni (upstream) and Tişăuţi (downstream) and air temperature in Suceava city area in 2019 (hourly data).
Figure 3. Comparison of water temperature at Mihoveni (upstream) and Tişăuţi (downstream) and air temperature in Suceava city area in 2019 (hourly data).
Water 12 01343 g003
Figure 4. Hourly differences in the diurnal regime of temperature in 2019 during summer (left column) and winter (right column) for: Suceava River downstream—Suceava River upstream (a,b), floodplain air—Suceava Weather Station (c,d), floodplain air—Suceava River downstream (e,f).
Figure 4. Hourly differences in the diurnal regime of temperature in 2019 during summer (left column) and winter (right column) for: Suceava River downstream—Suceava River upstream (a,b), floodplain air—Suceava Weather Station (c,d), floodplain air—Suceava River downstream (e,f).
Water 12 01343 g004
Figure 5. Comparisons of water and air temperatures in the Suceava city area in 2019: (a) annual regime of temperature, (b) average diurnal profiles.
Figure 5. Comparisons of water and air temperatures in the Suceava city area in 2019: (a) annual regime of temperature, (b) average diurnal profiles.
Water 12 01343 g005
Figure 6. Boxplots of water and air temperature time series (hourly data) for the following time intervals (the red line represents the median): (a) year 2019, (b) January 2019, (c) August 2019 (A1—Suceava Weather Station, A2—Suceava city in the floodplain, W1—Suceava River upstream, W2—Suceava River downstream).
Figure 6. Boxplots of water and air temperature time series (hourly data) for the following time intervals (the red line represents the median): (a) year 2019, (b) January 2019, (c) August 2019 (A1—Suceava Weather Station, A2—Suceava city in the floodplain, W1—Suceava River upstream, W2—Suceava River downstream).
Water 12 01343 g006
Figure 7. Surface plots of seasonally adjusted and normalised average hourly values of air and water temperature in Suceava city area in 2019: (a) air at Suceava Weather Station, (b) air in the floodplain, (c) Suceava River upstream and (d) Suceava River downstream.
Figure 7. Surface plots of seasonally adjusted and normalised average hourly values of air and water temperature in Suceava city area in 2019: (a) air at Suceava Weather Station, (b) air in the floodplain, (c) Suceava River upstream and (d) Suceava River downstream.
Water 12 01343 g007
Figure 8. Continuous wavelet transform analysis of hourly water temperature time series of Suceava River in 2019: (a) Mihoveni, (b) Tişăuţi.
Figure 8. Continuous wavelet transform analysis of hourly water temperature time series of Suceava River in 2019: (a) Mihoveni, (b) Tişăuţi.
Water 12 01343 g008
Figure 9. Wavelet coherence analysis of hourly water and air temperature time series of Suceava city in 2019: (a) Suceava River water temperature upstream and downstream, (b) air temperature at Suceava Weather Station and Suceava River water temperature upstream, (c) air temperature in the floodplain and Suceava River water temperature downstream.
Figure 9. Wavelet coherence analysis of hourly water and air temperature time series of Suceava city in 2019: (a) Suceava River water temperature upstream and downstream, (b) air temperature at Suceava Weather Station and Suceava River water temperature upstream, (c) air temperature in the floodplain and Suceava River water temperature downstream.
Water 12 01343 g009
Figure 10. Comparisons of water temperature of Suceava River and its urban tributaries during September 20–December 31, 2019: (a) temporal evolution of hourly data, (b) average diurnal profiles of water temperature.
Figure 10. Comparisons of water temperature of Suceava River and its urban tributaries during September 20–December 31, 2019: (a) temporal evolution of hourly data, (b) average diurnal profiles of water temperature.
Water 12 01343 g010
Figure 11. Links between rainfalls, wastewaters and the water temperature of Suceava River: (a) correlation between the rainfall amount and the wastewater treatment plant discharge, (b) average diurnal profile of Suceava River water temperature difference (downstream–upstream), (c) Suceava River water temperature difference (downstream–upstream) versus atmospheric precipitation at Suceava Weather Station during 3–4 August 2019 (48 h), (d) same as in (c), but during 14–15 July 2019 (24 h).
Figure 11. Links between rainfalls, wastewaters and the water temperature of Suceava River: (a) correlation between the rainfall amount and the wastewater treatment plant discharge, (b) average diurnal profile of Suceava River water temperature difference (downstream–upstream), (c) Suceava River water temperature difference (downstream–upstream) versus atmospheric precipitation at Suceava Weather Station during 3–4 August 2019 (48 h), (d) same as in (c), but during 14–15 July 2019 (24 h).
Water 12 01343 g011
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.
Water 12 01343 g012
Figure 13. Comparative temporal evolutions of various scaled correlations’ time series (smoothed with a span of 1 week/7 values) and water/air parameters of Suceava city area in 2019: (a) scaled correlations of upstream versus downstream time series, (b) downstream scaled correlations versus air temperature in the floodplain, (c) downstream scaled correlation versus water level at Tişăuţi.
Figure 13. Comparative temporal evolutions of various scaled correlations’ time series (smoothed with a span of 1 week/7 values) and water/air parameters of Suceava city area in 2019: (a) scaled correlations of upstream versus downstream time series, (b) downstream scaled correlations versus air temperature in the floodplain, (c) downstream scaled correlation versus water level at Tişăuţi.
Water 12 01343 g013
Table 1. Details of the monitoring stations used in this study.
Table 1. Details of the monitoring stations used in this study.
Site NameCoordinates (WGS84)Elevation
(m above Sea Level)
Riverbed-Banks/
Active Surface Type
Surrounding Land CoverData Type
Water monitoring stations
Suceava at Mihoveni47.681356°N,
26.199646°E
280concretecomplexown
Şcheia47.670982°N,
26.241660°E
276partly concretebuilt areaown
Cetăţii Creek47.647981°N,
26.274243°E
270partly concretebuilt areaown
Collector47.651907°N,
26.282219°E
269concretebuilt areaown
Suceava at Tişăuţi47.617615°N,
26.322629°E
260partly naturalpastureown
Air monitoring stations
Suceava Weather Station47.632349°N,
26.240747°E
359grasscultivated landpublic
EPA industrial type SV247.669097°N,
26.280840°E
289concretebuilt areapublic
Iţcani railway station area (TRN)47.678598°N,
26.229750°E
280partly grassbuilt areaown
AMBRO industrial area (AMB)47.661483°N,
26.270342°E
279grassbuilt areaown
Complex urban area (ITC)47.667320°N,
26.247910°E
276partly grasscomplexown
Water + air monitoring stations
Suceava at Iţcani47.671850°N,
26.246900°E
273natural/grassbuilt area/pasturepublic
Wastewater treatment plant (ACT)47.653000°N,
26.284566°E
271concrete/grasscomplexpublic/own
Table 2. Analysis of covariance (ANCOVA) of Suceava River water temperature upstream or downstream (dependent variable) versus other parameters (explanatory variables) - type III sum of squares analysis.
Table 2. Analysis of covariance (ANCOVA) of Suceava River water temperature upstream or downstream (dependent variable) versus other parameters (explanatory variables) - type III sum of squares analysis.
UPSTREAMDFSum of SquaresMean SquaresFPr > F
temperature WS110,074.59310,074.5931885.322<0.0001
precipitation WS146.71346.7138.7420.003
relative humidity WS1152.214152.21428.485<0.0001
solar radiation fl.175.03775.03714.0420.000
water level M1109.142109.14220.424<0.0001
DOWNSTREAMDFSum of SquaresMean SquaresFPr > F
temperature WS10.0000.0000.0020.964
precipitation WS10.0560.0560.4040.526
relative humidity WS10.0830.0830.5980.440
temperature fl.10.0090.0090.0620.803
solar radiation fl.10.9490.9496.8170.009
relative humidity fl.10.2260.2261.6220.204
precipitation I10.0030.0030.0230.879
water temperature M1611.040611.0404391.208<0.0001
water level M10.2250.2251.6180.204
water level T10.7570.7575.4410.020
discharge WTP10.0000.0000.0020.963
temperature WTP15.8955.89542.363<0.0001
Abbreviations: WS—Suceava Weather Station, M—Mihoveni, I—Iţcani, T—Tişăuţi, fl.—floodplain, WTP—wastewater treatment plant.

Share and Cite

MDPI and ACS Style

Briciu, A.-E.; Mihăilă, D.; Graur, A.; Oprea, D.I.; Prisăcariu, A.; Bistricean, P.I. Changes in the Water Temperature of Rivers Impacted by the Urban Heat Island: Case Study of Suceava City. Water 2020, 12, 1343. https://doi.org/10.3390/w12051343

AMA Style

Briciu A-E, Mihăilă D, Graur A, Oprea DI, Prisăcariu A, Bistricean PI. Changes in the Water Temperature of Rivers Impacted by the Urban Heat Island: Case Study of Suceava City. Water. 2020; 12(5):1343. https://doi.org/10.3390/w12051343

Chicago/Turabian Style

Briciu, Andrei-Emil, Dumitru Mihăilă, Adrian Graur, Dinu Iulian Oprea, Alin Prisăcariu, and Petruţ Ionel Bistricean. 2020. "Changes in the Water Temperature of Rivers Impacted by the Urban Heat Island: Case Study of Suceava City" Water 12, no. 5: 1343. https://doi.org/10.3390/w12051343

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

Article Metrics

Back to TopTop