Characteristics and Mechanism of the Environmental Capacity Changes in Haizhou Bay, Northern Jiangsu, China from 2006 to 2016

: Haizhou Bay is an open bay located in northern Jiangsu Province, China. This study analyzes the changes in the coastline, coastal development, and water quality of Haizhou Bay between 2006 and 2016. The box model method and numerical simulation are adopted to calculate the environmental capacities of Haizhou Bay in 2006 and 2016, analyze changes to environmental capacity features, and assess the inﬂuencing factors over this period. The scenario analysis method is used to discuss the inﬂuencing mechanism and degree of inﬂuence of factors (e.g., the water quality di ﬀ erence inside and outside the bay, and sea reclamation) on the environmental capacity and calculate the contribution of each inﬂuencing factor. The changes in terrain triggered by sea reclamation and water quality from 2006 to 2016 reduced the total environmental capacity of Haizhou Bay, with an inﬂuencing ratio of 0.198:0.802. In other words, poorer water quality inside the bay reduces the environmental capacity by a degree of 4.05 times that of sea reclamation. This study can o ﬀ er guidance on related future research aiming to protect the marine environment of Haizhou Bay and control the total amount of pollutants discharged into the sea.


Introduction
A bay is an area of sea that extends into the land, creating a curved coastline. As a transitional area between land and sea, bays are rich in resources such as marine life, ocean floor minerals, and marine oil or gas, and they are often used as ports for shipping and as tourist destinations because of their picturesque nature. However, coastal development inevitably harms the marine environment; for example, sea reclamation shortens the coastline. This leads to many changes to the coastal terrain, the characteristics of tidal movement, and the migration behavior of pollutants [1]. Sea reclamation also reduces the tidal prism (the tidal volume that a bay can carry) and the ability of seawater to remove pollutants, and it may even affect threatened bird habitats [2]. Moreover, bait and drugs for aquaculture may also pollute seawater. In the development and processing of coastal ports, sanitary sewage, ballast water, and accidental oil spills from ships at sea can degrade the quality of the marine environment [3,4]. In the process of expansion and development, urban civilizations will always have an impact on the surrounding hydrology and river landscape [5]. Land-sourced pollutants relating to upstream human activities flow into the sea through rivers, leading to grave threats to the marine environment of bays.
We analyzed the changes in the coastline, coastal development, and water quality of Haizhou Bay between 2006 and 2016. A numerical simulation was used to calculate the environmental capacities of the bay in 2006 and 2016. The scenario analysis method was employed to discuss the influencing mechanism and degree of several factors (e.g., the water quality difference inside and outside the bay and sea reclamation) on the environmental capacity and calculate the rates of contributions of these influencing factors. The exploration of changes in environmental capacity and quantification of the influencing factors in this study can offer guidance for related future research aiming to protect the marine environment of Haizhou Bay and control the total amount of pollutants discharged into the sea. Correspondingly, this research may promote the balanced and sustainable development of economy and society in the sea area of Haizhou Bay and adjacent land area. The scenario analysis method was employed to discuss the influencing mechanism and degree of several factors (e.g., the water quality difference inside and outside the bay and sea reclamation) on the environmental capacity and calculate the rates of contributions of these influencing factors. The exploration of changes in environmental capacity and quantification of the influencing factors in this study can offer guidance for related future research aiming to protect the marine environment of Haizhou Bay and control the total amount of pollutants discharged into the sea. Correspondingly, this research may promote the balanced and sustainable development of economy and society in the sea area of Haizhou Bay and adjacent land area.

Overview of Haizhou Bay
Lianyungang

