Vegetation Change and Its Response to Climate Change between 2000 and 2016 in Marshes of the Songnen Plain, Northeast China

: Songnen Plain is a representative semi-arid marshland in China. The Songnen Plain marshes have undergone obvious loss during the past decades. In order to protect and restore wetland vegetation, it is urgent to investigate the vegetation change and its response to climate change in the Songnen Plain marshes. Based on the normalized di ﬀ erence vegetation index (NDVI) and climate data, we investigated the spatiotemporal change of vegetation and its relationship with temperature and precipitation in the Songnen Plain marshes. During 2000–2016, the growing season mean NDVI of the Songnen Plain marshes signiﬁcantly ( p < 0.01) increased at a rate of 0.06 / decade. For the climate change e ﬀ ects on vegetation, the growing season precipitation had a signiﬁcant positive e ﬀ ect on the growing season NDVI of marshes. In addition, this study ﬁrst found asymmetric e ﬀ ects of daytime maximum temperature (T max ) and nighttime minimum temperature (T min ) on NDVI of the Songnen Plain marshes: The growing season NDVI correlated negatively with T max but positively with T min . Considering the global asymmetric warming of T max and T min , more attention should be paid to these asymmetric e ﬀ ects of T max and T min on the vegetation of marshes.


Introduction
Marshes are an important type of wetland ecosystem worldwide, and recent research suggests that marshes may have an important role in the global carbon cycle [1][2][3]. Due to extensive distribution, marshes may affect the global climate through biogeochemical (carbon uptake and releases) and biophysical (physical properties of the land surface) processes [2]. As a key component of wetland, vegetation serves as a natural interface linking together soil cycle, water cycle, and atmosphere cycle tightly [4,5]. Marsh vegetation plays an important role in maintaining marsh functions (e.g., wildlife habitat, biodiversity protection, water conservation, carbon sequestration) [6]. For example, marsh vegetation species such as Triarrhena spp. and Phragmites spp. can provide habitat and food for wintering migratory birds [7]. The structure and composition of marsh vegetation can have significant effects on biodiversity and carbon burial of the marsh ecosystem under a rapidly changing climate [8,9]. Due to the effects of human interference and climate change, global marshes are facing degradation and loss [2]. Under the background of global climate change, understanding the response of marsh vegetation to climate change can contribute to accurately predict future vegetation changes and thus provide scientific instructions for wetland conservation and restoration [10,11]. Because climate change has significant effects on vegetation in the world, exploring the impacts of climate change on vegetation has become a hot topic of current global change research [12,13]. It is reported that climate change, including changes in temperature and precipitation regimes, can have significant effects on wetland vegetation on global and regional scales [14]. For example, previous study found that the distribution of vegetation in tidal freshwater wetlands in the Patuxent River estuary of Maryland is mainly determined by flooding related to increasing precipitation [15]. In some cold mountain regions, however, the temperature is thought to have a greater impact on the growth of wetland vegetation than precipitation [16]. In the southeastern United States, climate change is reported to result in an expansion of black mangroves at the expense of some salt marsh vegetation [17]. In temperate climate regions, warming temperatures can lead to mangrove forest expansion and salt marsh ecosystems expense [18,19]. By contrast, some studies found a positive relationship between precipitation and vegetation abundance in tidal saline wetlands [20,21]. Studies focused on the mangrove-marsh ecotone in North America and the mangrove transition zones in New Zealand and eastern Australia have found that temperature regimes can have a large impact on tidal saline vegetation aboveground biomass and canopy height [17,[22][23][24]. Although many studies have analyzed the responses of vegetation to climate change in different regions of the world, few works focused on the responses of marsh vegetation to climate change [25,26]. Compared with other ecosystems, the marsh ecosystem has distinctive environmental conditions, which may result in distinct responses of vegetation to climate change [27]. Under the background of the global climate change, investigating the effects of climate change on marsh vegetation is crucial to further understand the relationships between climate change and vegetation on a global level.
Songnen Plain of China is a representative semi-arid region, which has a large area of marshes [28]. The climate in Songnen Plain is semi-arid, and vegetation in this region is very sensitive to climate change [29]. Many studies have investigated the responses of vegetation to climate change in the Songnen Plain, but most of them focused on the vegetation of the grassland ecosystem [30,31].  analyzed the impacts of temperature and precipitation on the vegetation of marshes in Northeast China, but they did not compare the seasonal responses of vegetation to temperature and precipitation changes. Due to climatic fluctuations and complicated physiological responses of plants to growth environments, the effects of climate change on vegetation are usually distinct in different seasons [32][33][34]. Therefore, the seasonal responses of vegetation in the Songnen Plain marshes to temperature and precipitation change are needed to further clarify. During the past decades, marshes in the Songnen Plain have undergone obvious loss and fragmentation [10,35]. In order to protect the ecological environment in the Songnen Plain of China, the Chinese government is taking massive action to protect and restore degraded wetlands in the Songnen Plain [36]. Investigating the vegetation change and its response to climate change in marshes of the Songnen Plain is important to provide valuable implications for the protection and restoration of wetland vegetation in this region. Based on daily maximum temperature (T max ) and minimum temperature (T min ) data, previous studies found asymmetric warming of global surface air temperatures: Nighttime minimum temperature increased faster than the daytime maximum temperature during the past decades [37,38]. Additionally, this asymmetric warming pattern would be more likely during day and night, instead of symmetric one [37,39]. It is interesting that some studies found asymmetric (i.e., different) effects of daytime and night-time temperatures on vegetation growth in the semi-arid region of China [37,40]. Shen et al. (2015) found that the night-time temperature had a positive effect, but the daytime temperature had a negative effect on the growth of temperate grasslands in China. It may be because warm daytime temperature could reduce water efficiency by enhancing evaporation, while enhancing leaf respiration caused by night-time warming may promote photosynthesis in the next day [37,40]. At present, it is still unclear whether there are asymmetric daytime and night-time temperatures on the growth of marsh vegetation in this region. Understanding the responses of marsh vegetation to day and night temperature can contribute to further explain the effects of asymmetric warming on ecosystem vegetation. In order to accurately predict future changes in the vegetation of Songnen Plain marshes, it is necessary to explore the separate effects of daytime and night-time temperatures on the growth of marsh vegetation in the Songnen Plain.
By using the normalized difference vegetation index (NDVI) data, this study analyzed the temporal and spatial changes of vegetation in the Songnen Plain marshes. Based on the climate data from 2000 to 2016, the relationships of the growing season NDVI with temperatures (including mean temperature, maximum temperature, and minimum temperature, referred to as T mean , T max, and T min hereafter) and precipitation were analyzed to explore the possible effects of climate change on marsh vegetation. In addition, we compared the responses of vegetation to temperature and precipitation changes in different months of the growing season. Based on previous related studies discussed above, we predict that marsh vegetation growth during 2000-2016 has been affected by changes in precipitation and temperature. The major scientific questions that this study seeks to address are: (1) Is the vegetation of marsh in Songnen Plain more affected by precipitation or temperature? (2) whether the maximum temperature and minimum temperature also have the asymmetric effects on the vegetation of marsh in Songnen Plain? As marshes in the Songnen Plain are a representative semi-arid marsh wetland in the world, this study can be expected to provide useful information about the effects of climate change on wetland vegetation in semi-arid regions, which will have implications for biodiversity, carbon cycling, and sustainability of wetland ecosystems. More specifically, it can contribute to the understanding of how temperature and precipitation regimes in recent years have affected wetland vegetation, and how the changes of these regimes associated with climate change could affect wetland functions in the future, a largely unresolved question in current wetland studies.

Study Area
The Songnen Plain is located in the central part of Northeast China, lying between 42 • 19 N-49 • 12 N and 121 • 36 E-126 • 36 E. This study area includes 35 counties of Heilongjiang province and 20 counties of Jilin province, with a total area of 11.51 × 10 4 km 2 ( Figure 1). The climate in the Songnen Plain is a semi-humid climate. The annual average temperature ranges from 4 • C to 6 • C, and total precipitation is from 350 to 650 mm, 80% of which is concentrated between May and September. Averaged over the period of 2000 to 2016, the mean annual temperature is about 4.61 • C and decreases from south to north. The average annual precipitation is about 475.92 mm in the Songnen Plain, increasing gradually from west to east ( Figure 2). The soil in the Songnen Plain mainly includes black soil, chernozem, meadow soil, Solonetz, Solonchak, and aeolian soil [18]. Natural vegetation in the study area is dominated by temperate steppe and temperate meadow. The Songnen Plain has a large area of marshes, which is an important distribution area of marsh wetlands in China [41]. The dominant native wetland plant generally includes Astragalus, Artemisia, Carex, Phragmites, Potamogeton, Polygonum, Salix, Adenophora, Allium, Chenopodium, Vicia, Corispermum, and Ranunculus [29].

Data
In this study, we used monthly average surface air temperature (T mean , T max, and T min ) and precipitation data during 2000-2016 from 19 weather stations distributed in the Songnen plain. These data were provided by the National Meteorological Information Center of China and were processed from daily temperature and precipitation dataset. The data quality of these data was guaranteed with a vigorous data assurance policy [38,42]. Normalized Difference Vegetation Index (NDVI) is defined as the difference between the red and near-infrared (NIR) reflectance divided by their sum [43,44]. As a good indicator of vegetation growth, NDVI correlates very well with vegetation coverage and primary productivity [45]. Amongst remotely sensed data from satellites, NDVI has been widely used to investigate the relationship between climate variation and vegetation growth because of its good quality [46][47][48]. To analyze vegetation change and its response to climate change, this study collected the Moderate Resolution Imaging Spectrometer (MODIS) NDVI data from 2000 to 2016, with a temporal resolution of 16 days and 250 × 250 m spatial resolution. Based on pixel reliability parameters, the initial quality control of these MODIS NDVI datasets was made to minimize their errors and biases [49]. In addition, two raster maps of marshes distribution were used in this study. These two maps of marsh distribution were provided by the Chinese Academy of Sciences, the accuracy of which was evaluated based on field survey [45]. By comparing these two periods of marsh distribution data in the Songnen Plain, we found that the loss in the area of marshes was 2300 km 2 , and the numbers of marsh patches increased by 45 during 2000-2016. The results are generally consistent with previous studies [10,50]. In order to reduce the possible impacts of land use change (e.g., marshes loss) on our results, this study only chose and analyzed the unchanged marshes in the Songnen Plain during the study period. However, it should also be acknowledged that it may not completely exclude the possible effects of wetland degradation and human activities, such as grazing, on our results, although we focused on the unchanged marshes.

Methods
According to previous studies, the growing season is defined as the period from May to September in the Songnen Plain [51]. Considering that the growth of vegetation mainly concentrates in the growing season [52,53], this study only investigates the spatiotemporal changes of the NDVI and climate factors during the growing season of 2000-2016. To further reduce the effects of atmosphere and clouds, solar angle on NDVI data, the Maximum Value Compositing (MVC) method was used to reconstruct the original NDVI data into the monthly NDVI data set [54,55]. According to the MVC method, the equation to calculate NDVI can be described as: where MNDVI j means the largest NDVI value for month j, j was the number of the study period (May to September), NDVI 1 and NDVI 2 mean NDVI value for the first and second half of a month, respectively. We used simple linear regression to calculate the trends of the NDVI and climate factors in this study [56]. The trend of NDVI at each pixel is described as follows: where Slope was the trend of NDVI, n was the number of studied years (in this study, n = 17 years), i was the ordinal number of the year, and NDVI i means the NDVI value for year i.
To investigate the possible effects of climate change on the NDVI, correlation coefficients between monthly NDVI and temperature (T mean , T max, and T min ), total precipitation was calculated according to the following formula: where R x and R y are the correlation coefficients between variable x and y, x i is the NDVI for year i, and y i is the temperature (T mean , T max, and T min ) or total precipitation for year i, x and y are the averaged NDVI and temperature (T mean , T max, and T min ) or total precipitation for all years. In this study, the average values of NDVI, precipitation, and temperature for marshes were obtained from the averages of all grid cells belonging to the marshes [28,34].
Considering the possible interactions between different climatic factors, partial least-squares (PLS) regression analysis [57] was also used to further test the impacts of precipitation and temperature on NDVI. A model coefficient (MC) and a Variable-Importance-in-the-Projection (VIP) value were calculated for each independent variable [45,57]. In PLS regression, we set the monthly average temperature and precipitation as the independent and monthly NDVI as the dependent variable for each pixel. According to previous studies, variables with VIP equal or greater than 0.8 are considered important for the model [45,57]. Compared with conventional regression approaches, PLS regression avoids the problem of autocorrelation between independent factors, such as the climatic factors of the individual months [45]. Due to its good performance, PLS regression has been widely used in analyzing the relationships between the dependent and independent variables [57].

Spatiotemporal Change of NDVI in the Songnen Plain Marshes
From 2000 to 2016, the growing season mean NDVI (referred to as the growing season NDVI hereafter) of marshes significantly increased by 0.06/decade over the Songnen Plain (p < 0.01). In different months of the growing season, monthly NDVI of marshes similarly increased significantly during the study period (Figure 3), with the largest increasing trends in May, June, and July (0.06/decade, 0.10/decade, 0.07/decade, respectively). Spatially, the changes of NDVI had apparent spatial heterogeneities in the Songnen Plain marshes from 2000 to 2016 (Figure 4). The growing season NDVI increased in 84% of the marshes in the study area, with the most obvious increase of NDVI in the central and south of Songnen Plain. By contrast, the reduced growing season NDVI accounted for only 13% of the marshes in the study area, which was mainly concentrated in the northeast of Songnen Plain (Figure 4). During the growing season, the spatial patterns of the monthly NDVI trend were generally similar to that of the growing season NDVI trend. However, the monthly average NDVI in August showed an obvious decreasing trend in the marshes of the central and south of Songnen Plain (Figure 4). Table 1 shows the correlations between the NDVI and temperature and precipitation for the whole marshes of the Songnen Plain. We found a significant positive relationship between the growing season precipitation with the growing season NDVI of the whole marshes in the Songnen Plain (p < 0.05). There was a moderately negative relationship between the growing season temperature and NDVI of the whole marshes in the Songnen Plain (p > 0.05). In addition, our results showed opposite effects of the growing season T max and T min on the growing season NDVI of the Songnen Plain marshes: The growing season NDVI had a negative correlation with T max, but had a positive correlation with T min (Table 1).  Table 1. Correlations between climate variables (precipitation, T mean , T max, and T min ) and NDVI of the Songnen Plain marshes, based on Pearson correlation analysis.

T mean T max T min Precipitation
Growing To further investigate the impacts of climate change on NDVI of the Songnen Plain marshes, we calculated the correlation between monthly NDVI and precipitation and temperature ( Table 1). The results showed that precipitation had a positive relationship with NDVI of marshes in all the months of the growing season, with the exception in August. Moreover, the positive correlation between precipitation and NDVI reached a significant level in May (p < 0.01). In August, however, we found a negative relationship of precipitation with NDVI of the Songnen Plain marshes. In terms of temperature impacts, we found that the monthly T mean had a positive correlation with NDVI in August, but had a negative correlation with NDVI in other months (Table 1). For the effects of daytime and night-time temperatures, we found asymmetric effects of T max and T min on NDVI of the Songnen Plain marshes. Similar to the growing season results, there were opposite effects of daytime T max and nighttime T min on NDVI of marshes for all the months of the growing seasons except August. By contrast, both T max and T min had a positive relationship with NDVI in August (Table 1).
Based on the PLS regression, we also found that the NDVI was strong positively correlated with precipitation in the growing season ( Figure 5). By contrast, the mean temperature had a weak negative relationship with the growing season NDVI. Consistent with correlation analysis results, the PLS regression also showed asymmetric effects of the growing season T max and T min on the growing season NDVI of marshes: The growing season NDVI was strong negatively associated with T max , but was strong positively associated with T min in the Songnen Plain marshes. For different months, all the PLS regression results also confirmed our findings about the effects of monthly temperature (including T mean , T max, and T min ) and precipitation on NDVI based on the correlation analysis results. The only slight difference between the results from PLS regression and correlation analysis was the different significance level for the effect of the climatic variable on NDVI, which is actually reasonable due to different statistical methods. Figure 5. Response of NDVI of the Songnen Plain marsh to temperature (including T mean , T max, and T min ) and precipitation change based on partial least-squares (PLS) regression. Model coefficients (MC) and the variable importance plots (VIP) values reflect the degree and the significance of the impact of climatic factors on the NDVI, respectively. GS refers to the whole growing season (May to September). We set the NDVI as the dependent variable and precipitation, T mean , T max and T min as the independent variables.
Spatially, the negative correlation between the growing season NDVI and precipitation was only concentrated in the north of Songnen Plain. The positive correlations between NDVI and temperature were also concentrated in the north of Songnen Plain. These opposite effects of T max and T min were more obvious in the central and north of the study area ( Figure 6). Similar to the growing season results, the spatial correlations between monthly NDVI and precipitation, T mean , T max, and T min have apparent spatial heterogeneities in the Songnen Plain marshes (Figure 7). With the exception of August, these were obviously asymmetric effects of T max and T min on NDVI in the central and north of the study area, which was consistent with the growing season results. In August, however, there are positive effects of both T max and T min on NDVI in most parts of the study area. Different from the results of the growing season and other months, precipitation in August had a significantly negative relationship with NDVI of marshes in the north of Songnen Plain.

The Correlations between Climate Variables and NDVI of the Songnen Plain Marshes
Over the whole marshes of the Songnen Plain, the growing season NDVI had a significant positive relationship with the growing season precipitation but had a weak negative correlation with temperature. These results suggest that, in recent years, the growing season NDVI has been mainly affected by precipitation and increasing precipitation could promote the growth of vegetation in the growing season. This finding was consistent with , who found that the vegetation of marshes is mainly affected by precipitation in arid and semi-arid regions of northeast China.
By further analyzing the relationships of NDVI with daytime and nighttime temperature, we first found asymmetric effects of T max and T min on NDVI of the Songnen Plain marshes. The growing season NDVI was negatively correlated with T max , but positively correlated with T min in the Songnen Plain marshes. This result indicates that daytime and night-time warming could prohibit and promote marsh vegetation growth in the Songnen Plain, respectively. In arid and semi-arid Songnen Plain, it is reported that warm daytime temperature could limit the growth of vegetation by enhancing evaporation and reducing soil water content [58]. By contrast, night-time warming may promote the growth of vegetation due to a rebound effect of compensatory stimulated photosynthesis in the next day [37,59]. These explanations may account for the asymmetric effects of daytime and night-time warming on NDVI of the Songnen Plain marshes.
In different months, we found that the precipitation had a negative relationship with NDVI in August, but had a positive correlation with NDVI in the other months. By contrast, the temperature had a positive relationship with NDVI of marshes in August but had a negative correlation with NDVI in the other months. The results of correlation analysis indicate that the NDVI was determined by precipitation and increasing precipitation could promote the growth of vegetation in most months of the growing season. In August, however, warm temperature could promote the growth of vegetation, and increasing precipitation had adverse impacts on vegetation growth. This may be because the precipitation in the Songnen Plain is the largest in August when the temperature becomes a limit factor affecting vegetation growth [16], and thus warm temperature is beneficial to the vegetation growth by enhancing photosynthesis [60]. In different months, we also found asymmetric effects of T max and T min on NDVI of marshes in May, June, July, and September. It indicates that warming daytime T max and nighttime T min will limit and promote the growth of vegetation in these months, respectively.
For the spatial correlations, the growing season precipitation generally had a positive relationship with the growing season NDVI, with the negative correlation only scattered in the north of Songnen Plain (p < 0.01). By contrast, there was a generally negative correlation between the growing season temperature and NDVI, with a positive correlation in the north of Songnen plain. The reason for this phenomenon may be because the temperature becomes a limit factor for vegetation growth in the north of Songnen Plain, where the temperature is relatively low due to high latitude. As a result, warm temperature promotes vegetation growth, but increasing precipitation is harmful to the growth of vegetation because of the low temperature caused by precipitation [34,40]. According to previous studies, the main vegetation types of marshes include Phragmites australis and Carex, which are mainly distributed in the north and south of Songnen Plain, respectively [61][62][63]. In the north of Songnen Plain, the water content is relatively abundant and temperature becomes the limit factor affecting the growth of Phragmites australis. It is reported that low temperatures caused by precipitation can limit the absorption of nutrition by plant roots and thus prohibits Phragmites australis growth of marshes [64]. In the south of Songnen Plain, the growth of Carex was mainly affected by the precipitation due to limited water content in the seasonal marsh regions. The different growth strategies of these two vegetation types [65] may partly explain the opposite effects of precipitation on NDVI in the north and south of marshes in the Songnen Plain.

Vegetation Change in the Songnen Plain Marshes
To further explain the temporal changes of NDVI in the Songnen Plain marshes, we analyzed the changes in temperature and precipitation during the study period. The results show that the growing season total precipitation significantly increased by 18.5 mm/decade (p < 0.05), but the growing season mean temperature showed a weak declining trend (−0.3 • C/decade) (Figure 8). In addition, there were opposite changes of the growing season T max and T min : The growing season T max showed a weak declining trend (−0.07 • C/decade), but the growing season T min showed a weak increasing trend   Based on the results of correlations between climate variables (precipitation, T mean , T max, and T min ) and NDVI (Table 2, Figure 3), we concluded that the increasing growing season precipitation accounts for the increasing growing season NDVI of the Songnen Plain marshes. For different months, the significantly increasing precipitation and decreasing T min may partly explain the NDVI increase in May, June, and September. In August, although the increasing temperature was likely to promote the growth of vegetation, the increasing precipitation could have exerted negative effects on the vegetation growth, which partly explains the weak increasing NDVI in August.
Spatially, the growing season NDVI of marshes in the Songnen Plain significantly increased during 2000-2016, with the largest increase of NDVI in the central and south of Songnen Plain (p < 0.01) ( Figure 4). Interestingly, we found the growing season precipitation showed the largest increase in the central and south of Songnen Plain (Figure 9). Therefore, the most obvious increase in the growing season precipitation in the central and south of Songnen Plain seems to explain the largest increase of NDVI of marshes in these regions. In August, the monthly precipitation of marshes showed an obvious decreasing trend in the central and south of Songnen Plain (Figure 10), which may partly account for the moderate decrease of NDVI in the marshes of these regions. Previous studies found that the marshes in the north of Songnen Plain have undergone obvious loss and degradation due to human activities [63]. In this study, we found that the reduced growing season NDVI was mainly concentrated in the northeast of Songnen Plain. Considering that the marsh loss and degradation can have obvious adverse effects on marsh vegetation [10], the marsh loss and degradation in the northeast of Songnen Plain may partly account for the decrease of NDVI in marshes of these regions. In order to restore the degraded marshes, some marsh restoration projects have been implemented by local governments in the central of Songnen Plain during recent years [10]. For example, a project of "diverting the Nenjiang River to the Momoge" was implemented in the Momoge region since 2002 to restore the degraded wetlands [66]. Through the implementation of a series of marshes restoration projects, the ecological environment and the vegetation coverage conditions have been greatly improved [67]. In this study, it is interesting to note that the largest increase of the growing season NDVI was concentrated in the central of Songnen Plain during 2000-2016. Therefore, it is likely that the ecological restoration project in the central of Songnen Plain may also partly explain the obvious increase of NDVI in this region [68].

Conclusions
Based on NDVI and climate data from 2000 to 2016, this study analyzed the spatiotemporal change in NDVI and the possible impact of precipitation and temperature (including T mean , T max, and T min ) on growth vegetation during the growing seasons. During 2000-2016, the growing season NDVI of marshes significantly increased at a rate of 0.06/decade over the whole Songnen Plain. Spatially, the changes of NDVI had apparent spatial heterogeneities in marshes, and the most obvious increase of vegetation was concentrated in the central and south of Songnen Plain. The results of correlation analysis indicated that the vegetation growth was determined by precipitation and the increased growing season precipitation could promote the growth of marsh vegetation in the Songnen Plain. By contrast, warmer growing season temperatures could inhibit vegetation growth, possibly due to warming-induced water deficiency. In August, however, we found that temperatures had a positive effect, but precipitation had a negative effect on the growth of marsh vegetation in the Songnen Plain. This phenomenon may be related to the fact that the amount of precipitation in the study area is the largest in August, when the precipitation is abundant for the vegetation growth, and warming temperatures are beneficial to vegetation growth by enhancing photosynthesis. In addition, this study first found asymmetric effects of daytime T max and nighttime T min on NDVI of the Songnen Plain marshes: The NDVI correlated negatively with T max but positively with T min in May, June, July, and September. It indicates that warming daytime and nighttime temperatures could have negative and positive effects on vegetation growth in the Songnen Plain marshes, respectively. Under the background of global asymmetric warming, we should pay more attention to these asymmetric effects of T max and T min on marsh vegetation, especially when it comes to modeling the changes in wetland vegetation caused by climate change in the Songnen Plain of Northeast China.
It should be noted that there may exist some limitations in the current study. First, the NDVI data with a resolution of 15 days could contain some uncertainties due to the effects of clouds, atmosphere, and solar angle. In addition, each NDVI pixel could only roughly reflect the conditions of vegetation within 250 × 250 m area. Second, the Maximum Value Compositing (MVC) and the interpolation method may contain some uncertainties. These uncertainties could have some effects on the results of vegetation change and its relationship with climate change in this study. Third, the correlations between NDVI and climate variables are actually more complicated because many other climatic and environmental factors may also have significant effects on NDVI. Further studies are still needed to investigate the possible impacts of other climatic and environmental factors as well as human activities on marsh vegetation in the Songnen Plain.