Using Modified Remote Sensing Imagery to Interpret Changes in Cultivated Land under Saline-Alkali Conditions

Managing the rapidly changing saline-alkali land under cultivation in the coastal areas of China is important not only for mitigating the negative impacts of such land on the environment, but also for ensuring long-term sustainability of agriculture. In this light, setting up rapid monitoring systems to assist decision-making in developing sustainable management plans is therefore an absolute necessity. In this study, we developed a new interpretation system where symbols are used to grade and classify saline-alkali lands in space and time, based on the characteristics of plant cover and features of remote sensing images. The system was used in combination with the maximum likelihood supervised classification to analyze the changes in cultivated lands under saline-alkali conditions in Huanghua City. The analysis revealed changes in the area and spatial distribution of cultivated under saline-alkali conditions in the region. The total area of saline-alkali land was 139,588.8 ha in 1992 and 134,477.5 ha in 2011. Compared with 1992, severely and moderately saline-alkali land areas decreased in 2011. However, non/slightly saline land areas increased over that in 1992. The results showed that the salinization rate of arable lands in Huanghua City decreased from 1992 to 2011. The moderately saline-alkali land southeast of the city transformed into non/slightly saline-alkaline. Then, severely saline-alkali land far from the coastal zone west of the city became moderately saline-alkaline. Spatial changes in cultivated saline-alkali lands in Huanghua City were such that the centers of gravity (CG) of severely and non/slightly saline-alkali land moved closer the coastline, while that of the moderately saline-alkali land moved from southwest coastal line to northwest. Factors influencing changes in cultivated lands in the saline-alkali ecosystem included climate, hydrology and human activity. Thus, studies are required to further explore these factors in order to build a better understanding into the relative contributions of the changes saline-alkali state on the functions of coastline ecosystems.


Introduction
Globally, some 83 million hectares of land is affected by salt, with a wide range of distribution in five continents and over 100 countries [1].Soil salinization is so much so extensive that it now spreads at the rate of two million hectares per year [2,3].This requires the development of effective preventive and management measures, including rapid reliable means of getting information on spatial and temporal distributions of saline-alkali lands across the globe.
Remote sensing techniques and data are widely used to detect and map salt-affected areas [4][5][6].Dwivedi and Rao [7] noted that the combination of Landsat TM bands 1, 3 and 5 was most suitable for the interpretation of salt-affected soils.Dwivedi and Sreenivas [8] believed that a substantial understanding of saline-alkali lands could be developed through principal component analysis, image differencing and band rationing.Metternichet and Zinck [9] used remote sensing images in combination with field survey and laboratory data to analyze the distribution of soil salinization.The choice of the time of remote sensing image is critical for the extraction of reliable soil salinization data.The spatial and temporal distributions of soil salinization in Jiefangzha Irrigation Sub-district were determined by analyzing remote sensing images acquired at the end of March.The irrigation sub-district is in the western part of Hetao Irrigation District in Inner Mongolia, China [10].Then, Wang [11] used remote sensing images for August and September to extract data on the salinization of the delta oasis of Weigan and Kuqa rivers.
In this study, we used remote sensing images acquired at the end of May or the beginning of June to determine the spatial and temporal variability of soil salinity in a coastal area in Huanghua City.At the time, soils with different levels of salts were under different vegetation covers of different reflectivity.As the features of each land use class in the composite TM images was clearly delineated in the standard false color, the images did not require complex transformation and index calculations in order to interpret the features.We developed an interpretation symbol system to grade and classify saline-alkali soils based on vegetation characteristics.The system was developed using data from site investigation, plant cover characteristics and remote sensing features.We used TM remote sensing image of standard false color composite (band 4, 3 and 2) without any complex image transformation and index calculation.Finally, the interpretation system was integrated with the maximum likelihood supervised classification method to derive relevant on saline-alkali soils.

Study Site
The west Bohai Sea in Huanghua City is part of a coastal alluvial plain in Hebei Province, China (Figure 1), covering a vast area of saline-alkali land in the region.The total area of the city is 220,000 ha, a large proportion (129,462.3ha) of which is cultivated.The main land cover includes water, marshland and residential land, with respective areas of 34,675.3ha, 25,752.2ha and 18,103.8ha.The other forms of land cover account for a very small fraction of the study area [12].The climate in the study area is warm temperate, semi-humid monsoon, with a mean annual precipitation of about 627 mm, occurring mainly between July and September.The mean annual potential evaporation is 1980 mm.The area, which is prone to summer drought and spring waterlogging, has shallow groundwater salinity of above 3 g/L.As there is lack of freshwater in the study area, it is mainly covered by salt-resistant vegetation.

Data Source
The remote sensing data, including two sets of Landsat (Row123/Path 32 of the WRS-2) and Thematic Mapper (TM) images (dated 27 May 1992 and 1 June 2011), were downloaded from the public domain [13].The Landsat scenes were acquired under clear sky conditions at the end of May or beginning of June.To eliminate bias in the comparison of the multi-temporal datasets, the data were referenced to a common coordinate grid.Image-to-image calibration was adopted and TM image with accurate geometric correction taken as reference image.In addition, data for 1992 and 2011 were geo-referenced to UTM coordinate system and the accuracy controlled within 0.5 grid (with 1 grid = 30 m ˆ30 m).The geographic referenced data were then clipped to the study area in ARCGIS.

