Next Article in Journal
Linking Land Tenure and Integrated Watershed Management—A Review
Previous Article in Journal
Public Environment Emotion Prediction Model Using LSTM Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Does China’s Municipal Solid Waste Source Separation Program Work? Evidence from the Spatial-Two-Stage-Least Squares Models

1
School of Economics, Zhejiang Gongshang University, 18 Xuezheng Road, Hangzhou 310018, China
2
Zhejiang Gongshang University Hangzhou College of Commerce, 66 Huancheng South Road, Hangzhou 311500, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(4), 1664; https://doi.org/10.3390/su12041664
Submission received: 25 December 2019 / Revised: 11 February 2020 / Accepted: 20 February 2020 / Published: 23 February 2020
(This article belongs to the Section Economic and Business Aspects of Sustainability)

Abstract

:
This paper evaluates the impact of the second municipal solid waste (MSW) source separation program on municipal solid waste generation (MSWG) in China. Without considering the spatial interactions between cities, the second MSW source separation program has a nonsignificant adverse effect on the per capita municipal solid waste generation (PMSWG). Relaxing the stable-unit-treatment-value assumption (SUTVA), which holds in most of the previous estimation literature on treatment effects, involving the spatial spillover effect among cities, as well as correcting the endogenous local policy has a significantly negative but not robust impact on the PMSWG. The estimation results of the generalized nesting spatial regression models (GNS) imply that the spatial interaction characteristics among Chinese prefecture-level cities may, if neglected, lead to underestimation of the reduction effects of the second MSW source separation policy on the absolute amount of PMSWG. More importantly, our study indicates that although not all the spatial econometric models support the significant reduction effect of source separation on the absolute amount of PMSWG, the source separation program significantly reduces the relative amount of PMSWG, and this result is robust in all spatial models.

1. Introduction

The problem of municipal solid waste (MSW) caused by economic growth and urbanization is expected to be a major worldwide challenge in the future, especially for developing countries. From the perspective of economic theory, the following three types of tax/subsidy schemes can be used to address the problem of MSW: Purchase-relevant instruments, discard-relevant instruments, and jointly-relevant instruments [1], which are for pricing solid waste collection and disposal services at the different stages of municipal solid waste management (MSWM). The economic research has found multiple combinations of these policies that can encourage an efficient MSWM [2,3,4,5]. However, all of these policy designs are flawed to varying degrees [6,7]. Due to the feasibility of various tax and subsidy policies in practice, the municipal authorities of developing countries have limited policy tools and mainly rely on traditional Command and Control Policies. Those policies involve imposing restrictions on the recycling proportion of all packaging waste [2], using a certain percentage of recycled materials [8], and implementing MSW source separation programs.
MSW source separation is one of the major MSW reduction approaches adopted worldwide. It has been proved to be an effective method to achieve waste minimization and circular economy in developed countries [9,10]. Similar attempts were under way in many developing countries, for example, the transformation of mandatory source separation under regulation in Malaysia [11], the encouragement of source separation through economic incentives in Ghana [12], the estimation of dynamic effects of diverse source separation rate on MSW collection and transportation in Thailand [13]. As the largest developing country and the largest MSW producer [14], China deserves more attention. Although there are distinct socio-political circumstances and traditions across developing countries and waste source separation program is not easily transferable, the lessons from waste classification practices, economic incentive measures, policy design, and implementation of waste separation in China may provide some novel understandings and future considerations relating to complex waste management system in other developing countries.
It is estimated that up to 70% of East Asian and Pacific MSW was produced in China [15]. In addition, as much as 1.32 to 3.52 million tons of plastic waste enters the ocean each year from China, the highest among the 20 countries assessed and more than three times as much as second-ranked Indonesia [16]. The same as other developing countries, Chinese cities are facing tremendous pressure in terms of MSWM [17]. In view of the critical situation of MSWM and the damage to the global environment and ecosystems caused by the rapid growth of MSW, China should hold the increase of MSW in check and improve the sustainability of MSWM. The governance experience of MSW in China is an important piece of the puzzle for global waste management.
China’s experience on MSW reduction has equal importance to other developing countries that are facing the same dilemma of both GDP per capita and MSW increasing. Currently, the local governments of China have two major policy instruments, i.e., pricing schemes and top-down source separation policies that are available to control the current quantities of solid waste. Since the pricing schemes have been verified to be a failure by previous studies in both China and other countries [18,19], MSW source separation is currently the dominating choice in China.
Unfortunately, although after more than a decade of attempts of MSW source separation, it remains unclear whether the MSW source separation program is working or not in China. This paper proposes a spatial-two-stage least squares (S2SLS) framework for evaluating the program impact on the per capita municipal solid waste generation (PMSWG) by involving the spatial spillover effects among cities as well as correcting the endogenous local policy.
The main contributions of this study are as follows. First and perhaps the most important, this paper allows for the interference of PMSWG among cities after confirming the pervasive “proximity” effects and that the independent assumption among adjacent units is no longer valid. So far, only a few papers have formally considered the presence of spatial spillovers from treated to untreated units, and this raises relevant econometrical estimation questions regarding the identification of treatment effects when interactions among units are plausible [20].
Second, the possibility of endogenous program choices leads to a biased estimation of the treatment effect. The omission of unobservable characteristics among cities in the econometric analysis, such as environmental awareness, might change the probability that a city implements “green” policies. Therefore, if the MSW source separation program is a purposeful action, identifying and controlling for the factors that lead the program to change is necessary [21].
As far as we know, few studies on MSW program evaluation consider these two potential problems above. The rest of the paper is organized as follows. Section 2 provides the background of China’s MSW source separation program and describes the methodology. Section 3 demonstrates the data and research design. Section 4 presents the empirical results and discussion. Section 5 concludes.

2. Background and Description of the Methodology

2.1. Background

Generally speaking, the MSW collection system in China was a mixed collection system before 2000, since there was no formal MSW separation policy at that time and the MSW was mixed by residents and discarded at fixed collection points in the community. Then, MSW was collected, transported and disposed of by the municipal authorities, which is called the formal collection system. Economically motivated, residents sort out some of their high commercial value MSW and sell it to recyclers. In addition, as in many other developing countries, the waste pickers, who act as the informal collection system [15], sort out the materials from the rest for sale, but most of the waste was mixed up and thrown away without separation.
The mixed collection system made it difficult for municipal departments to conduct secondary MSW separation, even some cities have waste separation facilities, but the efficiency is quite low. Therefore, a large amount of recyclable MSW is directly buried or incinerated. On the one hand, it decreases the efficiency of resource utilization, and on the other hand, terminal disposal pressure increases. It is estimated that only 20% of MSW is recycled and largely contributed by informal waste picking [14]. Therefore, to improve the recovery and separation efficiency of the formal collection system is critical to enhancing the overall efficiency of MSWM. Until now, the municipal charges for MSWM services have been levied mainly in two ways. The first is to levy charges attached to public utilities such as according to the water, electricity, or gas consumed by residents. The second is charging citizens through the residents’ committee or property management committee directly. Since the pricing schemes have been verified to be a failure by previous studies [18,19], the top-down encouragement of MSW source separation is currently the dominating choice in China.
So far, China has conducted three top-down source separation pilot programs. These pilots are exactly the policies that can affect the efficiency of the formal collection system. Because the policies are issued in the form of government regulations, the pilot cities are required to build a separation collection system and set up source classification facilities in communities. The difference between the three major programs is that the scope of policy implementation gradually expands from a small scope to more cities.
The first was launched in eight major cities, Beijing, Shanghai, Guangzhou, Shenzhen, Hangzhou, Nanjing, Xiamen, and Guilin, selected by the Ministry of Housing and Urban-Rural Development (MOHURD) in 2000 [22]. It did not define a clear category of MSW classification but advocated that eight pilot cities should carry out separation MSW collection as soon as possible to improve the performance of MSWM. In 2004, MOHURD promulgated the “Separation Method and Evaluation Criteria of Municipal Solid Waste” which divides MSW into six categories: recyclable, bulky, compostable, combustible, hazardous and other waste [23]. Cities were following this reference standard to enforce separation MSW collection.
After 14 years of preparation, each pilot city had formed a relatively systematic separation MSW collection scheme, such as different color bins were assigned at the community waste collection points, which were used to collect different kinds of MSW. The second started in 2015, MOHURD and the other five departments jointly issued a document that established 26 demonstration cities (districts) of MSW source separation [24]. The document stipulates that the separation coverage of MSW in residential communities should reach 90%, PMSWG should drop by 6%, and the recycling rate of MSW should reach 60%. Meanwhile, the second further emphasized the separation of food waste, since organic food waste is the dominating waste in developing counties such as China [25]. At the same time, the pilot cities have also begun to build food waste treatment facilities. However, this program actually started in March 2014, as the central government advocated an initiative that all local governments in China should participate actively [26].
The third was in 2017. The central government forwarded an announcement prepared by MOHURD in March requiring 46 cities to carry out mandatory MSW source separation policies [27], including all the provincial capital cities and public institutions, such as party, government and military organs, schools and hospitals. The purpose of the three top-down source separation programs was to improve the efficiency of waste recycling and reduction for the formal collection system. Meanwhile, from the perspective of residents, the participation rate will also affect the efficiency of the system. Mandatory programs are thought to have higher participation rates than voluntary programs [28], and some believe that source separation should be a key priority and that the program pilot in China should be extended to more cities [9].
Up to now, cities implementing source separation policies have formulated different separation standards according to their own characteristics. Most cities adopt the four-category approach, which divides the MSW into recyclable, food, hazardous, and other waste. However, the specific names of each city are slightly different. For example, food waste is named perishable waste in Hangzhou, as well as Shanghai calls food waste wet waste and other waste dry waste to make it easier for citizens to distinguish. In addition, the standard of Guangzhou even lists the categories of food waste in detail for the convenience of citizens to discard. When community residents discard waste, they need to put different waste into the corresponding color bin according to the regulations. Different garbage trucks would collect the waste according to the color of bins and transferred it to various waste facilities by local municipalities. Thus, it can increase the recycling rate of the formal collection system and also reduce the pressure of the landfill or incineration disposal facilities.
Diverting waste from landfills has many environmental benefits, such as the reduction in greenhouse gas emissions, the minimization of waste pretreatment operations, and energy recovery by composting [29]. Aphale et al. (2015) [30] verified that the recycling performances are positively related to the separation efficiencies. Thus, estimating the effectiveness of source separation is significant when considering the actual impact of the program. The government could benefit from it and improve the performance of MSWM. Since the implementation of the first pilot was 19 years ago, the public had little concern regarding environmental consciousness. Additionally, the municipal authorities did not provide enough support for it. Meanwhile, Han and Zhang (2017) [31] has proved the ineffectiveness of the pilot. The time interval between the third and the second pilot is relatively short, so the potential hysteresis of the second pilot may exit. Therefore, this paper mainly examined the causal effects of the second MSW source separation program on the PMSWG.

