The Impact of Pollution on China Marine Fishery Culture: An Econometric Analysis of Heterogeneous Growth

The damage caused by pollution is not equal across individuals within a given population, due to the differences in pollution level and coping ability. Identifying the different impacts of environmental pollution on agricultural economic growth and the source of the pollution is of great reference significance for realizing the organic integration of ecological protection and relative poverty alleviation. This study takes marine fishery culture as an example, with the cultured products data of 53 coastal cities in China from 1994 to 2018 and empirically captures the heterogeneous effects of environmental pollution on mariculture growth by constructing a theoretical analysis framework of heterogeneity in pollution damage. Results show that the deepening of marine environmental pollution will reduce the growth rate of marine cultured output, emerging significant regional differences. Nevertheless, under the background of the complex environmental situation currently, further damages mainly derive from the difference in the ability to deal with environmental pollution rather than the difference in pollution degree. In terms of different wastewater sources, the impact of industrial wastewater discharge on marine aquaculture production is not very significant. In contrast, the adverse shocks of urban domestic sewage are more prominent. Notably, the environmental pollution has a significant cumulative effect, which is not reversed. Therefore, a pollution shock in one period may leave the region permanently further behind. Based on this, in the follow-up process of pollution control, we should fully consider the heterogeneity in coping ability of individuals. While promoting industrial wastewater emission reduction, we should emphasize strengthening of the urban domestic sewage treatment strategies and give overall consideration to the ecological environment treatment and regional coordinated development.

The damage caused by pollution is not equal across individuals within a given population, due to the differences in pollution level and coping ability. Identifying the different impacts of environmental pollution on agricultural economic growth and the source of the pollution is of great reference significance for realizing the organic integration of ecological protection and relative poverty alleviation. This study takes marine fishery culture as an example, with the cultured products data of 53 coastal cities in China from 1994 to 2018 and empirically captures the heterogeneous effects of environmental pollution on mariculture growth by constructing a theoretical analysis framework of heterogeneity in pollution damage. Results show that the deepening of marine environmental pollution will reduce the growth rate of marine cultured output, emerging significant regional differences. Nevertheless, under the background of the complex environmental situation currently, further damages mainly derive from the difference in the ability to deal with environmental pollution rather than the difference in pollution degree. In terms of different wastewater sources, the impact of industrial wastewater discharge on marine aquaculture production is not very significant. In contrast, the adverse shocks of urban domestic sewage are more prominent. Notably, the environmental pollution has a significant cumulative effect, which is not reversed. Therefore, a pollution shock in one period may leave the region permanently further behind. Based on this, in the follow-up process of pollution control, we should fully consider the heterogeneity in coping ability of individuals. While promoting industrial wastewater emission reduction, we should emphasize strengthening of the urban domestic sewage treatment strategies and give overall consideration to the ecological environment treatment and regional coordinated development.
Keywords: environmental pollution, marginal damage amount, vulnerability, marine fishery culture, heterogeneous growth HIGHLIGHTS -Costs of pollution are unlikely to be evenly distributed across individuals. -Different marginal damages stem from differences in vulnerability.
-Urban sewage is more destructive than industrial wastewater. -Effect of pollution is cumulative and not reversed.

INTRODUCTION
Since the reform and opening-up policy, while the industry and agriculture enjoy rapid development, it is bringing about severe environmental problems. The aggravation of environmental pollution has not only restricted the high-quality development of the Chinese economy but also caused the loss of value of life of the Chinese residents at a trillion scale (Xie et al., 2019;Khan and Ozturk, 2020). Meanwhile, agriculture restricted by vulnerability and reliance on the natural environment has been subject to various constraints and negative influence of environmental pollution. The worsening of the eco-environment can not only impair the output of agricultural products, but also pose a serious threat to the quality and output of agricultural products and even the income of the farmers (Carter et al., 2017). In the field of natural sciences, a string of researches have already provided evidences for the influence of chemical-caused air pollution on agricultural production. The deepening of nitric oxide, ozone, and haze will impair the photosynthesis and disease resistance of agricultural products (Field et al., 2012;Huang et al., 2015;Rai, 2016). It is predicted that, by the year 2030, increase in the ozone concentration will have caused a drop of 12.1-16.4% in the global soybean output, a drop of 4.3-6.5% in the global corn output, and a drop of 10.6-15.6% in the global wheat output (Avnery et al., 2011). By resorting to the framework of economics and the econometric model, economic circles have also substantiated the negative impact of ozone pollution on grain output. They have also predicted that the losses caused by ozone pollution could be as high as 18 billion USD at a global scale, of which 6.45 billion USD would be from China (Wang and Mauzerall, 2004;Powell and Reinhard, 2016;Yi et al., 2016;YI et al., 2020).
As we all know, marine fishery plays an important role in the economic development of various countries, which is also affected by a set of external factors, such as climate change and environmental pollution (Kassa et al., 2018;Guo, 2020;Ameyaw et al., 2021;Kundu and Santhanam, 2021). The existing studies have shown that while environmental pollution affects fishery yield (Hai and Speelman, 2020), it also indirectly brings uncertainty to the growth of total factor productivity. As the saying goes, "Lucid waters and lush mountains are invaluable assets." With the growing attention paid by all walks of life toward environmental pollution, relevant subjects, such as economics and ecology, have vigorously examined the impact of environmental pollution on economic activities, particularly agricultural production, by an unprecedented depth and breadth (Shao, 2020;Shao et al., 2021). What's important, poverty-stricken rural areas in China are usually mired by environmental pollution, resource dependence, inequality, etc. (Shuai et al., 2019). Those short of capitals usually excessively rely on free natural resources and lack the ability to cope with environmental changes, thus aggravating environmental degradation (Ruch and Geyer, 2017). In additions, existing studies have almost unanimously agreed that individuals will take a serious of adaptive behaviors (e.g., diversification of production and adjustment of factors input) to diversify risks, in the face of changes in the external environment (Bezabih et al., 2010;Krüger et al., 2017;Maire et al., 2021).
However, there is a worth concern that the existing research studies may suffer from two sources of question. The first concerns that there is little observable research focusing on the effects of pollutions on marine production, especially fishery culture. On the one hand, against the background of consumption upgrading, the marine fishery culture as an important source of high quality proteins is a basic link for the Chinese nation to realize its generation-held dream of becoming a superpower in the field of marine development; on the other hand, the ecological status of marine is not optimistic, and a tremendous pressure of sewage is apparent (Bergland et al., 2020;Peng et al., 2021). These factors underscore the need for more detailed analyses and rigorous strategies to estimate the potential effects of pollution. What's important, the changes in environment may generate uneven distribution of costs across individuals within a given population (Hsiang et al., 2019). For example, agriculture usually takes up a high percentage in a the industrial development of a country and it lacks the ability to efficiently cope with changes, so impoverished countries are more vulnerable to global warming. As the new stage began, we will confront tougher challenges to address relative poverty. Therefore, at the heart of another question is that to capture the heterogeneity in the shocks of pollution and juggle mitigation environmental care with poverty alleviation.
To fill that gap, by taking marine fishery culture as an example, this study attempts to assess heterogeneity in hazards of environmental pollution and causes of that. It is hoped that theoretical bases and references can be provided to fuel the green, high-quality development of fishery production. Compared with the existing studies, the contributions to this strand of the literature can be summarized as follows: First, by recognizing heterogeneity in hazards of environmental pollution and causes, this study further enriches the research perspectives of the impact of environmental pollution on fishery culture production; second, for the purpose of increasing the universality and representativeness of the research findings, this study adopts sample data of 53 coastal cities over a long span of time for an empirical analysis.

Variation of Fishery Cultured Output and Cause of Difference From the Perspective of Environmental Pollution
There has been a consensus both in the field of natural sciences and social sciences on hazards of environmental pollution to the agricultural output. For example, the increase in nitric oxide, ozone, and haze in the air will impair the photosynthesis and disease resistance of crops, thus exerting a negative impact on the grain output. In cases where the maximum benefit is as a production goal, producers will adjust the allocation of production elements according to the current environmental status and technical conditions. In order to clarify the impact of environmental pollution on, this study builds the following analysis logic chart. In Figures 1A,B, the internal correlation of the output quantity, Q, and environmental pollution degree, w, with fishery cultured output losses, D, is presented, respectively. When the external environmental pollution reaches the degree of W L , then the optimal output for fishermen to maximize their profits should be Q 1 . In the short term, when the external environmental pollution keeps on worsening with the element price and consumer preference unchanged, namely when the pollution degree increases from W L to W H , its negative impact on the total output will be further deepened (as shown in Figure 1B), thus changing the correlation between the input and output of production elements, and reducing the marginal output value of elements (Saddique et al., 2020). Meanwhile, with the growing output losses, the marginal cost curve of products also increases from MC PL to MC PH . At the moment, fishermen will definitely reduce their fishery cultured output so as to maximize their profits, and thus lead to a decline of production ( Q=Q 1 − Q 2 ), which is equal to the change of damages ( D = D H − D L ) in Figure 1B.
Of special note is that the negative impact of environmental pollution on fishery cultured output is not linear. The negative impact might be varied with different degrees of pollution and different individual adaptability. First of all, the marginal hazards caused by environmental pollution toward seawater cultured output fishery output usually demonstrates an exponential distribution (Burke et al., 2015). Specifically, at the initial stage of pollution, increase in pollutants does not have an obvious impact on output, and the marginal hazards of the former to the latter usually turn out to be zero. Nevertheless, with the worsening of pollution, the marginal hazards will gradually increase. Therefore, to fishermen in different regions (with environmental pollution), worsening of environmental pollution might cause or even widen the gap of output among different regions. On the other hand, concerning individual adaptability toward environmental pollution, differences in individual producers can change the level of exposure toward pollutants and resistance against pollution. On that basis, the marginal hazards of pollutants will also be changed. For example, though increase in the ozone concentration can affect photosynthesis of crops, fishermen can not only change the ozone concentration of the crop environment through the establishment of greenhouses but also alleviate the impact of the ozone concentration on crop output by changing the input of elements, such as nitrogen fertilizers, in accordance with the crop growth. Furthermore, from the perspective of fishery cultivation, fishermen can change the aquatic environment by inputting algae, chemical purifiers, etc. This can help alleviate the hazards of water pollution. Therefore, as shown in Figure 1C, the marginal hazard curve of fisherman m with a stronger resistance and adaptability toward pollution is obviously lower than that of the fisherman n. At the moment, even if in cases where a lower level of external environmental pollution (W L ), due to heterogeneous adaptability, the differences in damages of environmental pollution across individuals who own the same initial fishery are apparent and thus will lead to an output gap D Lm -D Ln between fisherman n and fisherman m, as shown in Figure 1C.

Theoretical Model
With rich marine natural resources, while the marine fishery enjoys rapid development in China, it is reaching fullness. In cases where the maximum benefit is as a production goal, based on the theoretical model developed by Bond et al. (2010), and combining the impact path and mechanism of environmental pollution on fishery cultured output, this study builds the following production function model: (1) In Equation (1), Y denotes the total output; K denotes various production elements; A denotes the productivity of elements; P denotes the degree of environmental pollution; β and γ denote the impact of environmental pollution on the average fishery cultured output and productivity; and g i denotes the average output growth of fishery output in i region without the impact of environmental pollution. After taking logarithm of Equation (1), the cumulative effect of environmental pollution and capital accumulation on output should be included in Equation (1) to obtain the following dynamic growth model: In Equation (3), y it and k it denote the logarithm of the fishery cultured output and element input. Under the condition that the external market environment remains unchanged, the fishery economic growth rate, y it , of t period in i region can be obtained through the difference of Equation (3).
In order to facilitate the explanation and exploration of the impact of environmental pollution on fishery cultured output, it is assumed that output maintains at a steady status, namely, y it−j = y i . ρ is defined as the estimation coefficient of environmental pollution toward output; then, the impact of environmental pollution can be given by Equation (4).
If only the impact of the current environmental pollution is considered, that is, A it = g i + γ 0 P it and y it = β 0 P it + A it + k it + ε it , then Equation (5) can be simplified as below: In Equation (6), y it stands for economic growth. One can easily observe from Equation (6) that, when environmental pollution reaches a normal level, its average impact on fishery cultured output will reverse. For example, the aggravation of environmental pollution might cause a drop in output, but once environment returns to its average value, fishery cultured yields bounce back. Comparatively, the impact of environmental pollution on productivity will not reverse along with the recovery of environment. This means that environmental pollution in one period leaves the fishery culture economy permanently further behind.