Methods
A modified remote sensing imagery was used to interpret changes in cultivated lands under saline-alkali conditions in Huanghua City.The schematic representation of the modified method used to process the remote sensing images is depicted in Figure 2. Maximum likelihood Supervised classification was used to analyze land classification.Maximum likelihood classification (MLC) [14] is a method for determining a known class of distributions as the maximum for a given statistic.An assumption of normality is made for the training samples.The algorithm builds the probability density functions for each category.During classification, all unclassified pixels are assigned membership based on the relative likelihood (probability) of that pixel occurring within each category's probability density function [15].
Suppose there are G predefined categories, the unclassified image possesses m bands.According to Bayesian formula, the posterior probability of category k, P(Gk/x) is defined as:

Methods
A modified remote sensing imagery was used to interpret changes in cultivated lands under saline-alkali conditions in Huanghua City.The schematic representation of the modified method used to process the remote sensing images is depicted in Figure 2.

Methods
A modified remote sensing imagery was used to interpret changes in cultivated lands under saline-alkali conditions in Huanghua City.The schematic representation of the modified method used to process the remote sensing images is depicted in Figure 2. Maximum likelihood Supervised classification was used to analyze land classification.Maximum likelihood classification (MLC) [14] is a method for determining a known class of distributions as the maximum for a given statistic.An assumption of normality is made for the training samples.The algorithm builds the probability density functions for each category.During classification, all unclassified pixels are assigned membership based on the relative likelihood (probability) of that pixel occurring within each category's probability density function [15].
Suppose there are G predefined categories, the unclassified image possesses m bands.According to Bayesian formula, the posterior probability of category k, P(Gk/x) is defined as: Maximum likelihood Supervised classification was used to analyze land classification.Maximum likelihood classification (MLC) [14] is a method for determining a known class of distributions as the maximum for a given statistic.An assumption of normality is made for the training samples.The algorithm builds the probability density functions for each category.During classification, all unclassified pixels are assigned membership based on the relative likelihood (probability) of that pixel occurring within each category's probability density function [15].Suppose there are G predefined categories, the unclassified image possesses m bands.According to Bayesian formula, the posterior probability of category k, P(Gk/x) is defined as: where P(G k ) is the prior probability of category k; P(x/Gk) is conditional probability of observing x from G k (probability density function); and P(x) is the same for each category.If we have no knowledge about the prior distributions P(G k ), then we can assume that all categories are probable.Therefore, the likelihood function is determined by P(x/G k ), which is also called the likelihood of G k with respect to x.Given the logarithm, the discriminant function can be expressed as: where x = [x 1 , x 2 , . . ., x m ] T is the vector of a pixel; M k (x) is the likelihood function of x belonging to category k; and µ k and S k are the mean vector and covariance matrix of the Kth category.
A post-classification comparison was used to analyze land use changes.Here, the matrix tool in the image interpreter module of ERDAS IMAGINE software was used to get the map of spatial changes in a land use and transition matrix of different land use classes land for 1992-2011.
Then the center of gravity migration model was used to describe the overall intensity of change as well as the spatial variation characteristics of different degrees of the saline-alkali land.The spatial migration of the gravity center was demonstrated by the change in gravity center coordinate [16][17][18].The gravity center coordinate was calculated as: Y " where X and Y are the gravity center coordinate for cultivated land under saline-alkali conditions in the Huanghua City ( ˝), P j is the area of the jth patch (ha), X j and Y j are gravity center coordinate of the jth patch ( ˝), and n is the number of patches in the cultivated lands under saline-alkali conditions.

Saline-Alkali Land Classification
The land use in Huanghua City was divided into the following six classes based on the state land use classification standard [19].The classes included water body, salt pot, residential land, aquaculture and cultivated land.Cultivated land was further divided into three classes based on spectral signatures of land covers-severely saline-alkali land, moderately saline-alkali land and non/lightly saline-alkali land (Figure 3).

New Interpretation Symbol System
To determine the spatial and temporal changes in cultivated lands under saline-alkali conditions in the coastal region, a new interpretation symbol system consistent with the saline-alkali land classification was set up.Zhang [20] defined the classes of saline-alkali soils for Northeast China and developed an interpretation symbol system that combined image element (color) with spatial location.Dwivedi [21] set up an interpretation symbol system that combined image elements (color, texture, pattern, shape, size, association, etc.) and laboratory analysis of field soil samples data.The interpretation symbol system used in this study is based on site investigation, plant cover characteristics and remote sensing images features (Table 1).We used the spectral signature of characteristic plants in coastal cultivated land for the end of May or early June.It was driven by the vegetation characteristics at the site during the time of study.Severely saline-alkali lands were covered by salt-resistant vegetation such sa Suaeda glauca Bunge and Phragmites autralis (Figure 3), which was almost brown but had sporadic red.Moderately saline-alkali lands were covered by spring corn or cotton (Figure 3), with apparent gray and red image features in the ratio of about 50%.The image feature of winter wheat planted in non/slightly saline-alkali lands (Figure 3) was apparently red.