2.2. Description of the Methodology

In most empirical works, program evaluation with nonexperimental data needs to impose a strong assumption known as the stable-unit-treatment-value assumption, or SUTVA [32]. SUTVA requires that the treatment effects of units are independent of each other, that is, the spillover effects between the potential outcomes of units must be ruled out.
SUTVA might be plausible in many biomedical experiments and has been adopted by many previous program evaluation studies [33,34]. However, it is quite intuitive that SUTVA may not be true in many scenarios, such as in the research and development (R&D) activities of companies and epidemiology. It might be difficult to believe in the absence of “spillovers” among units [20]. Similarly, and a fortiori, the impact of implementing the second program is not only on the reduction of PMSWG but also, perhaps, on spillovers to the neighboring cities. So far, only a few papers have formally included the interaction between units in the literature on the estimation of program evaluation [20], let alone in the MSWM research.
Meanwhile, some believe that the SUTVA assumption is not testable in principle [34]. However, spatial factors actually play an important role in environmental and economic issues. Ignoring spatial factors may omit important variables, which can cause bias estimates for the model parameters [35]. According to the first law of geography [36], the sample data based on geographical location are often not independent of each other but have spatial dependence. Generally, the spatial dependence between adjacent units is stronger than that between distant units.
Emerging empirical literature has explored the spillover effects in the MSWM domain, such as the spillover effects of MSWM performance [37], MSW recycling outcomes [38], and PMWSG [39]. In terms of the spillover effect in PMSWG case, there may be several reasons why PMSWG in one city can be influenced by policies and practices adopted in neighboring cities.
Firstly, at the individual behavior level, the decrease in PMSWG represents the improvement of residents’ pro-environmental behaviors and indicators of sustainability. The psychological theory argues that behavior consistency effect and social identity effect may foster a positive spillover effect of pro-environmental behaviors [40]. Agovino et al. (2016) [41] adopted the waste collection habits as a proxy of pro-environmental behaviors, they found that provinces with good pro-environmental habits may positively influence neighboring ones in Italy. Crocita et al. (2016) [42] employed the separation collection of MSW as representing a pro-environmental attitude and confirmed the presence of positive proximity effects among virtuous neighbors on Italian provinces. In other words, if spatial units are surrounded by units with good waste separation habits, they are more likely to adopt the same behavior since they feel part of a virtuous circle [37]. Thus, a decrease in PMSWG in one city may lead to a decrease in PMSWG in nearby cities due to similar social and cultural environments. The spillover effect decreases as the distance increases, which is embodied in the specification of the spatial weight matrix described below.
Secondly, considering the regional differences, neighboring cities with similar characteristics may be easier to exchange and emulate each other’s experiences [38], and in this way, affect each other’s MSWM policies and, in turn, the PMSWG. Furthermore, due to the similarity of economic development level and consumption culture in adjoining regions, the implementation of MSWM policies, such as source separation program in pilot cities, is more likely to affect the MSW separation awareness and classification knowledge of residents in adjacent cities. These positive proximity effects were also observed by the pioneering work in Oskamp et al. (1991) [43]. Thus, the PMSWG of a city may be influenced by neighboring cities.
Thirdly, MSW is a part of human life and will be generated along with human economic activities. Moreover, the change in economic mode will inevitably affect the production and types of MSW. Economic convergence clubs within China have been verified [44], which implies that there will be different economic agglomeration blocks, and the interaction between adjacent cities within the blocks will be greater than that between cities farther away. If the economic growth rate is relatively high in one city, the adjoining cities may imitate its economic development model and industrial structural arrangement, as by-products of economic activity, thus the types and production of MSW in neighboring cities will be affected by the economic development and the MSWM policies of this city.
Logically speaking, the PMSWG of each city in China will demonstrate spatially related characteristics. The choropleth maps of PMSWG in China in 2002, 2016, and the average of 2002 to 2016 are presented in Figure 1. The darker areas indicate a higher PMSWG. The PMSWG appears to be spatially dependent since the cities with higher PMSWG are clustered together. That is, the stylized fact that PMSWG exhibits homogeneity in adjacent cities is observed. If the PMSWG of each city was randomly distributed with regard to location, we would not see the clusters. A formal test of spatial correlation will be conducted below.
In fact, in other fields, many scholars have also included the spatial interaction relationship in their research framework [45,46,47,48]. Therefore, if the pervasive spatial dependence is confirmed among cities in China, the spatial spillover effects should be considered in the program evaluation of MSW studies.
In addition, identifying the causal effect of the second program, T i , on PMSWG faces seemingly insurmountable problems in nonexperimental data, where T i = 1   iff city i   implemented the second program, otherwise T i = 0 . Generally, in the program evaluation literature, T i = 1 is called the treated group and T i = 0 refers to the control group. However, the program is usually not randomly assigned, and there may be various significant unobserved factors, including the alternatives available to cities, their environmental preferences, waste management achievement, and the potentially idiosyncratic political benefits for local bureaucrats in an interjurisdictional competition context [49]. More generally, the estimation of the casual effect of T i may be biased and inconsistent when T i is not randomly assigned, and this essential problem is called selection bias.
To control for the possibility of endogenous program choices, the decision of whether to adopt the second program or not was modeled as a function of observable exogenous variables. Then, the predicted probabilities for the program implementation were used to substitute the possible endogenous variable,   T i , in the PMSWG equation using two-stage least squares (2SLS). With data on ex-ante and ex-post program intervention, 2SLS allows a more flexible strategy to deal with unobserved heterogeneity and permits these characteristics to vary with the duration of the program.
Finally, we relaxed the SUTVA by including different spatial lag terms in 2SLS to control for the spatial interaction effects, and the treatment effect of the second program was estimated by setting different spatial weighting matrices. The data and detailed empirical strategy are introduced in the following section.

3. Data and Empirical Design

3.1. Data

The sample data include surveys both before and after the second source separation program intervention. Therefore, a panel data of 569 prefecture-level cities in mainland China between 2002 and 2016 are collected in this study. Using panel data could correct the bias caused by omitting variables and control for heteroskedasticity within different cities. We merged the panel data with the shapefile of prefecture-level cities [50], which has geographical information, to create spatial weighting matrices, and finally reduce our sample size to 288 prefecture-level cities.
All data are from Chinese official statistics, including the “China Statistical Yearbook”, the “China Statistical Yearbook of Urban Construction” and City Statistical Yearbook of each city [51,52]. The indicator variable SEP 2 = 1 if the city has implemented the second source separation program, and 0 otherwise. Let PROCAP = 1 if the city is the provincial capital city, and 0 otherwise. Some dummy variables of geographical location, EAST, MIDDLE, WEST, and NORTHEAST indicate that these cities are located in the east, middle, west, and northeast of China, respectively. As a summary, the descriptions of the variables are listed in Table 1.
The total and average MSWG in the treated cities are significantly larger than the average level of untreated cities. At the same time, the two groups have significant differences in terms of the GDP per capita, population density, and urbanization rate. The table illustrates that the data are unbalanced, and prefecture-level cities display considerable variation in MSWG, economic characteristics, and demographic variables, which may influence whether the cities implement the program or not, i.e., the program was not randomly assigned. This nontrivial difference may cause potential endogenous problems.

3.2. Specification with SUTVA

2SLS relaxes the exogeneity assumption of OLS and is robust to time-varying selection bias, which some identification strategies cannot control, such as the difference in differences estimation (DID). In the following section, we set the 2SLS econometric model with SUTVA, as in many previous program evaluation studies, and then relax the assumption by adding the spatial spillover effects in the model specification.

3.2.1. A Discrete Choice Model