Data Sources and Analytical Method
We collected American earth resources satellite (Landsat) images for 2006 and 2016 for an interpretation of remote sensing through ENVI (The Environment for Visualizing Images) to determine the coastline of the study area in the two years, so as to analyze the reclamation from 2006 to 2016. We collected 1:25,000 topographic maps of Haizhou Bay in 2006 and 1:20,000 topographic maps in 2016. According to the topographic maps of these two periods, the topographic erosion and sedimentation changes of Haizhou Bay were compared and analyzed. The ownership data of ocean development and utilization and ocean space planning in 2006 and 2016 were collected from the Ocean Management Department. The marine water quality data of Haizhou Bay in 2006 were from the marine water quality monitoring carried out by the Jiangsu Marine Environment Monitoring and Forecasting Center at Haizhou Bay. The marine water quality data of Haizhou Bay in 2016 were from the environmental quality report of Lianyungang Coastal Waters, which was issued by the Lianyungang Environmental Monitoring Center Station. ArcGIS was used for data processing and integration, information superimposition, and spatial analysis.

Data Sources and Analytical Method
We collected American earth resources satellite (Landsat) images for 2006 and 2016 for an interpretation of remote sensing through ENVI (The Environment for Visualizing Images) to determine the coastline of the study area in the two years, so as to analyze the reclamation from 2006 to 2016. We collected 1:25,000 topographic maps of Haizhou Bay in 2006 and 1:20,000 topographic maps in 2016. According to the topographic maps of these two periods, the topographic erosion and sedimentation changes of Haizhou Bay were compared and analyzed. The ownership data of ocean development and utilization and ocean space planning in 2006 and 2016 were collected from the Ocean Management Department. The marine water quality data of Haizhou Bay in 2006 were from the marine water quality monitoring carried out by the Jiangsu Marine Environment Monitoring and Forecasting Center at Haizhou Bay. The marine water quality data of Haizhou Bay in 2016 were from the environmental quality report of Lianyungang Coastal Waters, which was issued by the Lianyungang Environmental Monitoring Center Station. ArcGIS was used for data processing and integration, information superimposition, and spatial analysis.

The Box Model Method
The main methods to calculate the marine environmental capacity include the share rate method, optimization method, model trial method, box model method, and GIS (Geographic Information System) method. The box model method is widely used in the environmental evaluation of bays and estuaries. In this method, the water body is considered as a homogeneously mixed box with uniformly distributed properties such as temperature and concentration. Any local change in the box leads to global influence immediately, and the box will be homogeneous again at once. On such an assumption, the variable quantity of substances in the box per unit of time can be expressed by the amount of substance coming in minus that coming out plus the produced amount via the internal biochemical reaction [34]. That is, the seawater outside the bay carries pollutants into the bay during the tidal period. At the same time, the water in the bay carries pollutants away from the bay. Add to that the amount of variation in the biochemistry of the pollutants in the bay. We can calculate the total change in the amount of matter in the bay per unit time via the following equation: where V is the water body volume (m 3 ); C is the pollutant concentration of seawater in the box (mg/L); t is the tidal period (days); C in is the mean concentration of the seawater outside the box (mg/L); Q in and Q out represent the water volume flowing into and out of the box per unit time, respectively (m 3 ); β is the exchange fraction of internal seawater with external seawater; γ is the exchange fraction of external seawater with internal seawater; and S denotes the increase in substances produced via internal biochemical reactions per unit of time (kg). The pollutant studied in this paper is the chemical oxygen demand (COD). From a conservative point of view, the chemical and biological processes are not considered in the calculation of COD environmental capacity; instead, only the physical processes of convection and diffusion are considered. The marine environmental capacity can be obtained by a simplified calculation of the box model.
Calculation method: the water body is regarded as an ideal box of isoconcentration, and C 0 is the background concentration value of every internal point. Thus, pollutants will be homogeneously mixed the moment they enter the water body, and the marine environmental capacity can be calculated by where EC is the environmental capacity (tons); v is the water body volume (m 3 ); and C 0 and C S are the background and the water quality standard values (mg/L), respectively [35].

