Mapping Long-Term Spatiotemporal Dynamics of Pen Aquaculture in a Shallow Lake: Less Aquaculture Coming along Better Water Quality

: Pen aquaculture is the main form of aquaculture in some shallow lakes in eastern China. It is valuable to map the spatiotemporal changes of pen aquaculture in eutrophic lakes to assess its e ﬀ ect on water quality, thereby helping the relevant decision-making agencies to manage the water quality (WQ) of lakes. In this study, an automatic approach for extracting the pen aquaculture area was developed based on Landsat data. The approach integrates ﬁve algorithms, including grey transformation, discrete wavelet transform, fast Fourier transform, singular value decomposition and k-nearest neighbor classiﬁcation. It was successfully applied in the automatic mapping of the pen aquaculture areas in Lake Yangcheng from 1990 to 2016. The overall accuracies were greater than 92%. The result indicted that the practice of pen aquaculture experienced ﬁve stages, with the general area increasing in the beginning and decreasing by the end of the last stage. Meanwhile, the changes of nine WQ parameters observed from 2000 to 2016, such as ammonia (NH 3 -N), pH, total nitrogen (TN), total phosphorus (TP), chlorophyll a , biochemical oxygen demand (BOD), chemiluminescence detection of permanganate index (CODMn), Secchi disk depth (SDD) and dissolved oxygen (DO), were analyzed in the lake sectors of Lake Yangcheng, and then their relationships were explored with the percentage of pen aquaculture area. The result suggested that the percentage of pen aquaculture area exhibits signiﬁcantly positive correlations with NH 3 -N, TN, TP, chlorophyll a , BOD and CODMn, but signiﬁcantly negative correlations with SDD and DO. The experimental results may o ﬀ er an important implication for managing similar shallow lakes with pen aquaculture expansion and water pollution problems.


Introduction
China is the largest inland fishery producer in the world with a proportion of 54.5% of the global inland fishery production in 2012 (FAO, 2014). Lake aquaculture is an important food source and one of the most popular industries in China [1]. In China, the area of lake aquaculture accounts for 39.4% of inland aquaculture area based on a 5-year average area from 2006 to 2010 [2]. Pen aquaculture is the main form of aquaculture in shallow lakes [2]. It shifted from a traditional aquaculture with a natural-food based extensive culture in the 1980s and 1990s to a modern aquaculture with a high-density, Therefore, the main goal of this study is to develop an automatic approach for mapping time series pen aquaculture areas based on Landsat data, and then to assess their effect on water quality by combining the long-term and site-specific water quality observation data from Lake Yangcheng. The corresponding specific objectives include (i) developing an automatic and robust approach to map pen aquaculture areas in Lake Yangcheng from 1992 to 2016; (ii) exploring correlations between the percentage of pen aquaculture area and water quality parameters; and (iii) assessing the effectiveness of government measures implemented at various stages for controlling pen aquaculture activities to improve water quality.

Study Site
Lake Yangcheng is located in the northeast of the Lake Taihu basin, near the city of Suzhou, Jiangsu Province, China (31 • 21 -31 • 30 N and 120 • 39 -120 • 51 E) ( Figure 1). It is the third largest freshwater shallow lake in the Lake Taihu basin with a water surface area of approximately 108 km 2 and an average depth of 2.1 m. Two nearly parallel narrow strips of land divide the lake into three sectors with limited interconnection channels, named the Western Lake, Central Lake and Eastern Lake, with a water surface area of approximately 32.5 km 2 , 33.6 km 2 and 42.1 km 2 , respectively. In Lake Yangcheng, there are many inflowing and outflowing rivers ( Figure 1). According to multi-year observations, the inflowing water is generally from north to south and from west to east; its retention time is 0.38 year [38].
Remote Sens. 2020, 12, x FOR PEER REVIEW 3 of 21 by combining the long-term and site-specific water quality observation data from Lake Yangcheng. The corresponding specific objectives include (i) developing an automatic and robust approach to map pen aquaculture areas in Lake Yangcheng from 1992 to 2016; (ii) exploring correlations between the percentage of pen aquaculture area and water quality parameters; and (iii) assessing the effectiveness of government measures implemented at various stages for controlling pen aquaculture activities to improve water quality.