The first-stage regression is to construct a nonlinear binary response model to calculate the predicted probability of each city being selected into the second program, i.e., regressing the treatment dummy variable on the vector of exogenous covariates.
A discrete choice model is used to depict the decision of whether each local government implements the second program or not. The most prominent types are Logit and Probit, both of which can be derived under the assumption of the utility-maximizing behavior of the local government. The discrete choice model can be derived by random utility theory [53]. Each local government is assumed to have a tradeoff between the potential benefits and costs of implementing the program. A latent variable U i S E P 2 is defined as the net benefits for the city   i from adopting the program, as follows:
U i S E P 2 = γ Z i S E P 2 + ε i S E P 2             S E P 2 { 0 , 1 }
where U i S E P 2 is the net utility of whether city i chooses to implement the second source separation program ( S E P 2 ) or not. The vector Z i S E P 2 includes all the variables that may help determine the local government’s choices, such as the number of vehicles and equipment for municipal environmental sanitation, the number of Harmless Treatment Plants, and other exogenous characteristic variables of city i , and γ is the corresponding vector of coefficients to be estimated. ε i S E P 2 is the disturbance term. We cannot observe the net benefits U i S E P 2 from implementing the program. Instead, we only know whether a city has implemented such a program, i.e., as follows:
S E P 2 = { 1 ,     iff   U i S E P 2 > 0 0 ,     otherwise
Therefore, city i chooses to implement the program with the following probability:
P r ( S E P 2 = 1 | Z i S E P 2 ) = P r ( γ Z i S E P 2 + ε i S E P 2   > 0 | Z i S E P 2 ) = P r ( ε i S E P 2   > γ Z i S E P 2 | Z i S E P 2 ) = P r ( ε i S E P 2   < γ Z i S E P 2 | Z i S E P 2 )
If we assume that ε i S E P 2 follows a standard logistic distribution, the familiar conditional Logit model arises.
P r ( S E P 2 = 1 | Z i S E P 2 ) = ( 1 + exp ( γ Z i S E P 2 ) ) 1
We used this model to depict the choice of each city and to generate the predicted probability   P R ^ that each city will choose to adopt the program. P R ^ is used as the instrument variable (IV) to substitute the possible endogenous variable, S E P 2 , in the second stage regression.

3.2.2. The 2SLS Model

Conventionally, we use the following equation to estimate the treatment effect of the S E P 2 on PMSWG:
l n Y i t = α + β S E P 2 + θ X i t + μ i + φ t + ε i t
where Y i t denotes PMSWG, X i t is the vector of covariates that affect PMSWG and θ is the vector of coefficients on X i t . μ i is the city-level fixed effect, φ t represents the time effect, and ε i t is the disturbance component. β is the program impact that we want to estimate. The predicted probability P R ^ embodies only the exogenous variation in the treatment. Then, we use P R ^ as the exogenous IV for treatment placement, creating the following second-stage regression:
l n Y i t = α + β P R ^ + θ X i t + μ i + φ t + ε i t
In this study, taking the logarithmic form of some major continuous regressors can reduce the data volatility and the potential disturbance of the heteroscedasticity.

3.3. Spatial Econometric Specification

3.3.1. Spatial Dependence Tests

In the first step, the spatial dependence of sample data should be tested in advance. If there is no spatial dependence, the traditional regression model is sufficient. The Moran’s I index [54] and Moran scatter plot are the quantitative indexes that are usually used to measure the spatial dependence. Moran’s I is defined as follows:
I = N ( e W e ) S ( e e )
where N is the number of observational spatial units, e   is the residual vector of OLS, W is the spatial matrix, S is the normalized factor equal to the sum of all the elements in the weighting matrix. If we normalize the rows of the weighting matrix, Equation (7) reduces to the following:
I = e W e e e
In our study, W ( 1 )   is the distance threshold spatial weighting matrix, based on the specification that if city   i   and   j   are farther apart than a distance threshold   d , then w i j be set to 0, otherwise   w i j = 1 . To ensure that each city has at least one adjacent city, we set the distance threshold as the longest distance between city   i   and city   j   by using the method applied in spatial econometric practice [55]. Concretely, d   equals to 2109 km, which is the maximum arc distance between the two cities in the sample data. In accordance with the queen contiguity rule to specify   W ( 2 ) , if two cities, i   and   j , have a common border or vertex, then   w i j = 1 , otherwise   w i j = 0 . The specific forms of the two spatial matrices are as follows:
W ( 1 ) : w i j = { 1 ,         i f     d i j d   0 ,         i f     d i j > d   W ( 2 ) : w i j = { 1 ,       i f     i   j   c o n t i g u i t y 0 ,       i f     n o t

3.3.2. The S2SLS Models

If Moran’s I test and Moran scatter plot confirm the existence of spatial interaction, appropriate spatial econometric models should be set to control for spatial interference. We set the spatial correlation in the second stage of the above regression process for correct endogeneity as follows:
l n Y i t = α + β P R ^ + μ i + φ t + θ X i t + ρ W l n Y i t + δ W x j t + ( I σ W ) 1 ε i t
W   is the spatial weighting matrix, which characterizes the spatial relationships between the observable units and the elements w i j specify the potential spillover among interactive units   i and   j . W l n Y i t is the lag of the dependent variable, and ρ W   measures the spillover effects of nearby outcomes.   W x j t is the lag of the spatial correlated independent variable x j t , and δ W   measures the spillover of nearby covariates   x j t   to the outcome of unit   i .   σ   is the spatial correlation parameter of the residuals, which account for the spatial interactions among unobservable neighbor factors. Equation (10) is to estimate the absolute reduction effect of SEP2. At the same time, we will further examine the relative reduction impact of SEP2, that is, the treatment effect of SEP2 on the growth rate of PMSWG:
    d l n Y i t / L . l n Y i t = α + β P R ^ + μ i + φ t + θ X i t + ρ W l n Y i t + δ W x j t + ( I σ W ) 1 ε i t
  d l n Y i t   is the increase of PMSWG and L . l n Y i t   refers to the first lag of l n Y i t , i.e., l n Y i ( t 1 ) . Thus Equation (11) modifies the dependent variable for the growth rate of PMSWG in Equation (10).
The setting of Equation (10) is called the generalized nesting spatial model (GNS), which is the most general spatial econometric model. If δ = 0 ,   σ = 0 , the GNS will be degraded into a spatial lag model (SLM), also known as a spatial autoregressive model. If ρ = 0 ,   δ = 0 , the GNS will simply become a spatial error model (SEM), If σ = 0 , the GNS reduces to a spatial Durbin model (SDM). Thus, SEM, SLM, and SDM are special forms of GNS. In this study, the GNS is used to estimate the spatial regression model, to take full account of the spatial interaction.
In the spatial econometric model, we control the spatial dependence by setting three different spatial weight matrices. The specification of the contiguity spatial matrix W ( 2 )   (CON matrix) continued to be defined by the rule of queen contiguity previously described. The elements of the inverse-distance spatial matrix W ( 3 ) (IDIS matrix) equal to the reciprocal distance between cities, assuming spillover effects are proportional to the inverse of the distance between cities. The inverse-distance contiguity spatial matrix W ( 4 ) (IDIS-CON matrix) is a weighting matrix that contains the inverse distance for adjacent cities and is 0 otherwise. Three weighting matrices settings are shown below.
W ( 2 ) : w i j = { 1 ,       i f     i   j   c o n t i g u i t y   0 ,       i f     n o t W ( 3 ) : w i j = 1 / d i j W ( 4 ) : w i j = { 1 / d i j ,       i f     i   j   c o n t i g u i t y                           0 ,                       i f     n o t                                                                    

4. Results and Discussion

4.1. The 2SLS Evaluation

The 2SLS is adopted to control the self-selection problem. A Logit model is applied in the first stage to depict the decision of whether each local government implements the second separation program or not at the city level, i.e., regressing the treatment policy dummy variable (SEP2=1) on the vector of city-level exogenous covariates. Thus, in the first stage, we first include as explanatory variables a set of city-level exogenous variables correlated with T i but not correlated with the error term into Equation (1). Since the actual implementation of the program was in 2014, the probability of implementing the second MSW source separation program in each city was first estimated by using observable exogenous variables that were not affected by the policy before the implementation of the policy (2002–2013) and then, replacing the endogenous policy variables to test the program impact of the second MSW source separation. The first stage analysis is to generate the predicted probability P R ^ that each city will choose to execute the program. P R ^   is then used as the instrument variable (IV) to substitute the possible endogenous variable, SEP2, in the second stage regression. The reduction effects of the second source separation program on PMSWG can be estimated based on the correction of endogenous problems.
The regression results from the Logit model defined in Equation (4) are presented in Table 3. The fourth column in Table 3 displays the marginal effect of a change in the particular explanatory variable on the probability that a city adopts the second MSW source separation program.
After controlling for other covariates, the probability of implementing the program is estimated to increase by approximately 1.45% for an additional percentage of the number of vehicles and equipment for Municipal Environmental Sanitation (lnVEH) and to increase by 0.03% for a 1% increase in the percentage of Harmless Treatment Capacity per year (lnHLC). The number of Harmless Treatment Plants/Grounds (HLS) is also significantly positively related to the probability of enforcing this program. Perhaps the cities that have better performance regarding the infrastructure are more willing to apply for this top-down pilot program to acquire the potentially idiosyncratic political benefits.
The analysis also implies that the provincial capital cities are 8.89% more likely to implement the second MSW source separation program. In general, provincial capitals have more sufficient financial funds and developed economies, so they are more likely to adopt the program in response to this top-down political propaganda. From the marginal effects of other variables in Table 3, we can also infer the influence of other variables on the probability of cities choosing to implement the second source separation program. The results in Table 3 are instructive and interesting in their own right, but the primary purpose of estimating these discrete models is to create an exogenous prediction to substitute for the policy dummy variable S E P 2 in Equation (5).
In the second stage, the predicted probabilities of the second MSW source separation program ( P R ^ ) from the Logit model are used as an exogeneity predictor in Equation (6). To better approximate the exact relationship between environmental degradation and economic growth, Sobhee and K.Sanjeev (2004) [56] suggest that higher-order terms of major economic explanatory variables should be added to the regression. Therefore, we include GDP per capita (PGDP) in a nonlinear fashion in Equation (6), using the logarithmic form of PGDP and both its square and cubic form in the regression. The panel data estimation results of Equation (6) are reported in Table 4 when the SUTVA holds.
As shown in Table 4, the estimated coefficients for lnPGDP, (lnPGDP)2, and (lnPGDP)3 in all estimation approaches are significantly positive, negative, and positive, respectively, and the estimates are very close to each other. Therefore, it is appropriate to include higher-order terms of the economic variable in the regression analysis. The finding claims that the implementation of the second MSW source separation program is estimated to decrease PMSWG, except for the random effect model. However, neither estimate is significantly different from zero, which implies that after correcting the endogenous local policy program choice, the second MSW source separation program did not significantly reduce PMSWG. Thus, these estimation results cannot support the claims that China should vigorously promote the source separation program and extend it to more cities [9].
Some of the previous literature includes an impact evaluation of the demographic variables, population density, urban population, and urbanization rate on PMSWG [57,58,59]. We incorporated these covariates into the regression model and controlled the cross-product term of GDP per capita with the demographic variables and location dummy variables.
Since the location dummy variables do not change with time in the panel data, it is omitted in the individual fixed effect model and two-way fixed-effect model. At the same time, we also included the MSW pricing policy variables, i.e., the logarithmic form of the MSW disposal fee of each city (lnFEE), in the regression model. The estimation shows that the influences of the waste charging mechanism on the PMSWG seem to be negative in most approaches except for the random effect result, but the coefficients are insignificant in all models, suggesting the ineffective reduction effect of the waste charging policy on the PMSWG.
The null hypothesis that the constant terms are equal across the units of the individual fixed effect model is rejected. The statistic of the F test is 40.08, and the corresponding p-value is lower than 0.01 (with 287° of freedom [ d f ]), indicating that the pooled OLS would produce inconsistent estimates and the individual-specific heterogeneity should be controlled. In addition, the two-way fixed effect model adds the time effect to the individual fixed effect, and the time effects are jointly significant, at least at the 1% level ( F = 9.94 ,   d f = 14 ,   P < 0.01 ), suggesting that they should be included in the regression. The time effect is also controlled in the random effect model, and the c h i 2 test for the time effect also indicates that the null hypothesis of the insignificant time effect should be rejected ( c h i 2 = 238.65 ,   d f = 14 ,   P < 0.01 ) . In the Hausman-test framework for fixed effect and random effect model selection, the null hypothesis that the random estimator is consistent is soundly rejected ( c h i 2 = 51.06 ,   d f = 27 ,   P < 0.01 ). Therefore, the result of the random effect estimator is inconsistent, and the two-way fixed effect estimator is preferred.
In most previous works, it was generally assumed that there is no interaction among individuals, excluding spatial dependence, i.e., holding the SUTVA. However, if there is spatial interference among individuals, the estimation coefficient is biased. Therefore, the estimation shown in Table 4 that is obtained by ignoring spatial dependence is subject to bias and is not reliable. In the following part of this study, after confirming the existence of spatial dependence, spatial factors are controlled in the regression analysis to correct the estimation bias caused by ignoring spatial correlation.

4.2. Spatial Dependence Tests

We calculate the Moran’s I index and fit the Moran scatter plot by Geoda. Figure 2 shows the Moran scatter plots of the PMSWG of China’s 288 prefecture-level cities for 2002, 2016 and the average PMSWG from 2002 to 2016. The top three plots of Figure 1 use the setting of the spatial matrix W ( 1 ) , and the bottom three adopt W ( 2 ) . The low-low (lower left) and the high-high quadrants of each plot imply positive spatial interactions, and the low-high (upper left) and the high-low (lower right) quadrants indicate negative spatial dependence. As illustrated in Figure 2, most of the 288 cities appear in the upper right and lower left quadrants, 999 random permutation simulations were carried out for each Moran’s I index, and the pseudo p-values were all at least less than 0.01, which were significant at the level of 1%.
These results suggest that there is a nonnegligible positive global spatial dependence of PMSWG among cities. The appropriate spatial lag terms should be included to estimate the parameters of Equation (6). Otherwise, the regression results will be biased due to the omission of important variables.

4.3. The S2SLS Evaluation

In the last part of the empirical estimation process, we use Equation (10) to regress the PMSWG on the predicted probability of implementing the second MSW source separation program from the results of the Logit model in Table 3. Meanwhile, we also include three spatial spillover relationships in the following three ways:
  • include ρ W Y i t to allow adjacent outcomes to affect outcomes;
  • include δ W x j t to allow adjacent covariates to affect outcomes;
  • include ( I σ W ) 1 ε i t to allow adjacent unobservable errors to affect outcomes.
According to the above analysis, the two-way fixed effect model outperforms the other two models in Table 4. Therefore, we include the spatial weighting matrices,   W ( 2 ) , W ( 3 ) and W ( 4 ) , in the two-way fixed effect model, since the spatial econometric model is considering the spatial dependence on the basis of the traditional model with SUTVA. Then, the quasi-maximum likelihood (QML) estimator in Lee and Yu (2010) [60] is implemented to fit the GNS model.
It can be seen from Table 5 that the estimated coefficients of the spatial lag term of the dependent variable, economic variable (lnPGDP) and the error term, that is, ρ , δ , and σ , respectively, are all significantly positive at least at the level of 1% under the setting of three spatial weighting matrices. The results indicate that the PMSWG, the economic variable, and the unobservable factors that influence the PMSWG are not independent, that is, the SUTVA is invalid. The increase of the GDP per capita and the PMSWG in city i will have positive spillover effects on the PMSWG in city   j . Meanwhile, unobservable items that affect the PMSWG also have positive spillover effects between cities   i   and j. Thus, spatial effects should be considered in the regression model.
After controlling the spatial effects, the estimated coefficients of the nonlinear explanatory terms of the major economic explanatory variables, lnPGDP, (lnPGDP)2, and (lnPGDP)3, are all significant at the level of 1% under the setting of all the spatial weighting matrices, and their signs are positive, negative and positive, respectively, which are consistent with the estimated results of the conventional regression specification. However, after including spatial factors, ( P R ^ ) becomes significant at the level of 10% in the setting of the IDIS matrix, i.e., the implementation of the second MSW source separation program does decrease the PMSWG. However, this finding is not supported by the specification of the CON matrix and IDIS-CON matrix, which implies that we cannot infer the reduction effect of the source separation program safely.
At the same time, the estimated coefficient signs for demographic variables are also consistent with the previously estimated results without considering the spatial dependence relationship, indicating that the increase of the urban population will significantly reduce the PMSWG, which is consistent with the implication of the estimation results of the population density. This finding is similar to the observations that population agglomeration may increase the scale effect of MSWG [58], but it is contrary to the conclusion of Mazzanti and Zoboli (2008) [59]. Nonetheless, the increase of the urbanization rate does not necessarily represent the agglomeration of the population, but only the rise in the proportion of the urban population in a specific city administrative region. From this estimation, on the contrary, the increase of urbanization rate significantly increases the PMSWG, which confirms the research conclusion of Johnstone and Labonne (2004) [57] based on the utility maximization model.
In the spatial econometric model, the cross-product terms are still included in the regression analysis, and the signs of its estimation coefficients are still consistent with the traditional estimation method obtained from Table 4. Additionally, the estimated coefficient of MSW charging variable (lnFEE) is still not significantly different from zero in each spatial matrix specification, indicating that the waste pricing mechanism does not have the effect of inhibiting the PMSWG, which is consistent with the previous study [17].
From the above analysis considering spatial correlation, it can be seen that the implementation of the second MSW source separation program in spatial models significantly decreases the quantity of PMSWG to a certain degree in one of the spatial models. Thus, the previous estimation without concerning the spatial spillover effects may be biased downward, that is, the traditional regression with SUTVA may underestimate the impact of the program, in contrast to the spatial estimation.
Although the sample data are not able to support a significant and robust negative effect of the second MSW source separation program on PMSWG in all the spatial models, this negative estimation bias gives us important implications for further analysis.
Since the selection bias problem exists in our sample data, there are significant differences in the level of economic development between the cities that choose to implement the second MSW source separation and other cities. In the previous literature, empirical studies have been conducted on various linear and nonlinear connections of PMSWG with economic growth [61,62,63,64,65,66]. In many developing countries, the PMSWG has not been decoupled from economic growth, and even in some developed countries, absolute decoupling has not yet been achieved. That is, from the time trend perspective, the PMSWG has and will continue to grow for a long time.
Therefore, it makes more sense to test the impact of the source separation program on the relative growth of PMSWG. Next, we still adopt the same spatial weighting matrices above to control the spatial interdependence, and we replace the dependent variable with the growth rate of the PMSWG but only consider the linear relationship between it and the proxy variable of economic growth (lnPGDP), i.e., high-order terms of economic variables are excluded from the regression model. In this way, the effect of the second MSW source separation policy on the growth rate of the PMSWG is shown in Table 6.
The estimation results show that the estimated coefficients of the spatially lagged dependent variable ( ρ ), the economic variable ( δ ), and the error term ( σ ) are all significant at least at the level of 5%. The significant estimation coefficients of the spatial lag terms in Table 6 suggest that the GNS is appropriate to control the spatial correlations, but the growth rate of the PMSWG has a negative spillover effect on spatially related cities, which is in contrast to the estimation in Table 5. In addition, compared with the estimation in Table 5, the spatial spillover effect of the explanatory variable GDP per capita becomes significantly negative, namely, the increase of the GDP per capita in city i will curb the growth rate of the PMSWG in the spatial neighboring cities. Thus, obtaining wealth may be one of the best ways we can find to improve our environment today, the model fitting results are consistent with the theoretical derivation conclusion of Boucekkine and El Ouardighi (2016) [67].
The results also indicate that the higher the probability of implementing the second MSW source separation ( P R ^ ), the lower the growth rate of the PMSWG, and the estimates are very close to each other in the specification of all the three weighting matrices. By these estimates, when the change in the probability of the second MSW source separation is increased by 1%, the growth rate of the PMSWG will decrease by approximately 0.06% per person per year. However, the probability of implementing the second MSW source separation cannot be observed in reality and it is an instrumental variable for correcting endogenous dummy policy variable (SEP2=1, if treated, SEP2=0, otherwise) in the S2SLS regression context. We only observe whether the city has implemented the policy (SEP2=1) or not (SEP2=0), there is no middle ground in reality. In other words, compared with cities that did not implement the separation program, the enforcement of the second MSW source separation will result in a 5.79% decrease in the growth rate of the PMSWG, and it is a considerable value for the relative reduction effect of separation program. Thus, from the analysis, the second source separation program significantly reduced the growth rate of the PMSWG. Combined with the above analysis, the effect of the source separation program on the absolute amount of the PMSWG cannot be estimated steadily and consistently in all spatial models. However, it can still correct the negative bias due to ignoring the spatial dependences. Meanwhile, all three spatial models support that source separation policy significantly reducing the relative amount of PMSWG.
In the policy effect evaluation of source separation on the growth rate of PMSWG, the estimated coefficients of population density become no longer significant in all three models, and the estimated coefficients of the urbanization rate and urban population are still significantly positive and negative, respectively, at the level of 1%. The cross-product term of the GDP per capita and the urbanization rate, as well as the cross product term of the GDP per capita and urban population, are still significantly negative and positive at the level of at least 1% in the three spatial models, which indicates that the impact of urbanization rate and urban population on the growth rate of PMSWG depends on GDP per capita. In the case of the same rise of urbanization rate, the increase of GDP per capita will inhibit the growth rate of the PMSWG, while the identical rise in GDP per capita will lead to the increase of the growth rate of the PMSWG in the case of the same level of urban population increasing.
The reason may be the law of diminishing margins, which is the potential basis for the so-called Environmental Kuznets Curve [68], which refers to an inverted U-shaped relationship between environmental pollutants and economic growth indicators, following the pioneering work proposed by Kuznets (1955) [69]. Different cities may be at various stages of development in the Environmental Kuznets Curve [70]. In addition, with the same GDP per capita, the growth rate of the PMSWG in provincial cities is lower than that in nonprovincial cities. Furthermore, the effect of the MSW charging policy is still not significant in all models, which means the failure of the price policy of the MSWM in China.
Because the spatial econometric regression model explores the complex spatial dependency structure among spatial units, changing an explanatory variable or dependent variable of a particular spatial unit will affect the spatial unit itself, on the one hand, and all other spatial correlated units, on the other hand. This mutual spatial dependence will produce a feedback effect. Therefore, it will lead to a considerable bias when explaining the relationship between spatial units directly, according to the estimated parameters of the spatial regression model. Table 7 shows the average estimators of the direct effect, indirect effect, and total effect estimated by the Delta-Method under the specification of the IDIS matrix, and only the significant variables are listed.
The indirect effects of the predicted probability of the second MSW source separation program on the growth rate of the PMSWG are significantly negative under the three weighting matrices settings. This means that after correcting the endogenous local policy, the source separation policy has a negative spillover effect on the PMSWG. Take the contiguity spatial matrix as an example, when the probability of implementing the second MSW source separation program increases by 1%, the growth rate of the PMSWG in this city will decrease by 0.06%, and its spillover effect will lead to a decrease of 0.03% in the growth rate of the PMSWG in spatially related cities. Specifically, 31% of the total reduction effect is due to the spatial spillover effect.
In the effect evaluation of the source separation policy on the growth rate of the PMSWG, the estimated coefficients of the population density are no longer significant and are not reported in Table 7. Additionally, the direct effects and indirect effects of the GDP per capita, urban population and urbanization rate are still significant at the level of 1% in all the specifications of the spatial weighting matrices. Take the contiguity spatial weighting matrix, for example, among the total effects of one percent growth of GDP per capita, 55% of the negative effects are the further spillover effects on spatially dependent cities. Approximately 30% of the total effect of the increase of the urban population and urbanization rate are the negative and positive spillover effects, respectively.

5. Conclusions

Relaxing the SUTVA and correcting for endogenous local policy choices, we construct the GNS model to estimate the PMSWG reduction effect of China’s second MSW source separation program. This study contributes to the empirical literature by evaluating the impacts of endogenous policies with consideration for spatial correlations that cannot be ignored.
The analysis indicates that comparing with the results estimated from the models with SUTVA, the spatial estimation increases the significance of the impact of the second MSW source separation program under the specification of the IDIS matrix. That is, the reduction effect of the program on the absolute amount of PMSWG was underestimated with SUTVA. Although not all the spatial models support this absolute amount of PMSWG reduction, the impact on the relative amount of PMSWG reduction is significant and stable in alternative specifications of spatial weighting matrices.
Given that the absolute MSWG growth in China currently remains high, it is difficult to achieve the absolute decoupling between MSWG and economic growth in the short-run. Seeking feasible policies to curb the growth of PMSWG has great practical significance for China’s MSWM. According to the analysis of this study, the reduction effect of the source separation program on the growth rate of PMSWG is significant but still not satisfactory. Specifically, the estimates indicate a relative reduction rather than an absolute reduction effect of source separation program, and there is much room for policy improvement. In light of these findings, an urgent field for future research is to explore policy combinations to enhance the effectiveness of the source separation. It may be necessary to reinforce coordination in the policy implementation process and the whole chain of the MSWM. This study provides a quantitative factual basis for future policy revision and implementation.
With regard to the confirmation of nonnegligible positive spatial spillover effects between PMSWG at the city level in China, this study also encourages more researchers in developing counties to examine the existence of spatial dependence at different levels, such as the country level. If this interaction exists, the enhancement of the effectiveness of MSWM in a certain region not only can protect the environment of its own by reducing the PMSWG but also will have spillover effects on all other spatial-related regions. It also provides insights for further cooperation in MSWM among developing countries. However, the idiosyncratic characteristics of each county should be considered when transferring the conclusion of this paper.

Author Contributions

Conceptualization, L.Z. and J.Z.; methodology, J.Z.; software, J.Z.; validation, L.Z., J.Z. and Z.Z.; formal analysis, J.Z.; investigation, Z.Z.; resources, Z.Z.; data curation, Z.Z.; writing—original draft preparation, J.Z.; writing—review and editing, Z.Z.; visualization, J.Z.; supervision, L.Z.; project administration, L.Z.; funding acquisition, L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (Grant No. 71773114 and 71803175), the Zhejiang Social Sciences Federation (Grant No. 2019N32), the Natural Science Foundation of Zhejiang Province, China (Grant No. LY20G030002), and Education Department of Zhejiang Province, China (Grant No. Y201941978).

Acknowledgments

We are grateful for the generous help from Xiaojie Huang and Professor Xueyuan Wang. We owe much gratitude to Youzhen Lu and Fuxing Zou for their unselfish assistance and also appreciate the valuable comments by the two anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fenton, R.; Hanley, N. Economic Instruments and Waste Minimization: The Need for Discard-Relevant and Purchase-Relevant Instruments. Environ. Plan. A 1995, 27, 1317–1328. [Google Scholar] [CrossRef]
  2. Kinnaman, T.C. The economics of municipal solid waste management. Waste Manag. 2009, 29, 2615–2617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jenkins, R.R. The Economics of Solid Waste Reduction; Edward Elgar Publishing Limited: Hants, UK, 1993. [Google Scholar]
  4. Fullerton, D. Wenbo Wu Policies for Green Design. J. Environ. Econ. Manag. 1998, 36, 131–148. [Google Scholar] [CrossRef]
  5. Miranda, M.L.; Aldy, J.E. Unit pricing of residential municipal solid waste: Lessons from nine case study communities. J. Environ. Manag. 1998, 52, 79–93. [Google Scholar] [CrossRef]
  6. Fullerton, D.; Kinnaman, T.C. Household Responses to Pricing Garbage by the Bag. Am. Econ. Rev. 1996, 86, 971–984. [Google Scholar]
  7. Dinan, T.M. Economic efficiency effects of alternative policies for reducing waste disposal. J. Environ. Econ. Manag. 1993, 25, 242–256. [Google Scholar] [CrossRef]
  8. Palmer, K.; Walls, M. Optimal policies for solid waste disposal taxes, subsidies, and standards. J. Public Econ. 1997, 65, 193–205. [Google Scholar] [CrossRef]
  9. Tai, J.; Zhang, W.; Che, Y.; Feng, D. Municipal solid waste source-separated collection in China: A comparative analysis. Waste Manag. 2011, 31, 1673–1682. [Google Scholar] [CrossRef]
  10. Rousta, K.; Bolton, K.; Lundin, M.; Dahlén, L. Quantitative assessment of distance to collection point and improved sorting information on source separation of household waste. Waste Manag. 2015, 40, 22–30. [Google Scholar] [CrossRef]
  11. Moh, Y.C.; Abd Manaf, L. Solid waste management transformation and future challenges of source separation and recycling practice in Malaysia. Resour. Conserv. Recycl. 2017, 116, 1–14. [Google Scholar] [CrossRef]
  12. Owusu, V.; Adjei-Addo, E.; Sundberg, C. Do economic incentives affect attitudes to solid waste source separation? Evidence from Ghana. Resour. Conserv. Recycl. 2013, 78, 115–123. [Google Scholar] [CrossRef]
  13. Sukholthaman, P.; Sharp, A. A system dynamics model to evaluate effects of source separation of municipal solid waste management: A case of Bangkok, Thailand. Waste Manag. 2016, 52, 50–61. [Google Scholar] [CrossRef] [PubMed]
  14. Hoornweg, D.; Lam, P.; Chaudhry, M. Waste Management in China: Issues and Recommendations; World Bank: Washington, DC, USA, 2005. [Google Scholar]
  15. Mian, M.M.; Zeng, X.; Nasry, A.A.; Al-Hamadani, S.M. Municipal solid waste management in China: A comparative analysis. J. Mater. Cycles Waste Manag. 2017, 19, 1127–1135. [Google Scholar] [CrossRef]
  16. Jambeck, J.R.; Geyer, R.; Wilcox, C.; Siegler, T.R.; Perryman, M.; Andrady, A.; Narayan, R.; Law, K.L. Plastic waste inputs from land into the ocean. Science 2015, 347, 768–771. [Google Scholar] [CrossRef] [PubMed]
  17. Han, H.; Zhang, Z.; Xia, S. The crowding-out effects of garbage fees and voluntary source separation programs on waste reduction: Evidence from China. Sustainability 2016, 8, 678. [Google Scholar] [CrossRef] [Green Version]
  18. Fullterton, D.; Kinnaman, T.C. Garbage, Recycling, and Ilicit Burning or Dumping. J. Environ. Econ. Manag. 1995, 29, 78–91. [Google Scholar] [CrossRef] [Green Version]
  19. Han, H.; Zhang, Z.; Xia, S.; Li, H. The carrot or the stick: Individual adaption against varying institutional arrangements. J. Environ. Plan. Manag. 2018, 61, 568–596. [Google Scholar] [CrossRef]
  20. Cerulli, G. Econometric Evaluation of Socio-Economic Programs; Springer: Berlin/Heidelberg, Germany, 2015; Volume 49, ISBN 978-3-662-46404-5. [Google Scholar]
  21. Besley, T.; Case, A. Unnatural Experiments? Estimating the Incidence of Endogenous Policies. Econ. J. 2000, 110, 672–694. [Google Scholar] [CrossRef] [Green Version]
  22. Ministry of Housing and Urabn-Rural Develpment. Announcement of the Pilot Cities for the Source-Separation Collection of Municipal Solid Waste. 2000. Available online: http://www.mohurd.gov.cn/wjfb/200611/t20061101_156932.html (accessed on 21 February 2020).
  23. Ministry of Housing and Urban-Rural Development. Separation Method And Evaluation Criteria Of Municipal Solid Waste. 2004. Available online: http://www.mohurd.gov.cn/wjfb/200610/t20061031_155820.html (accessed on 21 February 2020).
  24. Ministry of Housing and Urban-Rural Development; National Development and Reform Commission; Ministry of Finance; Ministry of Environmental Protection; Ministry of Commerce. Announcement of the First Batch of Demonstration Cities (Districts) for Municipal Solid Waste Source Separation. 2015. Available online: https://ebook.chinabuilding.com.cn/zbooklib/book/detail/show?bookID=92441&SiteID=1 (accessed on 21 February 2020).
  25. Hoornweg, D.; Bhada-Tata, P. What a Waste: A Global Review of Solid Waste Management. 2012. Available online: http://documents.worldbank.org/curated/en/302341468126264791/What-a-waste-a-global-review-of-solid-waste-management (accessed on 21 February 2020).
  26. Ministry of Housing and Urban-Rural Development; National Development and Reform Commission; Ministry of Finance; Ministry of Environmental Protection; Ministry of Commerce. Announcement of the Ministry of Housing and Urban-Rural Development and Other Departments on the Application Work of Demonstration Cities (Districts) for Municipal Solid Waste Separation. 2014. Available online: http://www.mohurd.gov.cn/wjfb/201404/t20140425_217771.html (accessed on 21 February 2020).
  27. Ministry of Housing and Urban-Rural Development; National Development and Reform Commission. Implementation Scheme of Municipal Solid Waste Separation System. 2017. Available online: http://www.gov.cn/zhengce/content/2017-03/30/content_5182124.htm (accessed on 21 February 2020).
  28. Viscusi, W.K.; Huber, J.; Bell, J. Promoting Recycling: Private Values, Social Norms, and Economic Incentives. Am. Econ. Rev. 2011, 101, 65–70. [Google Scholar] [CrossRef] [Green Version]
  29. Wei, Y.; Li, J.; Shi, D.; Liu, G.; Zhao, Y.; Shimaoka, T. Environmental challenges impeding the composting of biodegradable municipal solid waste: A critical review. Resour. Conserv. Recycl. 2017, 122, 51–65. [Google Scholar] [CrossRef] [Green Version]
  30. Aphale, O.; Thyberg, K.L.; Tonjes, D.J. Differences in waste generation, waste composition, and source separation across three waste districts in a New York suburb. Resour. Conserv. Recycl. 2015, 99, 19–28. [Google Scholar] [CrossRef]
  31. Han, H.; Zhang, Z. The impact of the policy of municipal solid waste source-separated collection on waste reduction: A case study of China. J. Mater. Cycles Waste Manag. 2017, 19, 382–393. [Google Scholar] [CrossRef]
  32. Rubin, D.B. Bayesian Inference for Causal Effects. Ann. Stat. 1978, 6, 34–58. [Google Scholar] [CrossRef]
  33. Nichols, A. Causal inference with observational data. Stata J. 2007, 7, 507–541. [Google Scholar] [CrossRef] [Green Version]
  34. Fowlie, M.; Holland, S.P.; Mansur, E.T. What do emissions markets deliver and to whom? Evidence from Southern California’s NO X trading program. Am. Econ. Rev. 2012, 102, 965–993. [Google Scholar] [CrossRef] [Green Version]
  35. Anselin, L. Spatial effects in econometric practice in environmental and resource economics. Am. J. Agric. Econ. 2001, 83, 705–710. [Google Scholar] [CrossRef]
  36. Tobler, W.R. A Computer Movie Simulating Urban Growth in the Detroit Region. Econ. Geogr. 1970, 46, 234–240. [Google Scholar] [CrossRef]
  37. Agovino, M.; D’Uva, M.; Garofalo, A.; Marchesano, K. Waste management performance in Italian provinces: Efficiency and spatial effects of local governments and citizen action. Ecol. Indic. 2018, 89, 680–695. [Google Scholar] [CrossRef]
  38. Hage, O.; Sandberg, K.; Söderholm, P.; Berglund, C. The regional heterogeneity of household recycling: A spatial-econometric analysis of Swedish plastic packing waste. Lett. Spat. Resour. Sci. 2018, 11, 245–267. [Google Scholar] [CrossRef] [Green Version]
  39. Gui, S.; Zhao, L.; Zhang, Z. Does municipal solid waste generation in China support the Environmental Kuznets Curve? New evidence from spatial linkage analysis. Waste Manag. 2019, 84, 310–319. [Google Scholar] [CrossRef]
  40. Truelove, H.B.; Carrico, A.R.; Weber, E.U.; Raimi, K.T.; Vandenbergh, M.P. Positive and negative spillover of pro-environmental behavior: An integrative review and theoretical framework. Glob. Environ. Chang. 2014, 29, 127–138. [Google Scholar] [CrossRef]
  41. Agovino, M.; Crociata, A.; Sacco, P.L. Location matters for pro-environmental behavior: A spatial Markov Chains approach to proximity effects in differentiated waste collection. Ann. Reg. Sci. 2016, 56, 295–315. [Google Scholar] [CrossRef]
  42. Crociata, A.; Agovino, M.; Sacco, P.L. Neighborhood effects and pro-environmental behavior: The case of Italian separate waste collection. J. Clean. Prod. 2016, 135, 80–89. [Google Scholar] [CrossRef]
  43. Oskamp, S.; Harrington, M.J.; Edwards, T.C.; Sherwood, D.L.; Okuda, S.M.; Swanson, D.C. Factors influencing household recycling behavior. Environ. Behav. 1991, 23, 494–519. [Google Scholar] [CrossRef]
  44. Tian, L.; Wang, H.H.; Chen, Y. Spatial externalities in China regional economic growth. China Econ. Rev. 2010, 21, S20–S31. [Google Scholar] [CrossRef]
  45. Revelli, F. On spatial public finance empirics. Int. Tax Public Financ. 2005, 12, 475–492. [Google Scholar] [CrossRef]
  46. Elhorst, P.; Piras, G.; Arbia, G. Growth and convergence in a multiregional model with space-time dynamics. Geogr. Anal. 2010, 42, 338–355. [Google Scholar] [CrossRef]
  47. Moscone, F.; Tosetti, E.; Vittadini, G. Social interaction in patients’ hospital choice: Evidence from Italy. J. R. Stat. Soc. A 2012, 175, 453–472. [Google Scholar] [CrossRef] [Green Version]
  48. Moscone, F.; Knapp, M. Exploring the spatial pattern of mental health expenditure. J. Ment. Health Policy Econ. 2005, 8, 205–217. [Google Scholar]
  49. Tiebout, C.M. A Pure Theory of Local Expenditures. J. Polit. Econ. 1956, 64, 416–424. [Google Scholar] [CrossRef]
  50. Ministry of Natural Resources. Shapefile of Prefecture-Level Cities. Available online: http://bzdt.ch.mnr.gov.cn/ (accessed on 21 February 2020).
  51. National Bureau of Statistics of China. China Stiatistical Yearbook,2002–2016; China Statistics Press: Beijing, China, 2002–2016.
  52. National Bureau of Statistics of China. China Urban Construction Statistical Yearbook, 2002–2016; China Statistics Press: Beijing, China, 2002–2016.
  53. Train, K.E. Discrete Choice Methods with Simulation. 2009. Available online: https://eml.berkeley.edu/books/train1201.pdf (accessed on 21 February 2020).
  54. Moran, P.A.P. Notes on Continuous Stochastic Phenomena. Biometrika 1950, 37, 17–23. [Google Scholar] [CrossRef] [PubMed]
  55. Hao, Y.; Liu, Y.; Weng, J.-H.; Gao, Y. Does the Environmental Kuznets Curve for coal consumption in China exist? New evidence from spatial econometric analysis. Energy 2016, 114, 1214–1223. [Google Scholar] [CrossRef]
  56. Sobhee, K.S. The environmental Kuznets curve (EKC): A logistic curve? Appl. Econ. Lett. 2004, 11, 449–452. [Google Scholar] [CrossRef]
  57. Johnstone, N.; Labonne, J. Generation of Household Solid Waste in OECD Countries: An Empirical Analysis Using Macroeconomic Data. Land Econ. 2004, 80, 529–538. [Google Scholar] [CrossRef]
  58. Kinnaman, T.C.; Fullerton, D. Garbage and Recycling with Endogenous Local Policy. J. Urban Econ. 2000, 48, 419–442. [Google Scholar] [CrossRef]
  59. Mazzanti, M.; Zoboli, R. Waste generation, waste disposal and policy effectiveness: Evidence on decoupling from the European Union. Resour. Conserv. Recycl. 2008, 52, 1221–1234. [Google Scholar] [CrossRef]
  60. Lee, L.F.; Yu, J. Estimation of spatial autoregressive panel data models with fixed effects. J. Econom. 2010, 154, 165–185. [Google Scholar] [CrossRef]
  61. Berrens, R.P.; Bohara, A.K.; Gawande, K. Testing the inverted-U hypothesis for US hazardous waste. Econ. Lett. 1997, 55, 435–440. [Google Scholar] [CrossRef]
  62. Bloom, D.E.; Beede, D.N. The economics of municipal solid waste. World Bank Res. Obs. 1995, 10, 113–150. [Google Scholar]
  63. Cole, M.A.; Rayner, A.J.; Bates, J.M. The Environmental Kuznets Curve: An Empirical Analysis. Environ. Dev. Econ. 1997, 2, 401–416. [Google Scholar] [CrossRef]
  64. Fischer-Kowalski, M.; Amann, C. Beyond IPAT and Kuznets curves: Globalization as a vital factor in analysing the environmental impact of socio-economic metabolism. Popul. Environ. 2001, 23, 7–47. [Google Scholar] [CrossRef]
  65. Panayotou, T. Empirical Tests and Policy Analysis of Environmental Degradation at Different Stages of Economic Development. 1993. Available online: http://www.ilo.org/public/libdoc/ilo/1993/93B09_31_engl.pdf (accessed on 21 February 2020).
  66. Wang, P.; Bohara, A.K.; Berrens, R.P.; Gawande, K. A risk-based environmental Kuznets curve for US hazardous waste sites. Appl. Econ. Lett. 1998, 5, 761–763. [Google Scholar] [CrossRef]
  67. Boucekkine, R.; El Ouardighi, F. Optimal Growth with Polluting Waste and Recycling. In Dynamic Perspectives on Managerial Decision Making; Springer: Berlin/Heidelberg, Germany, 2016; pp. 109–126. [Google Scholar]
  68. Grossman, G.M.; Krueger, A.B. Economic Growth and the Environment. Q. J. Econ. 1995, 110, 353–377. [Google Scholar] [CrossRef] [Green Version]
  69. Kuznets, S. Economic growth and income inequality. Am. Econ. Rev. 1955, 45, 1–28. [Google Scholar]
  70. OECD. Indicators To Measure Decoupling of Environmental Pressure From Economic Growth. 2002. Available online: http://www.oecd.org/environment/indicators-modelling-outlooks/1933638.pdf (accessed on 21 February 2020).
Figure 1. Choropleth maps of per capita municipal solid waste generation (PMSWG) in China. (a) Shows the PMSWG of prefecture-level cities in China in 2002. (b) shows the year 2016, and (c) shows the average amount of PMSWG from 2002 to 2016.
Figure 1. Choropleth maps of per capita municipal solid waste generation (PMSWG) in China. (a) Shows the PMSWG of prefecture-level cities in China in 2002. (b) shows the year 2016, and (c) shows the average amount of PMSWG from 2002 to 2016.
Sustainability 12 01664 g001
Figure 2. Moran scatter plots of the PMSWG of the prefecture-level cities in China. (ac) show the Moran’s I indexes and the Moran scatter plots of the PMSWG of 288 prefecture-level cities in China in 2002, 2016 and the average PMSWG from 2002 to 2016, respectively, under the setting of spatial matrix W(1). The similar calculations and graphical representations under W(2) are illustrated in (df).
Figure 2. Moran scatter plots of the PMSWG of the prefecture-level cities in China. (ac) show the Moran’s I indexes and the Moran scatter plots of the PMSWG of 288 prefecture-level cities in China in 2002, 2016 and the average PMSWG from 2002 to 2016, respectively, under the setting of spatial matrix W(1). The similar calculations and graphical representations under W(2) are illustrated in (df).
Sustainability 12 01664 g002
Table 1. Descriptive statistics.
Table 1. Descriptive statistics.
VariableDefinition of VariableMeanSDMinMax
Dependent variable
PMSWGKilograms per capita per year of MSW270.741264.8481.4684818.248
City characteristic variables
PGDPGDP per capita (RMB yuan)30,793.4425,198.343430.049125,000
TOTFUrban maintenance and construction fund (RMB 10,000 yuan)241,000541,00021173,600,000
TOTEExpenditures for urban maintenance and construction fund (RMB 10,000 yuan)202,000430,00020452,800,000
SAINVFixed assets investment in city appearance and sanitation (RMB 10,000 yuan)6371.65614,954.3624108,000
BARThe ratio of built-up area in the total urban area (%)0.3650.2880.0031
VEHNumber of vehicles and equipment for municipal environmental sanitation (unit)258.491504.877103478
HLSNumber of Harmless Treatment Plants/ Grounds (unit)1.6691.808027
HLCHarmless Treatment Capacity per year (10,000 ton per year)41.796268.14080888.4465
Demographic variables
POPDThe number of persons per square kilometers3393.3322812.3632720,093
URBPUrban population (10000 person)231.214250.87491.7782420.46
URBRUrbanization rate (%)0.540.2830.0081
Policy variables
SEP21 if the city implements the second MSW source separation (0 otherwise)0.080.27101
FEEMSW disposal fee (RMB 10,000 yuan)1579.7733728.381926,524
Location variables
PROCAP1 if the city is the provincial capital (0 otherwise)0.1080.3101
EAST1 if the city is in the east of China (0 otherwise)0.3780.48501
MIDDLE1 if the city is in the middle of China (0 otherwise)0.2190.41301
WEST1 if the city is in the west of China (0 otherwise)0.2740.44601
NORTHEAST1 if the city is in the northeast of China (0 otherwise)0.1280.33501
The mean values of some characteristic variables between 2002 and 2016 in the treated and control cities are shown in Table 2.
Table 2. Comparison of city characteristic variables between the treated and the control group.
Table 2. Comparison of city characteristic variables between the treated and the control group.
SEP2lnMSWGlnPMSWGPGDPPOPDURBR
03.210.2229,380.633410.720.53
14.360.2447,071.413192.940.68
Total3.310.2230,793.433393.330.54
Table 3. Estimation of the probability of adopting the second MSW source separation program.
Table 3. Estimation of the probability of adopting the second MSW source separation program.
VariablesLogit
CoefficientMarginal Effects
lnVEH0.2780 **(0.1292)0.0145
lnHLC0.0059 ***(0.0017)0.0003
HLS0.1483 ***(0.0522)0.0077
lnSAINV−0.1097 **(0.0553)−0.0057
lnTOTF−0.5049 ***(0.1862)−0.2627
lnTOTE0.6144 ***(0.1894)0.0319
lnPGDP−0.0103(0.1443)−0.0005
BAR−0.9594 ***(0.2883)−0.0499
PROCAP1.7094 ***(0.2132)0.0889
CONS−4.6793 ***(1.1469)
Observations3456
Pseudo R20.2974
LR chi2572.3811
Notes: Standard errors are in parentheses. *, **, and *** indicate significance at the 0.10, 0.05, and 0.01 levels, respectively.
Table 4. Estimation results of the program’s impact on PMSWG with SUTVA.
Table 4. Estimation results of the program’s impact on PMSWG with SUTVA.
VariablesIndividual-Fixed EffectTwo-Way-Fixed EffectRandom Effect
P R ^ −0.1046(0.0837)−0.1275(0.0829)0.0400(0.0815)
lnPGDP6.2097 ***(1.7368)5.8375 ***(1.7520)5.9148 ***(1.7687)
(lnPGDP)2−0.7218 ***(0.1746)−0.6685 ***(0.1755)−0.6688 ***(0.1772)
(lnPGDP)30.0263 ***(0.0058)0.0243 ***(0.0059)0.0241 ***(0.0059)
lnPOPD−0.0377 ***(0.0065)−0.0236 ***(0.0067)−0.0155 **(0.0067)
URBR2.0456 ***(0.4011)1.5815 ***(0.4072)0.6766 *(0.4006)
lnURBP−1.7039 ***(0.0833)−1.5459 ***(0.0848)−1.5220 ***(0.0856)
lnFEE−0.0061(0.0046)−0.0048(0.0046)0.0035(0.0046)
lnPGDP×URBR−0.2445 ***(0.0394)−0.1949 ***(0.0398)−0.1194 ***(0.0396)
lnPGDP×lnURBP0.0839 ***(0.0087)0.0734 ***(0.0088)0.0810 ***(0.0088)
lnPGDP×PROCAP0.1036 ***(0.0269)0.1332 ***(0.0271)0.1046 ***(0.0274)
lnPGDP×EAST0.2232 ***(0.0235)0.2399 ***(0.0234)0.2491 ***(0.0234)
lnPGDP×MIDDLE0.1203 ***(0.0233)0.1150 ***(0.0230)0.1100 ***(0.0232)
lnPGDP×WEST0.2124 ***(0.0227)0.2066 ***(0.0225)0.1995 ***(0.0226)
PROCAP 0.0890(0.2898)
EAST −2.7992 ***(0.2451)
MIDDLE −1.5068 ***(0.2423)
WEST −2.4956 ***(0.2375)
CONS−7.9839(5.7293)−7.9787(5.8108)−6.9033(5.8653)
Observations432043204320
R-Square0.30780.31170.6880
Notes: Standard errors are in parentheses. *, **, and *** indicate significance at the 0.10, 0.05, and 0.01 levels, respectively.
Table 5. Estimation results of the program’s impact on PMSWG by GNS models.
Table 5. Estimation results of the program’s impact on PMSWG by GNS models.
VariablesCON-MatrixIDIS-MatrixIDIS-CON-Matrix
P R ^ −0.1309(0.0817)−0.1427 *(0.0812)−0.1148(0.0821)
lnPGDP6.7434 ***(1.6908)5.2242 ***(1.7507)6.5504 ***(1.6892)
(lnPGDP)2−0.7510 ***(0.1694)−0.5874 ***(0.1757)−0.7269 ***(0.1693)
(lnPGDP)30.0269 ***(0.0057)0.0214 ***(0.0059)0.0260 ***(0.0057)
lnPOPD−0.0212 ***(0.0068)−0.0230 ***(0.0067)−0.0200 ***(0.0068)
URBR1.0722 **(0.4231)0.9026 **(0.4073)0.8721 **(0.4162)
lnURBP−1.2972 ***(0.0954)−1.2747 ***(0.0924)−1.2589 ***(0.0927)
lnFEE−0.0048(0.0045)−0.0036(0.0045)−0.0043(0.0045)
lnPGDP × URBR−0.1361 ***(0.0417)−0.1192 ***(0.0400)−0.1169 ***(0.0409)
lnPGDP × lnURBP0.0457 ***(0.0102)0.0437 ***(0.0097)0.0421 ***(0.0098)
lnPGDP × PROCAP0.1661 ***(0.0267)0.1779 ***(0.0267)0.1657 ***(0.0267)
lnPGDP × EAST0.2495 ***(0.0262)0.1409 ***(0.0298)0.2347 ***(0.0256)
lnPGDP × MIDDLE0.0881 ***(0.0265)−0.0306(0.0326)0.0741 ***(0.0262)
lnPGDP × WEST0.1817 ***(0.0256)0.1277 ***(0.0266)0.1889 ***(0.0243)
CONS0.2479 ***(0.0028)0.2462 ***(0.0028)0.2475 ***(0.0028)
W
ρ 0.1271 ***(0.0261)0.4907 ***(0.0730)0.2280 ***(0.0324)
δ 0.0690 ***(0.0229)0.3163 ***(0.0648)0.0787 ***(0.0291)
σ 0.2192 ***(0.0376)0.8103 ***(0.0510)0.2553 ***(0.0456)
Observations432043204320
Pseudo R-Square0.10470.08210.1014
Notes: Standard errors are in parentheses. *, **, and *** indicate significance at the 0.10, 0.05, and 0.01 levels, respectively.
Table 6. Estimation results of the program’s impact on the growth rate of the PMSWG by GNS models.
Table 6. Estimation results of the program’s impact on the growth rate of the PMSWG by GNS models.
VariablesCON-MatrixIDIS-MatrixIDIS-CON-Matrix
P R ^ −0.0579 *(0.0333)−0.0583 *(0.0343)−0.0585 *(0.0341)
lnPGDP−0.0698 ***(0.0177)−0.0983 ***(0.0203)−0.0845 ***(0.0189)
lnPOPD−0.0017(0.0026)−0.0024(0.0028)−0.0017(0.0028)
URBR0.4894 ***(0.1407)0.5997 ***(0.1439)0.5852 ***(0.1431)
lnURBP−0.2985 ***(0.0388)−0.3811 ***(0.0394)−0.3630 ***(0.0386)
lnFEE0.0008(0.0018)0.0009(0.0019)0.0009(0.0019)
lnPGDP × URBR−0.0518 ***(0.0138)−0.0646 ***(0.0140)−0.0626 ***(0.0140)
lnPGDP × lnURBP0.0259 ***(0.0041)0.0344 ***(0.0041)0.0324 ***(0.0040)
lnPGDP × PROCAP−0.0301 ***(0.0113)−0.0406 ***(0.0114)−0.0372 ***(0.0114)
lnPGDP × EAST0.0035(0.0085)0.0195 *(0.0117)0.0069(0.0097)
lnPGDP × MIDDLE0.0040(0.0083)0.0236 *(0.0129)0.0093(0.0100)
lnPGDP × WEST−0.0004(0.0082)0.0060(0.0102)−0.0013(0.0094)
W
ρ −0.0042 *(0.0097)−0.0569 ***(0.0186)−0.0290 **(0.0124)
δ −0.3498 ***(0.0541)−0.3276 ***(0.0935)−0.1882 ***(0.0611)
σ 0.2498 ***(0.0727)0.0800 **(0.0351)0.0277 *(0.0697)
CONS0.0940 ***(0.0012)0.0952 ***(0.0011)0.0952 ***(0.0011)
Observations403240324032
Pseudo R-Square0.01820.00920.0169
Notes: Standard errors are in parentheses. *, **, and *** indicate significance at the 0.10, 0.05, and 0.01 levels, respectively.
Table 7. The direct, indirect and total effects under three spatial weighting matrices.
Table 7. The direct, indirect and total effects under three spatial weighting matrices.
VariablesCON MatrixIDIS MatrixIDIS-CON Matrix
DirectIndirectTotalDirectIndirectTotalDirectIndirectTotal
P R ^ −0.0584 *−0.0259 *−0.0842 *−0.0589 *−0.0217 *−0.0806 *−0.0587 *−0.0072 *−0.0659 *
(0.0343)(0.0189)(0.0510)(0.0339)(0.0135)(0.0468)(0.0342)(0.0050)(0.0386)
lnPGDP−0.0986 ***−0.1208 ***−0.2194 ***−0.0712 ***−0.0307 ***−0.1019 ***−0.0852 ***−0.0292 ***−0.1144 ***
(0.0204)(0.0431)(0.0545)(0.0180)(0.0143)(0.0276)(0.0190)(0.0103)(0.0242)
URBR0.6003 ***0.2663 ***0.8666 ***0.4981 ***0.1835 ***0.6816 ***0.5871 ***0.0717 ***0.6588 ***
(0.1440)(0.1293)(0.2367)(0.1432)(0.0667)(0.2008)(0.1436)(0.0317)(0.1637)
lnURBPOP−0.3814 ***−0.1692 ***−0.5506 ***−0.3038 ***−0.1119 ***−0.4158 ***−0.3642 ***−0.0445 ***−0.4087 ***
(0.0395)(0.0730)(0.0901)(0.0392)(0.0264)(0.0557)(0.0386)(0.0166)(0.0452)
Notes: Standard errors are in parentheses. *, **, and *** indicate significance at the 0.10, 0.05, and 0.01 levels, respectively.

Share and Cite

MDPI and ACS Style

Zhao, L.; Zou, J.; Zhang, Z. Does China’s Municipal Solid Waste Source Separation Program Work? Evidence from the Spatial-Two-Stage-Least Squares Models. Sustainability 2020, 12, 1664. https://doi.org/10.3390/su12041664

AMA Style

Zhao L, Zou J, Zhang Z. Does China’s Municipal Solid Waste Source Separation Program Work? Evidence from the Spatial-Two-Stage-Least Squares Models. Sustainability. 2020; 12(4):1664. https://doi.org/10.3390/su12041664

Chicago/Turabian Style

Zhao, Liange, Jianfeng Zou, and Zhijian Zhang. 2020. "Does China’s Municipal Solid Waste Source Separation Program Work? Evidence from the Spatial-Two-Stage-Least Squares Models" Sustainability 12, no. 4: 1664. https://doi.org/10.3390/su12041664

APA Style

Zhao, L., Zou, J., & Zhang, Z. (2020). Does China’s Municipal Solid Waste Source Separation Program Work? Evidence from the Spatial-Two-Stage-Least Squares Models. Sustainability, 12(4), 1664. https://doi.org/10.3390/su12041664

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop