Next Article in Journal
Sampling Uncertainties of Long-Term Remote-Sensing Suspended Sediments Monitoring over China’s Seas: Impacts of Cloud Coverage and Sediment Variations
Next Article in Special Issue
How Initial Forest Cover, Site Characteristics and Fire Severity Drive the Dynamics of the Southern Boreal Forest
Previous Article in Journal
Using ALS Data to Improve Co-Registration of Photogrammetry-Based Point Cloud Data in Urban Areas
Previous Article in Special Issue
Automatic Windthrow Detection Using Very-High-Resolution Satellite Imagery and Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landsat-Based Indices Reveal Consistent Recovery of Forested Stream Catchments from Acid Deposition

1
Global Change Research Institute of the Czech Academy of Sciences, Bělidla 986/4a, 603 00 Brno, Czech Republic
2
Laboratory on Geoinformatics and Cartography, Department of Geography, Faculty of Science, Masaryk University, Kotlářská 2, 61137 Brno, Czech Republic
3
Czech Geological Survey, Klárov 3, 118 21 Prague, Czech Republic
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(12), 1944; https://doi.org/10.3390/rs12121944
Submission received: 20 May 2020 / Revised: 8 June 2020 / Accepted: 14 June 2020 / Published: 16 June 2020
(This article belongs to the Special Issue Remote Sensing of Natural Forest Disturbances)

Abstract

:
Central European forests suffered from severe, large-scale atmospheric depositions of sulfur and nitrogen due to coal-based energy production during the 20th century. High deposition of acid compounds distorted soil chemistry and had negative effects on forest physiology and growth. Since 1994, continuous data on atmospheric deposition and stream runoff fluxes have provided evidence of ecosystem recovery from acidification. In this study, we combined for the first time mass budget data (sulfur deposition and total dissolved inorganic nitrogen (DIN) export) from the GEOMON monitoring network of headwater catchments with annual trajectories of vegetation indices derived from Landsat remote sensing observations. Time series of selected vegetation indices was constructed from Landsat 5, 7, and 8 using Google Earth Engine. Linear regression between the field data and vegetation indices was analyzed using R software. Biogeochemical responses of the forested catchment to declining acid deposition (driven by SO2 emission reduction) were consistent across all catchments covering various forest stands from different regions of the Czech Republic. Significant correlations were found with total sulfur depositions, suggesting that the forests are continuously and consistently prospering from reductions in acid deposition. Disturbance index (DI) was the only vegetation index that was well-related to changes in forest cover associated with salvage loggings (due to the forest decline) during the 1980s and 1990s. A significant relationship (R2 = 0.82) was found between the change in DI and DIN export in stream water. Regrowth of young forests in these highly affected areas tracks the most pronounced changes in total DIN export, suggesting a prominent role of vegetation in nitrogen retention. With the Landsat-derived DI, we could map decennial changes in forest disturbances beyond the small scale of the catchments to the regional level (demonstrated here for two protected landscape areas). This analysis showed the peak in forest disturbances to have occurred around the mid-1990s, followed by forest recovery and regrowth. Despite the improvement in forest ecosystem functioning over the past three decades in mountainous areas, emerging threats connected to changing climate will shape forest development in the near future.

Graphical Abstract

1. Introduction

Forests cover about 40% of the European Union (EU) territory [1], are regarded as the largest terrestrial ecosystem type, and accommodate a large proportion of European biodiversity. Forests provide a multitude of ecosystem services to humans, including climate regulation, timber production, water supply and regulation, clean air, and many others [2,3]. During the past three centuries, European forests have gone through fundamental changes in species composition, and most of the tree populations have been strongly affected by human activities [4]. In the Central European region, forests were largely converted to Norway spruce (Picea abies L., Karst) monocultures. Norway spruce is one of the economically important tree species, but, at the same time, it is vulnerable and susceptible to various biotic, abiotic, and anthropogenic disturbances when planted outside its natural habitats [5,6,7]. During the 20th century, Central European forests suffered from the severe, large-scale deposition of sulfur (S) and nitrogen (N) [8,9,10] that distorted soil chemistry and had negative effects on tree physiology and growth [11,12,13,14]. Among the highest deposition loads are those observed in the area of the Czech, Polish, and German borders, known as the “Black Triangle” [15] due to its high concentration of coal combustion power plants and topography that favored prolonged inversion events. Nevertheless, the Czech Republic is part of a region with the most pronounced decrease in regional SO2 emissions [10,16].
In Central Europe, the retreat of acid deposition started already in the mid-1980s and accelerated in the 1990s after the implementation of “clean air” legislation. The decline of regional SO2 emissions by 90% was followed by a consistent regional decline in sulfur (S) deposition [17]. Despite the modest reduction of N deposition in bulk precipitation, throughfall N fluxes remained generally unchanged. Nevertheless, acid-sensitive catchments started to recover from anthropogenic acidification, and nitrogen inorganic losses were substantially reduced [18], following the retreat of acid deposition, pointing towards better ecosystem immobilization of this limiting nutrient [19]—either by forest vegetation and/or by soil microbes.
In order to study the large-scale effects of reduced acid deposition on soil and forest recovery, a long-term observation network of 14 forested stream catchments, known as GEOMON (GEOchemical MONitoring), was established in the Czech Republic in 1994 [20]. Data from the monitoring program have served numerous studies to provide evidence for stream recovery from acidification [18,21,22] and reveal declines in stream eutrophication [23], increases in stream dissolved organic carbon [24], and effect of climate change on runoff patterns [25]. None of the studies to date have linked the long-term field observations with remote sensing (RS) data, which was the primary goal for this study.
Although the term “forest health” is frequently used, its definition is complex and might be ambiguous [26]. Despite this, many satellite-based RS methods for monitoring of forest health indicators have been developed (reviewed by [27,28]). According to Coppin et al. [29], those methods offer (1) consistent and repeatable procedures, (2) less expensive ways of obtaining information about land surface characteristics, and (3) possibilities for incorporating information from nonvisible parts of the electromagnetic spectrum. Landsat is one of the most used satellite missions for studying forest health properties [26], disturbance [30], and recovery [31]. Its main advantages are good coverage (global coverage with a revisit time of 16 days), sufficient spatial resolution (30 m pixel size), and, perhaps most importantly, its rich archive containing all images acquired since 1972 that allows for studying spatiotemporal changes in ecosystems for periods longer than 30 years [32,33]. Data from Landsat in this field have been incorporated successfully into studies dealing with forest health and structure on various scales ranging from global [34] to continental [35,36] and regional [31,37].
Many studies evaluating forest disturbances have been carried out using Landsat time series [38,39,40]. Usually, these take into account such large-scale ecological factors as bark beetle outbreaks [41], forest wildfires [36], windthrows [42], and harvesting [36,43]. Much less information can be found about the effects of disturbances caused by changes in biogeochemistry, such as soil acidification related to air pollution [44]. These continuous and subtler disturbances can be difficult to assess via RS because there are many environmental variables altering the RS product.
For studying forest disturbances, a disturbance index (DI) was proposed by Healey et al. [45]. It uses Tasseled Cap transformations of Landsat bands and their combinations and has been tested in both Russia and the USA, where it produced more accurate change classification than did original Landsat reflectance data. The DI was subsequently used in other studies located in regions with forest disturbances. Using linear regression, Eshleman et al. [30] found a relationship between the change in total dissolved nitrogen concentration and change in forest DI. A slightly different approach was carried out by Deel et al. [46], who introduced a cumulative version of DI for assessing the impacts of disturbance on the functioning of a forest. They found that this method produced significant relationships with percent canopy cover at plot level and with canopy percent N at catchment level.
The objective of the present study was to evaluate trends in Landsat RS-based indicators in relation to long-term field observations of total sulfur deposition (throughfall fluxes), as a regionally and temporary consistent measure of acid deposition and dissolved inorganic nitrogen (DIN) stream export, as a measure of nutrient retention capacity of forested catchments. We hypothesized that variation in stream water DIN leaching should track i) forest disturbance severity and ii) the extent of forest recovery from acidification.