Study Site
Lake Yangcheng is located in the northeast of the Lake Taihu basin, near the city of Suzhou, Jiangsu Province, China (31°21'-31°30' N and 120°39'-120°51' E) ( Figure 1). It is the third largest freshwater shallow lake in the Lake Taihu basin with a water surface area of approximately 108 km 2 and an average depth of 2.1 m. Two nearly parallel narrow strips of land divide the lake into three sectors with limited interconnection channels, named the Western Lake, Central Lake and Eastern Lake, with a water surface area of approximately 32.5 km 2 , 33.6 km 2 and 42.1 km 2 , respectively. In Lake Yangcheng, there are many inflowing and outflowing rivers ( Figure 1). According to multi-year observations, the inflowing water is generally from north to south and from west to east; its retention time is 0.38 year [38].  Lake Yangcheng is a representative lake with pen aquaculture of CMC in China. Each pen in the lake, made of 2-3 cm stretched mesh polythene netting tied to a bamboo framework, covers an area of about 1.3 ha with a length of about 130 m and a width of 100 m. The bamboo poles (8-10 cm in diameter and 4-5 m in length) of the framework are vertically inserted 60-100 cm into the bottom of the lake at intervals of 1.5-2.0 m [16]. According to field surveys and inquiries, CMC is the dominant cultured species in the pens due to its high demand and profit. The crabs are fed with artificial pelleted feeds supplemented with vegetable feeds (e.g., corn, wheat and pachyrhizus) and feeds of animal origin (e.g., trash fish) [11].