Classes of Saline-Alkali Land
Vegetation Feature Image Feature

Severely saline-alkali land Salt-resistant vegetation
Moderately saline-alkali land Spring corn, cotton Non/slightly saline-alkali land Winter wheat

New Interpretation Symbol System
To determine the spatial and temporal changes in cultivated lands under saline-alkali conditions in the coastal region, a new interpretation symbol system consistent with the saline-alkali land classification was set up.Zhang [20] defined the classes of saline-alkali soils for Northeast China and developed an interpretation symbol system that combined image element (color) with spatial location.Dwivedi [21] set up an interpretation symbol system that combined image elements (color, texture, pattern, shape, size, association, etc.) and laboratory analysis of field soil samples data.The interpretation symbol system used in this study is based on site investigation, plant cover characteristics and remote sensing images features (Table 1).We used the spectral signature of characteristic plants in coastal cultivated land for the end of May or early June.It was driven by the vegetation characteristics at the site during the time of study.Severely saline-alkali lands were covered by salt-resistant vegetation such as Suaeda glauca Bunge and Phragmites autralis (Figure 3), which was almost brown but had sporadic red.Moderately saline-alkali lands were covered by spring corn or cotton (Figure 3), with apparent gray and red image features in the ratio of about 50%.The image feature of winter wheat planted in non/slightly saline-alkali lands (Figure 3) was apparently red.

Classes of Saline-Alkali Land Vegetation Feature Image Feature
Severely saline-alkali land Salt-resistant vegetation

New Interpretation Symbol System
To determine the spatial and temporal changes in cultivated lands under saline-alkali conditions in the coastal region, a new interpretation symbol system consistent with the saline-alkali land classification was set up.Zhang [20] defined the classes of saline-alkali soils for Northeast China and developed an interpretation symbol system that combined image element (color) with spatial location.Dwivedi [21] set up an interpretation symbol system that combined image elements (color, texture, pattern, shape, size, association, etc.) and laboratory analysis of field soil samples data.The interpretation symbol system used in this study is based on site investigation, plant cover characteristics and remote sensing images features (Table 1).We used the spectral signature of characteristic plants in coastal cultivated land for the end of May or early June.It was driven by the vegetation characteristics at the site during the time of study.Severely saline-alkali lands were covered by salt-resistant vegetation such sa Suaeda glauca Bunge and Phragmites autralis (Figure 3), which was almost brown but had sporadic red.Moderately saline-alkali lands were covered by spring corn or cotton (Figure 3), with apparent gray and red image features in the ratio of about 50%.The image feature of winter wheat planted in non/slightly saline-alkali lands (Figure 3) was apparently red.

New Interpretation Symbol System
To determine the spatial and temporal changes in cultivated lands under saline-alkali conditions in the coastal region, a new interpretation symbol system consistent with the saline-alkali land classification was set up.Zhang [20] defined the classes of saline-alkali soils for Northeast China and developed an interpretation symbol system that combined image element (color) with spatial location.Dwivedi [21] set up an interpretation symbol system that combined image elements (color, texture, pattern, shape, size, association, etc.) and laboratory analysis of field soil samples data.The interpretation symbol system used in this study is based on site investigation, plant cover characteristics and remote sensing images features (Table 1).We used the spectral signature of characteristic plants in coastal cultivated land for the end of May or early June.It was driven by the vegetation characteristics at the site during the time of study.Severely saline-alkali lands were covered by salt-resistant vegetation such sa Suaeda glauca Bunge and Phragmites autralis (Figure 3), which was almost brown but had sporadic red.Moderately saline-alkali lands were covered by spring corn or cotton (Figure 3), with apparent gray and red image features in the ratio of about 50%.The image feature of winter wheat planted in non/slightly saline-alkali lands (Figure 3) was apparently red.

New Interpretation Symbol System
To determine the spatial and temporal changes in cultivated lands under saline-alkali conditions in the coastal region, a new interpretation symbol system consistent with the saline-alkali land classification was set up.Zhang [20] defined the classes of saline-alkali soils for Northeast China and developed an interpretation symbol system that combined image element (color) with spatial location.Dwivedi [21] set up an interpretation symbol system that combined image elements (color, texture, pattern, shape, size, association, etc.) and laboratory analysis of field soil samples data.The interpretation symbol system used in this study is based on site investigation, plant cover characteristics and remote sensing images features (Table 1).We used the spectral signature of characteristic plants in coastal cultivated land for the end of May or early June.It was driven by the vegetation characteristics at the site during the time of study.Severely saline-alkali lands were covered by salt-resistant vegetation such sa Suaeda glauca Bunge and Phragmites autralis (Figure 3), which was almost brown but had sporadic red.Moderately saline-alkali lands were covered by spring corn or cotton (Figure 3), with apparent gray and red image features in the ratio of about 50%.The image feature of winter wheat planted in non/slightly saline-alkali lands (Figure 3) was apparently red.

Classes of Saline-Alkali Land
Vegetation Feature Image Feature

Severely saline-alkali land Salt-resistant vegetation
Moderately saline-alkali land Spring corn, cotton Non/slightly saline-alkali land Winter wheat

Classification Templates
The training sets about water body, salt pot, residential land and aquaculture were selected based on the national land survey.The training sets in 1992 and 2011 were obtained based on the first and second national land survey, respectively.The training sets about cultivated land under saline-alkali condition using visual interpretation were selected using based on the new interpretation symbol system in 1992 and 2011.Training samples were selected for each class in order to establish the classification templates.Over ten training samples were merged for each final class.The contingency matrix tool was used to evaluate the classification templates, and the related error matrix given as Tables 2 and 3.In 1992, Classification accuracy of other classes is over 90%, except for marsh and aquaculture.In 2011, Classification accuracy of other classes is over 90%, except for marsh only.Classification accuracy of each class in templates was less than 15% in 1992 and 2011, so the classification templates were considered to be successful.

Saline-Alkali Land Information Extraction
Based on the classification templates, the land use classification maps were established using maximum likelihood supervised method in the ERDAS IMAGINE software (Intergraph Corporation, Huntsville, AL, USA) (Figure 4).The maps of spatial and temporal variations in the saline-alkali lands were obtained through post-classification comparison in ARCGIS environment (Figure 8).

Classification Templates
The training sets about water body, salt pot, residential land and aquaculture were selected based on the national land survey.The training sets in 1992 and 2011 were obtained based on the first and second national land survey, respectively.The training sets about cultivated land under saline-alkali condition using visual interpretation were selected using based on the new interpretation symbol system in 1992 and 2011.Training samples were selected for each class in order to establish the classification templates.Over ten training samples were merged for each final class.The contingency matrix tool was used to evaluate the classification templates, and the related error matrix given as Table 2 and Table 3.In 1992, Classification accuracy of other classes is over 90%, except for marsh and aquaculture.In 2011, Classification accuracy of other classes is over 90%, except for marsh only.Classification accuracy of each class in templates was less than 15% in 1992 and 2011, so the classification templates were considered to be successful.

Saline-Alkali Land Information Extraction
Based on the classification templates, the land use classification maps were established using maximum likelihood supervised method in the ERDAS IMAGINE software (Intergraph Corporation, Huntsville, AL, USA) (Figure 4).The maps of spatial and temporal variations in the saline-alkali lands were obtained through post-classification comparison in ARCGIS environment (Figure 8).

Validation of Remote Sensing Interpretation
We use statistical yearbook and previous research to verify accuracy of the land use classification maps of 1992 and field survey to verify accuracy of the land use classification maps of 2011.According to statistical yearbook and previous research [22], the area of cultivated land was about 130,000 ha.The area of cultivated land is 139,588.

Temporal Changes in Arable Saline-Alkali Land
The dynamics of soil salinization in the severely, moderately and non/slightly saline-alkali lands for the period from 1992 to 2011 are plotted in Figure 4. From the satellite images for 1992 and 2011, we the different classes of cultivated land areas and the related ratios (%) to the total area were statistically demarcated in ARCGIS environment.Figure 5 shows the different classes of cultivated areas as percentages of total cultivated area and Figure 6 shows the changes in area in the different classes for 1992-2011.As in Figure 2

Temporal Changes in Arable Saline-Alkali Land
The dynamics of soil salinization in the severely, moderately and non/slightly saline-alkali lands for the period from 1992 to 2011 are plotted in Figure 4. From the satellite images for 1992 and 2011, we the different classes of cultivated land areas and the related ratios (%) to the total area were statistically demarcated in ARCGIS environment.Figure 5 shows the different classes of cultivated areas as percentages of total cultivated area and Figure 6 shows the changes in area in the different classes for 1992-2011.As in Figure 2   For the period from 1992 to 2011, the total area of saline-alkali lands decreased by 3.66% from 139,588.8 ha to 134,477.5 ha, respectively.The area under severely saline-alkali lands changed the most (35.2%),decreasing from 31,425.8 ha in 1992 to 20,373.2ha in 2011.The area under moderately saline-alkali lands had a similar trend as that under severely saline-alkali lands, decreasing by 18.7% from 37,984.8 ha in 1992 to 30,888.2ha in 2011.However, non/lightly saline-alkali lands had the reverse trend, increasing by 18.6% from 70,178.1 ha to 83,216.1 ha.

Temporal Changes in Arable Saline-Alkali Land
The dynamics of soil salinization in the severely, moderately and non/slightly saline-alkali lands for the period from 1992 to 2011 are plotted in Figure 4. From the satellite images for 1992 and 2011, we the different classes of cultivated land areas and the related ratios (%) to the total area were statistically demarcated in ARCGIS environment.Figure 5 shows the different classes of cultivated areas as percentages of total cultivated area and Figure 6 shows the changes in area in the different classes for 1992-2011.As in Figure 2   For the period from 1992 to 2011, the total area of saline-alkali lands decreased by 3.66% from 139,588.8 ha to 134,477.5 ha, respectively.The area under severely saline-alkali lands changed the most (35.2%),decreasing from 31,425.8 ha in 1992 to 20,373.2ha in 2011.The area under moderately saline-alkali lands had a similar trend as that under severely saline-alkali lands, decreasing by 18.7% from 37,984.8 ha in 1992 to 30,888.2ha in 2011.However, non/lightly saline-alkali lands had the reverse trend, increasing by 18.6% from 70,178.1 ha to 83,216.1 ha.For the period from 1992 to 2011, the total area of saline-alkali lands decreased by 3.66% from 139,588.8 ha to 134,477.5 ha, respectively.The area under severely saline-alkali lands changed the most (35.2%),decreasing from 31,425.8 ha in 1992 to 20,373.2ha in 2011.The area under moderately saline-alkali lands had a similar trend as that under severely saline-alkali lands, decreasing by 18.7% from 37,984.8 ha in 1992 to 30,888.2ha in 2011.However, non/lightly saline-alkali lands had the reverse trend, increasing by 18.6% from 70,178.1 ha to 83,216.1 ha.
The changes in ratios of the areas of land under different saline-alkali land classes as percentages of total cultivated areas in 1992 and 2011 were also analyzed.The areas under severely and moderately saline-alkali lands as percentage of the total cultivated land area decreased from 22.5 to 15.1% and from 27.2% to 23.0%, respectively.The area under non/slightly saline-alkali lands increased by 11.6% for the period from 1992 to 2011 (Figure 5).It suggested that there was a drop in the degree of salinization of arable lands in Huanghua City, China.
As in Tables 2 and 3, the error in each class of the classification templates was less than 15%, which influenced the classification results.For example, the classifications accuracies of marshland were 87.7% and 88.56% for 1992 and 2011, respectively, so the area of marsh is less than the actual area.Although, other land use classes also had similar issues.However, this form of error did not compromise the overall trend of changes in land use the study area.

Spatial Changes in Saline-Alkali Arable
Figure 4 shows the spatial distribution of saline-alkali land in the study area in 1992 and 2011.It is clear that the severely saline-alkali lands were concentrated in the east of Huanghua City, near the coast.Most of the moderately saline-alkali lands were in the south of the city in 1992, largely shifting to the north of the city in 2011.The non/slightly saline-alkali lands occurred far from the coast.
The characteristics of the spatial changes in saline lands can be explained not only qualitatively (using land use maps), but also statistically (using the changes in different land use classes).The detected characteristics in land salinization in 1992-2011 are depicted in Figure 7 and the change matrix detailed in Table 4. From the change matrix (Table 4), it was clear that 24,770.6 ha of moderately saline-alkali lands and 7470.18ha of severely saline-alkali lands were converted into non/slightly saline-alkali lands in the study area.In addition, 3575.16ha of severely saline-alkali lands were converted into moderately saline-alkali lands in the region over the same period.Furthermore, 12,125.1 ha and 6540.3 ha of non/slightly saline lands were converted into moderately and heavily saline-alkali lands, respectively.Similarly, 2234.9 ha of moderately saline-alkali lands were converted into severely saline-alkali lands.
Sustainability 2016, 8, 619 9 of 14 The changes in ratios of the areas of land under different saline-alkali land classes as percentages of total cultivated areas in 1992 and 2011 were also analyzed.The areas under severely and moderately saline-alkali lands as percentage of the total cultivated land area decreased from 22.5 to 15.1% and from 27.2% to 23.0%, respectively.The area under non/slightly saline-alkali lands increased by 11.6% for the period from 1992 to 2011 (Figure 5).It suggested that there was a drop in the degree of salinization of arable lands in Huanghua City, China.
As in Tables 2 and 3, the error in each class of the classification templates was less than 15%, which influenced the classification results.For example, the classifications accuracies of marshland were 87.7% and 88.56% for 1992 and 2011, respectively, so the area of marsh is less than the actual area.Although, other land use classes also had similar issues.However, this form of error did not compromise the overall trend of changes in land use the study area.