Model Construction and Description of Variables
Based on the above theoretical model, this study constructs the following econometric regression model to estimate the impact of environmental pollution on fishery cultured output (Carroll and Weil, 1994;Bond et al., 2001): where, θ i denotes the regional fixed effect; θ rt denotes the time fixed effect, including the dummy variables and differences across regions; ε it denotes the error item; P denotes the environmental pollution of the lag item during L period; and ρ denotes the impact of environmental pollution on output growth. As major crops in China, field crops, such as wheat and paddy rice, are faced with not only a worsening atmospheric environment beset by increasing concentration of zone and haze, but also the aggravation of environmental pollution caused by the enrichment of waters and soil organic matters exceeding acceptable levels. In addition, the pollution degree, Poll totall , can be written as the total sewage discharge of the t year in the i region. It includes the industrial sewage discharge (Poll industry ) and the urban domestic sewage discharge (Poll living ). Different from the complex and various polluted environments facing field crops, environmental pollution confronting production of marine fishery mainly happens to waters. Besides, industrial sewage and domestic sewage of urban areas constitute primary sources of water pollution (Vennemo et al., 2009). Thereby, this study resorts to marine fishery production data of the coastal areas of China, and conducts an empirical analysis on the heterogeneity of hazards caused by environmental pollution. The explained variable g it can be written as the net output growth of seawater cultured products of the t year in i region.
In recognizing causes of the heterogeneity of environmental pollution, this study proceeds from the perspective of relative vulnerability and environmental pollution to construct relevant variables. The vulnerability of fishery culture depends on the sensitivity toward environmental pollution and adaptability measures adopted by humans. Compared with the former, the latter is by no means more referential for the design and formulation of relevant policies. A research study has revealed that worsening of environmental pollution has been a longlasting trend. With a high income and a large number of assets, fishermen can effectively alleviate the negative impact of external factors on production through pre-event and post-event defensive investment. Comparatively, impoverished fishermen with a lower income are more vulnerable toward external impact (Dell et al., 2012). On that basis, this study adopts the per capita disposable income of fishermen as a measure of vulnerability. The relatively vulnerable area is defined by dummy variables. In other words, areas where the per capita disposable income of current fishermen is lower than the average is defined as relatively vulnerable (poor). Assume that poor = 1. On the contrary, it is 0. Similarly, concerning the degree of environmental pollution, this study defines the area where the current total sewage discharge is higher than the average of samples as seriously polluted (worse). Assume that worse = 1. On the contrary, it is 0.
Based on the existing literatures, this study further constructs the relevant control variable matrix, Z, to cover the regional development characteristics, technological level, and circulation of agricultural products. As an important starting point for raising the level of education and promoting the progress of technology, the number of educational personnel is directly related to the improvement of productivity. Meanwhile, the auspice of government is an external force to be reckoned with in the technological advancements. In addition, as the income of the people keeps increasing, there has been an increasing requirement of the protein intake. Thus, with modern conditions of transportation, the growing demand for fresh seawater aquatic products could be satisfied in inner areas, providing an economic motivation of production. According to the above analysis, in order to control the social and impacts on marine fishery culture, this study employs following indicators: local public budget expenditure, expenditure for science and technology, educational fiscal expenditure, the number of educational personnel involved in regular higher education and research institutions, the local average annual population, freight traffic by highways, and area of paved roads, etc.
According to above theoretical analysis, to adopt Equation 6 for a quantitative analysis should meet the prerequisite that α i = 0. Therefore, this study further examines the linear correlation between the random disturbance item, ε it , in Equation (5) and the lag period, y it−1 , of growth, and in order to choose the prolonged lag period as an instrumental variable of growth, y it−1 , for a robustness test. Besides, in order to further identify the heterogeneity in hazards caused by pollution to fishery cultured output and causes of the heterogeneity, this study introduces the interaction item of environmental pollution, vulnerability (poor), and pollution degree (worse). Thereby, the model to estimate the heterogeneity in the impact of environmental pollution on output growth can be obtained, respectively: According to the theoretical model Equations 4, 5, the correlation between output growth during sample period should be verified before using the measurement model for estimation. This can help clarify whether the lag of output growth should be included into the regression model. Estimation results of the correlation coefficient are listed in Table 1. As shown in Table 1 1 , a weak correlation between the aquatic product output over different periods of time in different regions is not observed, which is different from the theoretical hypothesis above. Specifically, the current output growth and one lag are significantly and negatively correlated. This finding probably has a close bearing on the marine fishery production. As renewable biological resources, the stock of marine fishery resources is subject to not only the influence of biological factors, but also human production activities. This means that the resource stock after human exploitation can directly react upon the availability of future resources (Emery et al., 2017). Consequently, a significantly negative correlation is formed between the current output growth and the growth of Lag Phase 1. In light of the above conclusions, this study includes one lag of output growth into the panel estimation model of Equation (6).
Concerning potential problems, such as endogeneity caused by the first-order lag item of the explained item, this paper estimates the dynamic panel model Equations 7, 8 using the generalized method of moments (GMM). Compared with the differential GMM, the system GMM adopts the lags of horizontal value as the instrumental variable of the differential function. At the same time, the differential GMM and the horizontal GMM are combined. The lags of the differential variable are used as the instrumental variable of the horizontal equation. On that basis, the differential equation and the horizontal equation are adopted as an equation system for estimation. Therefore, the system GMM is more efficient in estimation. Moreover, the SE obtained by the model estimation can significantly address problems, such as weak instrumental variables of the horizontal lag item in the panel data model. Considering that, this study chooses the panel data of 53 coastal cities from 1994 to 2018, and uses the system GMM to estimate the impact of environmental pollution on growth of marine cultured output.