Landsat Image Data
Given the recognizable enclosure size, data consistency and traceability, we used Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and Operational Land Imager (OLI) with a long time-series record and a spatial resolution of 30 m to extract the pen aquaculture area and monitor its spatiotemporal change. In this study, a total of 26 scenes of cloud-free Landsat TM/ETM+/OLI images acquired from 1991 to 2016 were collected from the U.S. Geological Survey (USGS, https://earthexplorer.usgs.gov) ( Table A1). The collected images were acquired only between April and May because of the following reasons: (i) Since image spectral features of the pen facility are similar to the aquatic vegetation, especially floating-leaved vegetation, the image acquisition chosen in April or May can avoid such a confusion of spectral features between the pen facility and aquatic vegetation [39]. (ii) Crab seeds are transferred to aquaculture pens between mid-February and mid-April and are harvested at the end of October every year. Some pen facilities are removed or broken from November to the next February [10]. (iii) Some algae particles and budding floating-leaved vegetation float and aggregate around the pen facilities between April and May, which makes the spectral feature of the pen more significant from its background, and thus the information about the pen facility is easier to be extracted than that in the other months ( Figure 1b). In addition, since there were no available Landsat data in April and May in 2012 and 2014 (Table A1), Huangjing Charge Coupled Device (HJ CCD) images were used as the supplementary data to map the pen aquaculture in the two years. HJ-CCD images were acquired from the China Centre for Resources Satellite Data and Application (CRESDA). HJ-1A/1B CCD has a similar spectral band setting and spatial resolution to those of the first four bands of Landsat TM. Previous studies indicated that the data had a good consistency with Landsat data in application [40,41].
In order to develop an automatic method for extracting the pen aquaculture area, we chose a Landsat OLI image acquired on 22 May 2015, as a test image, and a high-resolution image called ZiYuan-3 (ZY-3) acquired on 28 April 2015, as a reference image. The ZY-3 satellite, launched in January 2012, is China's first civilian high-resolution optical satellite. It provides four multispectral (MS) bands with a spatial resolution of 5.8 m and a revisit cycle of 5 days. We visually interpreted the pen aquaculture area from the reference image to validate the result of the pen aquaculture area mapped from the test image.

Water Quality Data
We collected two datasets of water quality (WQ) data, including nine WQ parameters, namely, dissolved oxygen (DO), biochemical oxygen demand (BOD), pH, ammonia (NH 3 -N), total nitrogen (TN), total phosphorus (TP), chlorophyll a (Chla), Secchi disk depth (SDD) and chemiluminescence detection of permanganate index (COD Mn ). We defined the two datasets as Dataset A and Dataset B. Dataset A was provided by the Suzhou water environment monitoring department. Monthly WQ parameters were measured from seven regular stations from 2000 to 2016, except for 2001 and 2002. The seven regular stations, including three stations in Western Lake, two stations in Central Lake and two in Eastern Lake, were set according to the spatial variability of the WQ in the three lakes ( Figure 1). The average WQ at these stations in every lake could be roughly representative of the overall WQ state of the corresponding lakes. Dataset B was collected from 27 sampling sites in 2016, including eight sites in Western Lake, seven sites in Central Lake and 12 sites in Eastern Lake ( Figure 1). All the WQ parameters from the two datasets were measured using standard methods. SDD, pH and DO were directly measured in the field, while TN, NH 3 -N, TP, COD Mn and Chla were measured in the laboratory from samples collected from stations/sample sites. TP was determined by using ammonium molybdate spectrophotometry (GB11893-89), whereas TN was determined via alkaline potassium persulfate digestion and UV spectrophotometry (GB11894-89). CODMn was determined using the acidic potassium permanganate method (GB/T11892-89). NH 3 -N was determined using Nessler's reagent method (GB 7479-87), spectrophotometry and the molybdenum blue technique (GB 11893-89). In this study, based on Dataset A, we averaged the monthly WQ observations from Stations 1, 2 and 3 in Western Lake, Stations 4 and 5 in Central Lake, and Stations 6 and 7 in Eastern Lake, and then averaged them across 12 months to obtain an annual average of each lake sector for the correlation analysis with their corresponding pen aquaculture area. For the TN and TP observations from the seven stations in 2016, the spatial and seasonal variations in WQ parameters varied from the different stations ( Figure A1). Generally, there was the highest TN and TP concentration in Western Lake (Stations 1, 2 and 3), followed by Central Lake (Stations 4 and 5), and the lowest in Eastern Lake (Stations 6 and 7).
In addition, to explore the nutrient source in the three lake sectors, in 2016, monthly discharges of the inflow and outflow rivers surrounding the lakes were monitored in situ with an acoustic Doppler current profiler (Sontek, RiverSurveyor M9, SonTek, San Diego, CA, USA), and the corresponding TN and TP concentrations were also measured using standard lab measurement methods described above. Finally, the seasonal fluxes of TN and TP were calculated by multiplying the seasonal discharge and corresponding concentrations of TN and TP.

Feature Image Selection for Extracting Pen Aquaculture
In the Landsat test image, we first selected a sample line, including pixels of the pen aquaculture facility, water and aquatic vegetation, and then plotted the pixels' values of the blue, green, NIR and SWIR bands. Pixels of the pen facility have the greatest spectral difference from other objects in the NIR band (Figure 2), and thus the NIR image might be selected as the sensitive feature image from which to extract the pen aquaculture information.
Remote Sens. 2020, 12, x FOR PEER REVIEW 5 of 21 ammonium molybdate spectrophotometry (GB11893-89), whereas TN was determined via alkaline potassium persulfate digestion and UV spectrophotometry (GB11894-89). CODMn was determined using the acidic potassium permanganate method (GB/T11892-89). NH3-N was determined using Nessler's reagent method (GB 7479-87), spectrophotometry and the molybdenum blue technique (GB 11893-89). In this study, based on Dataset A, we averaged the monthly WQ observations from Stations 1, 2 and 3 in Western Lake, Stations 4 and 5 in Central Lake, and Stations 6 and 7 in Eastern Lake, and then averaged them across 12 months to obtain an annual average of each lake sector for the correlation analysis with their corresponding pen aquaculture area. For the TN and TP observations from the seven stations in 2016, the spatial and seasonal variations in WQ parameters varied from the different stations ( Figure A1). Generally, there was the highest TN and TP concentration in Western Lake (Stations 1, 2 and 3), followed by Central Lake (Stations 4 and 5), and the lowest in Eastern Lake (Stations 6 and 7).
In addition, to explore the nutrient source in the three lake sectors, in 2016, monthly discharges of the inflow and outflow rivers surrounding the lakes were monitored in situ with an acoustic Doppler current profiler (Sontek, RiverSurveyor M9, SonTek, San Diego, CA, USA), and the corresponding TN and TP concentrations were also measured using standard lab measurement methods described above. Finally, the seasonal fluxes of TN and TP were calculated by multiplying the seasonal discharge and corresponding concentrations of TN and TP.

Feature Image Selection for Extracting Pen Aquaculture
In the Landsat test image, we first selected a sample line, including pixels of the pen aquaculture facility, water and aquatic vegetation, and then plotted the pixels' values of the blue, green, NIR and SWIR bands. Pixels of the pen facility have the greatest spectral difference from other objects in the NIR band (Figure 2), and thus the NIR image might be selected as the sensitive feature image from which to extract the pen aquaculture information.

Extraction of Pen Aquaculture Information
The pen aquaculture area was extracted based on four steps, including image enhancement, image segmentation, feature extraction and target classification. Specifically, the procedures consist of five algorithms: grey transformation (GT), discrete wavelet transform (DWT), fast Fourier transform (FFT), singular value decomposition (SVD) and k-nearest neighbor (kNN) classification.
All algorithms were written and executed using MATLAB software, and some functions, including fft (), wavedec () and svd (), were directly called in MATLAB when conducting the DWT, FFT and SVD algorithms.

Extraction of Pen Aquaculture Information
The pen aquaculture area was extracted based on four steps, including image enhancement, image segmentation, feature extraction and target classification. Specifically, the procedures consist of five algorithms: grey transformation (GT), discrete wavelet transform (DWT), fast Fourier transform (FFT), singular value decomposition (SVD) and k-nearest neighbor (kNN) classification.
All algorithms were written and executed using MATLAB software, and some functions, including fft (), wavedec () and svd (), were directly called in MATLAB when conducting the DWT, FFT and SVD algorithms.
Detailed algorithms and steps are described as follows: 1. Image enhancement. Image enhancement is a process of adjusting the digital images to make a target (i.e., pen aquaculture area) easier to be identified. There are two steps, including image normalization and exponential GT, to enhance the image and highlight the pen facility [42]. The formula is as follows: where N(x) is the pixel value on a normalization image; F(x) is the pixel value on a feature image; F(m) is the maximum pixel value in the feature image; and G(r) is the transformed pixel value of the corresponding normalized image. 2.
Image segmentation. It was critical to automatically acquire an appropriate threshold for segmenting an image and identifying a target. Considering noise in the image after GT processing, it was relatively difficult to obtain an appropriate threshold for segmenting an image. In this study, a DWT was applied to remove noise and then to determine the corresponding thresholds [43,44]. The DWT transform can be expressed as follows: where f (x) is the histogram of the grey-level transformation image; a is the scaling factor determining the width of the wavelet; b is the shifting factor indicating the position of a mother wavelet; and ψ (a,b) (x) represents the wavelets that are transformed by scaling and shifting the mother wavelet ψ(x). In this study, a 0 = 2, b 0 = 1 and the mother wavelet ψ(x) was db4.
In this study, after the wavelet was reconstructed, a threshold for segmenting water and other objects could be identified and determined by finding out the first turning point to the left of the highest peak. Then the turning point was the threshold.

3.
Feature extraction. The pen facility has a similar spectral feature with aquatic vegetation, but it shows a regular shape and periodic spatial arrangement. In practice, an FFT proved to be an effective tool in extracting the time and spatial change frequency of the targets [20]. Therefore, in this study, we used a 2-dimensional discrete FFT (FFT2) to exact the features of the water, aquatic vegetation and different angles of a pen facility. The FFT2 formula is shown as follows: where F(u, v) is the frequency spectrum of the original image f (x, y). In this study, f (x, y) is a segmented image.
After FFT2, an SVD was executed to reduce the dimensionality and extract the features of different objects. The SVD provides a convenient way for breaking a matrix, which perhaps contains some data in which we are interested into simpler, meaningful pieces [45]. It is widely used in data compression, noise reduction, etc. It can be used to detect and extract weak signals from noisy data. For an SVD, we begin with an arbitrary really valued m × n matrix X.
where S is the m × n diagonal matrix; and V T is the n × n matrix. The columns of U and V are called the left and right singular vectors. In this study, X is the frequency-domain matrix derived from FFT2, with m = 15 and n = 15.

4.
Target identification. A kNN rule, one of the classic and top-performing classifiers, was used to identify the pen aquaculture area. It achieves a classification by calculating the similarity between the test sample (pixel) and all the training samples based on a discrimination function [46]. The similarity can be measured by Euclidean distances [47,48]. The discrimination function is as follows: where x i is the vector of class i in training set; x is the vector to be unclassified; and N is the number of classes. In this study, there existed water, aquatic vegetation and the pen grids with different angles, which can be generalized as 45 • , 75 • and 90 • .

Validation of Extracting Result
Given the lack of actual area data of pen aquaculture, we used the aquaculture distribution area created by visually interpreting the high resolution ZY-3 images as a reference to assess the extracted result from the Landsat test image. We defined the classification map created with the automatic algorithms as Map A and the visual interpretation map as Map B. By overlay analysis, we could get the area in agreement (i.e., the area shown in both Maps A and B), the omission area (i.e., the aquaculture area identified as non-aquaculture area) and the commission area (i.e., the non-aquaculture area identified as aquaculture area) to conduct a spatial accuracy assessment. Further, the overall accuracy as well as the errors of commission and omission could be calculated using the following formulas: Overall

Extraction of Pen Aquaculture Area and Validation
Pixel grey values of the normalized blue, green, red, NIR, SWIR1 and SWRI2 bands in a sample line, including the typical pen, aquatic vegetation and water pixels ( Figure 2). It was significant that only the NIR band exhibited significant and periodic peaks of the pen facility compared with all the other five bands. Thus, NIR could be selected as a feature band to extract the pen area.
The NIR image, enhancement image, segmented image and extracted resultant map of the pen aquaculture area are shown in Figure 3. Obviously, the pen grid feature was enhanced and was identified more easily (Figure 3b). However, some noise, including cloud and aquatic vegetation, were also enhanced (e.g., noise within the red circles in Figure 3b). A threshold of 231 was obtained by the DWT and wavelet reconstruction to segment the image ( Figure A2). The segmentation image indicated that some noise, such as light cloud (e.g., some noises within the red circles in Figure 3b), were removed by the wavelet transformation, but some noise, such as aquatic vegetation and light cloud (e.g., some noise within the red circles in Figure 3c), were not removed. Finally, the distribution map of pen aquaculture in Lake Yangcheng after FFT ( Figure A3) and kNN were obtained (Figure 3d). Obviously, due to FFT2 processing, some vegetation regions without pen facilities but with similar spectral feature as them had been excluded (e.g., vegetated areas within the red circles in Figure 3c).  (1) and (2), segmented image based on Equation (3) (c) and extracted resultant map of pen aquaculture based on equation (6) (d). The white areas within the red circles are noise (e.g., aquatic vegetation and light cloud).
The extracted result was transferred into a shapefile ( Figure 4a) and compared with the visual interpretation result (Figure 4b) in GIS mapping software. Spatially, most of the omission and commission areas were located along edges of the pen aquaculture. Most of misclassified areas were actually aquatic vegetation with a high coverage (e.g., the area within the yellow circle in Figure 4c). Overall accuracy of the exacted result was 92.95% with a commission error of 16.72% and an omission The extracted result was transferred into a shapefile ( Figure 4a) and compared with the visual interpretation result (Figure 4b) in GIS mapping software. Spatially, most of the omission and commission areas were located along edges of the pen aquaculture. Most of misclassified areas were actually aquatic vegetation with a high coverage (e.g., the area within the yellow circle in Figure 4c). Overall accuracy of the exacted result was 92.95% with a commission error of 16.72% and an omission error of 7.05% ( Table 1). The relative difference is −2.46~4.18% between the monitoring area of pen aquaculture in 2000, 2001, 2002 and 2015 using the proposed approach and the area recorded in the corresponding reference years ( Table 2).
Remote Sens. 2020, 12, x FOR PEER REVIEW 9 of 21  and are the omission area (i.e., the aquaculture area identified as nonaquaculture) and commission area (i.e., the non-aquaculture area identified as aquaculture area), respectively; is the area in agreement in both Maps A and B. and are the omission error and commission error, respectively.

Spatiotemporal Changes in Pen Aquaculture from 1992 to 2016
The spatiotemporal changes in pen aquaculture were different among the three lake sectors from 1992 to 2016 (Figures 5 and 6). In Eastern Lake, the percentage of pen aquaculture area to the whole lake area reached 90% in 1993, decreased to 55.80% in 1994, then increased to 93.54%, and finally decreased to 40.38% in 2016. In Central Lake, the percentage increased from 57.42% in 1993 to 89.95% in 2014, and then decreased to 37.03% in 2016. In Western Lake, it increased from 46.08% in 1993 to 87.64% in 1998, and then kept unchanged until 2008 when it was removed completely. In general, the change process for the entire lake can be categorized into five stages. Stage I was before 1993 when there was not any pen facility there. Stage II was from 1992 to 2002. Pen aquaculture started in 1993 with a percentage of 66.81%, and then increased to 87.71% in 2002. Stage III was between 2002 and  Overall accuracy = 92.95%; E C = 16.72%; E O = 7.05%.
Note: S A and S B are the area by the approach proposed and the visual interpretation method, respectively; S O and S C are the omission area (i.e., the aquaculture area identified as non-aquaculture) and commission area (i.e., the non-aquaculture area identified as aquaculture area), respectively; S x is the area in agreement in both Maps A and B. E O and E C are the omission error and commission error, respectively.

Spatiotemporal Changes in Pen Aquaculture from 1992 to 2016
The spatiotemporal changes in pen aquaculture were different among the three lake sectors from 1992 to 2016 (Figures 5 and 6). In Eastern Lake, the percentage of pen aquaculture area to the whole lake area reached 90% in 1993, decreased to 55.80% in 1994, then increased to 93.54%, and finally decreased to 40.38% in 2016. In Central Lake, the percentage increased from 57.42% in 1993 to 89.95% in 2014, and then decreased to 37.03% in 2016. In Western Lake, it increased from 46.08% in 1993 to 87.64% in 1998, and then kept unchanged until 2008 when it was removed completely. In general, the change process for the entire lake can be categorized into five stages. Stage I was before 1993 when there was not any pen facility there.

Long-Term Trends in Water Quality and Correlations with the Percentage of Pen Aquaculture
From 2000 to 2016, WQ was improved in the three lake sectors with decreasing NH3-N, TN, Chla, BOD, CODMn and TP and increasing SDD and DO, especially after 2017 (Figure 7). The trends were more significant in the Central and Eastern Lakes than Western Lake, and more obvious at Stage

Long-Term Trends in Water Quality and Correlations with the Percentage of Pen Aquaculture
From 2000 to 2016, WQ was improved in the three lake sectors with decreasing NH 3 -N, TN, Chla, BOD, COD Mn and TP and increasing SDD and DO, especially after 2017 (Figure 7). The trends were more significant in the Central and Eastern Lakes than Western Lake, and more obvious at Stage IV. WQ was worst in Western Lake, and it was better in Central Lake than Eastern Lake.

Advantages and Uncertainty of the Proposed Approach
Satellite remote sensing has been proved to be an effective tool to map aquaculture distribution at a large scale [19,26,28,29,31,52]. Sensitive image features of the aquaculture facilities (e.g., cages and pens) and classification methods are key in mapping an aquaculture area. In this study, the NIR band was proved to be more sensitive to pen facilities than other bands (Figure 2), and thus effective in mapping the pen facilities. This is because of the amount of floating-leave vegetation and algae particles that aggregate around the pen facilities (due to wind and waves) (Figure 1b) during the image acquisition time (i.e., between April and May), which make the pen show a significant vegetation spectral feature (i.e., the high reflection in the NIR band and strong absorption in the Red band). In addition, although there exist a few methods for mapping aquaculture facilities (e.g., cages and pens), they have disadvantages and limitations. For example, an object-oriented image analysis is a very commonly used method for extracting aquaculture facilities due to their regular shapes and special locations [33][34][35][36][37]53,54]. However, the determination of thresholds of the feature parameters relies on interested image features and manual checking. The threshold determined for one time period or a specific image could be inadequate when applied to other dates or images due to different influences, including water properties, atmospheric conditions, sun-view geometry, etc. [55]. Therefore, when using time-series Landsat data to extract the change information of pen aquaculture,

Advantages and Uncertainty of the Proposed Approach
Satellite remote sensing has been proved to be an effective tool to map aquaculture distribution at a large scale [19,26,28,29,31,49]. Sensitive image features of the aquaculture facilities (e.g., cages and pens) and classification methods are key in mapping an aquaculture area. In this study, the NIR band was proved to be more sensitive to pen facilities than other bands (Figure 2), and thus effective in mapping the pen facilities. This is because of the amount of floating-leave vegetation and algae particles that aggregate around the pen facilities (due to wind and waves) (Figure 1b) during the image acquisition time (i.e., between April and May), which make the pen show a significant vegetation spectral feature (i.e., the high reflection in the NIR band and strong absorption in the Red band). In addition, although there exist a few methods for mapping aquaculture facilities (e.g., cages and pens), they have disadvantages and limitations. For example, an object-oriented image analysis is a very commonly used method for extracting aquaculture facilities due to their regular shapes and special locations [33][34][35][36][37]50,51]. However, the determination of thresholds of the feature parameters relies on interested image features and manual checking. The threshold determined for one time period or a specific image could be inadequate when applied to other dates or images due to different influences, including water properties, atmospheric conditions, sun-view geometry, etc. [52]. Therefore, when using time-series Landsat data to extract the change information of pen aquaculture, the existing methods are inefficient because of the low degree of automation. Compared with the existing methods, our approach could offer an effective solution for finding out the appropriate thresholds for different images acquired in different dates to automatically map pen aquaculture, and thus it has a large potential in tracking the spatiotemporal changes of pen aquaculture area using time-series satellite images without manual intervention. The proposed approach included (i) automatically determining the thresholds for segmenting water and other objects by using DWT; (ii) extracting the identification features of the water, aquatic vegetation and pen facility according to their intrinsic attributes and shape characteristics using FFT and SVD; and (iii) automatically mapping pen aquaculture area by using the kNN rule based on these identification features.
When executing the approach to other lakes or images, the following few points and limitations should be noted: (i) the approach is only suitable for mapping pen aquaculture area using images acquired from April to May (see the detailed reasons in Section 2.2.1); (ii) for images with different spatial resolutions, some parameters for the algorithms need to be adjusted and modified, such as an appropriate matrix size involving the number of pixels containing a regular pen input in the algorithms of FFT2 and SVD; (iii) the algorithms included in our approach are not unique, and others with the same functions might be more appropriate due to different images and sites condition. For example, the GT method for highlighting a pen facility in visual can be replaced by other image enhancement methods, such as a gradient transform [53]. A Gaussian fitting also may be used to replace WT to remove the noise and find a threshold [54].