Determination of the Calculation Boundary
The main function of the box calculation model is the determination of the calculation boundary. The ECG calculation result varies according to the scope enclosed by boundaries. Based on the analysis of the Haizhou Bay landform and terrain and their change trends, tidal dynamic field, water flux, and the present coastal development, the boundary that best separates the areas inside and outside the bay was delineated based on the following assumptions. First, the boundary cannot be disturbed by internal factors of the bay, and the submarine topography should be relatively stable. Second, the boundary should be perpendicular or parallel to the direction of the tidal hydrodynamic force here. Third, coastal development should be included.
Through a previous analysis of terrain change and the comparison between the submarine topography in 2006 and 2016, a few changes were found in the terrain outside the 16 m isobath, indicating that the terrain of areas deeper than 16 m are rarely disturbed by these influencing factors. Thus, the calculation boundary should be drawn deeper than the 16 m isobath.
Haizhou Bay is controlled by a rotary tidal wave that moves from north to south but exhibits a rectilinear current along the shore. When the tide rises, the open-side tide flows into Haizhou Bay from northeast to southwest (NE-SW). When the tide recedes, it runs off the bay from southwest to northeast (SW-NE). Except in the flank, there is a sharp incline between the tide and isobath or coastline, indicating that the tide mainly tends to flow back and forth between both sides. In this regard, the calculation boundary should be perpendicular or parallel to the direction of the tidal hydrodynamic force, namely in the NW-SE or NE-SW direction.
The coastal development of Haizhou Bay in 2016 was mainly carried out in the sea area 30 km from the coast and shallower than the 16 m isobath, 91.5% of which is used for aquaculture.
After consideration of factors such as terrain, tidal power, and coastal development, the perpendicular cross-sections were determined as the outer boundary. Cross-sections A and B are in the SE-NW and SW-NE directions, respectively. The boundary and the present coastline jointly define the scope of the ECG calculation ( Figure 2).
Water 2020, 12, x FOR PEER REVIEW 5 of 18 Haizhou Bay is controlled by a rotary tidal wave that moves from north to south but exhibits a rectilinear current along the shore. When the tide rises, the open-side tide flows into Haizhou Bay from northeast to southwest (NE-SW). When the tide recedes, it runs off the bay from southwest to northeast (SW-NE). Except in the flank, there is a sharp incline between the tide and isobath or coastline, indicating that the tide mainly tends to flow back and forth between both sides. In this regard, the calculation boundary should be perpendicular or parallel to the direction of the tidal hydrodynamic force, namely in the NW-SE or NE-SW direction.
The coastal development of Haizhou Bay in 2016 was mainly carried out in the sea area 30 km from the coast and shallower than the 16 m isobath, 91.5% of which is used for aquaculture.
After consideration of factors such as terrain, tidal power, and coastal development, the perpendicular cross-sections were determined as the outer boundary. Cross-sections A and B are in the SE-NW and SW-NE directions, respectively. The boundary and the present coastline jointly define the scope of the ECG calculation ( Figure 2).
The cross-sectional tidal flux is the product of the cross-sectional area and the flow velocity perpendicular to the cross-section. The tidal fluxes of cross-sections A and B were summed via Delft 3D and considered to represent the exchanged water body volume.

Calculation of Marine Environmental Capacity
The environmental capacity of the gulf (ECG) of Haizhou Bay comprises the static environmental capacity (ECS) and dynamic environmental capacity (ECD). The ECS is the difference between the standard and background values of still water quality, and the ECD is the change in tidal prism (water purification ability) from the tide rising and falling: The cross-sectional tidal flux is the product of the cross-sectional area and the flow velocity perpendicular to the cross-section. The tidal fluxes of cross-sections A and B were summed via Delft 3D and considered to represent the exchanged water body volume.