Data Sources
As a main source of high-quality proteins in China, the yield of seawater cultured products reached 20.65 million tons in 2019, accounting for about 31.8% of the total aquatic products. Meanwhile, the proportion of seawater fishing products has decreased to about 18.7%. The cultured products have become the type of seawater aquatic products with the largest cultured output in China. It is important to note that little adaptation targeted at offsetting the shock of pollution is observable, although the seawater fishing industry is facing increasing environmental pressures. Therefore, with growing proportion of seawater fishing products, choosing marine fishery culture as a sample for the study can reflect to a certain extent the impact of pollution on the marine fishery. According to the production, acreage, and fisher household income levels of the different main producing areas in 2019, we picked 53 prefecture-level cities of 11 provinces 2 , autonomous regions, and municipalities excluding the Hong Kong, Macau, and Taiwan, as our sampling regions, by using system sampling and K-means cluster. Data selected span from 1991 to 2018. In 2018, the output of seawater cultured products of prefecture-level city samples takes up 95.7% of the national total, which suggests that the sample regions selected are representative enough. Data about the output of seawater cultured products are collected from the "China's Rural Statistical Yearbook" and statistical yearbooks of each cities and the provinces where they belong. Variables, such as the industrial sewage discharge, local average annual pollution, local public budget, spending on science and technology, spending on education, highway freight volume, and actual urban road area, are from "China's Marine Statistical Yearbook, " "China's Marine Yearbook, " and "China's Urban Statistical Yearbook." Discharge of urban domestic sewage is measured by the total domestic water of residents. Regarding other missing data, this study replaces them with the average of sample data at the early and late stage.
Descriptive statistics of variables required by this study are presented in Table 2. As one observes in output across regions can also provide evidence for this point. At the same time, the industrial sewage discharge, domestic sewage discharge, and per capita disposable income of fishermen demonstrate obvious regional and temporal differences across coastal regions. The aforesaid differences can not only suggest a high representativeness of samples but also lay a solid data basis for the follow-up empirical analysis.

Stationary Test of Panel Data
In order to ensure the accuracy of model regression and avoid the appearance of "spurious regression, " this study carries out stationary test of all variables. The Levin Lin and Chu T statistical magnitude (LLC), which are used to judge the existence of the common unit root, and the Fisher ADF statistical magnitude (ADF), which is used to judge the existence of the individual unit root, are both adopted. The lag order is judged by the Akaike Information Criterion (AIC) and Schwarz Criterion (SC). The results of the test are shown in Table 3 below. Results suggest that all variables have denied the existence of the unit root and eliminated the possibility of "spurious regression."