Main Factors Affecting Water Quality in Lake Yangcheng
Generally, excess nitrogen (N) and phosphorus (P) loading from the point and non-point sources in a watershed are considered as one of the main factors damaging the ecological quality of lakes [38,40,55]. In Lake Yangcheng, besides the external nutrient loadings such as industrial effluents and agriculture non-point sources, intensification through increasing use of pelleted feeds and the expansion of pen aquaculture area was another main factor inducing lake eutrophication and deteriorating of WQ [2]. The amount of food fed to the crabs depends on the growth stages of the crab, water quality, temperature and so on. Approximately 0.8 kg of food per acre is put into the pen at the beginning of crab growth, and it will be increased at 10% every three days. About 1/3 of the feed nutrients are consumed by fish and the CMC, with the remaining being released into the immediate environment [5,8]. The excessive pelleted feed and aquaculture wastes, including residual solid feed, fish (or crab) feces and soluble excretions enriching nitrogen and phosphorus, could result in or aggravate lake eutrophication [15,17].
Exchange of nutrients between a lake and watershed is performed through inflowing and outflowing rivers. In Lake Yangcheng, the inflowing rivers are distributed in the west and north lines, and outflowing rivers in the south and east lines ( Figure A4a). According to the field surveys and measurements in 2016, about 4765.12 tons TN and 331.08 tons TP were inflowed into Western Lake from the west line, and the loading amounts of TN and TP were about 142.24 ton/km 2 and 9.88 ton/km 2 , respectively. In addition, there was no pen aquaculture in the lake. Consequently, in Western Lake, the poor WQ and high nutrients might be mainly caused by nutrient loadings from the corresponding watershed, which might explain why Western Lake had the worst WQ although the pens had been removed since 2008 (Figure 7). Water quality of the Western Lake has been improved since 2008, even though it was still the worst among the three lake sectors (Figures 5 and 7). Besides removing the pen aquaculture in the lake, in the corresponding watershed, a series of other measures have been also taken in order to reduce the nutrient loading into the lake, such as removing or closing hundreds of the heavy polluting enterprises and restaurants along with the lake, shrinking the aquaculture pond area, intensifying and improving sewage treatment, etc. These measures made contributions to nutrient reductions in the inflow and water quality improvement of the lakes. In Central and Eastern Lakes, about 640.17 tons TN and 56.091 tons TP inflowed into the lakes from the north line, and the loading amounts of TN and TP were about 8.45 ton/km 2 and 0.91 ton/km 2 , respectively. They only account for 1/17 of the loading amounts of TN and 1/10 of the TP in Western Lake. However, given its retention time of only 0.38 year and the interconnection of the three lake sectors, it is difficult to quantify and separate nutrient contributions from pen aquaculture and from the surrounding area in the two lakes based on the limitations of the existing data we collected. However, pen aquaculture activities certainly could have a negative effect on water quality in Lake Yangcheng because the TN and TP inside the pens were always higher than those outside the pens in both lake sectors ( Figure A5). In the future, based on multi-source data, including survey data of watershed pollution sources, detailed nutrient measurements and simulations inside and outside lakes, as well as the pen areas being mapped by our approach, etc., it is possible to develop an integrated model for quantitatively estimating the contribution rate of each source.