Spatial Changes in Saline-Alkali Arable
Figure 4 shows the spatial distribution of saline-alkali land in the study area in 1992 and 2011.It is clear that the severely saline-alkali lands were concentrated in the east of Huanghua City, near the coast.Most of the moderately saline-alkali lands were in the south of the city in 1992, largely shifting to the north of the city in 2011.The non/slightly saline-alkali lands occurred far from the coast.
The characteristics of the spatial changes in saline lands can be explained not only qualitatively (using land use maps), but also statistically (using the changes in different land use classes).The detected characteristics in land salinization in 1992-2011 are depicted in Figure 7 and the change matrix detailed in Table 4. From the change matrix (Table 4), it was clear that 24,770.6 ha of moderately saline-alkali lands and 7470.18ha of severely saline-alkali lands were converted into non/slightly saline-alkali lands in the study area.In addition, 3575.16ha of severely saline-alkali lands were converted into moderately saline-alkali lands in the region over the same period.Furthermore, 12,125.1 ha and 6540.3 ha of non/slightly saline lands were converted into moderately and heavily saline-alkali lands, respectively.Similarly, 2234.9 ha of moderately saline-alkali lands were converted into severely saline-alkali lands.Generally, land conversion was mainly between moderately and non/slightly saline-alkali lands.Moreover, land conversion among the different classes of saline-alkali lands resulted in an Generally, land conversion was mainly between moderately and non/slightly saline-alkali lands.Moreover, land conversion among the different classes of saline-alkali lands resulted in an overall decline in the area of land under moderately and severely saline-alkali conduction.In other words, land salinization in Huanghua City decreased from 1992 to 2011.As seen in Table 4, 3704.22,562.87 and 3645.99 ha of marshland was converted into non/slightly, moderately and heavily saline-alkali lands, respectively.Meanwhile, some cultivated land under non/lightly, moderately and severely saline-alkali conditions were converted into marshland.Some cultivated lands were also converted into residential lands.Although errors exist in the values reported above, such errors are not large enough to affect the overall dynamics of land use change in the study area.
The changes in the different land use classes were also analyzed using centroid shift for the period from 1992 to 2011.The moderately saline-alkali lands in southeast Huanghua City were converted into non/slightly saline-alkaline lands, and severely saline-alkali lands far from the coast (west of the city) converted into moderately saline-alkaline lands (Figure 7).The changes in the center of gravity for the different classes cultivated lands under saline-alkali condition (Figure 8).For 1992-2011, the gravity center of cultivated land under severely and non/slightly saline-alkali conditions shifted 3.47 km and 9.55 km towards the southeast, respectively.Cultivated lands under moderately saline-alkali conditions shifted 31.89km towards the northwest.The center of gravity of severely and non/slightly saline-alkali land shifted closer towards the coastline, while that of moderately saline-alkali lands shifted from the southwest coastline towards the northwest.overall decline in the area of land under moderately and severely saline-alkali conduction.In other words, land salinization in Huanghua City decreased from 1992 to 2011.As seen in Table 4, 3704.22,562.87 and 3645.99 ha of marshland was converted into non/slightly, moderately and heavily saline-alkali lands, respectively.Meanwhile, some cultivated land under non/lightly, moderately and severely saline-alkali conditions were converted into marshland.Some cultivated lands were also converted into residential lands.Although errors exist in the values reported above, such errors are not large enough to affect the overall dynamics of land use change in the study area.
The changes in the different land use classes were also analyzed using centroid shift for the period from 1992 to 2011.The moderately saline-alkali lands in southeast Huanghua City were converted into non/slightly saline-alkaline lands, and severely saline-alkali lands far from the coast (west of the city) converted into moderately saline-alkaline lands (Figure 7).The changes in the center of gravity for the different classes cultivated lands under saline-alkali condition (Figure 8).For 1992-2011, the gravity center of cultivated land under severely and non/slightly saline-alkali conditions shifted 3.47 km and 9.55 km towards the southeast, respectively.Cultivated lands under moderately saline-alkali conditions shifted 31.89km towards the northwest.The center of gravity of severely and non/slightly saline-alkali land shifted closer towards the coastline, while that of moderately saline-alkali lands shifted from the southwest coastline towards the northwest.Note that: A1 is center of gravity of cultivated lands under non or light saline-alkali conditions in 1992; A2 is centers of gravity of cultivated lands under moderately saline-alkali conditions in 1992; A3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 1992; B1 is centers of gravity of cultivated lands under non or lightly saline-alkali conditions in 2011; B2 is centers of cultivated lands under gravity of moderately saline-alkali conditions in 2011; and B3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 2011.

Method Feasibility
Largely, image acquisition period decides image features.The specific seasons of image acquisition highlight research objectives and procedures.Guan and Liu [23] noted that images acquired in May were critical for distinguishing different classes of land use in the Yellow River Delta.In the Huanghua City study area, the specific season of images acquisition is end of May or early June.During this time, different classes of cultivated land have distinct vegetation characteristic.The high contrast in image features during this time makes the classification process very successful.
In this study, the classification accuracy of the templates exceeded 90%.The exceptions only were marsh and salt pot, which were developed based on the national land survey and the new interpretation symbol system (Tables 3 and 4).By combining the classification templates with the maximum likelihood supervised classification method, information was derived on cultivated lands under saline-alkali conditions in the coastal ecosystem of Huanghua City.This information was used to analyze the area and spatial distribution of cultivated land under saline-alkali conditions in 1992-2011.A total of 11 field sites were identified in the land use map for Huanghua City in 2011.The land use classes of the 11 field sites matched well with interpretation results.
The total area of cultivated land under saline-alkali conditions was 91,100 ha.The proportions of non/slightly (salinity < 0.3%), moderately (0.3% < salinity < 0.5%) and severely (salinity > 0.5%) saline-alkali lands in 2006 were 81.67%, 12.55% and 5.77%, respectively [22].The total area of cultivated land under saline-alkali lands land were 70,987 ha, with proportions of non/slightly, moderately and severely saline-alkali lands in 2011 of 61.9%, 23% and 15.1%, respectively.The order of the proportion of each class of land use was non/slightly > moderately > severely.The accuracy of all the classes was high and therefore reasonable.The results showed that salinization of cultivated lands in the Huanghua city study area dropped during the period from 1992 to 2011 [24].Thus, the new interpretation symbol system and its classification system were suitable for analysis of spatial and temporal variations in saline-alkali lands in coastal regions.The modified remote sensing imagery characteristics were dynamic in real-time and suitable for an accurate monitoring of large areas.The method was also labor, time and cost efficient compared with conventional field measurements [3,25,26].
Vegetation classes as indicators for salinity had strong correlation with electrical conductivity.In this study, climate was not significant factor driving vegetation growth and distribution in the study area.Because this was surprising, it requires further research to more accurately determine the role of climate in the growth and distribution of vegetation in cultivated lands under saline-alkali conditions.