Analysis of Benchmark Estimation Results
The estimation results of the impact of environmental pollution on output growth are shown in Table 4. Model (1) demonstrates the impact of environmental pollution with the total sewage discharge as the measure; Model (2) further clarifies the impact of environmental pollution on output growth of aquatic products. The P-value of AR (1) and AR (2) test 3 in Table 4 suggests that the residual series of the differential equation is found with the  Table 4 show that the core explaining variable, the total sewage discharge (Poll totall ), has a significantly negative impact on the output growth of aquatic products. This to some extent, suggests that the aggravation of environmental pollution might hinder the development of the regional marine fishery culture. After considering different types of pollution sources, this study finds that the domestic sewage discharge (Poll living ) also has a significantly negative impact on the output growth of aquatic products, and that the negative impact is more obvious than the negative impact of the total sewage discharge. Of special note is that the industrial discharge (Poll industry ) has a positive impact on the output growth of aquatic products, but the finding fails to pass the significance test. As stated above, the pollution can impair the nutrition and quality of waters, thus restricting the growth of fishes. At the same time, increase in harmful material content, such as water nitrogen and phosphorus, will not only cause variations of fish diseases but also increase the probability of infection among fish. This will significantly restrict the breeding and reproduction of fish. As to other control variables, the total population has a significantly positive impact on the fishery aquatic cultivation growth (Hasan et al., 2019). series correlation test (AR) and Sargan test on estimation results. The Arellano-Bond series correlation test makes the null hypothesis that the residual sequence of the differential equation is not found with the second-order series correlation.
As important sources of high-quality proteins, fishery cultivation plays an indispensable role in the daily diets of humans. As the income of the people keeps increasing, there has been an increasing requirement of the protein intake. This explains why the growing population will drive the growth of the fishery cultivation industry. The public fiscal budget expenditure, spending on science and technology, and the number of faculties in institutions of higher learning all have a significantly positive impact on the growth of marine fishery production. Fiscal spending of government on science and technology and increase in top-notch technological talents are important propellers of scientific and technological development. They can, in addition to improving the regional total-factor productivity, provide effective countermeasures for the fishery to cope with climate and environmental changes and drive the steady growth of fishery production. It should be noted that the educational fiscal expenditure has a significantly negative impact on the marine fishery production growth. This can, to some extent, indicate that increase in educational fiscal spending is not conducive to the development of marine fishery culture. According to labor marked theory, in cases where all workers face the same labor market opportunities, educational fiscal spending could effectively improve the years of education for individuals, especially low-income populations, leading to a parallel upward shift in the wage frontier and a redistribution in labors. In other words, under the conditions that the wage gap across sectors is opening, more years of education could improve the probability for the fishers to obtain a non-agricultural employment with a relatively higher wages. Therefore, the improvement of educational fiscal expenditure provides an economic motivation of transfer of the labor from fishery industry to non-agricultural sectors (e.g., manufacturing, hotels, and catering services), thus leading to a decrease in the growth of marine fishery culture.
On that basis, the amount of environmental pollution is interacted with pollution degree and vulnerability degree in Model (3)-(5), which are used to further verify the sources of heterogeneity in hazards of environmental pollution. Regression results indicate that, with the pollution level and vulnerability degree introduced, the heterogeneity in hazards of environmental pollution across development zones is more obvious. Besides, the backward regions are suffering to more serious damage. In terms of the causes of heterogeneity, the coefficient on the interaction between total sewage discharge and pollution level, although negative, did not pass the significance test, suggesting that the amount of marginal damage caused by pollution did not differ according to the difference in pollution level, i.e., the difference in pollution level was not the main cause of the heterogeneity in marginal damage from pollution at this stage. As a center of economy, the eastern coastal regions are densely populated with industrial enterprises and its rapid economic development has led to more serious environmental pollution, which in turn has led to a relatively high level of pollution in the eastern region as a whole, thus reducing the difference in marginal damage from pollution to some extent. On the other hand, as a typical public good, the quality of the marine environment has obvious externalities and non-exclusivity. An increase in environmental pollution in some areas can not only directly affect the place where the pollution occurs but can also have cascading negative externalities. At the same time, the estimated coefficient of the interaction between vulnerability and total sewage discharge is significantly positive at the 1% confidence level, indicating that increased coping capacity can effectively mitigate the negative shocks of pollution. In addition, for different sources of pollution, the estimated coefficients of the interaction coefficients of industrial sewage discharge, urban domestic sewage discharge, and pollution level and vulnerability level are basically the same, i.e., the impact of environmental pollution on marine fishery culture in relatively vulnerable areas is more pronounced. Due to the differences in factor structure and economic development patterns, there is a clear bias and agglomeration in technological progress across regions (Inyinbor Adejumoke et al., 2018;Jayaswal et al., 2018). Developed regions have a clear advantage in terms of material base and technological availability compared to lagging regions, and therefore, have a wider range of effective tools to deal with the negative impacts of environmental pollution.