Effects of Pen Aqauculture Govermance
The pen aquaculture area in Lake Yangcheng experienced five change stages, as illustrated in Figure 7. Change trends of TP and TN showed a similar trend with those of pen aquaculture area percent in Central and Eastern Lakes (Figure 9). Their changes were closely related to local government measures from Stages III to V. At Stage I, there was no pen facility in the lakes and traditional aquaculture with natural food, such as grass carp and aquatic vegetation, and the corresponding water quality state was not clear due to the lack of corresponding WQ observation data. At Stage II, pen aquaculture started and then fast expanded without any measure and reached its peak (Figure 9). At this stage, because of the high production from pen aquaculture, the economic profitability was paid too much attention instead of considering the deterioration of the WQ. Consequently, in 2000, the mean TN level reached the surface WQ criteria of Class V (TN ≥ 2.0 mg/L) and the TP level was almost Class IV (TP ≥ 0.1). At Stage III, the WQ parameters were at the worst levels, and the TP and TN levels kept at Class IV and Class V, respectively. The deterioration of the WQ induced a considerable concern by the Suzhou Fisheries Bureau, and they proposed that pen aquaculture in Lake Yangcheng should be shrunken to 40% of the water surface area, and it was approved by the Standing Committee of the People's Congress of Jiangsu Province in 2004. From Figures 6 and 9, we see that pen aquaculture experienced only a slight decrease in the Central and Eastern Lakes from 2004 to 2007. Stage IV was located at an important transition of pen aquaculture area. Due to a major pollution event, an algae bloom outbreak in Lake Taihu in 2007, the Jiangsu Province Water Resources Department and Ocean and Fisheries Bureau jointly issued a proposal about controlling water pollution in the Taihu basin involving Lake Yangcheng and undertook unprecedented efforts to compulsively remove pen aquaculture in some lakes in the Taihu basin. For example, the extent of pen culture area in East Taihu Lake had been reduced from 11,267 ha to 3000 ha in 2009. Obviously, those efforts made a significant effect on improving the WQ (Figures 7 and 9). The total percentage of pen aquaculture in the Central and Eastern Lakes declined fast, from 80% in 2007 to 32% in 2010, and all pens in Western Lake were removed in 2007 due to the worst WQ (Figures 6 and 7). Consequently, the TN and TP decreased significantly but were still at the Class IV and V levels (Figure 9). At Stage V, the relevant governments continued to appeal that the pen aquaculture area was to be reduced to make a way for the health requirement of a good water ecology in Lake Yangcheng. Meanwhile, some measures that could control pollution and improve the WQ were also taken, such as controlling the amount of feeds, establishing ecological dams and limiting the pollutant discharge from surrounding watersheds [1]. As a result, the WQ was further improved from Class IV to Class III for TP and from Class V to Class IV for TN.
Therefore, local government measures were important driving factors of pen aquaculture decrease from Stage III to Stage V. Meanwhile, the TN and TP levels had significantly similar change patterns as the percentage of pen areas in the Central and Eastern Lakes from Stage III to Stage V ( Figure 9). As shown in Figures 8 and A5, the government measures regarding removing and controlling pen aquaculture had an important positive effect on water quality in the Central and Eastern Lakes. Therefore, local government measures were important driving factors of pen aquaculture decrease from Stage III to Stage V. Meanwhile, the TN and TP levels had significantly similar change patterns as the percentage of pen areas in the Central and Eastern Lakes from Stage III to Stage V ( Figure 9). As shown in Figure 8 and Figure A5, the government measures regarding removing and controlling pen aquaculture had an important positive effect on water quality in the Central and Eastern Lakes.