Calculation of Marine Environmental Capacity
The environmental capacity of the gulf (ECG) of Haizhou Bay comprises the static environmental capacity (ECS) and dynamic environmental capacity (ECD). The ECS is the difference between the standard and background values of still water quality, and the ECD is the change in tidal prism (water purification ability) from the tide rising and falling: Static Environmental Capacity (ECS) The ECS shows the pollutant-carrying ability of seawater enclosed by the land and ocean boundaries in the research area. It can be obtained via calculation of the theoretical maximum environmental capacity, the used environmental capacity, and the remaining environmental capacity: where ECS b is the theoretical maximum ECG of Haizhou Bay under a given standard of water quality, and ECS d is the used ECG in the background condition. ECS b is the product of the standard water quality of marine functional zoning and the water body volume, which is calculated as follows: where C b is the standard water quality of marine functional zoning-in this paper, it refers to the second-class seawater quality (mg/L); V x is the water body volume in the box (m 3 ); S is the sea area enclosed by the box boundary (m 2 ); and L is the mean depth of seawater in the box (m). The water quality grade is based on the "Sea Water Quality Standard" (GB3097-1997). This is the National Standard of the People's Republic of China, which contains the standards of various pollutants such as inorganic nitrogen, active phosphate, chemical oxygen demand (COD), and so on. We selected the COD parameter as a tracer of marine pollution in this study. ECS d is the product of the present pollutant concentration in Haizhou Bay and the water body volume in the box, which is calculated as follows: where C x is the present pollutant concentration of Haizhou Bay (mg/L). Then, the remaining environmental capacity (ECS) can be expressed by Equation (4).

Dynamic Environmental Capacity (ECD)
ECD refers to the change in tidal prism from the rising and falling of the tide, namely the change in the exchange capacity of the water body. The exchanged seawater will leave Haizhou Bay through the tidal movements, carrying polluted water to somewhere far away from the land, where the water quality is good. In the next tidal rising, external clean seawater flows into Haizhou Bay. This circulation demonstrates the physical cleaning ability of seawater.
ECD is the product of the present pollutant concentrations inside and outside Haizhou Bay and the exchanged seawater volume, which are calculated as follows: Water 2020, 12, 2990 where C w is the present pollutant concentration outside Haizhou Bay (mg/L), and V n is the mean exchanged water body volume of Haizhou Bay (m 3 ). The basic assumption of theoretical maximum dynamic environmental capacity (ECD x ) is that the water in the bay is already the upper limit of the control water quality standard, while the water outside the bay is clean and does not contain pollutants. The calculation method is the product of the pollutant concentration of the control water quality standard and the net tidal capacity, which is calculated as follows: where C b is the standard water quality of marine functional zoning-in this paper, it refers to the second-class seawater quality (mg/L); and V n is the mean exchanged water body volume of Haizhou Bay (m 3 ).

Numerical Simulation Method
Based on the data of exchanged water body volume, the box model method was adopted to calculate the environmental capacity in accordance with the physical change of pollutants in the water body and indicate the physical diffusion ability of the water body. Delft 3D (numerical simulation software) was employed to create a hydrodynamic force model of the tide and simulate the tidal field of Haizhou Bay [36][37][38][39].

Model Grid and Parameter Setting
The Lianyungang area, spanning 113 km from north to south and 120 km from east to west, is taken as the scope for model calculation in this paper. Within this scope, the northwest and southwest boundaries are the fixed boundaries of Lianyungang. The northeast and southeast boundaries are the open boundaries of the sea area. The control method of boundary conditions is the tidal level control. Orthogonal grids are used to draw the model. There are 387,325 grid cells in the model. In the process of grid drawing, the sea area of Haizhou Bay is encrypted. The smallest grid cell is 82 m × 69 m. The Manning coefficient is 0.02. The water temperature is 5 • C. The annual average wind speed is 5.5 m/s, and the normal wind direction is east (Figure 3). For more details, please refer to [40]. where Cb is the standard water quality of marine functional zoning-in this paper, it refers to the second-class seawater quality (mg/L); and Vn is the mean exchanged water body volume of Haizhou Bay (m 3 ).

Numerical Simulation Method
Based on the data of exchanged water body volume, the box model method was adopted to calculate the environmental capacity in accordance with the physical change of pollutants in the water body and indicate the physical diffusion ability of the water body. Delft 3D (numerical simulation software) was employed to create a hydrodynamic force model of the tide and simulate the tidal field of Haizhou Bay [36][37][38][39].