2. Materials and Methods

2.1. Study Sites

This study focused on small stream catchments distributed within the GEOMON network [18,20], established in 1994 to study changes in precipitation and stream chemistry due to reductions in SO2 and NOx + NH3 emissions and their impacts on forest ecosystem functioning. Currently, data from 23 years of continuous monitoring are available to study ecosystem biogeochemical changes [18]. The network consists of 15 catchments, representing a variety of environmental conditions in the Czech Republic. Two catchments (Litávka (LIT) and Spálenec (SPA)) were excluded from this study due to limited temporal data records. Therefore, 13 catchments were considered. Their names (and three-letter abbreviations) are presented, and characteristics described in Table 1, and their locations are shown in Figure 1. The catchments are predominantly covered by Norway spruce (Picea abies L. Karst) forests. The only catchments with substantial coverage of broadleaf forest (European beech, Fagus sylvatica L.) are JEZ and LES. Minor representation of broadleaf trees can be found in the POM catchment. The catchments’ sizes range from 22 ha (PLB) to 262 ha (MOD), and they are situated mainly in montane and submontane areas (often being parts of protected landscape areas). Most of the catchments are climatically and lithologically predisposed to be sensitive to acidification, with the exception that PLB is well buffered by its ultrabasic bedrock [47]. The three catchments UDL, UHL, and JEZ were most exposed to direct air pollution in the second half of the 20th century and suffered from large-scale deforestation that reached its maximum in the mid of 1990s [44,48].

2.2. Field Data

Input fluxes (atmospheric deposition) and stream water runoff fluxes (output) were regularly monitored on a monthly basis using standardized field and analytical methodology in the GEOMON network. Regular monitoring at all catchments started since 1994; however, on some sites (JEZ, LYS, PLB, SAL), measurements started even earlier. Both open area (bulk) and throughfall precipitation were collected using polyethylene collectors. Stream water runoffs were computed by means of V-notch weirs or flumes with established rating curves. Although water samples were analyzed for detailed chemical composition (summarized in Oulehle et al. [18]), for the purposes of this study, we focused solely on total sulfur deposition derived from the throughfall precipitation fluxes (further denoted as S-input), collected bellow tree canopies, and DIN stream export (further denoted DIN-output). Ammonia concentration contributed minor to the inorganic nitrogen fluxes in stream water; therefore, DIN export equaled NO3 export. Sulfate and nitrate concentrations were measured by ion chromatography (Knauer 2000). Monthly observations were aggregated into annual fluxes (kg ha−1 year−1) covering the period 1994–2016, by multiplying the volume-weighted concentration of solutes by water fluxes.

2.3. Landsat Data and Annual Trajectories

The processing workflow is summarized in Figure 2. The time series of Landsat images for 13 studied stream catchments was accessed and analyzed using Google Earth Engine (GEE, [49]; https://earthengine.google.com). We used existing GEE image collections of surface reflectance products from Landsat-5 thematic mapper (TM), Landsat-7 enhanced thematic mapper plus (ETM+), and Landsat-8 operational Landsat imager (OLI). Landsat 5 had started acquiring images in 1984, preceding our field observations by 10 years. We used this time period to assess the prior behavior of the forest ecosystems before the field observations within the GEOMON network were begun and to analyze decennial changes in regional forest conditions starting from 1984.
Landsat images were first filtered according to the following criteria: i) images only from the peak vegetation season June to September were included, ii) cloudy or partially cloudy images were removed, iii) images with more than 25% of no data were removed from the study (this applied to Landsat 7 data only because of the scan line corrector failure). The numbers of valid images varied substantially between the catchments (Table 1). We had only 69 images for MOD, which is located in the top parts of the Giant Mountains (Krkonoše) and has frequent cloud cover, but as many as 336 images for ANE, which is a small catchment located in the central part of the Czech Republic having less cloud cover, and additionally is covered by two Landsat footprints (path/row: 191/25 and 191/26). There was a considerable gap in available Landsat images over the Czech Republic between 1996 and 1999 caused by data transfer regulations that affected many regions of the world [50]. All available images per catchment were used to compute zonal statistics for each year, consisting of median and standard deviation for each band and selected vegetation indices (Table 2). Median was used instead of mean because it was not susceptible to outliers, as could be the case when dealing with large amounts of RS data. In this way, annual trajectories for the period 1994–2016 were composed of each catchment from spectral reflectance and vegetation indices that were selected to be sensitive to forest health changes.
One of those vegetation indices (VIs) evaluated was the DI proposed by Healey et al. [45]. It was computed from a Tasseled Cap transformation, consisting of a linear combination of six Landsat bands, into the three components—brightness (B), greenness (G), and wetness (W). The DI was calculated as a linear combination of those three components after they had been first rescaled (Br, Gr, Wr) using their mean (Bμ, Gμ, Wμ) and standard deviation (Bσ, Gσ, Wσ) values computed from reference populations of pixels (Equations (1)–(4)).
B r = ( B B μ ) / B σ
G r = ( G G μ ) / G σ
W r = ( W W μ ) / W σ
D I = B r ( G r + W r )
The reference populations in this step consisted of pixels for all catchments in 1984. This was done so that the DI could be compared throughout all the catchments with the forest condition in the first year of images acquisition. Although on a larger scale, DI should be computed only over forested pixels, in our study, we regarded all pixels as forested because all catchments with the exception of MOD are almost completely covered by tree canopy.
Other vegetation indices (Table 2) were chosen for their sensitivity to different canopy characteristics. These included the widely used normalized difference vegetation index (NDVI), which is sensitive to canopy greenness; enhanced vegetation index (EVI), which was designed to be more sensitive to changes in areas having high biomass and to reduce the influence of atmospheric conditions; plant senescence reflectance index (PSRI), which is sensitive to leaf senescence by estimating the carotenoids and chlorophyll ratio; transformed chlorophyll absorption reflectance index (TCARI), which is sensitive to chlorophyll pigments.

2.4. Data Analysis

First, annual trajectories of VIs were plotted as normalized z-scores and visually compared with trajectories of the field-measured S-input and DIN-output at catchment level. Second, the Pearson correlation coefficient was computed for each catchment between field measurements and each individual vegetation index time series data (p > 0.05). Third, adopting the method of Eshleman et al. [30], we computed the change in DI and field-measured S-input and DIN-output as the difference between the linear trend line values at the end and beginning of the observation period. Using an ordinary least-squares regression, we related changes in forest conditions as measured by the difference in DI and measured sulfur deposition and stream water DIN. The linear trend of those scatterplots and R2, indicating the strength of the given relationship, were calculated. All these statistical analyses were done in R software [55].
Furthermore, we used DI to analyze decennial changes in forest condition since the beginning of Landsat observations at the regional level. For this purpose, we selected two protected landscape areas: Eagle Mountains (Orlické Hory) hory, where the UDL catchment is located, and Slavkov Forest (Slavkovský les), where the PLB and LYS catchments are located. The two protected landscape areas were selected because UDL and PLB were identified (see Results) to be the most contrasting catchments in terms of forest condition changes. Cloud-free images during peak vegetation season were chosen in 1984, 1994, 2003, 2014, and 2019 for Eagle Mountains and in 1985, 1995, 2003, 2015, and 2019 for Slavkov Forest. In just one case, we were unable to find a completely cloud-free scene (for Slavkov Forest from around the year 2013), so we chose the least cloudy image from that year. Unlike how we handled DI for small catchments, here we selected only forested pixels by applying a mask derived from forest index. Forest index introduced by Ye et al. [56] was computed from 2004 (Eagle Mountains) and 2005 (Slavkov forest) Landsat images. The mask was created by applying an empirical threshold of the forest index (3.6 for Eagle Mountains, 4.0 for Slavkov forest). By doing this, we were able to directly compare every pixel and, furthermore, we were able to distinguish pixels that had been deforested in the past (direct pixel comparison would not have been possible if we would have used a different forest mask for each individual image because the deforested pixels would in that case not have been included).
We calculated percentages of growth/recovery, no change, and disturbed pixels—both for the entire protected landscape areas and for smaller catchments located within them. Pixels in images were classified into three classes using mean (μ) and standard deviation (σ): class 1—growth/recovery (> μ + 0.5σ), class 2—no change (between μ + 0.5σ and μ − 0.5σ), and class 3—disturbed (< μ − 0.5σ) [57].

3. Results

3.1. Long-Term Trends in Measured Sulfur Deposition and Stream Water Dissolved Inorganic Nitrogen Export

Standard scores for bulk sulfur deposition (S-input) and stream water nitrogen export (DIN-output) across all GEOMON catchments are shown in Figure 3. It showed a strong decline in S-input during the 1990s. The most pronounced decline of sulfur deposition was observed at the JEZ, UDL, and UHL, which were those most exposed to acid deposition. The maximum deposition rate for JEZ was 90 kg S ha−1 year−1 measured in 1995, and, in 20 years, until 2016, it decreased to 13.6 kg S ha−1 year−1. Similarly, the maximum deposition for UDL was 84 kg S ha−1 year−1 measured in 1998, and the lowest was 9.9 kg S ha−1 year−1 in 2012. At UHL, the maximum was 72 kg S ha−1 year−1 in 1995, and a minimum was 6.1 kg S ha−1 year−1 in 2011. This declining trend was consistently observed also at all other sites. In less polluted regions of the central part of the Czech Republic, the maximum value of 23.7 kg S ha−1 year−1 was observed at ANE in 1996, and it decreased to 6.8 kg S ha−1 year−1 in 2015.
In general, the stream water DIN-output also exhibited a decreasing trend, but the variability between the years and between catchments was larger than for S-input (as indicated by the wider shaded area in Figure 3b than in 3a). Especially at the UHL, UDL, and MOD, the maximum DIN-output values (maximum 13.2 kg N ha−1 year−1 measured in 1994 at MOD, and a maximum of 23.7 kg N ha−1 year−1 measured in 1997 at UDL) were observed in the late 1990s, which well corresponded to forest deterioration in the top parts of the northern mountain ranges [58,59]. Some of the catchments (such as ANE, LKV, and SAL, located in the Bohemian-Moravian Highlands), however, exhibited almost no trend and variation between years in stream water DIN-output.

3.2. Long-Term Trends in Remote Sensing Indicators

Correlation coefficients between field measurements and Landsat VIs are presented in Figure 4 (crosses indicate non-significant correlations (p > 0.05)). The correlation coefficients were calculated for the period since 1994 (when the field measurements began). With the RS data, however, we were able to track forest changes since 1984.
Results showed that more significant correlations were found between S-input and VIs than between DIN-output and VIs. The S-input showed a significant correlation with DI (6 catchments), NDVI (10 catchments), and TCARI (10 catchments). Only DI showed a strong correlation (Pearson’s r > 0.75) on three catchments (LYS, UDL, and UHL). Fewer correlations were found for EVI and PSRI across the catchments. The most noteworthy catchments were, on the one hand, MOD and ANE, where almost no correlations with the VIs were observed, and, on the other hand, UDL and UHL, where correlations were observed with almost all VIs.
In general, almost no significant correlations were observed between VIs and DIN-output. The significant correlations were mainly identified for DI (LYS, POM, UDL, and UHL).
Based on the results of the correlation analysis, the most promising VIs for tracing forest changes in relation to S-input were NDVI and TCARI. Nevertheless, when inspecting annual trajectories of the individual VIs (Figure 5), NDVI showed a positive increasing trend for the 1984–2016 observation period for all the catchments, regardless of known deforestation at UDL and UHL catchments in the early 1990s. Figure 5 shows trend lines for standard scores between field measurements and Landsat VIs. For the sake of brevity, only seven catchments were shown here, as they well represented the geographic areas in the Czech Republic and overall variety between the catchments. Vegetation indices were shown from 1984, meaning that they preceded field observations of S-input and DIN-output by 10 years. The overall steeper slope could be seen in S-input, which also explained the stronger correlations for S-input seen in Figure 4.
Among all the VIs, we chose, for several reasons, to use DI for further analysis. For DIN-output, none of the indices were particularly successful, but DI showed a more significant correlation than did any of the other indices. What is more important, only DI was able to trace forest decline and partial deforestation at catchments most affected by acid deposition (namely JEZ, UDL, and UHL) that occurred before the uniform field measurements started (Figure 5). NDVI and TCARI were not sensitive to these changes. PSRI showed some sensitivity, but, overall, DI was the most sensitive index. Furthermore, by using DI, we were able to directly compare our results with those from other studies [30,44], also using this index for assessing forest disturbances. Focusing on the trend change within the entire observation period (1995–2016), we observed a positive linear relationship between the S-input and DI with R2 = 0.58 (Figure 6a). A similar relationship was found also between DIN-output and DI, but the regression coefficient was higher at R2 = 0.82 (Figure 6b).
Decennial changes in DI over the two protected landscape areas are presented in Figure 7 (Eagle Mountains with UDL) and Figure 8 (Slavkov Forest with PLB). In Eagle Mountains, we could observe a severe increase in DI that corresponded to larger forest disturbances located mainly at the mountain ridges, where UDL is located, between the years 1984 and 1994. On the contrary, no dramatic changes in forest disturbances were detected in the region of Slavkov Forest, except likely small-scale logging activities. Increased DI in 2013 in the northwest area was attributed to imperfections in Landsat image data, such as cloud cover. We classified the entire area into three categories and analyzed the proportion of disturbed, no change, and growth/recovery areas (summarized in Figure 9). We could observe the increase of disturbed pixels from 26.2% to 34.4% across the area of Eagle Mountains as a whole and from 27.7% to 94.5% in UDL between 1984 and 1994. In 2004, only 2.3% of pixels belonged to the growth/recovery category. By 2013, the proportion of those pixels had risen to 64.3%. Finally, in 2019, the growth/recovery category occupied 80.7% of the UDL catchment area.

4. Discussion

4.1. Acidification Retreat across GEOMON Stream Catchments

Precipitation acidity reversal was driven almost exclusively by the reduction in SO2 emission, and it was consistent across all monitored stream catchments. The decline in S deposition of almost 90% over the past 25 years is noteworthy and paralleled the consistent declines in stream DIN export and increase in precipitation pH [18]. Despite reductions in N deposition over this period, the magnitude of decreases in DIN stream concentration was larger, indicating increasing immobilization of N in GEOMON catchments [18]. The temporarily increasing N retention could have been caused by i) greater immobilization in forest biomass and/or ii) greater immobilization in the soil environment. If the former of these was the pivotal mechanism, then individual indices derived from RS data should be able to detect sink in forest biomass, as was demonstrated using DI for UDL and UHL catchments. Nevertheless, improvement to soil through root exudates from tree labile carbon supply might have a direct impact on nutrient immobilization in soil by stimulating soil organic matter processing and microbial N retention. The deterioration of tree growth and rapid improvement after a decline in an acid deposition have been recently revealed by extensive dendrochronological research undertaken in areas affected by acid rain [13,60,61].