Conclusions
In this study, an approach for automatically extracting pen aquaculture was first developed, and the accuracy of mapping pen aquaculture was higher than 92%. Then, the spatial and area changes of pen aquaculture in Lake Yangcheng were mapped from 1992 to 2016. Through this study, several conclusions could be derived from our experimental results as follows: The analysis results offer practical guidance for the local government to manage pen aquaculture in Lake Yangcheng. Meanwhile, they could also present an important implication for managing pen

Conclusions
In this study, an approach for automatically extracting pen aquaculture was first developed, and the accuracy of mapping pen aquaculture was higher than 92%. Then, the spatial and area changes of pen aquaculture in Lake Yangcheng were mapped from 1992 to 2016. Through this study, several conclusions could be derived from our experimental results as follows: The analysis results offer practical guidance for the local government to manage pen aquaculture in Lake Yangcheng. Meanwhile, they could also present an important implication for managing pen aquaculture and promoting WQ for similar shallow lakes that have pen aquaculture expansion and water pollution problems, especially for those in eastern China.  Acknowledgments: Acknowledgement for the data support from Lake-Watershed Sub Center, National Earth System Science Data Center, National Science & Technology Infrastructure of China (http://lake.geodata.cn).

Conflicts of Interest:
The authors declare no conflict of interest.    Figure A2. Histograms of grey level image (a) and wavelet reconstruction (b). The threshold was 231 according to Figure A2b, and all pixels with a value greater than 231 were counted as water in the DWT image. Figure A3. One-dimensional projections of FFT2 for 45°, 75° and 90° pens, aquatic vegetation and water. Note that the 45°, 75° and 90° enclosure represent the angles between pen facility and horizontal line. A pen facility area had distinct and periodical features compared with water and aquatic vegetation. Figure A2. Histograms of grey level image (a) and wavelet reconstruction (b). The threshold was 231 according to Figure A2b, and all pixels with a value greater than 231 were counted as water in the DWT image. Figure A2. Histograms of grey level image (a) and wavelet reconstruction (b). The threshold was 231 according to Figure A2b, and all pixels with a value greater than 231 were counted as water in the DWT image. Figure A3. One-dimensional projections of FFT2 for 45°, 75° and 90° pens, aquatic vegetation and water. Note that the 45°, 75° and 90° enclosure represent the angles between pen facility and horizontal line. A pen facility area had distinct and periodical features compared with water and aquatic vegetation.

Conflicts of Interest:
The authors declare no conflict of interest.