Assessment of Carbon Flux and Soil Moisture in Wetlands Applying Sentinel-1 Data

The objectives of the study were to determine the spatial rate of CO2 flux (Net Ecosystem Exchange) and soil moisture in a wetland ecosystem applying Sentinel-1 IW (Interferometric Wide) data of VH (Vertical Transmit/Horizontal Receive—cross polarization) and VV (Vertical Transmit/Vertical Receive—like polarization) polarization. In-situ measurements of carbon flux, soil moisture, and LAI (Leaf Area Index) were carried out over the Biebrza Wetland in north-eastern Poland. The impact of soil moisture and LAI on backscattering coefficient (σ°) calculated from Sentinel-1 data showed that LAI dominates the influence on σ° when soil moisture is low. The models for soil moisture have been derived for wetland vegetation habitat types applying VH polarization (R2 = 0.70 to 0.76). The vegetation habitats: reeds, sedge-moss, sedges, grass-herbs, and grass were classified using combined one Landsat 8 OLI (Operational Land Imager) and three TerraSAR-X (TSX) ScanSAR VV data. The model for the assessment of Net Ecosystem Exchange (NEE) has been developed based on the assumption that soil moisture and biomass represented by LAI have an influence on it. The σ° VH and σ° VV describe soil moisture and LAI, and have been the input to the NEE model. The model, created for classified habitats, is as follows: NEE = f (σ° Sentinel-1 VH, σ° Sentinel-1 VV). Reasonably good predictions of NEE have been achieved for classified habitats (R2 = 0.51 to 0.58). The developed model has been used for mapping spatial and temporal distribution of NEE over Biebrza wetland habitat types. Eventually, emissions of CO2 to the atmosphere (NEE positive) has been noted when soil moisture (SM) and biomass were low. This study demonstrates the importance of the capability of Sentinel-1 microwave data to calculate soil moisture and estimate NEE with all-weather acquisition conditions, offering an important advantage for frequent wetlands monitoring.


Introduction
Wetlands are very sensitive ecosystems, functioning as a habitat for many unique plants and animals.Wetlands play a significant role in many ecological functions including: carbon balance, regulation of the hydrological cycle, flood control, and the improvement and maintenance of water quality [1].These areas remove atmospheric CO 2 through sequestration and subsequent peat accumulation, which can help to mitigate climate change effects.Those climate variations may also

•
Classification of the wetland vegetation habitat types based on optical and microwave satellite images.

•
Developing models for soil moisture assessment applying σ • calculated from Sentinel-1 data.
• Examination of influence of LAI and soil moisture on σ • calculated from Sentinel-1 data.

•
Developing NEE models applying σ • calculated from Sentinel-1 VH and VV data.

•
Mapping spatial distribution of NEE over a test site area.
The investigation was carried out over the Biebrza Wetland located in North-East Poland from April to October in 2013-2015.For the year 2013 in-situ data for the classification of vegetation habitats have been carried out, and for the years 2014-2015 ground measurements of soil-vegetation parameters and carbon fluxes were done.

Test Site
The study was carried out at Biebrza Wetland test site, which is located in North-East Poland (UL coordinates: 53 • 50 N, 22 • 10 E; LR coordinates: 53 • 05 N, 23 • 10 E). Figure 1 presents the localization of the test site in Poland on the orthophotomap (Figure 1a) and the Biebrza Wetland area on Landsat 8 OLI (Operational Land Imager) RGB (Red, Green, Blue spectral bands 4,3,2) image registered on 3 August 2014 (Figure 1b).The Biebrza Wetland was designated as a wetland site of international importance in the NATURA 2000 Network and since 1995 has been under the protection of the RAMSAR Convention on Wetlands.The test site contains Biebrza National Park (BNP) established in 1993 with a total area of 59.233 ha which includes 15.547 ha of forests, 18.182 ha of agricultural land, and 25.494 ha of wetlands-the most valuable habitats of the park.It is a flat area with an average altitude of about 105 m above sea level, situated close to the Belarus border.The main river is Biebrza, which flows out near the eastern border of Poland and whose recipient is the Narew River.
Remote Sens. 2016, 8, 756 3 of 22 The investigation was carried out over the Biebrza Wetland located in North-East Poland from April to October in 2013-2015.For the year 2013 in-situ data for the classification of vegetation habitats have been carried out, and for the years 2014-2015 ground measurements of soil-vegetation parameters and carbon fluxes were done.