4.2. Landsat Data Suitability

The opportunity to obtain cloud-free Landsat images for such relatively small study areas was generally good, although the number of images available varied greatly among the sites (Table 1). For some catchments located in mountainous areas with frequent cloud cover (such as MOD in the Giant Mountains), fewer suitable images were found. There was also a considerable gap in available Landsat images over the Czech Republic between 1996 and 1999 caused by data transfer regulations that affected many regions of the world [50]. For assessing forest disturbance on a larger spatial scale, however, and as we demonstrated here for two protected landscape areas, it was difficult to find a completely cloud-free scene for the entire area during the main growing season. Although there are other satellite sensors, such as MODIS (with a temporal frequency of 1–2 days), that allow for seasonal dynamics to be evaluated, data availability since 2000, together with the coarse spatial resolution of 250 m, made the MODIS product unsuitable for assessment in this study. The recently launched Sentinel-2 offers great potential, with enhanced spatial, spectral, and temporal resolutions compared to Landsat. Its potential for long-term trend studies is still to be realized because Sentinel-2’s consistent data acquisition began only in 2015. The incompatibility between field measurements and RS image data and frequent gaps throughout the year due to increased probability of cloud cover during the autumn and winter periods made it impossible to study both the trend between years and the seasonal variations as suggested, for instance, by Verbesselt et al. [62].
Computation of annual VI trajectories was fully automated in GEE and in a manner such that we were able to implement several measures to ensure a high degree of consistency between the images, cloud filtering, filtering of the Landsat 7 striping effect, and selection of median value that was less susceptible to outliers than mean would have been. Even all these measures, however, could not eliminate variability between the years. This likely could be attributed to inappropriate cloud mask and, thus, the potential inclusion of areas affected by artifacts of atmospheric corrections. Other reasons could be temporal inconsistency (even though we considered only images from the vegetation season (June–September)) or topographic effects [63]. Landsat data collections in the GEE area were processed to level-2 atmospherically corrected surface reflectance (Tier 1). Although Landsat 5–7 is atmospherically processed using the LEDAPS algorithm [64] and Landsat 8 data using the LaSRC algorithm [65], we did not observe any noticeable offset in yearly median values after 2013, for which Landsat 8 data are available.

4.3. Trends in Remote Sensing Indicators

Significant correlations were found between total S deposition and VIs. This was especially the case for NDVI, which combined joint temporal development of S deposition with concurrent aging of the forests and forest management interventions (namely large-scale harvests in selected catchments). The decline in S deposition was observed consistently for all catchments, and so was a consistent increase in NDVI, which was attributed to growth in forest biomass. NDVI trends extended by 10 years of RS observations prior to the field monitoring did not detect any known forest decline in the top parts of the Jizera (UHL) and Eagle (UDL) Mountains. We observed neither expected decrease in NDVI nor visible changes in EVI and TCARI, and this could be explained by the rapid growth of herbaceous understory. The general trend of NDVI nevertheless inversely mirrored the trend in S deposition (Figure 5), suggesting the likely positive role of declining precipitation acidity on forest growth. This view is further supported by a dendrochronological study done in the Giant Mountains [13]. Contrarily, DI trends exhibited a rapid increase (similarly as observable for the PSRI) around the mid-1990s and then a steady decrease corresponding to forest regrowth (Figure S1 in the online Supplementary Materials). DI was successfully applied by Kupková et al. [44] and Mišurec et al. [57] to study forest recovery in the Ore Mountains (Krušné hory) within the Czech Republic.
When observing a linear relationship between the change of DI and S-input or DIN-output (Figure 6), we could infer that the majority of catchments are situated in the middle cluster, with the exception of UDL and UHL, for which there was a significant decrease of DI in the studied time period (the DI was high there in 1995 because of the deforestation, which is documented in Figure 7 for UDL). Therefore, the relationships were driven by factors reliant primarily on forest management decisions (e.g., harvest) rather than on continuous improvement in forest status. In contrast, PLB is experiencing a subtle increase in DI due to recent selective logging activities connected to bark beetle infestation.
Our simple linear model, examining the relationship between change in DI and change in stream water DIN, was consistent with the results of Eshleman et al. [30], but it was created from a longer data time series (22 years instead of the 6 years studied by Eshleman). Although it showed an even stronger relationship, in drawing conclusions about the strength of the relationship, we must admit that this was mainly due to the influence of data for a few catchments situated outside the big cluster of catchments, showing no change in DI and stream water DIN export. The small number and area (<200 ha) of catchments did not allow for examining the effect of the spatial distribution of forest disturbances. This effect was assessed by Cowles et al. [66], who analyzed data from a large catchment with many embedded sub-catchments. It would be interesting to analyze the trends while including more catchments experiencing major positive or negative changes in forest conditions. This is likely to happen soon due to a rapidly spreading bark beetle outbreak that is destroying productive spruce forests in the Czech Republic [67]. Therefore, we could expect that catchments at middle altitudes (below 800 m) would sooner or later be largely deforested due to salvage fellings. Catchments located in the Moravian–Bohemian highlands (ANE, LKV, and SAL) were already been influenced during 2019. Unfortunately, field data on their S-input and N-output were not yet available at the time of this study. As stated above (in 4.1), the DIN export revealed more year-to-year variability compared to stream nitrate concentration due to the variable hydrological conditions. Oulehle et al. [18] showed a consistent decline of stream nitrate concentration across GEOMON catchments. Such a general decline that is steeper than the decline in a deposition is consistent with an overall improvement of N immobilization in forested ecosystems, either in vegetation (consistent with NDVI and TCARI observations) or in soil by higher N microbial immobilization [68].
Large-scale analysis of DI illustrated a similar trend of forest regrowth in the Eagle Mountains protected landscape area after the 1990s, as was observed by Kupková et al. [44] for the eastern part of the Ore Mountains. On those maps (Figure 7 and Figure 8), DI was also sensitive to the spatial resolution of the satellite system at the edges between the studied forest and other land cover types, such as bare soil, and, therefore, it is important to carefully select only forest pixels or pixels that had been forested in the past to compare those between the years. In future years, even better results will be obtained by using satellite systems with better spatial resolutions, like Sentinel-2.
An increase in the proportion of growth/recovery class, seen in Figure 9b, could be associated with stricter pollution control policies and desulfurization of coal combustion power plants. Both of these have made reforestation possible in the most affected areas. By contrast, an increase in disturbed pixels, seen in Figure 9d, could be attributed to sparse logging activities that occurred at our study site.

5. Conclusions

This study combined more than 20 years of field observations of total sulfur deposition and stream water dissolved inorganic nitrogen export from small forested catchments in the Czech Republic with vegetation indices derived from the Landsat time-series data.
Biogeochemical responses of forested catchments to declining acid deposition (driven by a reduction in SO2 emission) were consistent across catchments. The reduction of DIN export was not exclusively driven by a reduction in N deposition, but higher forest ecosystem retention was detected. Consequently, RS indices of forest cover properties showed marked changes over the past three decades. DI was well-related to changes in forest cover associated with forest management activities (i.e., salvage loggings during the 1980s and 1990s). Regrowth of young forests in these highly affected areas tracked the most pronounced changes in DIN export, thus suggesting the prominent role of vegetation in N retention. The observed temporal changes in DI on particular catchments were mirrored on the larger, landscape-scale, thereby pointing to the extent of forested areas historically affected by acid rain. Contrary to DI, the NDVI and TCARI indices linearly and consistently changed across a majority of catchments, independently of the mean age of forests within individual catchments. This might suggest that the forests are continuously and consistently prospering due to reduced acid deposition.
A combination of independent approaches using mass budget data and remote sensing indices complemented by recent dendrochronological studies (e.g., [13,60,69]) have revealed consistent improvement in forest status, including increased productivity and consequent N ecosystem retention, as a response to a reversal in acid deposition.
Despite the improvement in forest ecosystem functioning over the past three decades in mountainous areas, emerging threats connected to changing weather patterns will shape forest development in the near future. Our study demonstrated that RS is a vital approach for detecting adverse impacts of environmental changes on forest functioning at large (landscape) and small (headwater catchment) scales. Satisfactory detection of change magnitude is dependent on the proper selection of RS indices, with DI being most effective for changes in forest cover, whereas NDVI and TCARI are more suitable for detecting continuous changes associated with forest growth and/or productivity.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-4292/12/12/1944/s1, Figure S1: Time-series and trends of VIs, DIN-output and S-input for all catchments.

Author Contributions

Conceptualization, M.Š. and L.H.; formal analysis, M.Š.; funding acquisition, L.H.; methodology, M.Š.; project administration, L.H.; resources, F.O. and P.K.; software, M.Š.; supervision, L.H.; visualization, M.Š.; writing – original draft, M.Š., F.O., and L.H.; writing – review and editing, F.O., P.K., R.J., and K.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Czech Science Foundation, grant number 17-05743S and by the Grant Agency of the Masaryk University (grant number MUNI/A/1356/2019).

Acknowledgments

We would like to thank the anonymous reviewers whose comments have improved this manuscript.

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. Eurostat. Available online: https://ec.europa.eu/eurostat/web/products-eurostat-news/-/EDN-20180321-1 (accessed on 3 January 2020).
  2. Barredo, J.I.; Bastrup-Birk, A.; Teller, A.; Onaindia, M.; Fernández De Manuel, B.; Madariaga, I.; Rodríguez-Loinaz, G.; Pinho, P.; Nunes, A.; Ramos, A.; et al. Mapping and Assessment of Forest Ecosystems and their Services–Applications and Guidance for Decision Making in the Framework of MAES; Publications Office of the European Union: Brussels, Belgium, 2015. [Google Scholar]
  3. Brockerhoff, E.G.; Barbaro, L.; Castagneyrol, B.; Forrester, D.I.; Gardiner, B.; González-Olabarria, J.R.; Lyver, P.O.; Meurisse, N.; Oxbrough, A.; Taki, H.; et al. Forest biodiversity, ecosystem functioning and the provision of ecosystem services. Biodivers. Conserv. 2017, 26, 3005–3035. [Google Scholar] [CrossRef] [Green Version]
  4. Jansen, S.; Konrad, H.; Geburek, T. The extent of historic translocation of Norway spruce forest reproductive material in Europe. Ann. For. Sci. 2017, 74, 56. [Google Scholar] [CrossRef]
  5. Hlásny, T.; Turčáni, M. Persisting bark beetle outbreak indicates the unsustainability of secondary Norway spruce forests: case study from Central Europe. Ann. For. Sci. 2013, 70, 481–491. [Google Scholar] [CrossRef] [Green Version]
  6. Mezei, P.; Grodzki, W.; Blaženec, M.; Jakuš, R. Factors influencing the wind–bark beetles’ disturbance system in the course of an Ips typographus outbreak in the Tatra Mountains. For. Ecol. Manag. 2014, 312, 67–77. [Google Scholar] [CrossRef]
  7. Schelhaas, M.-J.; Nabuurs, G.-J.; Schuck, A. Natural disturbances in the European forests in the 19th and 20th centuries. Glob. Chang. Biol. 2003, 9, 1620–1633. [Google Scholar] [CrossRef]
  8. Rydval, M.; Wilson, R. The Impact of Industrial SO2 Pollution on North Bohemia Conifers. Water Air Soil Pollut. 2012, 223, 5727–5744. [Google Scholar] [CrossRef]
  9. Stoddard, J.L.; Jeffries, D.S.; Lukewille, A.; Clair, T.A.; Dillon, P.J.; Driscoll, C.T.; Forsius, M.; Johannessen, M.; Kahl, J.S.; Kellogg, J.H.; et al. Regional trends in aquatic recovery from acidification in North America and Europe. Nature 1999, 401, 4. [Google Scholar] [CrossRef]
  10. Vestreng, V.; Myhre, G.; Fagerli, H.; Reis, S.; Tarrasón, L. Twenty-five years of continuous sulphur dioxide emission reduction in Europe. Atmos. Chem. Phys. 2007, 7, 3663–3681. [Google Scholar] [CrossRef] [Green Version]
  11. Cienciala, E.; Altman, J.; Doležal, J.; Kopáček, J.; Štěpánek, P.; Ståhl, G.; Tumajer, J. Increased spruce tree growth in Central Europe since 1960s. Sci. Total Environ. 2018, 619–620, 1637–1647. [Google Scholar] [CrossRef]
  12. Hruška, J.; Cienciala, E. Long-term acidification and nutrition degradation of forest soils—A limiting factor of the present forestry. Czech Geol. Surv. 2005. [Google Scholar]
  13. Kolář, T.; Čermák, P.; Oulehle, F.; Trnka, M.; Štěpánek, P.; Cudlín, P.; Hruška, J.; Büntgen, U.; Rybníček, M. Pollution control enhanced spruce growth in the “Black Triangle” near the Czech–Polish border. Sci. Total Environ. 2015, 538, 703–711. [Google Scholar] [CrossRef] [PubMed]
  14. Santrůcková, H.; Santrůcek, J.; Setlík, J.; Svoboda, M.; Kopácek, J. Carbon isotopes in tree rings of Norway spruce exposed to atmospheric pollution. Environ. Sci. Technol. 2007, 41, 5778–5782. [Google Scholar] [CrossRef] [PubMed]
  15. Grubler, A. Trends in global emissions: Carbon, sulfur, and nitrogen. In Encyclopedia of Global Environmental Change; Munn, T., Ed.; Causes and Consequences of Global Environmental Change; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2001; Volume 3, ISBN 978-0-471-97796-4. [Google Scholar]
  16. Kopáček, J.; Posch, M.; Hejzlar, J.; Oulehle, F.; Volková, A. An elevation-based regional model for interpolating sulphur and nitrogen deposition. Atmos. Environ. 2012, 50, 287–296. [Google Scholar] [CrossRef]
  17. Oulehle, F.; Kopáček, J.; Chuman, T.; Černohous, V.; Hůnová, I.; Hruška, J.; Krám, P.; Lachmanová, Z.; Navrátil, T.; Štěpánek, P.; et al. Predicting sulphur and nitrogen deposition using a simple statistical method. Atmos. Environ. 2016, 140, 456–468. [Google Scholar] [CrossRef] [Green Version]
  18. Oulehle, F.; Chuman, T.; Hruška, J.; Krám, P.; McDowell, W.H.; Myška, O.; Navrátil, T.; Tesař, M. Recovery from acidification alters concentrations and fluxes of solutes from Czech catchments. Biogeochemistry 2017, 132, 251–272. [Google Scholar] [CrossRef]
  19. LeBauer, D.S.; Treseder, K.K. Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology 2008, 89, 371–379. [Google Scholar] [CrossRef] [Green Version]
  20. Fottová, D.; Skořepová, I. Changes in Mass Element Fluxes and their Importance for Critical Loads: Geomon Network, Czech Republic. Waterairsoil Pollut. 1998, 105, 365–376. [Google Scholar] [CrossRef]
  21. Krám, P.; Hruška, J.; Shanley, J.B. Streamwater chemistry in three contrasting monolithologic Czech catchments. Appl. Geochem. 2012, 27, 1854–1863. [Google Scholar] [CrossRef]
  22. Navrátil, T.; Kurz, D.; Krám, P.; Hofmeister, J.; Hruška, J. Acidification and recovery of soil at a heavily impacted forest catchment (Lysina, Czech Republic)—SAFE modeling and field results. Ecol. Model. 2007, 205, 464–474. [Google Scholar] [CrossRef]
  23. Oulehle, F.; McDowell, W.H.; Aitkenhead-Peterson, J.A.; Krám, P.; Hruška, J.; Navrátil, T.; Buzek, F.; Fottová, D. Long-Term Trends in Stream Nitrate Concentrations and Losses Across Watersheds Undergoing Recovery from Acidification in the Czech Republic. Ecosystems 2008, 11, 410–425. [Google Scholar] [CrossRef]
  24. Hruška, J.; Krám, P.; McDowell, W.H.; Oulehle, F. Increased Dissolved Organic Carbon (DOC) in Central European Streams is Driven by Reductions in Ionic Strength Rather than Climate Change or Decreasing Acidity. Environ. Sci. Technol. 2009, 43, 4320–4326. [Google Scholar] [CrossRef]
  25. Lamacova, A.; Hruska, J.; Kram, P.; Stuchlik, E.; Farda, A.; Chuman, T.; Fottova, D. Runoff Trends Analysis and Future Projections of Hydrological Patterns in Small Forested Catchments. Soil Water Res. 2014, 9, 169–181. [Google Scholar] [CrossRef] [Green Version]
  26. Lausch, A.; Erasmi, S.; King, D.J.; Magdon, P.; Heurich, M. Understanding Forest Health with Remote Sensing -Part I—A Review of Spectral Traits, Processes and Remote-Sensing Characteristics. Remote Sens. 2016, 8, 1029. [Google Scholar] [CrossRef] [Green Version]
  27. Pause, M.; Schweitzer, C.; Rosenthal, M.; Keuck, V.; Bumberger, J.; Dietrich, P.; Heurich, M.; Jung, A.; Lausch, A. In Situ/Remote Sensing Integration to Assess Forest Health—A Review. Remote Sens. 2016, 8, 471. [Google Scholar] [CrossRef] [Green Version]
  28. Lausch, A.; Erasmi, S.; King, D.J.; Magdon, P.; Heurich, M. Understanding Forest Health with Remote Sensing-Part II—A Review of Approaches and Data Models. Remote Sens. 2017, 9, 129. [Google Scholar] [CrossRef] [Green Version]
  29. Coppin, P.; Jonckheere, I.; Nackaerts, K.; Muys, B.; Lambin, E. Review ArticleDigital change detection methods in ecosystem monitoring: A review. Int. J. Remote Sens. 2004, 25, 1565–1596. [Google Scholar] [CrossRef]
  30. Eshleman, K.N.; McNeil, B.E.; Townsend, P.A. Validation of a remote sensing based index of forest disturbance using streamwater nitrogen data. Ecol. Indic. 2009, 9, 476–484. [Google Scholar] [CrossRef]
  31. Griffiths, P.; Kuemmerle, T.; Baumann, M.; Radeloff, V.C.; Abrudan, I.V.; Lieskovsky, J.; Munteanu, C.; Ostapowicz, K.; Hostert, P. Forest disturbances, forest recovery, and changes in forest types across the Carpathian ecoregion from 1985 to 2010 based on Landsat image composites. Remote Sens. Environ. 2014, 151, 72–88. [Google Scholar] [CrossRef]
  32. Cohen, W.B.; Goward, S.N. Landsat’s Role in Ecological Applications of Remote Sensing. BioScience 2004, 54, 535–545. [Google Scholar] [CrossRef]
  33. Wulder, M.A.; Masek, J.G.; Cohen, W.B.; Loveland, T.R.; Woodcock, C.E. Opening the archive: How free data has enabled the science and monitoring promise of Landsat. Remote Sens. Environ. 2012, 122, 2–10. [Google Scholar] [CrossRef]
  34. Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A.; Thau, D.; Stehman, S.V.; Goetz, S.J.; Loveland, T.R.; et al. High-Resolution Global Maps of 21st-Century Forest Cover Change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lehmann, E.A.; Wallace, J.F.; Caccetta, P.A.; Furby, S.L.; Zdunic, K. Forest cover trends from time series Landsat data for the Australian continent. Int. J. Appl. Earth Obs. Geoinf. 2013, 21, 453–462. [Google Scholar] [CrossRef]
  36. White, J.C.; Wulder, M.A.; Hermosilla, T.; Coops, N.C.; Hobart, G.W. A nationwide annual characterization of 25years of forest disturbance and recovery for Canada using Landsat time series. Remote Sens. Environ. 2017, 194, 303–321. [Google Scholar] [CrossRef]
  37. Potapov, P.V.; Turubanova, S.A.; Tyukavina, A.; Krylov, A.M.; McCarty, J.L.; Radeloff, V.C.; Hansen, M.C. Eastern Europe’s forest cover dynamics from 1985 to 2012 quantified from the full Landsat archive. Remote Sens. Environ. 2015, 159, 28–43. [Google Scholar] [CrossRef]
  38. Cohen, W.B.; Yang, Z.; Kennedy, R. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 2. TimeSync—Tools for calibration and validation. Remote Sens. Environ. 2010, 114, 2911–2924. [Google Scholar] [CrossRef]
  39. Pickell, P.D.; Hermosilla, T.; Frazier, R.J.; Coops, N.C.; Wulder, M.A. Forest recovery trends derived from Landsat time series for North American boreal forests. Int. J. Remote Sens. 2016, 37, 138–149. [Google Scholar] [CrossRef]
  40. Liang, L.; Hawbaker, T.J.; Zhu, Z.; Li, X.; Gong, P. Forest disturbance interactions and successional pathways in the Southern Rocky Mountains. For. Ecol. Manag. 2016, 375, 35–45. [Google Scholar] [CrossRef] [Green Version]
  41. Senf, C.; Pflugmacher, D.; Wulder, M.A.; Hostert, P. Characterizing spectral–temporal patterns of defoliator and bark beetle disturbances using Landsat time series. Remote Sens. Environ. 2015, 170, 166–177. [Google Scholar] [CrossRef]
  42. Oeser, J.; Pflugmacher, D.; Senf, C.; Heurich, M.; Hostert, P. Using Intra-Annual Landsat Time Series for Attributing Forest Disturbance Agents in Central Europe. Forests 2017, 8, 251. [Google Scholar] [CrossRef]
  43. Senf, C.; Pflugmacher, D.; Hostert, P.; Seidl, R. Using Landsat time series for characterizing forest disturbance dynamics in the coupled human and natural systems of Central Europe. ISPRS J. Photogramm. Remote Sens. 2017, 130, 453–463. [Google Scholar] [CrossRef]
  44. Kupková, L.; Potůčková, M.; Lhotáková, Z.; Albrechtová, J. Forest cover and disturbance changes, and their driving forces: A case study in the Ore Mountains, Czechia, heavily affected by anthropogenic acidic pollution in the second half of the 20th century. Environ. Res. Lett. 2018, 13, 095008. [Google Scholar] [CrossRef]
  45. Healey, S.; Cohen, W.; Zhiqiang, Y.; Krankina, O. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection. Remote Sens. Environ. 2005, 97, 301–310. [Google Scholar] [CrossRef]
  46. Deel, L.N.; McNeil, B.E.; Curtis, P.G.; Serbin, S.P.; Singh, A.; Eshleman, K.N.; Townsend, P.A. Relationship of a Landsat cumulative disturbance index to canopy nitrogen and forest structure. Remote Sens. Environ. 2012, 118, 40–49. [Google Scholar] [CrossRef]
  47. Krám, P.; Oulehle, F.; Štědrá, V.; Hruška, J.; Shanley, J.B.; Minocha, R.; Traister, E. Geoecology of a Forest Watershed Underlain by Serpentine in Central Europe. Northeast. Nat. 2009, 16, 309–328. [Google Scholar] [CrossRef]
  48. Vávrová, E.; Cudlín, O.; Vavříček, D.; Cudlín, P. Ground vegetation dynamics in mountain spruce (Picea abies (L.) Karsten) forests recovering after air pollution stress impact. Plant Ecol. 2009, 205, 305–321. [Google Scholar] [CrossRef]
  49. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  50. Goward, S.; Arvidson, T.; Williams, D.; Faundeen, J.; Irons, J.; Franks, S. Historical Record of Landsat Global Coverage: Mission Operations, NSLRSDA, and International Cooperator Stations. Photogramm. Eng. 2006, 15. [Google Scholar]
  51. Liu, H.Q.; Huete, A. A feedback based modification of the NDVI to minimize canopy background and atmospheric noise. IEEE Trans. Geosci. Remote Sens. 1995, 33, 457–465. [Google Scholar] [CrossRef]
  52. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the Great Plains with ERTS. In Proceedings of the 3rd ERTS Symposium, Washington, DC, USA, 10–14 December 1973; NASA: Washington, DC, USA, 1974. [Google Scholar]
  53. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Rakitin, V.Y. Non-destructive optical detection of pigment changes during leaf senescence and fruit ripening. Physiol. Plant. 1999, 106, 135–141. [Google Scholar] [CrossRef] [Green Version]
  54. Haboudane, D.; Miller, J.R.; Tremblay, N.; Zarco-Tejada, P.J.; Dextraze, L. Integrated narrow-band vegetation indices for prediction of crop chlorophyll content for application to precision agriculture. Remote Sens. Environ. 2002, 81, 416–426. [Google Scholar] [CrossRef]
  55. R Core Team. R: A language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 3 June 2020).
  56. Ye, W.; Li, X.; Chen, X.; Zhang, G. A Spectral index for Highlighting Forest cover from Remotely Sensed Imagery; Jackson, T.J., Chen, J.M., Gong, P., Liang, S., Eds.; SPIE Asia-Pacific Remote Sensing: Beijing, China, 2014; p. 92601L. [Google Scholar]
  57. Mišurec, J.; Kopačková, V.; Lhotáková, Z.; Campbell, P.; Albrechtová, J. Detection of Spatio-Temporal Changes of Norway Spruce Forest Stands in Ore Mountains Using Landsat Time Series and Airborne Hyperspectral Imagery. Remote Sens. 2016, 8, 92. [Google Scholar] [CrossRef] [Green Version]
  58. Kooijman, A.M.; Emmer, I.M.; Fanta, J.; Sevink, J. Natural regeneration potential of the degraded Krkonoše forests. Land Degrad. Dev. 2000, 11, 459–473. [Google Scholar] [CrossRef]
  59. Vacek, S.; Hůnová, I.; Vacek, Z.; Hejcmanová, P.; Podrázský, V.; Král, J.; Putalová, T.; Moser, W.K. Effects of air pollution and climatic factors on Norway spruce forests in the Orlické hory Mts. (Czech Republic), 1979–2014. Eur. J. For. Res. 2015, 134, 1127–1142. [Google Scholar] [CrossRef]
  60. Ponocná, T.; Chuman, T.; Rydval, M.; Urban, G.; Migaɬa, K.; Treml, V. Deviations of treeline Norway spruce radial growth from summer temperatures in East-Central Europe. Agric. For. Meteorol. 2018, 253–254, 62–70. [Google Scholar] [CrossRef]
  61. Treml, V.; Ponocná, T.; Büntgen, U. Growth trends and temperature responses of treeline Norway spruce in the Czech-Polish Sudetes Mountains. Clim. Res. 2012, 55, 91–103. [Google Scholar] [CrossRef] [Green Version]
  62. Verbesselt, J.; Hyndman, R.; Zeileis, A.; Culvenor, D. Phenological change detection while accounting for abrupt and gradual trends in satellite image time series. Remote Sens. Environ. 2010, 114, 2970–2980. [Google Scholar] [CrossRef] [Green Version]
  63. Qiu, S.; Lin, Y.; Shang, R.; Zhang, J.; Ma, L.; Zhu, Z. Making Landsat Time Series Consistent: Evaluating and Improving Landsat Analysis Ready Data. Remote Sens. 2019, 11, 51. [Google Scholar] [CrossRef] [Green Version]
  64. Masek, J.G.; Vermote, E.F.; Saleous, N.E.; Wolfe, R.; Hall, F.G.; Huemmrich, K.F.; Gao, F.; Kutler, J.; Lim, T.-K. A Landsat Surface Reflectance Dataset for North America, 1990 2000. IEEE Geosci. Remote Sens. Lett. 2006, 3, 68–72. [Google Scholar] [CrossRef]
  65. Vermote, E.; Justice, C.; Claverie, M.; Franch, B. Preliminary analysis of the performance of the Landsat 8/OLI land surface reflectance product. Remote Sens. Environ. 2016, 185, 46–56. [Google Scholar] [CrossRef]
  66. Cowles, T.R.; McNeil, B.E.; Eshleman, K.N.; Deel, L.N.; Townsend, P.A. Does the spatial arrangement of disturbance within forested watersheds affect loadings of nitrogen to stream waters? A test using Landsat and synoptic stream water data. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 80–87. [Google Scholar] [CrossRef]
  67. Zahradník, P.; Zahradníková, M. Salvage felling in the Czech Republic‘s forests during the last twenty years. Cent. Eur. For. J. 2019, 65, 12–20. [Google Scholar] [CrossRef]
  68. Tahovská, K.; Kaňa, J.; Bárta, J.; Oulehle, F.; Richter, A.; Šantrůčková, H. Microbial N immobilization is of great importance in acidified mountain spruce forest soils. Soil Biol. Biochem. 2013, 59, 58–71. [Google Scholar] [CrossRef]
  69. Altman, J.; Fibich, P.; Santruckova, H.; Dolezal, J.; Stepanek, P.; Kopacek, J.; Hunova, I.; Oulehle, F.; Tumajer, J.; Cienciala, E. Environmental factors exert strong control over the climate-growth relationships of Picea abies in Central Europe. Sci. Total Environ. 2017, 609, 506–516. [Google Scholar] [CrossRef]
Figure 1. Locations of stream catchments within the GEOMON monitoring network in the Czech Republic.
Figure 1. Locations of stream catchments within the GEOMON monitoring network in the Czech Republic.
Remotesensing 12 01944 g001
Figure 2. Field and image data processing workflow (VI – vegetation index, S – sulfur, DIN – dissolved inorganic nitrogen).
Figure 2. Field and image data processing workflow (VI – vegetation index, S – sulfur, DIN – dissolved inorganic nitrogen).
Remotesensing 12 01944 g002
Figure 3. Long-term trends observed in (a) sulfur deposition and (b) stream water total dissolved inorganic nitrogen export as normalized z-scores computed from all catchments. Bold lines show median values for all sites, gray-shaded areas show variations between catchments (10th and 90th percentiles), and dashed lines show the linear trends.
Figure 3. Long-term trends observed in (a) sulfur deposition and (b) stream water total dissolved inorganic nitrogen export as normalized z-scores computed from all catchments. Bold lines show median values for all sites, gray-shaded areas show variations between catchments (10th and 90th percentiles), and dashed lines show the linear trends.
Remotesensing 12 01944 g003
Figure 4. Pearson correlation coefficients computed between measured sulfur deposition (S-input) or stream water dissolved inorganic nitrogen export (DIN-output) and selected vegetation indices from the Landsat satellite data. Crosses indicate non-significant correlations (p > 0.05).
Figure 4. Pearson correlation coefficients computed between measured sulfur deposition (S-input) or stream water dissolved inorganic nitrogen export (DIN-output) and selected vegetation indices from the Landsat satellite data. Crosses indicate non-significant correlations (p > 0.05).
Remotesensing 12 01944 g004
Figure 5. Trends (dashed lines) in long-term standard scores for measured sulfur deposition (S-input) and stream water dissolved inorganic nitrogen export (DIN-output) during the observation period 1994–2016 and trends in vegetation indices extracted from Landsat satellite data for the observation period 1984–2016. Only 7 out of 13 stream catchments were shown here, as they well represented the geographic areas in the Czech Republic and overall variety among the catchments. Trends for all of the catchments are presented in Figure S1, available online in the Supplementary Materials.
Figure 5. Trends (dashed lines) in long-term standard scores for measured sulfur deposition (S-input) and stream water dissolved inorganic nitrogen export (DIN-output) during the observation period 1994–2016 and trends in vegetation indices extracted from Landsat satellite data for the observation period 1984–2016. Only 7 out of 13 stream catchments were shown here, as they well represented the geographic areas in the Czech Republic and overall variety among the catchments. Trends for all of the catchments are presented in Figure S1, available online in the Supplementary Materials.
Remotesensing 12 01944 g005
Figure 6. Linear relationships between changes from 1995 to 2016 in forest disturbance index (DI) and (a) total sulfur deposition (S-input), and (b) stream water dissolved inorganic nitrogen export (DIN-output). Bold lines show the linear trends, grey-shaded areas show 95% confidence level intervals. The coefficient of determination and p-values are included in the graph.
Figure 6. Linear relationships between changes from 1995 to 2016 in forest disturbance index (DI) and (a) total sulfur deposition (S-input), and (b) stream water dissolved inorganic nitrogen export (DIN-output). Bold lines show the linear trends, grey-shaded areas show 95% confidence level intervals. The coefficient of determination and p-values are included in the graph.
Remotesensing 12 01944 g006
Figure 7. Decennial changes in the Landsat-derived disturbance index computed for the Eagle Mountains protected landscape area (protected since 1969), with detailed insight into the UDL catchment.
Figure 7. Decennial changes in the Landsat-derived disturbance index computed for the Eagle Mountains protected landscape area (protected since 1969), with detailed insight into the UDL catchment.
Remotesensing 12 01944 g007
Figure 8. Decennial changes in Landsat-derived disturbance index computed for the Slavkov Forest protected landscape area (protected since 1974), with detailed insight into the PLB catchment. Note that the increased value of DI in the northwest area of the Landsat image from 2013 was due to cloud cover (there was not a single cloud-free Landsat image available for the entire area around 2013).
Figure 8. Decennial changes in Landsat-derived disturbance index computed for the Slavkov Forest protected landscape area (protected since 1974), with detailed insight into the PLB catchment. Note that the increased value of DI in the northwest area of the Landsat image from 2013 was due to cloud cover (there was not a single cloud-free Landsat image available for the entire area around 2013).
Remotesensing 12 01944 g008
Figure 9. Percentages of disturbed (DI < μ − 0.5σ), no change (DI between μ − 0.5σ and μ + 0.5σ), and growth/recovery (DI > μ + 0.5σ) pixels of (a) Eagle Mountains protected landscape area, (b) zoom to the UDL catchment, (c) Slavkov Forest protected landscape area, and (d) zoom to the PLB catchment. Statistics were computed from Figure 7 and Figure 8.
Figure 9. Percentages of disturbed (DI < μ − 0.5σ), no change (DI between μ − 0.5σ and μ + 0.5σ), and growth/recovery (DI > μ + 0.5σ) pixels of (a) Eagle Mountains protected landscape area, (b) zoom to the UDL catchment, (c) Slavkov Forest protected landscape area, and (d) zoom to the PLB catchment. Statistics were computed from Figure 7 and Figure 8.
Remotesensing 12 01944 g009
Table 1. Characteristics of GEOMON stream catchments in this study and the number of cloud-free Landsat images for the period 1984-2016 that were used to compute annual trajectories.
Table 1. Characteristics of GEOMON stream catchments in this study and the number of cloud-free Landsat images for the period 1984-2016 that were used to compute annual trajectories.
CodeCatchment NameArea
(ha)
Elevation
(m a.s.l.)
Mean Temp. (°C)Mean Precip. (mm)Bedrock TypeForest Age No. of Images
ANEAnenský potok27480–5408.0673Gneiss40–60336
CERČervík185640–8906.01212Sandstone, Claystone40–60213
JEZJezeří261475–9256.0815Gneiss60–80168
LESLesní potok70400–5008.0615Granite40–60259
LIZLiz99830–10255.5892Gneiss60–100180
LKVLoukov66470–6607.5754Granite40–80261
LYSLysina27830–9505.01005Granite40–60192
MODModrý potok2621010–15503.01727Phyllite, Mica schist100–12069
PLBPluhův bor22690–8006.0861Serpentinite60–100228
POMPolomka69510–6407.0697Gneiss60–80231
SALSalačova Lhota168560–7457.0675Gneiss60–80206
UDLU dvou louček33880–9505.01502Gneiss20–40246
UHLUhlířská187780–8705.51250Granite, Granodiorite20–40183
Table 2. Summary of vegetation indices computed from Landsat image data in this study (B – blue, G – green, R – red, and NIR – near-infrared band).
Table 2. Summary of vegetation indices computed from Landsat image data in this study (B – blue, G – green, R – red, and NIR – near-infrared band).
Vegetation IndexEquationReference
DI – Disturbance indexsee eq. 1–4Healey et al. [45]
EVI – Enhanced vegetation index2.5*[(NIR-R)/(NIR+6*R-7.5*B+1)]Liu and Huete [51]
NDVI – Normalized difference vegetation index(NIR-R)/(NIR+R)Rouse J. W. et al. [52]
PSRI – Plant senescence reflectance index(R-G)/NIRMerzlyak et al. [53]
TCARI – Transformed chlorophyll absorption reflectance index3*[(NIR-R)-0.2(NIR-G)*(NIR/R)]Haboudane et al. [54]

Share and Cite

MDPI and ACS Style

Švik, M.; Oulehle, F.; Krám, P.; Janoutová, R.; Tajovská, K.; Homolová, L. Landsat-Based Indices Reveal Consistent Recovery of Forested Stream Catchments from Acid Deposition. Remote Sens. 2020, 12, 1944. https://doi.org/10.3390/rs12121944

AMA Style

Švik M, Oulehle F, Krám P, Janoutová R, Tajovská K, Homolová L. Landsat-Based Indices Reveal Consistent Recovery of Forested Stream Catchments from Acid Deposition. Remote Sensing. 2020; 12(12):1944. https://doi.org/10.3390/rs12121944

Chicago/Turabian Style

Švik, Marian, Filip Oulehle, Pavel Krám, Růžena Janoutová, Kateřina Tajovská, and Lucie Homolová. 2020. "Landsat-Based Indices Reveal Consistent Recovery of Forested Stream Catchments from Acid Deposition" Remote Sensing 12, no. 12: 1944. https://doi.org/10.3390/rs12121944

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