Changes in Cultivated Lands
Soil texture and irrigation water quality influence the soil salinity.Climate, terrain features and groundwater table also cause soil salinization.Then, human activity aggravates soil salinization to different levels [27][28][29].Soil salinization has obvious seasonal characteristics, as dramatic increase in evaporation (due to high spring temperatures) results in salt deposition.Soil desalination occurs in summer due to soil eluviation driven by summer rains.The low rainfall in autumn coupled with low temperatures results in low vaporization and salt deposition.When soils become frozen in wheat fields and groundwater depth steadily drops, the result is low salt accumulation.
Groundwater depth sits at about 0.5-4 m below the land surface, with salinity above 5 g/L in most regions.In the study area, groundwater supplements soil water and soil salt because of its shallow depth and high salinity.Since the 1960s, authorities in Huanghua City have embarked on the construction of drainage systems in cultivated lands to improve saline-alkali condition.The drainage systems consist of networks of ditches, which convey water to and from the field.There is 8 m width ditch in every 50 m wide field.However, the drainage system northwest of Huanghua City has a series of problems including deposition, blocking, falling depth, reduced flow, etc.Most frustrating among the problems is the rapid rate at which the drains have been filled by sedimentation.Drainage and its related problems have favored increase in groundwater level.Decreasing groundwater depth enhances soil salinization.However, the drainage system southeast of Huanghua City has remained intact since 1992.This has meant a decreasing groundwater level and low soil salinization.

Conclusions
In this paper, a modified method was developed for application in interpreting changes in remote sensing images for cultivated land under saline-alkali conditions in coastal regions.The modified method was a new interpretation symbol system for grading and classifying cultivated land under saline-alkali conditions.The system was developed based on site investigation, plant cover characteristics and remote sensing images features.Then, classification templates were established based on the interpretation symbol system and the maximum likelihood supervised classification used to generate land use maps for Huanghua City in Hebei Province China.
The total area of saline-alkali land decreased during the period from 1992 to 2011 in Huanghua City, hence a drop in the degree of salinization of arable land in the city.The spatial distribution of saline-alkali land changed with time and the center of gravity of severely, moderately and non/slightly saline lands changed as well.The center of gravity of severely and non/slightly saline-alkali lands shifted closer towards the coastline, whereas that of moderately saline-alkali lands shifted from the southwest coastline towards the northwest.

Figure 1 .
Figure 1.Map of China depicting the location depicting Hebei Province (top left), that of the province depicting the study area (right) and an expanded map of the study area (bottom left).

Figure 2 .
Figure 2. A schematic representation of the method of modification of the remote sensing imagery used in the study.

Figure 1 .
Figure 1.Map of China depicting the location depicting Hebei Province (top left), that of the province depicting the study area (right) and an expanded map of the study area (bottom left).

Figure 1 .
Figure 1.Map of China depicting the location depicting Hebei Province (top left), that of the province depicting the study area (right) and an expanded map of the study area (bottom left).

Figure 2 .
Figure 2. A schematic representation of the method of modification of the remote sensing imagery used in the study.

Figure 2 .
Figure 2. A schematic representation of the method of modification of the remote sensing imagery used in the study.

Figure 3 .
Figure 3. Photos of different classes of saline-alkali lands taken in the field in 2011.

Figure 3 .
Figure 3. Photos of different classes of saline-alkali lands taken in the field in 2011.

Figure 3 .
Figure 3. Photos of different classes of saline-alkali lands taken in the field in 2011.

Table 1 .
Classes of saline-alkali lands in Huanghua City and the interpretation symbols of TM (standard false color composite).Classes of Saline-Alkali Land Vegetation Feature Image Feature Severely saline-alkali land Salt-resistant vegetation Moderately saline-alkali land Spring corn, cotton Non/slightly saline-alkali land Winter wheat Moderately saline-alkali land Spring corn, cotton

Figure 3 .
Figure 3. Photos of different classes of saline-alkali lands taken in the field in 2011.

Table 1 .Figure 3 .
Figure 3. Photos of different classes of saline-alkali lands taken in the field in 2011.

Figure 4 .
Figure 4. Land use maps for the period 1992 and 2011 in Huanghua City, Hebei Province.
8 ha by analyzing the land use classification map of 1992.They are about consistent.We choose 11 sites in the land use classification map of 2011 and did field investigation through Global Position System (GPS) in 2011.The land use classes of the 11 field sites matched well with interpretation results.

Figure 4 .
Figure 4. Land use maps for the period 1992 and 2011 in Huanghua City, Hebei Province.

3. 5 .
Validation of Remote Sensing Interpretation We use statistical yearbook and previous research to verify accuracy of the land use classification maps of 1992 and field survey to verify accuracy of the land use classification maps of 2011.According to statistical yearbook and previous research [22], the area of cultivated land was about 130,000 ha.The area of cultivated land is 139,588.8ha by analyzing the land use classification map of 1992.They are about consistent.We choose 11 sites in the land use classification map of 2011 and did field investigation through Global Position System (GPS) in 2011.The land use classes of the 11 field sites matched well with interpretation results.
, there were changes in the other classes of land use in 1992-2011 in the study area.Residential land area was increased by 27.86% from 25,521.6 ha to 32,631.7 ha during 1992-2011.Aquaculture area shrunk, decreasing from 23,065.5 ha in 1992 to 13,468.5 ha in 2011.The area of salt spot increased from 6512.8 ha to 7684.4 ha.Area under water bodies significantly increased from 4362.8 ha to 17,070.8ha over the 20-years period.Marshland area also decreased from 18,998.2 ha in 1992 to 11,716.6 ha in 2011 in Huanghua City in Hebei Province, China.
, there were changes in the other classes of land use in 1992-2011 in the study area.Residential land area was increased by 27.86% from 25,521.6 ha to 32,631.7 ha during 1992-2011.Aquaculture area shrunk, decreasing from 23,065.5 ha in 1992 to 13,468.5 ha in 2011.The area of salt spot increased from 6512.8 ha to 7684.4 ha.Area under water bodies significantly increased from 4362.8 ha to 17,070.8ha over the 20-years period.Marshland area also decreased from 18,998.2 ha in 1992 to 11,716.6 ha in 2011 in Huanghua City in Hebei Province, China.

Figure 5 .
Figure 5. Ratio of the areas of land under different saline-alkali land classes as percentages of total cultivated areas in 1992 and 2011 in Huanghua City in Hebei Province, China.

Figure 6 .
Figure 6.Temporal changes in the areas of land under different saline-alkali land classes from the period from 1992 to 2011 in Huanghua City in Hebei Province, China.

Figure 5 .
Figure 5. Ratio of the areas of land under different saline-alkali land classes as percentages of total cultivated areas in 1992 and 2011 in Huanghua City in Hebei Province, China.
, there were changes in the other classes of land use in 1992-2011 in the study area.Residential land area was increased by 27.86% from 25,521.6 ha to 32,631.7 ha during 1992-2011.Aquaculture area shrunk, decreasing from 23,065.5 ha in 1992 to 13,468.5 ha in 2011.The area of salt spot increased from 6512.8 ha to 7684.4 ha.Area under water bodies significantly increased from 4362.8 ha to 17,070.8ha over the 20-years period.Marshland area also decreased from 18,998.2 ha in 1992 to 11,716.6 ha in 2011 in Huanghua City in Hebei Province, China.

Figure 5 .
Figure 5. Ratio of the areas of land under different saline-alkali land classes as percentages of total cultivated areas in 1992 and 2011 in Huanghua City in Hebei Province, China.

Figure 6 .
Figure 6.Temporal changes in the areas of land under different saline-alkali land classes from the period from 1992 to 2011 in Huanghua City in Hebei Province, China.

Figure 6 .
Figure 6.Temporal changes in the areas of land under different saline-alkali land classes from the period from 1992 to 2011 in Huanghua City in Hebei Province, China.

Figure 7 .
Figure 7. Spatial changes in land use for the period 1992-2011 in Huanghua City, China.

Figure 7 .
Figure 7. Spatial changes in land use for the period 1992-2011 in Huanghua City, China.

Figure 8 .
Figure 8. Centers of gravity of cultivated lands under saline-alkali conditions in Huanghua City of Hebei Province, China.Note that: A1 is center of gravity of cultivated lands under non or light saline-alkali conditions in 1992; A2 is centers of gravity of cultivated lands under moderately saline-alkali conditions in 1992; A3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 1992; B1 is centers of gravity of cultivated lands under non or lightly saline-alkali conditions in 2011; B2 is centers of cultivated lands under gravity of moderately saline-alkali conditions in 2011; and B3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 2011.

Figure 8 .
Figure 8. Centers of gravity of cultivated lands under saline-alkali conditions in Huanghua City of Hebei Province, China.Note that: A1 is center of gravity of cultivated lands under non or light saline-alkali conditions in 1992; A2 is centers of gravity of cultivated lands under moderately saline-alkali conditions in 1992; A3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 1992; B1 is centers of gravity of cultivated lands under non or lightly saline-alkali conditions in 2011; B2 is centers of cultivated lands under gravity of moderately saline-alkali conditions in 2011; and B3 is centers of gravity of cultivated lands under severely saline-alkali conditions in 2011.

Table 1 .
Classes of saline-alkali lands in Huanghua City and the interpretation symbols of TM (standard false color composite).

Table 2 .
Error matrix for the classifications of land use in 1992 in Huanghu City, China (unit: pixel).

Table 3 .
Error matrix for the classifications of land use in 2011 in Huanghu City, China (unit: pixel).