Analysis of Results From Models With Lags
Based on the former panel regression model, this paper adds the industrial sewage discharge with five lags at most to better understand and estimate the impact of environmental pollution and heterogeneity in the impact. Table 5 presents results from the model with no lags, one lag, three lags, or five lags.
According to regression results, under the condition that the lags of total sewage discharge are introduced, the negative impact of environmental pollution on marine fishery culture is still significant. If the lag is not introduced, the increase in the unit sewage discharge can lead to a decrease of 0.249 in the output growth of the backward region. With one lag of total sewage discharge added, the impact of the sewage discharge on the output growth increase to 0.309. After more lags are introduced, cumulative effects of environmental pollution are generally similar. The unit sewage discharge increase reduces growth by 0.215-0.766. It can be seen that although the estimated coefficients of the lags are interleaved, the sum of the estimated coefficients of all lags is not significantly zero. This suggests to some extent that there is a significant negative cumulative effect of environmental pollution on regional marine fish farming production, with environmental pollution not only affecting current average output levels but also affecting productivity through lagged effects. This may be due to the fact that environmental pollution increases the probability of the occurrence and spread of fish fry diseases, which in turn leads to a reduction in cultured output. At this point, producers, as rational economic agents, will inevitably reduce the amount of inputs of factors of production based on the goal of profit maximization, resulting in lower productivity that has a lagged effect. When the capital stock is adjusted to a new steady state, the permanent impact on productivity may affect subsequent capital accumulation, resulting in a lagged growth effect. As can  be seen, there is not only a horizontal but also a growth effect from the negative effects of environmental pollution, both of which may have an impact on the output growth of marine fishery culture. The difference is that when pollution returns to normal levels, the effect of environmental pollution on average output can be restored to normal levels, i.e., its horizontal effect on output will be reversed. In contrast, the lagged growth effect on productivity is not reversed by an improvement in environmental pollution, and the presence of this effect means that the environmental declines in one period leaves the fishery culture economy permanently further behind. According to above theoretical analysis, environmental pollution not only affects the average output, but also affects the fishery productivity. The difference is that the impact on output disappears as pollution degree returns to its prior state. However, the impact on the fishery growth rate will not reserve, along with the alleviation of the pollution degree. Though the estimation coefficient of the lag is an interleaving positive and negative number, the sum of estimation coefficients on lags is obviously not significantly zero. This means that the impact of environmental pollution on the output growth is significantly long term. In other words, its impact on the productivity is obviously higher than its impact on the average output. This also means that the aggravation of pollution over one period might cause the backwardness of this region in the long term.
In addition to examining the impact of total sewage discharge, this study also takes into account differences in the impact of different pollution sources. According to analysis results, this study adds the lags of the industrial sewage discharge and the domestic sewage discharge. The dynamic panel model is built for a better understanding and estimation of the impact of environmental pollution and the heterogeneity of the impact. Tables 6, 7 presents results from the model with no lags, one lag, three lags, or five lags. As shown in Table 6, as the lags increases, the impact of the core explaining variable, industrial sewage discharge on the output growth gradually shifts from positive to negative, but it fails to pass the significance test. The result is also reflected in the cumulative effect of the industrial sewage discharge. By calculating the cumulative effects of different lags, this study easily observes that, though the above result fails to pass the significance test, the impact of industrial sewage discharge on marine fishery production growth becomes increasingly obvious, with the passage of time. As the central economic regions of China, industrial development in the eastern coastal regions of China is obviously clustered. This allows the industry to capitalize on the efficiency of labor distribution among enterprises to significantly promote the concentrated governance of industrial pollution. After that, the negative impact of water pollution can be maximally alleviated through collective discharge and concentrated treatment (Wang et al., 2021).
As to urban domestic sewage, estimated results of dynamic panel model in Table 7 shows that the impact of the current domestic sewage discharge on output growth is significantly negative in the models with lags. This can, to some extent, reveal that the increase in domestic sewage discharge can significantly hinder the development of marine fishery culture industry. Of special note is that the domestic sewage discharge   exhibits significantly negative cumulative effects. With three lags of sewage discharge introduced, the marginal impact of the domestic sewage discharge on output growth increases from 0.302 to 0.933. Summing up the above analysis results, this study points out that the impact of the urban domestic sewage on the marine fishery production should not be ignored, compared with the impact of industrial sewage discharge. Though the Chinese government has adopted a series of severe environmental regulations, most are mainly targeted at industrial enterprises. Seldom are there any regulations intended for the treatment of urban domestic sewage. Meanwhile, as the resource consumption and pollutants increase along with the increasing population, domestic sewage has been a critical part of the water pollution treatment (Zuo et al., 2019).

Comparison of Estimated Results During Different Periods of Time
In 2017, the Chinese government launched the Second National Agricultural Census, which led to the adjustment of statistical standards of the marine fishery cultured output in different regions, such as redefining the standards to compute shell products (oysters, scallops, etc.). In order to prevent differences of data standards from impairing the precision of estimated results, and also considering differences in recalculating the output of marine fishery cultured product s across regions, this study adopts 2007 as the demarcation and divides samples into two periods, namely, from 1994 to 2006 and 2007 to 2018. The current effects of the total sewage discharge and the cumulative effect of one, three, and five lags are analyzed.

Comparison of Estimated Results in Different Regions
Statistical errors can be easily found with statistical modes, standards, and accuracy in underdeveloped regions (Ciccone and Jarociński, 2010;Deaton and Heston, 2010). The credibility of estimated results heavily depends on the accuracy of the survey data. Therefore, according to the regional differences of economic development in China, this study divides the samples into the northern region and southern region with Yangtze River as the demarcation for a robustness test 4 . The test results are shown in Table 9. Similarly, considering the direction and significance of the estimation coefficient of the model, this study observes that the impact of the regional total sewage discharge on the output growth has not yet undergone obvious changes, and that the significance has also improved. Combining the above robustness test results, it is self-evident that the research findings obtained by estimation in this study are consistent, robust, reliable, and persuasive.

CONCLUSION AND RECOMMENDATIONS
This study takes marine fishery culture as an example, based on the seawater cultured products data of 53 coastal cities in China from 1994 to 2018. From theoretical analysis to empirical discrimination, this study captures the heterogeneous impact of environmental pollution on output growth. It is found that the weakening effect of the environmental pollution on the output growth of seawater products has apparent regional heterogeneity, and this damage difference mainly comes from the different ability to deal with the environmental pollution among different regions, rather than the difference of pollution degree. In terms of different wastewater sources, the impact of industrial wastewater discharge on marine fishery culture is not very significant. In contrast, the adverse shocks of urban domestic sewage are more prominent. Notably, the impact of environmental pollution has a significant cumulative effect, and this negative impact is just irreversible. Therefore, the aggravation of environmental Huizhou, Shanwei, Yangjiang, Dongguan, Zhongshan, Chaozhou, Jieyang, Beihai, Fangchenggang, Qinzhou, Haikou and Sanya.
pollution will further widen the gap in marine cultured output between regions. Under the background of the green sustainable development of marine fishery in China, a full understanding of the possible negative impact of pollution is necessary to actively take measures to respond to the challenge. Generally speaking, there are two ways to deal with environmental pollution: "adaptation" and "mitigation." First in terms of "adaptation, " an emphasis and increase in research, development, and the promotion of relevant technologies (e.g., the warning mechanism of environmental pollution, the development on new varieties, the purification equipment research, and development) are needed. More importantly, we should pay attention to the universality of technological progress. By constructing a reasonable technology subsidy policy, with the help of knowledge popularization and technology promotion, we can guide the transfer of advanced technology to backward areas, which improve the coping ability of fishers there effectively and urge them to take adaptive behavior to alleviate the negative impact of environmental pollution and realize the balanced development among regions. Secondly, in terms of "mitigation, " on the one hand, innovation should be actively promoted and new techniques for urban wastewater should be explored to address the environmental damage of pollutions; on the other hand, improvements are needed in the chemical articles for daily use, such as liquid detergent and toilet cleaners, etc., and more pollution-free articles should be developed. Meanwhile, given the spatial spillover effect and relevance of marine pollution, through the identification of interest demands of different subjects, we shall establish more diversified governance incentives and compensation mechanisms, scientifically delimit the cross-provincial and cross-municipal boundaries of marine pollution control, and increase the public responsibility awareness and supervision efficiency by further building a communication platform between the government and non-governmental organizations.
Although our study confirms the role of pollution in the production of marine fishery culture, further investigation is necessary to improve the accuracy of pollution reduction strategies. One of the main premises for the application of the model in this study is that the marine fishery culture development has been saturated, whereas marine aquaculture has not reached fullness in the early years, and some areas have not even been developed due to environmental protection and other factors. Another limitation of our study is that the industry production conditions are poorly treated. For instance, biological, and technological conditions, and economic environment for the producers, e.g., product markets and input markets (labor market) in fishery production, have not been considered. Therefore, a more substantial data should be collected, which is important for future research to compare the shocks of pollutions and evaluate the treated effectiveness of different related strategies.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.