Societal Implications of Forest and Water Body Area Evolution in Czechia and Selected Regions

: Land cover evolution is an environmental factor that can be used to characterize forest ecosystem services (FES). This study aims to analyze the change in forest cover and water bodies between 1990 and 2018 in the whole Czech Republic, and in the Central Bohemian and South Moravian regions, and its effects on freshwater provision. Additionally, we attempt to understand the societal implications of water quality, public perception, and environmental investment on natural ecosystems. Forest cover and water body data were obtained from the Corine land cover database, while water quality and investment were compiled from the Czech Statistical Ofﬁce. Public perceptions on the Czech FES were collected from a national survey. Between 1990 and 2018, forest cover has increased by 3.94% and water bodies by 7.65%; however, from 2014 to 2018, severe droughts were reported that compromised the availability of surface water, presumably on artiﬁcial structures, causing an increase in the occupied area. Regarding public perception, respondents with less education, and the older population, obtained an assessment of the low performance of the FES, while the water quality and investment indicate that environmental funding has contributed to improving the quality of outﬂow water from the wastewater treatment plants, fulﬁlling all the allowed limits of the urban wastewater treatment directive. Hence, a multidisciplinary approach can help decision makers promote policies that integrate environmental management measures, investment protection, and contribute to sustainable development.


Introduction
Ecosystem services are the benefits that humans obtain from functional ecosystems [1,2], with water provision being among the essential roles of forest ecosystem services (FES) [3,4], while trees are expected to protect and provide clean drinking water [5][6][7]. For human wellbeing, the role of water stretches from consumption to essential recreation elements [8,9]. Since the early 1990s, ecosystem services have been extensively studied to understand interactions and pressures that affect them [6,10], while evaluating the FES has become essential information to policymakers [11]. Ecosystems related to water have been studied by analyzing inland and marine waters, and are classified as hydrologic ecosystem services [12]. Some studies have reported that the results of the perception and evaluation of water, for instance, the optics and perception of water quality towards pollution and aesthetics [13,14], create an evaluation and verbal description of water landscape [9], and of the perception of drinking water quality [15][16][17].
Furthermore, land cover is defined as the bio-or physical cover on the surface of the Earth, which is categorized into different land classes, e.g., grassland, forest area, etc. [18][19][20]. The water body is one of these categories. It can be of temporary or permanent sphere with an area of 78,870 km 2 , a population of 10,693,939, and is divided into 14 regions [76], of which the following two are part of the study: Central Bohemian (Středočeský kraj), and South Moravian (Jihomoravský kraj), see Figure 1.
The Central Bohemian region is the largest region of Czechia and is in the central part of the historical region of Bohemia. It has an area of 10,928 km 2 and 1,385,141 inhabitants [76]. The main sources of water pollution are industry, agriculture, and pollution from small municipalities, which lack sewer connections and wastewater treatment plants (WWTPs). In 2018, agricultural land occupied 60.2% of the territory of the region, forest land 27.5%, and surface water 1.9%. Since 2000, the conversion of arable land into permanent pasture has been reported in built-up areas and courtyards [77].
The South Moravian region reports 1,191,989 inhabitants in an area of 7188 km 2 [76]. In 2018, water quality was reported as polluted to very heavily polluted, due to intensive agricultural management and the textile industry. Agricultural land occupies 58.8% of the region, forests 28.1%, and water areas 2.2%. Since 2005, the cultivated area of agricultural land has decreased and the area of permanent pasture has increased for green urban areas, industrial or commercial units, and mining areas [78]. The Central Bohemian region is the largest region of Czechia and is in the central part of the historical region of Bohemia. It has an area of 10,928 km 2 and 1,385,141 inhabitants [76]. The main sources of water pollution are industry, agriculture, and pollution from small municipalities, which lack sewer connections and wastewater treatment plants (WWTPs). In 2018, agricultural land occupied 60.2% of the territory of the region, forest land 27.5%, and surface water 1.9%. Since 2000, the conversion of arable land into permanent pasture has been reported in built-up areas and courtyards [77].
The South Moravian region reports 1,191,989 inhabitants in an area of 7188 km 2 [76]. In 2018, water quality was reported as polluted to very heavily polluted, due to intensive agricultural management and the textile industry. Agricultural land occupies 58.8% of the region, forests 28.1%, and water areas 2.2%. Since 2005, the cultivated area of agricultural land has decreased and the area of permanent pasture has increased for green urban areas, industrial or commercial units, and mining areas [78].

Data Selection, Acquisition and Preprocessing
The selection of factors that affect the hydrological ecosystem services in the Czech Republic was carried out considering the conditions of the study area and the availability of reliable information, to have enough data to allow the variables to be correlated to identify patterns or trends. National reports on the state of the water resource were searched, finding interesting topics to compare their behavior, such as floods, droughts, and water quality; however, when incorporating other comparison factors such as environmental investment and public perception to make the analysis more meaningful, incompatibilities were identified to correlate variables such as the period of events in the case of floods and droughts. In this way, the quality of water and environmental investment in the WWTPs, to compare their operation with the perception of the public and their connection with the evolution of plant cover, were the most relevant.
To answer the objectives of the paper, research data were collected from different sources. Forest cover and water body evolution data were collected from the Corine land cover (CLC) database. For the societal implication indicators, water quality, instruments on water pollution abatement and control, and the public perception and evaluation on the Czech forest water provisioning services were compiled from different sources, as can be seen in Figure 2.

Forest Cover and Water Bodies
The land cover of Czechia (in ha) was extracted from the Corine land cover (CLC) database from the Copernicus Land Monitoring Service [79] and processed with ArcGIS Desktop version 10.8 [80]. The database was developed by interpreting satellite images,  Landsat, SPOT, IRS P6 LISS, RapidEye, and Sentinel-2; each inventory and update was  developed according to the technology and information available in the elaboration period.  The CLC database contains an inventory and four updates from the years 1990, 2000, 2006, 2012 and 2018. It has evolved according to the progress in the instruments and methods of interpretation and image processing, with a scale of 1: 100,000 and a minimum mapping unit (MMU) of 25 hectares, the geometric precision of 1990 is ≤50 m, while for 2000, 2006 and 2012 it is ≤25 m and for 2018 it is ≤10 m with the Sentinel-2 satellite data. It is a free access product that offers to monitor land cover over time to have reliable information about the phenomena that occur on Earth [81]. For further information on the methodology and instructions for use, please consult the following link: https://land.copernicus.eu/ user-corner/technical-library/clc-product-user-manual/view (accessed on 1 June 2021).
The CLC database includes the following five categories: (1) artificial surfaces; (2) agricultural areas; (3) forest and semi-natural areas; (4) wetlands; and (5) water bodies and 44 classes of land cover. The CLC in Czechia consists of 29 classes (Figure 1). Categories of forest and semi-natural areas, and inland waters are further described in Table 1. The forest areas from the CLC analyzed in this paper were under the category of broad-leaved, coniferous, and mixed forests, which were compared to the information collected from the Ministry of Agriculture (MoA) of the Czech Republic [82,83]. Meanwhile, water body data were compiled from the watercourses and water bodies of the CLC, which were compared with the national reports collected from the yearbooks of Czech Statistical Office Section 3.1 land use or land use balance every six years from 2000 to 2018 [84]. This was due to the fact that annual reports on water resource management do not contain information about surface area occupied by water bodies. The earliest national forestry report was in 1998 and 2000 for the water body land balance, so in this paper, the year 1990 was not included in the comparison analysis of the CLC and the national reports, while the area analysis included the periods reported by the CLC.

Water Quality
Raw data of physicochemical water quality parameters, i.e., biological oxygen demand (BOD), chemical oxygen demand (COD), undissolved substances, nitrogen totals (N total), phosphorus totals (P total), and wastewater treatment plants characteristics, in the period of 2009-2019 were extracted from the annual reports of the Czech Statistical Office on water supply systems, sewerage and watercourses [84]. The units of the measured parameters for both inflow and outflow wastewater treatment plants are expressed in tons/year. The raw data indicate the analysis method for biological oxygen demand using a standard incubation test period of 5 days at 20 • C and chemical oxygen demand using the dichromate method.
Removal efficiency (RE) or load reduction among parameters was also compared based on the required minimum threshold level of load reduction suggested by the urban wastewater treatment directive (UWWTD) [85,86].

Instruments on Water Pollution Abatement and Control
Economic instruments (EIs) are designed to achieve environmentally sustainable objectives [87] using economic or market incentives to generate desired behavior between both producers and consumers [88]. Since there are several classifications of EIs, this research used the categories proposed by the Organisation for Economic Co-operation and Development (OECD) [89]. In the Czech Republic, there are a wide series of instruments reported to address water pollution, such as (1) taxes for permitted discharge of wastewater into groundwater, or surface water; (2) subsidies supporting de-sludging of ponds, flood damage mitigation, environmental education, and enlightenment; and (3) voluntary approaches such as eco-labeling system, EMAS-eco-management and audit scheme. This study used the environmentally motivated subsidies in the form of environmental protection expenditure accounts [89] to provide information on the economic resources devoted to the prevention, reduction, and elimination of water pollution [90] in order to compare Remote Sens. 2021, 13, 4019 6 of 27 the investment with the quality of the water in the wastewater treatment plants. In 2019, more than CZK 10.7 billion was invested in support for water management [91]. According to the methodology of environmental accounts, environmental protection investment refers to expenses for the acquisition of assets for environmental protection and expenses related to environmental protection activities; while environmental noninvestment expenditure refers to operational activities including wage costs, rent payments, energy, and other materials, supplies, and payments for services of the companies that manage the activities. The wastewater management budget includes the construction of water treatment plants and sewerage systems to control water quality and other activities and facilities the prevention of pollution generation [92].
Environmental protection investment and environmental non-investment expenditure of Czechia, in CZK thousands current prices, with inflation, average exchange rate in 2019: EUR 1 = CZK 25.36 [92], on wastewater management for the period of 2005-2019, were collected from the annual reports of the Czech Statistical Office on environmental accounts [93].

Public Perception
Forest perception is any opinions or thoughts about the features of FES in water provisioning [5]. The performance in this paper is defined as the respondent's perceived level of the selected forest products or services based on their knowledge and experiences [94].
A computer-assisted web interviewing (CAWI) technique in November 2020 was performed in collaboration with an external market research company, STEM/MARK, a.s. (Prague, Czech Republic). The recruitment was performed by sending the questionnaires to the potential respondents on the company's list based on age, sex, education level, region, and village size (total recruited respondents aged 16 to 65 years (N = 1509) in the whole Czech Republic). All returned questionnaires were included in the analysis (100%). The questionnaire was designed as a closed-ended tool. The four degrees of the respondent's expectations of the Czech forest services related to provision and regulation of water were used, from one (1 = strongly disagree) to four (4 = strongly agree). One category was added for those who could not judge. Five degrees of scores, whereby 1 (one) was not very important and 5 (five) was very important, were applied for the perceived performance level of the Czech forest services related to the provision of water. The education levels of the respondents were categorized as a high school graduate without the secondary school leaving certificate (maturita), or with the certificate and higher education (university level).

Data Analysis
The CLC was processed by clipping the land cover layers with the region layer of the Czech Republic downloaded from the Nomenclature of Territorial Units for Statistics (NUTS) 2021 [95]. Afterwards, the spatial joint and the area calculation of the different coverage areas present in the Czech Republic were used, then the attributes of the tables were exported to Excel and finally, the graphical outputs of land cover location, losses and gains were generated.
Distribution of the non-categorical data was first tested for normality using the Kolmogorov-Smirnov test. The forest land and water bodies in Czechia were presented in mean and standard deviation. Pairwise comparisons in the observed years (1990, 2000, 2006, 2012, and 2018) were analyzed using the Wilcoxon signed-rank nonparametric test (for non-categorical and not normally distributed data). Similarly, the forest and water body area changes were analyzed using the Wilcoxon signed-rank test for the years 1990-2006, 2006-2018, and 1990-2018. Gain, losses, the rate of forest area, and water body change [33,96] were calculated using the following Formula (1): where area (a) was the forested landscape and water bodies (in ha) at the earliest reported year in the paper (1990), while area (b) referred to the final year (2018). Removal efficiency or load reduction [41,97] for each of the water quality parameters measured per year was calculated using the following Formula (2):

RE or Load reduction (%) =
In f low − Out f low In f luent × 100, The comparison between environmental protection investment, environmental noninvestment expenditure, and water pollution removal efficiency was analyzed using Spearman's nonparametric correlation for the period 2005-2019.
For comparing means of the age of respondents (non-categorical and normally distributed data) between regions, the data were analyzed using ANOVA, while a chi-square test was used to analyze the categorical data (gender of the respondents, education level, city size, and frequency of forest visit). Binary logistic regression was applied to define the significant predictors of the public expectations, and the perceived performance level of the Czech FES related to the provision of freshwater (1 = low score of public expectation/perceived performance level). The variables included in the analysis were age, sex (1 = female), education level (1 = maturita and higher), frequency of forest visit (1 = rarely, twice a year, and less), size of the city, and region.
The statistical significance in all analyses was designated as a p-value of less than 0.05. Statistical analysis was performed using IBM SPSS Statistics version 25 (IBM Corp., Armonk, NY, USA).

Forest Cover and Water Bodies Area at National Level
The forested landscape of the country occupies about 33.7% of the country's total area. From the national report in 2018, the forests consisted of 71.5% coniferous and 27.3% broad-leaved trees, respectively, and about 1.2% was forest land without trees. Of the total coniferous trees, Norway spruce (Picea abies) accounted for 50.0%, followed by 16.2% pine, beech (8.6%), oak (7.3%), birch (2.8%), and others (8.7%) [98]. Of the total Czech forested area, 75% was dedicated to the production forests [99,100]. Additionally, more than 50% of the forests in the country are managed by the state, in which the primary direction of the management (more than 90%) is as a forest-based enterprise [98]. From the CLC, high forest cover was found in South and Central Bohemia, followed by the Plzeň region.
More than half of the area of Czech forests are owned by the state (54.09%). As many as 19.27% of the forests are owned by individual forest owners, followed by municipal and urban forests (17.17%), church and other religious organizations (5%), legal entities (3.27%), and cooperative and associations (1.19%). Furthermore, the rest of the land (0.01%) is other types of forest ownership. The majority share of the state ownership is managed by the Forests of the Czech Republic (state forestry company), and there are also other organizations. Therefore, organizations established by the state have a significant impact on the forest management implementation that is related to the total forest land in the country [101]. Given that the state can pursue uniform management goals, while the goals of the same groups of owners may be different, and due to their share in the ownership structure, they have less opportunity to influence the state forestry policy. The Czech hydrological system consists of the following three main watersheds: the basins of the Elbe, the Odra, and the Morava rivers [102]. In 2018, the water bodies reported by the CLC occupied about 0.74% of the total area of the country; a high water body area was found in South Bohemia and Central Bohemia, followed by the South Moravian region. For the same year, the national report in the category of surface water bodies reported an area of 2.11%. Figure 3 depicts the geospatial representation of forest cover and water bodies in Czechia in 1990-2006, 1990-2018, and 2006-2018. The CLC database generated about 1.5% more forest areas (range: 0.72 to 1.90) in the selected years, as compared to the national report by the Czech MoA. The low proportions of additional areas in the CLC were similar to the forested landscape without trees [98,103,104]. The CLC provided information on coniferous, broad-leaved, and mixed-forest regions; notwithstanding, the Czech forestry report only presented the areas of coniferous and broad-leaved forests.   Table 2 shows the Czech forest areas' evolution based on the CLC classification by the type of forest trees. Ordinarily, based on the total forested landscape and the areas by forest type (coniferous, broad-leaved, and mixed forest), significant increases have been observed in 1990, compared to the years 2006, 2012, and 2018. No significant increase in the forest areas was found in the broad-leaved and mixed forests between 1990 and 2018. The coniferous forest area slightly decreased after 2012, and even more in 2018, compared to 2006.
Regarding the evolution of the Czech water body area by the type of inland water, there were no significant changes in the area concerning watercourses and bodies in the different periods. However, the total water bodies showed a significant increase between 1990 and 2018. Similarly, as reported in the national land use balance, surface water increased its area in the period 2000 to 2019. Between 2000, 2006, and 2012, there was no significant increase (Table 2).  Figure 3 indicates, in red, that to the north in the Liberec and Hradec Králové regions, to the northeast in the Moravian-Silesian and Olomouc regions, and the west in the Plzeň region, there was a reduction in the forest cover in the studied areas (coniferous, broad-leaved, and mixed forest), which became transitional woodland-shrubbery, that is, transitional areas with shrubbery or herbaceous vegetation with occasional scattered trees; however, they remain in the category of forest and semi-natural areas. This is consistent with the results of the study conducted in 2018, where, according to the scale resolution (20 m), the land cover analysis identified a reduction in coniferous forests due to disturbances of the bark beetle and the Kyrill windstorm (2007), while deciduous and mixed forests persisted during the study period [105]. At the same time, Figure 3 shows a violet gain area, which occurs in small patches distributed throughout the territory, and, specifically, in the regions of Karlovy Vary, Ústí nad Labem, Liberec, Vysočina, South Bohemia, and Plzeň. Table 2 shows the Czech forest areas' evolution based on the CLC classification by the type of forest trees. Ordinarily, based on the total forested landscape and the areas by forest type (coniferous, broad-leaved, and mixed forest), significant increases have been observed in 1990, compared to the years 2006, 2012, and 2018. No significant increase in the forest areas was found in the broad-leaved and mixed forests between 1990 and 2018. The coniferous forest area slightly decreased after 2012, and even more in 2018, compared to 2006. Regarding the evolution of the Czech water body area by the type of inland water, there were no significant changes in the area concerning watercourses and bodies in the different periods. However, the total water bodies showed a significant increase between 1990 and 2018. Similarly, as reported in the national land use balance, surface water increased its area in the period 2000 to 2019. Between 2000, 2006, and 2012, there was no significant increase (Table 2).
Forest cover and watershed losses and gains in the country were also evaluated, as depicted in Table 3. The results presented a slight decrease in the total Czech forest area (−0.98%) from 2006 to 2018, especially in the coniferous forest-around −3.40%. In contrast, the broad-leaved and mixed forests consistently increased over the years. Furthermore, the water body in 1990-2018 showed a steady incremental gain.

Forest Cover and Water Bodies in the Case Study Areas
In 1990, the forests in the South Moravian region occupied about 18.8% of the total Czech forest area, which decreased to 16.7% in 2018, while the proportion of forested landscape in the Central Bohemian region increased slightly, from 11.1% in 1990 to 12.0% in 2018 ( Table 2). The results of the CLC evolution in the selected years showed that the forest cover in the case study areas agreed with the national data. Forested landscapes in Central Bohemia and South Moravia had the highest increment in 2006, compared to 1990 (+5.90%, +2.42%, respectively). Meanwhile, a decrease in the total forest area was also found in these regions from 2006 to 2018, by −0.20% and −1.72%, respectively. The Central Bohemian region consistently showed the highest increase in total forest area, and the lowest loss of forest areas as well.
The water bodies in Central Bohemia increased by about 6.8% in 1990-2018, while they decreased by about 1.4% in South Moravia in the same period. On the other hand, the watercourses in Central Bohemia reported significant differences based on the Wilcoxon signed-rank test when comparing 1990 with 2000, 2012, and 2018, with a slight increase of around 2.0% between 1990 and 2018. Conversely, South Moravia reported a reduction of 10.6% in the same period (Table 2).   ; an increase in all the water quality parameters is shown. In 2019, the WWTPs in Czechia had fulfilled the threshold level of removal efficiency among the measured parameters suggested by the UWWTD (BOD, COD, undissolved substance, total N, and total P), by more than 75%, 70-90%, 50%, 75%, and 75%, respectively.  ; an increase in all the water quality parameters is shown. In 2019, the WWTPs in Czechia had fulfilled the threshold level of removal efficiency among the measured parameters suggested by the UWWTD (BOD, COD, undissolved substance, total N, and total P), by more than 75%, 70-90%, 50%, 75%, and 75%, respectively.