Model Grid and Parameter Setting
The Lianyungang area, spanning 113 km from north to south and 120 km from east to west, is taken as the scope for model calculation in this paper. Within this scope, the northwest and southwest boundaries are the fixed boundaries of Lianyungang. The northeast and southeast boundaries are the open boundaries of the sea area. The control method of boundary conditions is the tidal level control. Orthogonal grids are used to draw the model. There are 387,325 grid cells in the model. In the process of grid drawing, the sea area of Haizhou Bay is encrypted. The smallest grid cell is 82 m × 69 m. The Manning coefficient is 0.02. The water temperature is 5 °C. The annual average wind speed is 5.5 m/s, and the normal wind direction is east (Figure 3). For more details, please refer to [40].

Model Validation
The research focus of this paper is the change mechanism of Marine environmental capacity in Haizhou Bay from 2006 to 2016.Therefore, the calculation of the model and the validation of the  Compared with 2016, due to the marine development activities in Haizhou Bay and the input of pollutants into the sea, the topography and seawater quality have changed, the tidal volume of Haizhou Bay has changed, and the marine environment capacity will also change.

The Influencing Factors and Changes of Bay Environmental Capacity
Compared with 2016, due to the marine development activities in Haizhou Bay and the input of pollutants into the sea, the topography and seawater quality have changed, the tidal volume of Haizhou Bay has changed, and the marine environment capacity will also change.
We hope to study the impact of topography and water quality on the Marine environmental capacity of Haizhou Bay by means of quantitative calculation. So we assumed two computational boundaries and scenarios.
Scenario After that, comparing the actual marine environmental capacity in 2006 and 2016 with the hypothetical scenario A and scenario B, the impact of topography or water quality on the environmental capacity could be quantified.
The difference between Scenario A and the 2006 environmental capacity reflects the water quality on the environmental capacity; the difference between the 2016 environmental capacity and Scenario A reflects the impact of topographic changes on the environmental capacity.
The comparison shows that the difference between environmental capacity in 2006 and Scenario B reflects the impact of topographic changes on environmental capacity, while the difference between Scenario B and the environmental capacity in 2016 reflects the impact of water quality on the environmental capacity. The changes in isobath show that the sea area above the 16 m isobath was still in the silting state, and the sea area below the 16 m isobath was basically stable. increased by 1.5 times, and the sea area use types increased from two to five. The area of sea reclamation over these 10 years was 34.9 km 2 .

Changes in Natural Conditions
From 2006 to 2016, great changes occurred in the coastline, submarine topography, seawater quality, and coastal development of Haizhou Bay, threatening the environment and potentially strongly affecting the ECG.

ECG Calculation Results for 2006 and 2016
According to the data for Haizhou Bay's terrain in 2006, L was 10.367 m, S was 101.034 million m 2 , and Vx was 10.474 billion m 3 . According to the data in the "Current Conditions and Assessment In 2006, the mean concentrations of COD, labile phosphate, and inorganic nitrogen of Haizhou Bay met the standards of second-class seawater quality (Table 1). In 2016, the mean concentrations of COD and labile phosphate still met the standard, whereas that of inorganic nitrogen exceeded it. The mean concentrations of all three variables were all raised, and that of inorganic nitrogen had increased by 284.5%, signaling an overall tendency of the water quality of Haizhou Bay to worsen.  (Figure 6a). In 2016, 414 km 2 of sea area was developed in the study area: 91.5% for open-type marine culture, 1.6% for sea-enclosure marine culture, 0.1% for tourism and entertainment, 0.9% for communications and transportation, and 5.9% for urban construction (Figure 6b). The area of marine development and utilization in the study area increased by 1.5 times, and the sea area use types increased from two to five. The area of sea reclamation over these 10 years was 34.9 km 2 .