Test Site
The study was carried out at Biebrza Wetland test site, which is located in North-East Poland (UL coordinates: 53°50′′N, 22°10′′E; LR coordinates: 53°05′′N, 23°10′′E).Figure 1 presents the localization of the test site in Poland on the orthophotomap (Figure 1a) and the Biebrza Wetland area on Landsat 8 OLI (Operational Land Imager) RGB (Red, Green, Blue spectral bands 4,3,2) image registered on 3 August 2014 (Figure 1b).The Biebrza Wetland was designated as a wetland site of international importance in the NATURA 2000 Network and since 1995 has been under the protection of the RAMSAR Convention on Wetlands.The test site contains Biebrza National Park (BNP) established in 1993 with a total area of 59.233 ha which includes 15.547 ha of forests, 18.182 ha of agricultural land, and 25.494 ha of wetlands-the most valuable habitats of the park.It is a flat area with an average altitude of about 105 m above sea level, situated close to the Belarus border.The main river is Biebrza, which flows out near the eastern border of Poland and whose recipient is the Narew River.
The soils are mainly hydrogenic, with a dominant majority of peat soils in various stages of mouldering.The peatland complex covers about 90.000 ha and is the largest and the least disturbed of its kind in Central Europe.Wetlands floods occur annually in the spring and, besides precipitation, are the main supply of moisture into the peat soil.The site is characterized by large homogeneous fields of the following non-forest habitats: reeds, sedges, sedge-moss, grass, grass-herbs.
The weather is the coldest in Poland, with long winters (up to 120 days below zero) and a short growing season.The length of the summer ranges from 77 to 85 days, with daytime temperatures averaging 15 °C.The level of precipitation is one of the lowest in Poland-around 510 mm.The dry year was 2015, and normal years were 2013 and 2014.The soils are mainly hydrogenic, with a dominant majority of peat soils in various stages of mouldering.The peatland complex covers about 90.000 ha and is the largest and the least disturbed of its kind in Central Europe.Wetlands floods occur annually in the spring and, besides precipitation, are the main supply of moisture into the peat soil.The site is characterized by large homogeneous fields of the following non-forest habitats: reeds, sedges, sedge-moss, grass, grass-herbs.
The weather is the coldest in Poland, with long winters (up to 120 days below zero) and a short growing season.The length of the summer ranges from 77 to 85 days, with daytime temperatures averaging 15