Public Perception
Of the total recruited respondents, 1338 were included in the analysis because they were able to judge the expectations of the water provisioning services of the Czech forests and their performance level. Furthermore, 170 and 147 respondents from the Central Bohemian and South Moravian regions, respectively, were also selected for the analysis. Table 4 presents the general characteristics of the respondents in each region, and compares the results to the total of the Czech respondents selected. The Central Bohemian respondents were significantly older than the others, most of them possessing the national secondary school leaving certificate (maturita). Most of the respondents also lived in a town with less than 20,000 inhabitants, although the area was the second most populated region after the capital city of Prague [106]. Meanwhile, the highest proportion of respondents with higher education levels (66.9%), living in the town with more than 20,000 inhabitants, were found in the South Moravian region (38.3%).

Public Perception
Of the total recruited respondents, 1338 were included in the analysis because they were able to judge the expectations of the water provisioning services of the Czech forests and their performance level. Furthermore, 170 and 147 respondents from the Central Bohemian and South Moravian regions, respectively, were also selected for the analysis. Table 4 presents the general characteristics of the respondents in each region, and compares the results to the total of the Czech respondents selected. The Central Bohemian respondents were significantly older than the others, most of them possessing the national secondary school leaving certificate (maturita). Most of the respondents also lived in a town with less than 20,000 inhabitants, although the area was the second most populated region after the capital city of Prague [106]. Meanwhile, the highest proportion of respondents with higher education levels (66.9%), living in the town with more than 20,000 inhabitants, were found in the South Moravian region (38.3%). Generally, no statistically significant differences were found between the selected regions compared with the entire respondents. The total Czech respondent's expectation was confirmed by the perceived performance level based on the results of Spearman correlation (r = 0.442), a similar trend was found among the respondents in the South Moravian; in contrast, the Central Bohemian respondents had less expectation on the forests in providing the water source. The mean of forests' expectation was above three and the mean perceived performance level in providing water sources was 4.3 in the three study areas (Figure 8).   (33) 22.0 (294) ta are presented using mean ± sd, or % (n); for comparing means between regions, non-categorical data were analyzed using OVA, while a chi-square test was used to analyze the categorical data; * p-value < 0.05; ** p-value < 0.001.
Generally, no statistically significant differences were found between the selected regions compared with the entire respondents. The total Czech respondentʹs expectation was confirmed by the perceived performance level based on the results of Spearman correlation (r = 0.442), a similar trend was found among the respondents in the South Moravian; in contrast, the Central Bohemian respondents had less expectation on the forests in providing the water source. The mean of forestsʹ expectation was above three and the mean perceived performance level in providing water sources was 4.3 in the three study areas (Figure 8).

Discussion
In order to answer the research questions, the studied factors were grouped to present the interactions and the effects identified throughout the periods analyzed, with respect to the hydrologic ecosystem service.

Discussion
In order to answer the research questions, the studied factors were grouped to present the interactions and the effects identified throughout the periods analyzed, with respect to the hydrologic ecosystem service.

Forest and Water Area Evolution in the Czech Republic
Since 1995, when the Forest Act came into force, Czechia has declared the intention to protect forest areas that provide non-wood functions to the forest, such as the protection of water, soil, and the preservation of biological diversity, among others. Those areas are called special purpose forests and are recognized by the state forest administration, or by the initiative of the forest owner [107]. Similarly, the legal framework of the forest sector in the Czech Republic, following European policies, has allowed the sustainable development of the forest, through various documents to identify forest problems and propose solution measures, as well as monitoring, harvesting, disease control, research, grants, and incentives. The evolution of the National Forest Programs has allowed ministerial cooperation and the establishment of governmental entities, such as the Forests of the Czech Republic State Enterprise, Forest Management Institute, and Forestry and Game Management Research Institute [108,109]. Moreover, the analysis of the relationship between water and forests is aligned with the EU forest strategy, by recognizing the role of forests and their multifunctionality in the provision of socio-economic services [110].
A significant increase in the Czech forest area, since the 1990s to the selected years, have been observed. The highest incremental percentage (+4.87%) was found in 1990-2006 (differences: 127,494.39 ha). After 2006, the forest areas slightly decreased, which was the most pronounced in 2018 (−0.98%), cf. Table 2. Until 2017, the net annual increment was consistently higher than the total felling. However, from 2018 onwards, the total felling exceeded the net annual increment, due to a high demand for timber in a small region [98,100], exacerbated by massive salvage logging that was caused by bark beetle infestation of the coniferous trees [83,111]. These problems have caused the country to experience economic, social, and environmental losses. During the last few years, bark beetles have been reported to have exponentially increased in Europe, due to climate change [112][113][114]. Furthermore, the state owns more than half (54.09%) of the Czech forests, and most of the forests are managed by the state forest enterprises. The remaining forest land is owned by private entities (19.27%), municipalities (17.17%), communes and associations (5.00%), as well as legal entities (3.27%) [98]. The state forests applied a more bureaucratic system to forest management, and were the role model for the private forests. Since the bark beetle infestation has caused more economic damage for the private than the state forests, the private owners need to take action earlier, e.g., by diversifying their products. The country also provided subsidies, compensation, and recovery payments of EUR 259 million in 2018-2019 [115]. According to the national accounts in 2019, private owners received CZK 8,046,000 for landscape and biodiversity protection, including, for example, the restoration and preservation of habitat and species, natural and semi-natural types of landscape, and environmental stability elements [56].
The larger area of forest (range: 0.72 to 1.90) reported in the CLC in the selected years, compared to the national report by the Czech MoA, and the smaller area of water bodies reported in the CLC (58,294 ha), compared with the land balance report (167,000 ha), can be attributed to the following two reasons: (1) the minimum mapping unit of 25 ha in the CLC, and (2) the differences in classification categories. The significant disparities in the total areas might have also been caused by the fact that the mixed-forest areas in the CLC were reported as either coniferous or broad-leaved forests in the national report. The land use balance only includes surface water bodies, while the CLC includes watercourses as well. Complementary data, with a higher level of detail, are required in different periods to be able to compare the changes in the water bodies with greater certainty.
The slightly decreased forest areas after 2012, and even more in 2018, compared to the year 2006, can be attributed to the bark beetle calamity experienced by the country, especially after 2017, which resulted in a massive clear-cutting of coniferous trees [116].
In regards to the increasing water body area in both the CLC and the official land balance report, this phenomenon is surprising compared to the global trend of the reduction in freshwater area due to natural phenomena [117], and the pressure exerted on freshwater worldwide, since almost 70% is affected by the industrial sector, through its use and pollution [118]. Additionally, it is widely known that Czechia is highly dependent on rainfall and runoff conditions [66]. Between 2014 and 2018, the reported surface water levels were affected by the accumulated drought of previous years; in 2018, the total precipitation was the lowest since 1961 [111]. Thus, it is presumed that the water bodies included in the CLC and the national report include artificial structures that have been used to improve the management of waterways since 1991, as water channels for irrigation or drainage fields [65,111], which distorts the mapping of the area occupied by water bodies.

Forest-Water Interaction and Public Perception
Forest-water interactions have been reported to provide essential and complex pathways for optimizing land use practices and water availability. The water system is susceptible to climate change; as the Earth's surface warmed, freshwater experienced, for example, a high rate of evaporation, a change in precipitation quantity and quality, and modifications of soil strength and biodiversity, which altered the quality and quantity of the groundwater. The frequency of flooding has also changed [119][120][121]. Forests are closely related to water, and also to the changes in land cover locations, such as riparian areas, alluvial plains, and wetlands. Any modifications in the quality of the structure and functioning of the landscape reduces the diversity and ecological processes taking place in an ecosystem [35,60,61,121]. Forests contribute to the regulation of the water cycle. Excessive rainwater is channeled along branches and trunks, reducing erosion of the soil, and the concentration of leaves in the soil facilitates the infiltration of water into the water table [35,122]. Water quality improves in the forest because there is less discharge of pollutants compared to industrial, settlement, and agricultural areas. The forest root network acts as a filter, reducing the levels of pollutants [121,123]. Additionally, the composition and state of the forest can vary the effects on water, e.g., heterogeneous forests have higher interception and redistribution of rainfall [122]. In cases such as this, sustainable forest management is expected to restore and improve groundwater recharge, global surface cooling, and water flow.
The low score of the FES performance evaluation was significantly associated with low education levels, and negatively associated with the age of the respondents (p = 0.001). A higher education level was found to have a link with positive behavior towards forest conservation [124]. The results of the survey implied the importance of forest education, especially on the role of the forest ecosystem with water provisioning services, to various audiences with different education levels, for supporting the forest policy to protect these areas, so they can function properly. The results indicate the importance of promoting the role of the Czech FES in water provisioning services among the older population and lower education levels, by improving the knowledge and awareness of the role of the Czech FES. The social protection of the forests will also increase, e.g., willingness to protect the trees from extended logging activities, and providing support to mixed forests to increase the health and vitality of the forests from insect attack [35,39].