ECG Calculation Results for 2006 and 2016
According to the data for Haizhou Bay's terrain in 2006, L was 10.367 m, S was 101.034 million m 2 , and Vx was 10.474 billion m 3 . According to the data in the "Current Conditions and Assessment  (2005-2010), most of the studied areas were agricultural and fishery areas, which should meet the class II seawater quality standard. According to the "Sea water quality standard (GB 3097-1997)", the concentration limit of COD class II water quality standard is 3 mg/L. Delft 3D software was used to calculate the average water exchange capacity of the study area. The model calculation time was from January 22nd to January 30th, in 2006. The spring tide and neap tide were included in the simulation period, and the daily average water exchange volume of the spring tide and neap tide was obtained. The mean exchanged water body volume of Haizhou Bay (V n ) was 5 (Table 2).
According to the data for Haizhou Bay's terrain in 2016, L was 10.413 m, S was 97.541 million m 2 , and V x was 10.157 billion m 3 . According to the data in the "Report on Offshore Environmental Quality of Lianyungang of Jiangsu Province", 2016, C x and C w were 1.63 and 1.40 mg/L, respectively. According to the Marine functional zoning of Jiangsu Province (2011-2020), most of the studied areas were agricultural and fishery areas, which should meet the class II seawater quality standard. The concentration limit of COD class II water quality standard is 3 mg/L. The mean exchanged water body volume of Haizhou Bay (V n ) was 5

Calculation Results for the Hypothetical Scenarios
Similarly, using the box model and numerical simulation method, we calculate the COD environmental capacity of Scenario A and Scenario B (Table 3).  In 2016, compared with 2006, the static environmental capacity decreased by 2948.44 tons; the theoretical maximum static environmental capacity decreased by 951.75 tons, and the used static environmental capacity increased by 1996.69 tons. ECS d is influenced by seawater quality and water body volume. We assumed that the water body in Haizhou Bay was not polluted, namely a mean COD of 0 mg/L, and that its change was only influenced by sea reclamation projects. Under these circumstances, ECS b decreased by 951.75 tons. The environmental capacity of the bay decreased by 1996.69 tons due to the change of water quality between 2006 and 2016.
In 2016, compared with 2006, the theoretical maximum dynamic environmental capacity decreased by 1082.25 tons, and the dynamic environmental capacity increased by 183.35 tons. ECD is influenced by water quality changes inside and outside of the bay and the change in the net tidal prism. We assumed that the water quality change inside the bay was the poorest scenario of ideal ECD reduction (just meeting the standard of the second-class water quality), while that outside the bay was the best scenario (clean water), which represents the theoretical maximum of water pollution diffusion of unit net tidal prism). The mean COD was 3 mg/L inside the bay and 0 mg/L outside the bay, with a 3 mg/L difference, and the ECG reduction that only considered net tidal prism decline was 10.8225 million tons. However, water quality inside and outside the bay, and especially the difference between them, cannot be ignored because a greater difference leads to more diffuse pollutants in total. The actual changes in water quality inside and outside the bay, their difference, and the coastal change due to sea reclamation increased ECD by 183.35 tons from 2006 to 2016. In other words, although sea reclamation reduced the net tidal prism of the water body exchange by 360.75 million m 3 /d and the theoretical ECD by 10.8225 million tons, the difference of water quality inside and outside the bay cannot be considered negligible. The COD difference inside and outside the bay increased from 0.185 mg/L in 2006 to 0.23 mg/L in 2016, signaling a growth in the ability of the unit water body to exchange and transport pollutants. Hence, the decline in net tidal prism between 2006 and 2016 raised the ECD.
According to this analysis, ECD was influenced by the net tidal prism and the difference of water quality inside and outside the bay. The theoretical ECD merely shows the environment quality change (ideal maximum change) triggered by the changes in water depth, terrain, and the net tidal prism owing to factors such as sea reclamation; thus, it is limited to a theoretical discussion instead of a practical calculation of ECG. However, variations in the water quality difference can control ECG change and increase ECD, even under the overall degradation of water quality. Seawater deterioration may lead to the continuous increase in ECG instead of having no influence, because the degradation decreases ECS b , and the polluted water can enter the bay through water body exchange. In the long term, the water body quality outside the bay will worsen, diminishing the difference of water quality inside and outside the bay and further reducing the ECD.

ECG-Influencing Factors and Ratios
To further discuss the influence of changes in terrain and water quality on ECG in different factor combinations, their influencing ratios were analyzed in combinations #1 (terrain combined with water quality) and #2 (water quality combined with terrain). Based on different boundary conditions, how the changes in terrain caused by sea reclamation, water quality because of water quality variation, and differences in water quality inside and outside the bay influence the environmental capacity were examined according to #1 and #2 (Table 4, Figure 7). Considering that the theoretical maximum static environmental capacity and theoretical maximum dynamic environmental capacity are only affected by the change in net tidal capacity, and that there is a big gap between the theoretical maximum static environmental capacity and the actual environmental capacity, this part mainly discusses the topographic change caused by reclamation: the water quality change between inside and outside the bay on the static residual environmental capacity and dynamic environmental capacity; meanwhile, the environmental capacity of the bay is discussed according to paths 1 and 2.   The ratios of terrain to water quality on ECG reduction, ECSb decline, and ECD change were 0.198:0.802, 0.16:0.84, and −0.408:1.408, respectively. As a result, the rates of contributions of the two factors on ECG reduction, ECSb decline, and ECD change were 19.8% and 80.2%, 16% and 84%, and −40.8% and 140.8%, respectively. Poorer water quality led to an ECG reduction 4.05 times that of sea reclamation, ECSb decline 5.25 times that of sea reclamation, and ECD 3.45 times that of sea reclamation in absolute value, since a larger water quality difference increases ECG, whereas sea reclamation decreases ECG by diminishing the net tidal prism.

Conclusions
In view of the influencing factors of the sea area capacity of Haizhou Bay, in order to effectively protect the marine environment of Haizhou Bay, it is suggested that research be carried out on the control line of regional reclamation, with attention paid to the superimposed influence of reclamation on the marine environmental capacity and control of the reclamation scale. It is necessary to strengthen the management and control of the marine environment in Haizhou Bay, including the supervision of marine development activities and sewage discharge of rivers entering the sea. The total amount of pollutants discharged into Haizhou Bay should be controlled, and discharge standards should be strictly implemented to meet the discharge standards. We should strengthen the supervision of pollutant discharge from marine development activities. We should take necessary environmental protection measures to prohibit the direct discharge of pollutants in Haizhou Bay. We should control the pollution of mariculture wastewater, control the scale of aquaculture, and adopt ecological breeding methods to reduce the production and discharge of aquaculture wastewater. In order to improve the marine environment of Haizhou Bay and promote the social and economic sustainable development of Haizhou Bay, measures such as the control of the marine environment and the total amount of pollutants discharged from the sea must be taken.
The main innovation of this study is to expand the spatial and temporal boundaries of marine environmental capacity research. The marine environmental capacity of Haizhou Bay in 2006 and The difference between Scenario B and the actual environmental capacity in 2006 reflects the impact of topographical changes on environmental capacity. The results show that the topographic change caused by reclamation resulted in a total reduction of 577.51 tons in the sea area compared with 2006, including 510.77 tons of static residual environmental capacity and 66.74 tons of dynamic environmental capacity. This shows that the reduction of the bay area and the weakening of water exchange capacity caused by reclamation activities lead to the reduction of static and dynamic environmental capacity. The difference between Scenario A and the actual environmental capacity in 2006 reflects the impact of water quality changes on environmental capacity. The calculation results show that the static residual environmental capacity decreased by 2513.8 tons due to the deterioration of water quality in the bay, while the dynamic environmental capacity increased by 266.32 tons per day due to the increase of water quality difference between inside and outside the bay. However, in terms of the total environmental capacity of the bay per day, there was still a decrease of 2247.48 tons, which shows that the deterioration of water environment quality was still the main reason for the decrease in environmental capacity.
The ECG changes of #1 and #2 conformed to the actual ECG changes from 2006 to 2016, indicating that both #1 and #2 can be adopted to precisely simulate the joint effect of changes in terrain and water quality on ECG and thus can be used for analysis.
In #1, the influencing ratios of terrain to water quality on ECG reduction, ECS b decline, and ECD change were 0.209:0.791, 0.173:0.827, and −0.364:1.364, respectively. In #2, the influencing ratios of terrain to water quality on ECG reduction, ECS b decline, and ECD change were 0.187:0.813, 0.147:0.853, and −0.453:1.453, respectively. This means that a calculation based on terrain changes (#1) will overstate the impact of terrain changes, while a calculation based on water quality changes (#2) will overstate the impact of changes in the water quality. Therefore, to comprehensively consider both factors, the mean of #1 and #2 was adopted in this calculation, and the ratio of terrain to water quality on ECG reduction, ECS b decline, and ECD change were 0.198:0.802, 0.16:0.84, and −0.408:1.408, respectively. As a result, the poorer water quality led to an ECG reduction 4.05 times that of sea reclamation, an ECS b decline by 5.25 times, and an ECD decline by 3.45 times that of sea reclamation in terms of absolute value.
In this study, only physical diffusion was considered, while biological and chemical factors were ignored. Thus, the true ECG is likely to be higher. Future studies should also include the combined influence of physical diffusion and biochemical decomposition on the environmental capacity. The ratios of terrain to water quality on ECG reduction, ECS b decline, and ECD change were 0.198:0.802, 0.16:0.84, and −0.408:1.408, respectively. As a result, the rates of contributions of the two factors on ECG reduction, ECS b decline, and ECD change were 19.8% and 80.2%, 16% and 84%, and −40.8% and 140.8%, respectively. Poorer water quality led to an ECG reduction 4.05 times that of sea reclamation, ECS b decline 5.25 times that of sea reclamation, and ECD 3.45 times that of sea reclamation in absolute value, since a larger water quality difference increases ECG, whereas sea reclamation decreases ECG by diminishing the net tidal prism.

Conclusions
In view of the influencing factors of the sea area capacity of Haizhou Bay, in order to effectively protect the marine environment of Haizhou Bay, it is suggested that research be carried out on the control line of regional reclamation, with attention paid to the superimposed influence of reclamation on the marine environmental capacity and control of the reclamation scale. It is necessary to strengthen the management and control of the marine environment in Haizhou Bay, including the supervision of marine development activities and sewage discharge of rivers entering the sea. The total amount of pollutants discharged into Haizhou Bay should be controlled, and discharge standards should be strictly implemented to meet the discharge standards. We should strengthen the supervision of pollutant discharge from marine development activities. We should take necessary environmental protection measures to prohibit the direct discharge of pollutants in Haizhou Bay. We should control the pollution of mariculture wastewater, control the scale of aquaculture, and adopt ecological breeding methods to reduce the production and discharge of aquaculture wastewater. In order to improve the marine environment of Haizhou Bay and promote the social and economic sustainable development of Haizhou Bay, measures such as the control of the marine environment and the total amount of pollutants discharged from the sea must be taken.
The main innovation of this study is to expand the spatial and temporal boundaries of marine environmental capacity research. The marine environmental capacity of Haizhou Bay in 2006 and 2016 is calculated by a box model. The influencing factors of marine environmental capacity are analyzed, and the contribution rate of various influencing factors is calculated. This study explores the change and mechanism of marine environmental capacity and can provide a reference for other similar studies.
Author Contributions: L.S. and J.W. carried out the numerical simulation and cartography and analyzed the data; H.Z. provided support to the model setup and validation; L.S., J.W., H.Z., and M.X. discussed the results; the paper was written and edited by L.S. and J.W. the paper edited by L.S. and J.W. All authors have read and agreed to the published version of the manuscript.
Funding: This research was funded by a Jiangsu Provincial Marine Science and Technology Innovation Project, grant number HY2018-3.

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