Data Sets
The datasets used for the study were: in-situ measurements, TerraSAR-X (TSX) and Landsat 8 OLI images used for vegetation habitat types classification, and Sentinel-1 (S1) data used for examining the relationship of backscatter to soil moisture, LAI, and NEE.Table 1 presents the list of microwave satellite data taken into the research: TerraSAR-X (X-band, wavelength 31 mm, frequency 9.65 GHz, ENL-1 (Effective Number of Looks), and Sentinel-1 (C-band, wavelength 5.55 cm, frequency 5.405 GHz, ENL-5).TerraSAR-X images were registered in ScanSAR imaging mode and single polarization-VV.The TSX acquired product type was Geocoded Ellipsoid Corrected (GEC) in UTM projection.Sentinel-1 images were registered in Interferometric Wide (IW) mode, incidence angle 35 • and dual polarization-VV and VH.The S1 acquired product type was Ground Range Detected (GRD) in WGS84 ellipsoid.The optical data-Landsat 8 OLI (Operational Land Imager) image registered on 19 May 2013 was applied in combination with TSX images in the classification procedure.This image was registered in nine spectral bands from 0.43 to 1.38 µm with a spatial resolution of 30 m for Bands 1 to 7 and 9, and 15 m for Band 8-panchromatic.
The in-situ measurements have been done from April to October of the years 2013-2015 at nearly 25 points distributed in various sites of the Biebrza Wetland in order to assess the distribution of soil moisture, vegetation parameters, and CO 2 flux over various types of vegetation habitats.During the field campaigns, each of the sites was treated three times in three different localizations around a measurement point for the calculation of the mean value.The dates of in-situ measurements corresponded to the dates of Sentinel-1 overpasses.The positions of the measurement plots were determined using GPS.Instantaneous flow of CO 2 between the ground surface and the atmosphere was measured using a static chamber method [40].The concentration of CO 2 (ppm) and air temperature ( • C) inside the plastic transparent chamber (40 × 40 × 35 cm) was measured by a portable non-dispersive infrared sensor (SenseAir ® , SenseAir SE) in one minute intervals.With this set-up, carbon flux (NEE) was first measured under light conditions immediately followed by a measurement of ecosystem respiration (RESP) after darkening of the chamber.The ground measurements also involved simultaneously measured significant vegetation and soil parameters such as: soil moisture in the layer of 0-15 cm (there is strong relationship between 0-5 cm and 0-15 cm for the peat soil) using the TDR method (TRIME-FM, IMKO Micromodultechnik GmbH, Ettlingen, Germany), LAI (LAI 2000, LiCor Inc., Lincoln, NE, USA), wet and dry biomass (in the laboratory from samples harvested inside the whole surface covered by the chamber following the flux measurement).All in-situ and satellite data have been stored in ACCESS and EXCEL databases.

Methods
The positions of the measurement plots, has been determined by GPS and were further used for preparing the measurement points layer which was needed for extracting satellite data.
Sentinel-1 IW and TSX ScanSAR microwave images have been processed using ESA Sentinel-1 Toolbox (S1TBX, developed for ESA by Array Systems Computing-Toronto, ON, Canada, in partnership with DLR-National Aeronautics and Space Research Centre, Germany, Brockmann Consult-Hamburg, Germany, and OceanDataLab-Locmaria Plouzane, France) in SNAP (Sentinel Application Platform, and is a common architecture for all Sentinel Toolboxes being jointly developed by Brockmann Consult and Array Systems Computing) to derive a backscatter coefficient (σ • ) (dB).The processing included: speckle filtering applying Gamma-MAP filter 7 × 7, radiometric calibration and data conversion to backscattering coefficient.Then, scenes were geometrically registered to local projection PUWG1992 and S1 σ • values which corresponded to measurement points were extracted (5 × 5 pixels-within this area three in-situ measurements have been carried out to calculate mean value) using ERDAS software (Hexagon Geospatial/Intergraph ® , Norcross, GA, USA).Landsat 8 OLI images have been processed using ERDAS software.Then, the scene acquired in UTM map projection and WGS84 ellipsoid was geometrically registered to local projection PUWG1992.Atmospheric Correction of TOA reflectance (Top of Atmosphere) was not applied.It was not necessary for the classification process, which uses only one cloud-free image, as the training data and the data to be classified are in the same relative scale.[41].
The wetland vegetation habitats have been classified using combined one Landsat 8 OLI image and three TerraSAR-X ScanSAR VV images, all registered in the year 2013.These images have been chosen for classification, because of their high spatial resolution and the period of the growing season (before the first cut).Moreover, X-wave is dedicated to vegetation studies.Also, there was no possibility to acquire adequate optical and Sentinel-1 images for the study period.Landsat 8 OLI spectral bands have been resampled to pixel size of TSX images i.e., 18 m.The pixel-based supervised maximum likelihood classification method was applied using the ERDAS software.This type of classification has been used due to availability of adequate ground truth data collected in 2013 during field campaigns.The representative, homogeneous ground reference data were used for training and testing purposes.There were no changes in wetland vegetation type and cover from 2013 to 2015.
The collected time series of CO 2 flux concentration were validated in terms of temporal linearity.Fluxes with a correlation coefficient of the slope R 2 > 0.95 were accepted for further analysis.The CO 2 flux rate (NEE and RESP) was then calculated using the following equation described by Mosier and Mack in [42]: where: E-CO 2 flux rate (mg m −2 •h −1 ); V-chamber volume (0.056 m 3 ); A-collar area (m 2 ); ∆C ∆t -average rate of CO 2 concentration change with time; T-air temperature in chamber ( • C).NEE and RESP values were then predicted to µmol CO 2 m −2 •s −1 assuming that 1 µmol CO 2 weights 44 µg.
In-situ measured soil-vegetation parameters and satellite microwave data were analysed for wetland vegetation habitats.Simple and multiple regression analyses were used to create the best models-first to examine the relationship between soil moisture and the backscattering coefficient of Sentinel-1, as well as the relationship between LAI and σ • of S1.Data were checked for normal distribution and significance prior to all analysis.Statistical analyses for modeling NEE applying Sentinel-1 VV and VH data were done in STATISTICA and EXCEL software.

Classification of Wetland Vegetation Habitats
To estimate the relationship between microwave satellite data and soil-vegetation parameters precisely, there is a need to perform a classification of wetland habitats, as the relationships vary by the vegetation type.For the classification procedure, the following Landsat 8 OLI channels acquired on 19 May 2013 were chosen: 2 (0.45-0.51), 3 (0.53-0.59 µm), 4 (0.64-0.67 µm), 5 (0.85-0.88 µm), 6 (1.36-1.38 µm), and 7 (1.57-1.65 µm) along with three TerraSAR-X ScanSAR VV images acquired on 21 May 2013, 7 June 2013 and 23 June 2013.The validation and accuracy of the classification were performed using ground truth observations.As the results of wetland habitats classification, water and seven vegetation classes were distinguished: two forest habitats (coniferous and deciduous) and five representing most dominant non-forest wetland habitats: reeds, sedge-moss, sedges, grass-herbs, and grass.Table 2 presents the error matrix for classified data in %.The overall accuracy equals to 87.27% with 0.8333 of overall Kappa Statistics.Table 3 presents accuracy totals for each of the classes.  2 presents the error matrix for classified data in %.The overall accuracy equals to 87.27% with 0.8333 of overall Kappa Statistics.Table 3 presents accuracy totals for each of the classes.Figure 2 presents a map of classified wetland habitat types.

Soil Moisture
In-situ soil moisture (SM) values (ranged from 10% to 100% volumetric.) from the test site have been related to the backscattering coefficient (σ • ) calculated from Sentinel-1 IW VV and VH data, applying various regression models.The best results have been found applying VH polarization and a linear model performed for particular vegetation habitat types separately (Figure 2).These results are presented in Table 4.The relationships between SM and σ • at the 95% confidence level is statistically significant (p < 0.05).The highest correlation belongs to σ • calculated from Sentinel-1 IW VH data for sedge-moss and sedges, the most typical wetland habitats.Slightly better results have been found applying an exponential model.Derived equations (Table 4) were used to map soil moisture distribution over the Biebrza Wetland.Figures 3 and 4 present two maps composed using Sentinel-1 VH data acquired on 11 June 2015 and 10 July 2015, respectively, as an example.It can be noted that much lower SM value in July compared to June was caused due to severe drought, which started in May.Figures 5 and 6 present relationships between soil moisture measured in-situ and calculated from Sentinel-1 IW VV and VH data (respectively) applying a statistically derived algorithm (Table 5).The analyses were performed for all wetland habitats (for a representative sample).The coefficient of determination (R 2 ) between SM and σ • VH is comparable to the R 2 obtained between SM and σ • VV.The R 2 are slightly lower from that obtained for particular habitat classes separately (Table 4).Table 5 presents results of these statistical analyses which could be applied for mapping soil moisture without the need for classification of the vegetation habitat types.There is a relationship between leaf area index (LAI) and wet biomass (Bw) (both measured insitu) independent of vegetation habitat types.Figure 7 presents the best fit which was a power model    There is a relationship between leaf area index (LAI) and wet biomass (Bw) (both measured in-  There is a relationship between leaf area index (LAI) and wet biomass (Bw) (both measured in-situ) independent of vegetation habitat types.Figure 7 presents the best fit which was a power model at a 95% confidence level.The data also include measurements taken from the previous studies.A large increase in R 2 was achieved when particular habitat types (Figure 2) were included in the analysis separately.Table 6 presents results of statistical analyses performed for each of the classified non-forest wetland habitats.The highest R 2 has been found for sedges, the lowest for reeds, very often water logged.at a 95% confidence level.The data also include measurements taken from the previous studies.A large increase in R 2 was achieved when particular habitat types (Figure 2) were included in the analysis separately.Table 6 presents results of statistical analyses performed for each of the classified non-forest wetland habitats.The highest R 2 has been found for sedges, the lowest for reeds, very often water logged.LAI can be estimated by applying optical or microwave satellite data [31,43,44].In the present study the analysis between LAI and the backscattering coefficient calculated from Sentinel-1A data has been investigated by means of regression models.Vegetation covering the soil plays a certain role on the registered microwave signal.In the present study, it has been examined how vegetation (represented by LAI) affects the S1 microwave signal.At fixed LAI values, soil moisture is in simple linear correlation with σ° VH, but otherwise the dependence of LAI on backscatter has a polynomial character.The linear multiply regression model of σ° VH includes these two components as the independent variables.The resulted regression equation was as follows: σ° VH = −22.45+ 0.05 × SM + 0.55 × LAI − 0.086 × LAI 2  (2) where R 2 = 0.75, p < 0.05. Figure 8    LAI can be estimated by applying optical or microwave satellite data [31,43,44].In the present study the analysis between LAI and the backscattering coefficient calculated from Sentinel-1A data has been investigated by means of regression models.Vegetation covering the soil plays a certain role on the registered microwave signal.In the present study, it has been examined how vegetation (represented by LAI) affects the S1 microwave signal.At fixed LAI values, soil moisture is in simple linear correlation with σ • VH, but otherwise the dependence of LAI on backscatter has a polynomial character.The linear multiply regression model of σ • VH includes these two components as the independent variables.The resulted regression equation was as follows: where R 2 = 0.75, p < 0.05. Figure 8   The partial correlations for the particular components have been examined.For the model presented by Equation (2) the R was equal to 0.85 for SM, and for LAI-R was equal to 0.22 (expression of LAI in Equation ( 2) was treated as one component).It means that the influence of SM for the explanation of the variance of σ° VH, after establishing the influence of LAI, is several times stronger than the influence of LAI when SM is constant.Figure 9 illustrates this case in 3D space.A full range of soil moisture under fixed LAI caused the variation of σ° VH up to 5 dB.The case of LAI under fixed SM caused a variation of σ° VH of only 1.2 dB.It can be concluded that the influence of LAI on σ° is slightly over the statistical error (0.89 dB).The partial correlations for the particular components have been examined.For the model presented by Equation (2) the R was equal to 0.85 for SM, and for LAI-R was equal to 0.22 (expression of LAI in Equation ( 2) was treated as one component).It means that the influence of SM for the explanation of the variance of σ • VH, after establishing the influence of LAI, is several times stronger than the influence of LAI when SM is constant.Figure 9 illustrates this case in 3D space.A full range of soil moisture under fixed LAI caused the variation of σ • VH up to 5 dB.The case of LAI under fixed SM caused a variation of σ • VH of only 1.2 dB.It can be concluded that the influence of LAI on σ • is slightly over the statistical error (0.89 dB).
The influence of the two components, LAI and SM, on σ • VH values when the soil moisture is low was examined.The model described by Equation (3) presents the relationship of LAI and SM for the conditions when soil moisture was lower than 40%.σ • VH = −11.6 + 0.046 × SM + 1.12 × LAI − 0.24 × LAI 2  (3) Figure 10 presents the relationship between observed and modelled σ • VH (dB) for SM <40, where R 2 = 0.6, p < 0.05, number of observations-52, RMSE = 0.73 dB, MBE = 0.1 dB.When LAI is constant, the partial correlation for SM explaining σ • VH variation gave R equal to 0.61, and R was equal to 0.48 for LAI explaining σ • VH variation when SM is constant.
for the explanation of the variance of σ° VH, after establishing the influence of LAI, is several times stronger than the influence of LAI when SM is constant.Figure 9 illustrates this case in 3D space.A full range of soil moisture under fixed LAI caused the variation of σ° VH up to 5 dB.The case of LAI under fixed SM caused a variation of σ° VH of only 1.2 dB.It can be concluded that the influence of LAI on σ° is slightly over the statistical error (0.89 dB).The influence of the two components, LAI and SM, on σ° VH values when the soil moisture is low was examined.The model described by Equation (3) presents the relationship of LAI and SM for the conditions when soil moisture was lower than 40%.
σ° VH = −11.6 + 0.046 × SM + 1.12 × LAI − 0.24 × LAI 2  (3) Figure 10 presents the relationship between observed and modelled σ° VH (dB) for SM <40, where R 2 = 0.6, p < 0.05, number of observations-52, RMSE = 0.73 dB, MBE = 0.1 dB.When LAI is constant, the partial correlation for SM explaining σ° VH variation gave R equal to 0.61, and R was equal to 0.48 for LAI explaining σ° VH variation when SM is constant.With fixed low soil moisture values, the range of backscatter variance depended on LAI and could increase up to 3.2 dB.When LAI is on a fixed level, the variance of backscatter caused by SM gets up to 2.4 dB, (Figure 11).Contribution of LAI to backscatter VH was higher when SM values were low.High difference of σ • VH between two points covered by vegetation does not mean that the soil moisture is really different.In the model (3) σ • VH >−19 dB represents a high value of LAI for SM of about 30%-40%.The same values of σ • VH, i.e., >−19 dB in the model ( 2) occurred for high SM (>70%) independently of LAI (Figures 9 and 11).

NEE Modeling
NEE values measured in-situ (Section 2.2) were related to various in-situ and satellite data to derive prediction models.The analyses were performed for particular vegetation habitat types (Figure 2) separately.The best models were found by applying Sentinel-1 to both VH and VV polarization (Table 7).These were multiple linear regression models.Since the p-value was less than 0.05, there is a statistically significant relationship between the variables at the 95.0% confidence level.The R 2 indicates that the model explains no more than 51%-58% of the variability in NEE.Thus, these preliminary results have to be validated and the model has to be corrected to improve statistical correlation.Equations from Table 7 have been applied for the NEE mapping over the Biebrza Wetland.Figures 12 and 13

NEE Modeling
NEE values measured in-situ (Section 2.2) were related to various in-situ and satellite data to derive prediction models.The analyses were performed for particular vegetation habitat types (Figure 2) separately.The best models were found by applying Sentinel-1 to both VH and VV polarization (Table 7).These were multiple linear regression models.Since the p-value was less than 0.05, there is a statistically significant relationship between the variables at the 95.0% confidence level.The R 2 indicates that the model explains no more than 51%-58% of the variability in NEE.Thus, these preliminary results have to be validated and the model has to be corrected to improve statistical correlation.Equations from Table 7 have been applied for the NEE mapping over the Biebrza Wetland.

Classification of Wetland Vegetation Habitat Types
Classification of wetland vegetation habitat types using remote sensing data has been undertaken many times.The methods, most often used, concerned the application of optical data [44][45][46][47][48]. Dabrowska-Zielinska et al. [44] derived six wetland habitats from the classification of MERIS (Medium Spectral Resolution Imaging Spectrometer) images acquired on 21 April 2003: anthropogenic, pastures, sedge-mursh, rush-sedge, reeds, and scrubs.Wright and Gallant [48] distinguished five classes: aquatic bed (vegetation on or below the water surface), emergent (herbaceous hydrophytes), forested, scrub-shrub (height < 6 m), and unconsolidated shore, applying Landsat TM imagery.Recently, the application of microwave data for the classification of vegetation habitats has been expanded by Evans et al. [49,50].Santoro et al. [51] who studied ALOS PALSAR (Advanced Land Observing Satellite Phased Array type L-band Synthetic Aperture Radar) data found that better classification results were obtained from HV than from HH polarization.Dabrowska-Zielinska et al. [43] and Conforth et al. [52] also presented wetland classification using ALOS PALSAR data acquired in HV polarization.Flores De Santiago et al. [53] found that a combination of HH and HV polarization modes from ALOS PALSAR data was better at separating forested wetlands using an object-based classification approach.Bwangoy et al. [54] applied optical (Landsat TM (Thematic Mapper) and ETM+ (Enhanced Thematic Mapper Plus) and JERS-1 (Japanese Earth Resources Satellite-1) data to classify the wetland and non-wetland classes of the Congo Basin.Ozesmi and Bauer [55] reported that the combination of radar and optical data will improve wetland classification.
There is a lack of a detailed classification of European wetland vegetation.In this study we decided to apply microwave and optical data for the classification of wetland vegetation habitats.Compared to the classification results undertaken using ALOS PALSAR HV [43], in the present study it was possible to distinguish wetland habitat types more precisely.The obtained classification results of wetland habitat types (two forest and five non-forest habitats) presented in Figure 2 were not found in the literature before and are promising in wetland changes monitoring.

Soil Moisture
Soil moisture retrieval methods from single polarisation, active, spaceborne, microwave observations were described by Barrett et al. in [56].Baghdadi et al. [57] applied ENVISAT ASAR data for soil moisture retrieval over bare soils.The research on application of radar data for soil moisture assessment has been published for inland wetlands at High Latitudes applying C-band SAR (Synthetic Aperture Radar) data-Reschke et al. in [58], Bartsch et al. in [59], and Kasischke et al. in [60].They stated that C-band backscatter may be enhanced by inundation or high levels of soil moisture, what has been presented in our study.Paloscia et al. applied ENVISAT ASAR [35] and Sentinel-1 [61] images for soil moisture assessment over vegetated areas.Kornelsen and Coulibaly [62] conducted a comprehensive literature review of SAR soil moisture retrieval methodology, its limitations and potential.Hornáček et al. [63] presented the potential for surface soil moisture retrieval via the Change Detection method using Sentinel-1 data.Dorigo et al. [64] presented the evaluation of CCI (Climate Change Initiative) soil moisture product using ground based observations.Soil moisture for the Biebrza Wetland has been investigated by authors in a previous study using microwave data acquired by ERS-1/2 SAR VV, ENVISAT ASAR IS2 HH, and ALOS PALSAR HH [25,43,44,65,66].In the present study Sentinel-1 IW VV and VH data have been applied for the assessment of soil moisture (there was no S1 acquisition in HH polarisation).The best results have been found applying VH polarization.In general, the results are consistent with the literature, however most of the research presented in the literature was done for other than wetlands areas.The models developed by the authors will be validated and corrected using S1 and in-situ data acquired in the year 2016.
For soil moisture change detection, the image differencing method was used.Differencing involved the subtraction of backscattering coefficient values between two different date images.The assessment of LAI using ENVISAT MERIS and ASAR data was presented by Dabrowska-Zielinska et al. in [44,62,63].This was done for the wetland habitats distinguished from the classification of MERIS image.LAI prediction, applying ALOS PALSAR HV data, was presented by Dabrowska-Zielinska et al. in [43,65,66].The obtained results showed that with the increase of LAI there was an increase in the backscattering coefficient, which indicates the low signal attenuation by vegetation in L-band, contrary to results obtained in C-and X-band.For LAI assessment, the ASAR mode with higher incidence angles gives better results, which have been examined by authors in other studies [67], applying VV IS6 data for the agriculture area.In the present study, the influence of LAI and SM on σ° calculated from Sentinel-1 data has been examined.Two linear multiply regression models of σ° VH that include these two components as the independent variables have been developed.Comparing the results of model (2) and model (3), it has to be noted that the  The assessment of LAI using ENVISAT MERIS and ASAR data was presented by Dabrowska-Zielinska et al. in [44,62,63].This was done for the wetland habitats distinguished from the classification of MERIS image.LAI prediction, applying ALOS PALSAR HV data, was presented by Dabrowska-Zielinska et al. in [43,65,66].The obtained results showed that with the increase of LAI there was an increase in the backscattering coefficient, which indicates the low signal attenuation by vegetation in L-band, contrary to results obtained in C-and X-band.For LAI assessment, the ASAR mode with higher incidence angles gives better results, which have been examined by authors in other studies [67], applying VV IS6 data for the agriculture area.In the present study, the influence of LAI and SM on σ • calculated from Sentinel-1 data has been examined.Two linear multiply regression models of σ • VH that include these two components as the independent variables have been developed.Comparing the results of model (2) and model (3), it has to be noted that the influence of LAI on σ • VH value is much stronger when soil moisture values are low (<40%).The obtained results are promising and will be validated in the future work using new satellite and in-situ data.

NEE Modeling
The studies on carbon exchange in grassland areas applying in-situ and satellite MODIS (Moderate-Resolution Imaging Spectroradiometer) and TerraSAR-X StripMap VV data have been presented by authors in [68].The results showed significant correlations between NEE and Normalised Difference Vegetation Index -NDVI (R 2 = 0.56), and between surface temperature (Ts) and RESP (R 2 = 0.54).Also, a significant relationship has been found between NEE and σ • calculated from TSX data (R 2 = 0.53)-the higher CO 2 absorption corresponded to higher NDVI and σ • values.The proposed model for Net Ecosystem Exchange is based on a simple relationship with Sentinel-1 VV and VH data, which are actually delivered with high temporal resolution.This gives the possibility of frequent monitoring of the wetlands ecosystem, especially important under high soil-vegetation dynamics.However, these preliminary results have to be validated and corrected to improve statistical correlation.The carbon flux depends mainly on vegetation and humidity conditions.The study explained why the microwave data of Sentinel-1 VH and VV give promising results to the NEE estimates.Changes in soil moisture and vegetation cover (cuts or grazing) of the wetlands vegetation causes the carbon balance disturbance.Since wetland ecosystems are very important in carbon cycling, they have to be monitored regularly, which is possible by applying Sentinel-1 data.

Conclusions
At the Biebrza wetland study site it has been demonstrated that: (1) Wetland vegetation habitats have been classified using a combination of one optical (i.e., Landsat 8 OLI) and three microwave (i.e., TerraSar-X VV) images.The remote sensing based classification distinguished several wetlands' non-forest classes that have not been noticed before, which is novel.(2) Soil moisture could be assessed using Sentinel-1 data acquired in VH polarisation.However, these preliminary results have to be validated and models corrected using new acquisitions.Also, if available in the future, S1 HH polarisation will be included for SM modelling.
(3) There is the influence of soil moisture and vegetation biomass represented by LAI on σ • VH, however soil moisture impact on backscatter is stronger.(4) Comparing the results of model (2) and model (3), it has to be noted that the influence of LAI on σ • VH value is much stronger when soil moisture is low.(5) NEE in-situ measurements had positive values during dry soil conditions and low biomass means that CO 2 was released into the atmosphere.Also, when dry vegetation was over the water table in the flooded area, NEE values were positive.For the wet soil moisture conditions and high biomass, CO 2 absorption dominated (NEE negative values).( 6) NEE could be assessed by applying combined Sentinel-1 VV and VH data.However, to obtain better accuracy, statistical correlation should be improved.These preliminary results will be validated and the model corrected by applying new satellite acquisitions and in-situ data.
The research has been done for wetlands where moisture of the peat soil cannot be compared to the moisture of mineral soils in agriculture areas.The peat soils were constantly covered by vegetation in the vegetation growth cycle.That's why the study on the application of radar data for moisture variation of peatland is difficult and in a sense unique.The results of the study could be useful for the CCI soil moisture product presented by Dorigo et al. [64].

Figure 1 .
Figure 1.Localisation of the test site in Poland on the orthophotomap (a); and the Biebrza Wetland area on the Landsat 8 OLI RGB (4,3,2) image registered on 3 August 2014 (b).

Figure 1 .
Figure 1.Localisation of the test site in Poland on the orthophotomap (a); and the Biebrza Wetland area on the Landsat 8 OLI RGB (4,3,2) image registered on 3 August 2014 (b).

Figure 2 .
Figure 2. Map of Biebrza Wetland vegetation habitat types.Figure 2. Map of Biebrza Wetland vegetation habitat types.

Figure 2 .
Figure 2. Map of Biebrza Wetland vegetation habitat types.Figure 2. Map of Biebrza Wetland vegetation habitat types.

Figure 3 .
Figure 3. Map of soil moisture based on Sentinel-1 VH data acquired on June 2015.

Figure 4 .
Figure 4. Map of soil moisture based on Sentinel-1 VH data acquired on July 2015.

Figure 4 .
Figure 4. Map of soil moisture based on Sentinel-1 VH data acquired on July 2015.

Figure 4 .
Figure 4. Map of soil moisture based on Sentinel-1 VH data acquired on 10 July 2015.

Figure 5 .
Figure 5. Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VV data through statistically derived algorithm.

Figure 6 .
Figure 6.Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VH data through statistically derived algorithm.

Figure 5 .
Figure 5. Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VV data through statistically derived algorithm.

Figure 5 .
Figure 5. Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VV data through statistically derived algorithm.

Figure 6 .
Figure 6.Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VH data through statistically derived algorithm.

Figure 6 .
Figure 6.Relationship between observed soil moisture (SM) and calculated from Sentinel-1 VH data through statistically derived algorithm.

Figure 7 .
Figure 7. Relationship between LAI and Bw (wet biomass) for all considered habitats.

Figure 7 .
Figure 7. Relationship between LAI and Bw (wet biomass) for all considered habitats.

Figure 10 .
Figure 10.Relationship between σ° VH observed and predicted by applying Equation (3).With fixed low soil moisture values, the range of backscatter variance depended on LAI and could increase up to 3.2 dB.When LAI is on a fixed level, the variance of backscatter caused by SM gets up to 2.4 dB, (Figure 11).Contribution of LAI to backscatter VH was higher when SM values were low.High difference of σ° VH between two points covered by vegetation does not mean that the soil moisture is really different.In the model (3) σ° VH >−19 dB represents a high value of LAI for SM of about 30%-40%.The same values of σ° VH, i.e., >−19 dB in the model (2) occurred for high SM (>70%) independently of LAI (Figures 9 and 11).
present maps of NEE distribution over classified wetland habitats for two dates: 11 June 2015 and 10 July 2015.The NEE values varied from 4 to −10 μmol CO2 m −2 •s −1 .The prevailing NEE positive values (CO2 emission to atmosphere, mean from 1 to 4 μmol CO2 m −2 •s −1 ) at the date 11 June 2015 belong to grass, grass-herbs, and sedges habitats, the biomass values of which were low at that time (harvested or grazed).For the date 10 July 2015 CO2 absorption dominated (NEE negative values, mean from −1 to −10 μmol CO2 m −2 •s −1 ) in the most moist habitats, i.e., reeds and sedge-moss, while emission was noted for the lower soil moisture values (NEE values were positive, mean: 0-2 μmol CO2 m −2 •s −1 ).
Figures 12 and 13 present maps of NEE distribution over classified wetland habitats for two dates: 11 June 2015 and 10 July 2015.The NEE values varied from 4 to −10 µmol CO 2 m −2 •s −1 .The prevailing NEE positive values (CO 2 emission to atmosphere, mean from 1 to 4 µmol CO 2 m −2 •s −1 ) at the date 11 June 2015 belong to grass, grass-herbs, and sedges habitats, the biomass values of which were low at that time (harvested or grazed).For the date 10 July 2015 CO 2 absorption dominated (NEE negative values, mean from −1 to −10 µmol CO 2 m −2 •s −1 ) in the most moist habitats, i.e., reeds and sedge-moss, while emission was noted for the lower soil moisture values (NEE values were positive, mean: 0-2 µmol CO 2 m −2 •s −1 ).

Figure 12 .
Figure 12.Map of NEE distribution over Biebrza wetland on 11 June 2015.

Figure 14 22 Figure 14
Figure 14 presents a soil moisture difference map created by applying data from the dates: 10 July 2015 and 11 June 2015.The advantage of this technique is that, in cases where vegetation surface roughness doesn't change in time, the difference in backscatter between two dates can be related solely to a change in the soil moisture content.Soil moisture and vegetation represented by LAI were observed as having influence on the value of σ • calculated from Sentinel-1 VH data.Soil moisture is in simple correlation with σ • , but the dependence of LAI has a polynomial character.The dependence of LAI on backscatter is associated with low or high soil moisture.That's why the research on the influence of soil moisture and vegetation cover should be continued under different soil-vegetation conditions.Remote Sens. 2016, 8, 756 16 of 22 Figure 14 presents a soil moisture difference map created by applying data from the dates: 10 July 2015 and 11 June 2015.The advantage of this technique is that, in cases where vegetation surface roughness doesn't change in time, the difference in backscatter between two dates can be related solely to a change in the soil moisture content.Soil moisture and vegetation represented by LAI were observed as having influence on the value of σ° calculated from Sentinel-1 VH data.Soil moisture is in simple correlation with σ°, but the dependence of LAI has a polynomial character.The dependence of LAI on backscatter is associated with low or high soil moisture.That's why the research on the influence of soil moisture and vegetation cover should be continued under different soil-vegetation conditions.

Figure 14 .
Figure 14.Soil moisture difference map created by applying dates: 10 July 2015 and 11 June 2015.

Table 1 .
Microwave data acquisitions over the Biebrza Wetland test site.

Table 2 .
Error matrix for classified data in %.

Table 3 .
Accuracy totals for each class.

Table 4 .
Results of relationships at a 95% confidence level between measured soil moisture (SM) and σ • calculated from Sentinel-1A VH data.

Table 5 .
Results of the relationships at the 95% confidence level between measured SM and calculated SM from Sentinel-1 VV and VH data.Influence of LAI and Soil Moisture on σ • Calculated from Sentinel-1 Data.

Table 6 .
Results of the relationships at the 95% confidence level between measured in-situ LAI and Bw.

Table 6 .
Results of the relationships at the 95% confidence level between measured in-situ LAI and Bw.

Habitat R 2 St. Error of Est. Number of Data p-Value Equation
presents the scatterplot of this model.The root mean square error (RMSE) and the mean bias error (MBE) were used as the model performance indicators for the training and validation datasets.The RMSE was equal to 0.89 dB and MBE was equal to 0.002 dB for the training dataset (49 points).RMSE was equal to 1.8 dB and MBE was equal to −0.37 dB for the test containing 63 observations.

Table . 7
. NEE (Net Ecosystem Exchange) models at 95% confidence level derived by applying in-situ and Sentinel-1 satellite data.

Table 7 .
NEE (Net Ecosystem Exchange) models at 95% confidence level derived by applying in-situ and Sentinel-1 satellite data.
• S1VH − 2.72 × σ • S1VV TerraSAR-X (TSX) a radar Earth observation satellite, a joint venture being carried out under a public-private-partnership between the German Aerospace Center (DLR) and