Water Quality-Instruments Interaction and Public Perception
Biological oxygen demand (BOD) and chemical oxygen demand (COD) have been globally acknowledged as indicators for organic load in wastewater [125,126]. Nitrogen totals (N total) and phosphorus totals (P total) may suggest nutrient stability and indicate eutrophication [127][128][129]. N total and P total are typical representative parameters of domestic wastewater [130] related to the human excreta [67], as well as household chemicals, such as detergents [131], and substances related to agricultural practices [132], while undissolved substances may be closely related to the existence of suspensions in water [133].
High organic loads in South Moravia (BOD 26,180 tons/year in 2019) may be caused by the high production of wastewater related to both demographic and industrial issues. Figure 3a reveals that the drop in BOD load occurred in 2010 for WWTPs in South Moravia, but there was only a slight decrease in COD load in the same year. This is expected to be due to the presence of potential substances that are not simply degraded by biological activities, since BOD indicates the level of biodegradable fraction in wastewater [134]. Further regional observation and confirmation regarding this assumption are required to identify certain relevant practices that this load might have potentially originated from. Among the studied regions, South Moravia represented the highest population density in 2019 [76] (165.2 km 2 ), indicating a higher production of domestic wastewater. Meanwhile, industrial activities may generate a high load of COD (61.014 tons/year in 2019) containing compounds that are not readily biodegradable and inhibit the activity of aerobic bacteria. In view of the fact that industry constitutes the most essential sector in South Moravia, there is a tendency to enhance the higher value of production, due to shifts in structures and businesses, resulting in a higher intensity of industrial activities. Moreover, foreign investors predominantly allocate their investments in Brno and Modřice-the locations with the most extensive industrial terminals in South Moravia [135].
Generally, South Moravia tends to take a lead in terms of the load values of undissolved substances during the studied period, which predictably occurred due to the fact that more inhabitants lived in residential areas connected to the sewerage system from 2009 (87.9%) to 2019 (90.3%) than in Central Bohemia (66.5% in 2009 and 74.4% in 2019) [91,136]. If a larger number of inhabitants linked to the sewerage system is associated with a presumption of more wastewater discharged to the WWTPs inflow, it leads to an assumption that the WWTPs examined in South Moravia received a heavier load of undissolved substances compared to Central Bohemia. However, for further observation, we would require the information related to the sampling season frequency in WWTP per year, as the presence of undissolved substances is also influenced by dry or rainy weather [137,138]. Concerning nitrogen, South Moravia showed higher values, similar to the other parameters, possibly due to a greater number of inhabitants living in sewage system-equipped municipalities (1,073,942) [139], leading to the higher production of discharged nutrient loads.
Meanwhile, total phosphorus is frequently appraised as an indicator of the nutrients transported by runoff (dissolution or erosion), mainly from agricultural activities [132,140]. To the best of our knowledge, the assumption regarding the pattern of P total load in Figure 4f is associated with the mode of agricultural management and practices, since it can influence the fluctuation of the phosphorus state in the soil [140]. In addition, the increase in total P load in WWTPs inflow, starting from 2011 to 2019, is also predicted to be due to the increase in agricultural area in both the regions between 2006 and 2018 (based on our data analysis from CLC database), indicating that more agricultural practices were intensified. A specific evaluation of agricultural management, practices, policy, and sampling method from both the regions should be taken into account to support our prediction regarding the pattern of P total load. The analysis of water quality parameters allows the visualization that the water treated in the plants can reduce its pollutant load; however, the quality of the other water bodies that cover the Czech Republic is unknown, and that can affect the provision of freshwater. According to the annual water reports, regions with polluted to highly polluted water quality are identified [77,78,141], especially near urban settlements and industrial zones, which can compromise the provision of freshwater and other services [25,30,35].
The removal efficiency average of all the parameters that are comparatively examined for Central Bohemia, South Moravia, and the nationwide Czech Republic encountered high average values, proven by the ability to fit removal threshold requirements for all the measured parameters, adjusted by UWWTD. Despite the concern of South Moravia possessing a low level of water quality, as observed by several reports [48,78], the data analysis in our study pointed out that this region occupied the highest value of RE average over a period of time, compared to Central Bohemia, concerning its performance in achieving a higher RE for all the measured parameters. A study by Janosova [48] reported that a WWTP in Kyjov, South Moravia, employed physical and biological processes for the pretreatment of collected sewage and combinative systems, such as a circulating activation tank, secondary settling tanks, and a pumping installation used for recirculated and excess sludge. This set of techniques can be associated with the removal efficiency analysis in South Moravia.
However, further examination is necessary to identify the detailed specification of the treatment appliances of the WWTPs applied on a regional scale, corresponding to the raw data, as well as their retention time processing and the state of the sampling season regime.
According to previous results, the significant improvement of water treatment systems and the increase in the number of WWTPs can be attributed to a series of environmental policies and economic instruments implemented to conserve natural resources and protect the landscape [70,142,143]. Specifically, compliance with the UWWTD, to modernize and build new wastewater treatment plants before 2010, was a requirement that was established upon the entry into the EU. Since 1997, large financial investments have been required, mainly supported by the Instrument for Pre-Accession Structural Assistance, and structural and cohesion funds, before and after entering the European Union, respectively [143,144], to comply with the EU Directives on water.
Despite not finding a significant correlation between the environmental protection investment and removal efficiency of wastewater treatment in WWTPs, it is evident that the investment used in the infrastructure increases the capacity of WWTPs in the country, as well as the level of compliance with UWWTD, as confirmed by the OECD in its report [71]. However, it was identified that the environmental investment reports lack specific information on the costs and expenses of the specific activities that are developed for wastewater management. More studies are required that use detailed information on investments to verify, with greater certainty, the effectiveness of these investments.
The challenges due to climate change require investments to continue to implement more and better measures of mitigation and adaptation, as structures for water storage, innovation for the use of ground and surface water, capture and reuse of wastewater, protection of sources from which the water supply originates, and alignment with the Sustainable Development Goals on clean water, sanitation, and water life [117].
Environmental investment in processes that improve, treat, or remedy the physicochemical state of water is considered essential if properly developed [145]. Additionally, it is vital to develop land use and water policies that integrate people's perceptions, maintaining a balance between ecosystem services, users, and economic development [61]. In addition, the general public and two case study areas (the Central Bohemian and South Moravian regions) perceived that the Czech forests could still provide and deliver a good clean water source (range of mean was 3.3-3.4 from the total four scores). The predictor analysis, using binary logistic regression, revealed that the respondents' age had a negative and significant correlation with the low score of perception on the Czech forest ecosystem's ability to provide a clean water source. Community perception could provide current insights and firsthand experiences of the people [51,52], as our preliminary results of public perception showed that the younger generation seemed to have an unfavorable view of the forest ecosystem's ability to provide clean water. Future research is needed to investigate the reasons for their opinions.

Conclusions
The change in forest cover, based on the CLC database, in 1990, 2000, 2006, 2012, and 2018 presented the evolution of forests in the Czech Republic, and pointed to the importance of addressing the decline in forest areas in 2018, which was confirmed by the current national forest report. Two of the main causes were the bark beetle incursion and weather events, such as the Kyrill windstorm of 2007. The decline in forested landscapes has been reported in the case study regions as well.
The water body area showed an increase between 1990 and 2018. It is presumed that the CLC database and the land use balance report included artificial structures that distorted the mapped area. However, the comparison of water reports allowed us to identify that between 2014 and 2018, the levels of surface water were affected by the accumulated drought from the previous years. The relatively stable forest areas in the Central Bohemian and South Moravian regions were positively linked to the water bodies.
The expectations and perceived performance of the community were used to understand the relationship between the FES and the provision of freshwater. It was identified that respondents with a low educational level and an older population score provided a low-performance valuation of the FES, which calls attention to the importance of integrating the forest ecosystem with the water supply services to encourage various audiences, with different education levels, to support forest policies to protect forest areas, so that they can function properly.
The Czech Republic and the study regions met the removal efficiency or load reduction threshold among the measured parameters suggested by the UWWTD (BOD, COD, undissolved substance, N total, and P total), by more than 75%, 70-90%, 50%, 75%, and 75%, respectively, in wastewater treatment plants. Also, the number of WWTPs increased from 2158 in 2009 to 2759 in 2020.
The significant improvement in effluence quality and the increased number of WWTPs can be attributed to environmental protection expenditures that were created to implement environmental policies in wastewater management. Therefore, it is suggested to continue investing in controlling water pollution, building new wastewater treatment plants, and improving the existing ones to reduce the harmful risks for hydrological ecosystem services.
The environmental policies can significantly influence the services of hydrological ecosystems, according to land use, investment programs and the development of sustainable activities, directives to control pollution, and initiatives that promote knowledge transfer to generate awareness and commitment of society. In this sense, it is expected that the framework of European and Czech environmental policies will be maintained, where short-and long-term strategies establish an increase in forest cover, public participation, and investment in technological improvement, to reduce pollution in water and other resources.
It was identified that the hydrological ecosystems of the Czech Republic and the study areas present natural and anthropic pressures, such as the bark beetle infestation, climatic events that affect the water regime, and the availability of the resource. Similarly, it was identified that water resources present a high load of pollutants before being treated in the WWTPs, which modifies the conditions and availability of the water bodies that do not enter the water treatment system. Finally, it was identified that some groups in society lack knowledge about the benefits and effects of the forest in providing services, which could eventually put pressure on the resource, by not valuing its ecosystem importance.
The study provides information on multi-pronged analyses that can be developed to investigate the interactions between the components of nature, society, and the economy, allowing decision makers to use their findings to promote management and investment measures according to the needs of hydro-ecosystem services. Similarly, the findings can be used to promote inter-ministerial cooperation, to develop activities, spend budgets more effectively, and improve the current state of resources that provide ecosystem services to society. Additionally, this kind of analysis opens up possibilities for future research that promotes understanding the complexity of ecosystems, and identifying the changing needs of forests and water resources for sustainable development. Data Availability Statement: CLC database is available on https://land.copernicus.eu/pan-european/ corine-land-cover (accessed on 1 June 2021); water quality and environmental accounts are available on https://www.czso.cz/ (accessed on 15 June 2021).