Introduction

A river basin is a complex system characterized by its integrity, regional interconnectivity, and the interactions between its various components [1]. As a typical public good that spans multiple administrative regions, the water resources within these basins present serious risks to human health and the environment when polluted due to their spatial spillover effects and negative externalities [2, 3]. Furthermore, in these regions, industrial activities and local government governance measures are closely intertwined, making the management of river basin water pollution a matter that encompasses the interests of multiple regions and stakeholders, thus creating a complex dynamic interaction [4]. Furthermore, the fragmentation of regional jurisdictions often leads to challenges such as “free-riding” and the “tragedy of the commons” [5-8]. Consequently, promoting inter-regional collaborative governance is widely regarded as a critical element of comprehensive watershed pollution control. However, differences in geographical location, industrial structure, and levels of economic development across cities result in variations in the costs and benefits of pollution control, generating conflicts of interest that hinder cooperation in water pollution management [9]. Therefore, balancing the costs and benefits among regions and stakeholders remains a significant challenge in the process of cross-regional collaborative water pollution governance.

Existing research primarily focuses on the governance interactions between upstream and downstream regions or between the mainstream and tributaries. For example, Yuan et al. [10] first compared the optimal emission levels of multiple pollutants in upstream and downstream regions under both noncooperative and cooperative game scenarios. Chen et al. [11] demonstrated that basin-wide collaborative governance can foster cooperation and resource sharing between upstream and downstream regions, optimize water pollution control, and promote the sustainable economic development of the basin. Song et al. [12] introduced a water quality credit market and developed a stochastic differential game model of upstream and downstream efforts, discussing the optimal feedback equilibrium of watershed environmental quality under three different scenarios. These studies demonstrate that the governance measures of upstream local governments affect both the pollution burden on downstream local governments and the overall water quality. However, such studies typically focus on the strategic interactions within a vertical governance structure, often overlooking the unique interactions that may exist between the local governments on both sides of the river during the governance process. Joint prevention and control on both sides of trans - provincial river basins is a key and difficult issue in water pollution management [13]. The two major challenges driving this governance complexity are “high costs and insufficient internal motivation” and “unclear pollution liability” [14, 15], both of which hinder the effectiveness of transfer payments and reward–punishment mechanisms while weakening external enforcement. Nevertheless, the existing literature provides limited insights into this issue. Therefore, enhancing governance efficiency requires addressing cross-bank cooperation and competition through institutional innovation, cross-regional collaboration, and interest-balancing mechanisms [16, 17].

Scholars have primarily used the evolutionary game model to research the water pollution management problems between government and enterprises. Yang et al. [18] utilizes an evolutionary game model to compare scenarios with and without pollution rights trading and employs stability analysis and sensitivity testing using phase diagram techniques. Wang et al. [19] developed an evolutionary game model involving local governments, enterprises, and the public to analyze the stable strategies for cooperative water pollution control in the Yangtze River Basin. However, due to information asymmetry and the limited rationality of the participants, it is difficult for stakeholders to achieve equilibrium through a one-time decision. In contrast, the differential game model allows for the dynamic representation of the continuous decision-making process of both the government and enterprises at different time points, better capturing the timeliness and policy flexibility inherent in real-world governance [20, 21]. More importantly, transboundary water pollution management involves not only local governments but also industrial enterprises and the central government. Industrial enterprises directly impact pollution levels through their emission reduction efforts [22], while the central government shapes policies, provides subsidies, and sets standards that guide and constrain the actions of local governments and enterprises [23]. However, existing research often overlooks the roles of enterprises and the central government, limiting a comprehensive understanding of the dynamics of water pollution governance.

Existing research on upstream and downstream water pollution and air pollution control provides important policy insights for addressing cross-border water pollution, especially in the context of climate change [24]. The former focuses on the stability of collaborative alliances [25] and the distribution of costs and benefits [26], while the latter examines the effectiveness of policies such as central regulation [27] and goal assessment [28] in regional pollution management. However, the applicability and implementation pathways of transferring these regulatory policies remain to be explored. The two main approaches to ecological compensation are fiscal transfers, where the government compensates affected regions [29, 30], and market-based mechanisms [31, 32], like carbon trading, which incentivize ecological preservation through market transactions. Liu et al. [33] showed that the Ecological Compensation Policy in river basins significantly optimizes industrial structures, while its impact on economic growth and technological advancement remains inconspicuous. Xu et al. [34] concluded that bilateral ecological compensation between upstream and downstream regions most effectively reduces pollutant emissions, making it the optimal form of ecological compensation. However, compared to the ecological compensation mechanism, the cost-sharing approach offers greater adaptability and flexibility. Yang et al. [35] developed a Stackelberg game model based on a cost-sharing mechanism for a Y-shaped river basin. Through numerical analysis, they demonstrated that this mechanism can enhance the effectiveness of pollution elimination, achieving Pareto improvements in both environmental and economic domains.

The main contributions of this paper are as follows. (1) We shift the focus from the traditional upstream–downstream governance model to the largely underexplored cross-bank governance relationship in transboundary river basins. (2) We extend the game-theoretic framework to incorporate not only local governments on both banks but also industrial enterprises and the central government – stakeholders that have often been overlooked in prior studies.

(3) The competition factor between local governments on both sides of the river is considered, capturing how economic development pressures lead to reduced pollution control efforts. (4) The differential game model is used to analyze the pollution control efforts of governments and enterprises from a dynamic perspective.

Materials and Methods

Establishment and Solution of Differential Game Models

Model Assumptions

Assumption 1: Consider a trans-jurisdictional river basin where the left and right banks (GL and GR, respectively) are governed by adjacent but independent local governments, each responsible for water pollution control within their jurisdictions. Meanwhile, the central government (GC) formulates macro-level governance policies for the entire river basin. Industrial production activities occur on both sides of the river, with industrial enterprises (E) discharging wastewater into the river, thereby acting as the primary sources of pollution.

Assumption 2: Industrial production activities on both the left and right banks discharge wastewater containing organic matter and various pollutants into the river. Let the total pollutant emissions from industrial enterprises on the left bank be denoted as qL, qL>0, and those from the right bank as qR, qR>0. The central government may impose an environmental protection tax ω>0 on local governments based on pollutant emissions, aiming to guide their pollution control decisions.

Assumption 3: Let the pollution control efforts of GL, GR, and GC be denoted as uL, uR and uC, respectively, where 0<uL, uR, uC<1, which reflects the policies, personnel, and financial investments of these governments. Additionally, the pollution control efforts of E are represented by uE, 0<uE<1, indicating efforts to reduce pollutant emissions through technological upgrades, end-of-pipe treatment, and clean production methods. Following the established literature (e.g., Gao et al. [36-38]), the cost of pollution control for each stakeholder is represented as a quadratic function of their respective effort levels. Therefore, the cost functions of each participant in the wastewater treatment process are represented as shown in Equation (1).

(1)
C1(t)=12uL2,C2(t)=12kuR2,C3(t)=12cuC2,C4(t)=12euE2

Where k>0 is a coefficient capturing cost discrepancies between the two regions and c>0, e>0 represent the cost coefficients for GC and E, respectively.

Assumption 4: E's utility function at time t can be represented as follows:

(2)
LE(t)=ψuE+δq(t)

Where ψ>0 is the social reputation impact coefficient of wastewater treatment and δ(δ>0) is the external influence coefficient of water pollution control on the firm.

Assumption 5: The reduction of pollutants within the basin results from the collective efforts of GL and GR, GC and E, assuming that their unit pollution control efforts eliminate the following amounts of pollutants, respectively: ξ, η, θ and φ. Formally, the differential equation can be described as follows:

(3)
q˙(t)=ξuL+ηuR+θuC+ϕuE−αq(t)

Where q(t) represents the effectiveness of pollution control, while α denotes the marginal diminishing efficiency of pollution treatment.

Assumption 6: Under economic development pressures, the two local governments on the river’s banks exhibit reluctance to independently undertake high-cost pollution control measures. To capture this effect, we define a competition coefficient β∈[0,1]. A higher value of β indicates a stronger tendency for both bank governments to reduce their pollution abatement efforts in response to economic constraints, reflecting heightened competitive pressures that favor development goals over environmental quality.

Without competition or central government supervision, the baseline pollution control efforts of GL and GR are uL0 and uR0. When competition and oversight are introduced, their efforts become:

(4)
uL=uL0−βuR+d1uC,uR=uR0−βuL+d2uC

Where β∈[0,1] captures the influence of economic pressures on reducing their efforts, and d1, d2>0 represent the central government’s regulatory impact on their actions. For the industrial enterprise, its abatement effort uE depends on local government pressure and policies, and is modeled as:

(5)
uE=uE0+d3(uL+uR)

Similarly, d3>0 represent the extent to which the regulatory intensity of GL and GR influences the industrial enterprise’s behavior, indicating that stronger oversight and more stringent policies lead the enterprise to invest more in pollution reduction.

Assumption 7: Pollution control contributes to overall social welfare as captured by the relation:

(6)
S(t)=S0+λq(t)

Where S0 denotes the initial welfare level, and λ>0 represents the incremental coefficient reflecting how improvements in pollution abatement translate into higher social welfare.

Assumption 8: The conversion of changes in social welfare into respective benefits varies among GL, GR and E. Let πL, πR, πE>0 denote their distinct benefit conversion coefficients, reflecting differing economic development levels and regional conditions. Additionally, all stakeholders employ a common discount rate ρ>0, indicating a uniform temporal valuation of future payoffs across the central government, local governments, and the enterprise.

Non-Cost Sharing Scenario

In the Nash noncooperative game scenario, all parties make independent decisions.

The objective function of GL is JGLN:

(7)
JGLN=∫0∞[ξuLr1−12uL2+πLS(t)−(qL−ξuL)ω+μ1ηuR+χϕuE]e−ρtdt

Here, the parameters μi(i = 1, 2) represent the economic benefits each local government obtains from the other’s pollution-control efforts. χ represents the spillover effect coefficient of enterprise governance. And ri(i = 1, 2, 3, 4) respectively represent the governance marginal utility coefficients of GL, GR, GC and E.

The objective function of GR is JGRN:

(8)
JGRN=∫0∞[ηuRr2−k2uL2+πRS(t)−(qR−ηuR)ω+μ2ξuL+χϕuE]e−ρtdt

The objective function of GC is JGCN:

(9)
JGCN=∫0∞[θuCr3−c2uC2+S(t)+(qL+qR−ξuL−ηuR)ω]e−ρtdt

The objective function of E is JEN:

(10)
JEN=∫0∞[−e2ue2+πES(t)+r4LE(t)]e−ρtdt

Cost-Sharing Scenario

In this case, the central government funds pollution control and allocates cost-sharing ratios (m1, m2, 0<m1 + m2<1) to the left-bank and right-bank regional governments [39, 40], while the regional governments provide subsidies (n1, n2, 0<n1 + n2<1) to incentivize industrial enterprises’ pollution control efforts. At this point, the objective functions of each game participant can be expressed as Equations (11) to (14):

(11)
JGLS=∫0∞[ξuLr1−1−m12uL2−en12uE2+πLS(t)−(qL−ξuL)ω+μ1ηuR+χϕuE]e−ρtdt
(12)
JGRS=∫0∞[ηuRr2−k(1−m2)2uL2−en22uE2+πRS(t)−(qR−ηuR)ω+μ2ξuL+χϕuE]e−ρtdt
(13)
JGCS=∫0∞[θuCr3−c2uC2−12m1uL2−k2m2uR2+S(t)+(qL+qR−ξuL−ηuR)ω]e−ρtdt
(14)
JES=∫0∞[−(1−n1−n2)e2uE2+πES(t)+r4LE(t)]e−ρtdt

Results and Discussion

Non-Cost Sharing Scenario

To obtain the Markov Perfect Equilibrium solution for the non-cost sharing game, a set of bounded, continuous, and differentiable value functions VGLN(q(t)), VGRN(q(t)), VGCN(q(t)) and VEN(q(t)) are constructed. These functions satisfy the Hamilton-Jacobi-Bellman (HJB) equations with respect to the variable q(t), expressed as:

(15)
ρVGLN(q(t))=maxuL≥0{ξuLr1−12uL2+πLS(t)−(qL−ξuL)ω+μ1ηuR+χϕuE+VGLN′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(16)
ρVGRN(q(t))=maxuR≥0{ηuRr2−k2uR2+πRS(t)−(qR−ηuR)ω+μ2ξuL+χϕuE+VGRN′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(17)
ρVGCN(q(t))=maxuC≥0{θuCr3−c2uC2+S(t)+(qL+qR−ξuL−ηuR)ω+VGRN′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(18)
ρVEN(q(t))=maxuE≥0{−e2uE2+πES(t)+r4[ψuE+δq(t)]+VEN′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}

According to the first-order necessary conditions, the first-order partial derivatives of the right-hand side of Equations (15) to (18) with respect to uL, uR, uC and uE are set to zero. The equilibrium solutions of the game can then be derived as follows:

(19)
uLN=ξ(r1+ω)−βμ1η+χϕd3+VGLN′(q(t))(ξ−ηβ+ϕd3)uRN=η(r2+ω)−βμ2ξ+χϕd3+VGRN′(q(t))(η−ξβ+ϕd3)kuCN=(θr3−ω(ξd1−d2β1−β2+ηd2−d1β1−β2)+VGCN′(q(t))(θ+ξd1−d2β1−β2+ηd2−d1β1−β2+ϕd3(d1−d2β1−β2+d2−d1β1−β2)))cuEN=r4ψ+ϕVEN′(q(t))e

The expression for VGLN(q(t)), VGRN(q(t)), VGCN(q(t)), and VEN(q(t)) is given as follows:

(20)
VGLN(q(t))=AGLq(t)+BGL,VGRN(q(t))=AGRq(t)+BGRVGCN(q(t))=AGCq(t)+BGC,VEN(q(t))=AEq(t)+BE

By the substitution method, substituting Equation (21) into Equations (15)-(18) yields the following expressions:

(21)
AGLN=πLλρ+α,AGRN=πRλρ+α,AGCN=λρ+α,AEN=λπE+δr4ρ+α
(22)
BGLN=1ρ[ξuLN(r1+ω)−12(uLN)2+πLS0−qLω+μ1ηuRN+χϕuEN+πLλρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]BGRN=1ρ[ηuRN(r2+ω)−k2(uRN)2+πRS0−qRω+μ2ξuLN+χϕuEN+πRλρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]BGCN=1ρ[θuCNr3−c2(uCN)2+S0+(qL+qR−ξuLN−ηuRN)ω+λρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]BEN=1ρ[−e2(uEN)2+πES0+r4ψuEN+λπE+δr4ρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]

By substituting Equation (21) into Equation (19), the optimal control strategies for each participant can be obtained as follows:

(23)
uLN=ξ(r1+ω)−βμ1η+χϕd3+πLλρ+α(ξ−ηβ+ϕd3)
(24)
uRN=η(r2+ω)−βμ2ξ+χϕd3+πRλρ+α(η−ξβ+ϕd3)k
(25)
uCN=(θr3−ω(ξd1−d2β1−β2+ηd2−d1β1−β2)+λρ+α(θ+ξd1−d2β1−β2+ηd2−d1β1−β2+ϕd3(d1−d2β1−β2+d2−d1β1−β2)))c
(26)
uEN=r4ψ+ϕλπE+δr4ρ+αe

By substituting Equation (21) and (22) into Equation (20), the optimal payoff functions for each participant are obtained as follows:

(27)
VGLN(q(t))=πLλρ+αq(t)+1ρ[ξuLN(r1+ω)−12(uLN)2+πLS0−qLω+μ1ηuRN+χϕuEN+πLλρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]
(28)
VGRN(q(t))=πRλρ+αq(t)+1ρ[ηuRN(r2+ω)−k2(uRN)2+πRS0−qRω+μ2ξuLN+χϕuEN+πRλρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]
(29)
VGCN(q(t))=λρ+αq(t)+1ρ[θuCNr3−c2(uCN)2+S0+(qL+qR−ξuLN−ηuRN)ω+λρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]
(30)
VEN(q(t))=λπE+δr4ρ+αq(t)+1ρ[−e2(uEN)2+πES0+r4ψuEN+λπE+δr4ρ+α(ξuLN+ηuRN+θuCN+ϕuEN)]

At this point, the optimal value of social welfare is:

(31)
S(t)N=S0+λ(ξuLN+ηuRN+θuCN+ϕuEN−αq(t))

Cost-Sharing Scenario

In this game model, the central government shares a certain proportion of the pollution control costs with the left-bank and right-bank regional governments. The regional governments, in turn, provide incentive subsidies to the enterprises within their respective jurisdictions. The amount of subsidies provided by the regional governments is linked to the pollution control efforts of the enterprises. The objective functions of the participants can be expressed as:

(32)
ρVGLS(q(t))=maxuL≥0{ξuLr1−1−m12uL2−en12uE2+πLS(t)−(qL−ξuL)ω+μ1ηuR+χϕuE+VGLS′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(33)
ρVGRS(q(t))=maxuR≥0{ηuRr2−k(1−m2)2uR2−en22uE2+πRS(t)−(qR−ηuR)ω+μ2ξuL+χϕuE+VGRS′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(34)
ρVGCS(q(t))=maxuC≥0{θuCr3−c2uC2−12m1uL2+(qL+qR−ξuL−ηuR)ω−k2m2uR2+S(t)+VGRS′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}
(35)
ρVES(q(t))=maxuE≥0{−(1−n1−n2)e2uE2+πES(t)+r4[ψuE+δq(t)]+VES′(q(t))[ξuL+ηuR+θuC+ϕuE−αq(t)]}

The equilibrium solutions of the game in the cost-sharing mechanism can be obtained as Equations (36).

(36)
uLS=ξ(r1+ω)−βμ1η+χϕd3−d3en1uES+VGLN′(q(t))(ξ−ηβ+ϕd3)1−m1*uRS=η(r2+ω)−βμ2ξ+χϕd3−d3en2uES+VGRS′(q(t))(η−ξβ+ϕd3)k(1−m2*)uCS=1c(θr3−ω(ξd1−d2β1−β2+ηd2−d1β1−β2)−m1*uLSd1−d2β1−β2−km2*uRSd2−d1β1−β2+VGCS′(q(t))(θ+ξd1−d2β1−β2+ηd2−d1β1−β2+ϕd3(d1−d2β1−β2+d2−d1β1−β2)))uES=r4ψ+VES′(q(t))ϕ(1−n1*−n2*)e

Similarly, the value function is a linear function of q(t). Thus, using the method of substitution of coefficients, the following equations are obtained:

(37)
VGLS′(q(t))=πLλρ+α,VGRS′(q(t))=πRλρ+α,VGCS′(q(t))=λρ+α,VES′(q(t))=λπE+δr4ρ+α

Substituting Equation (37) into Equation (36), the following result is obtained:

(38)
uLS=ξ(r1+ω)−βμ1η+χϕd3−d3en1uES+πLλρ+α(ξ−ηβ+ϕd3)1−m1*
(39)
uRS=η(r2+ω)−βμ2ξ+χϕd3−d3en2uES+πRλρ+α(η−ξβ+ϕd3)k(1−m2*)
(40)
uCS=1c(θr3−ω(ξd1−d2β1−β2+ηd2−d1β1−β2)−m1*uLSd1−d2β1−β2−km2*uRSd2−d1β1−β2+λρ+α(θ+ξd1−d2β1−β2+ηd2−d1β1−β2+ϕd3(d1−d2β1−β2+d2−d1β1−β2)))
(41)
uES=r4ψ+λπE+δr4ρ+αϕ(1−n1*−n2*)e

Then, substituting Equations (38)-(41) into the central government’s HJB and setting the first-order partial derivatives of the right-hand side with respect to mi, ni(i = 1, 2) to zero, we obtain Equations (42):

(42)
m1*=2(λρ+α−ω)ξ−ξ(r1+ω)+βμ1η−χϕd3+d3en1uES−πLλρ+α(ξ−ηβ+ϕd3)2(λρ+α−ω)ξ+ξ(r1+ω)−βμ1η+χϕd3−d3en1uES+πLλρ+α(ξ−ηβ+ϕd3)m2*=2(λρ+α−ω)η−η(r2+ω)+βμ2ξ−χϕd3+d3en2uRS−πRλρ+α(η−ξβ+ϕd3)2(λρ+α−ω)η+η(r2+ω)−βμ2ξ+χϕd3−d3en2uRS+πRλρ+α(η−ξβ+ϕd3)

By substituting Equation (38)-(41) into the objective functions of both the left-bank and right-bank local governments, and applying the first-order necessary conditions by setting the partial derivatives with respect to n1 and n2 to zero, the following results are obtained:

(43)
n1*=ϕd3e(πLλρ+α+χ)ξ(r1+ω+πLλρ+α)+c+πLλρ+αθm1cd1−d2β1−β2n2*=ϕd3e(πRλρ+α+χ)η(r2+ω+πRλρ+α)+c+πRλρ+αθkm2cd2−d1β1−β2

The optimal benefit functions for all participants are presented as shown in Equation (44) to (47):

(44)
VGLS(q(t))=πLλρ+αq(t)+1ρ[ξuLS(r1+ω)−1−m12uLS2−en12uES2+πLS0−qLω+μ1ηuRS+χϕuES+πLλρ+α(ξuLS+ηuRS+θuCS+ϕuES)]
(45)
VGRS(q(t))=πRλρ+αq(t)+1ρ[ηuRS(r2+ω)−k(1−m2)2uRS2−en22uES2+πRS0−qRω+μ2ξuLS+χϕuES+πRλρ+α(ξuLS+ηuRS+θuCS+ϕuES)]
(46)
VGCS(q(t))=λρ+αq(t)+1ρ[θuCSr3−c2uCS2−12m1uLS2−k2m2uRS2+S0+(qL+qR−ξuLS−ηuRS)ω+λρ+α(ξuLS+ηuRS+θuCS+ϕuES)]
(47)
VES(q(t))=λπE+δr4ρ+αq(t)+1ρ[−(1−n1−n2)e2uES2+πES0+r4ψuES+λπE+δr4ρ+α(ξuLS+ηuRS+θuCS+ϕuES)]

Under the cost-sharing mechanism, the maximum social welfare is shown in Equation (48).

(48)
S(t)S=S0+λ(ξuLS+ηuRS+θuCS+ϕuES−αq(t))

The following steps involve deriving and comparing the optimal control strategies of each participant under both scenarios.

According to Equations (23) and (38), we obtain:

(49)
uES−uEN=ψr4+(δr4+λπE)ϕα+ρe·(n1+n21−n1−n2)>0uLS−uLN=n1(δr4(η+ξ)+λπE(η+ξ)+d3r4(αψ+δϕ+λϕπE+ψρ))(α+ρ)(1−m1)(n1+n2−1)<0uRS−uRN=n2(δr4(η+ξ)+λπE(η+ξ)+d3r4(αψ+δϕ+λϕπE+ψρ))k(α+ρ)(1−m2)(n1+n2−1)<0uCS−uCN=−1c[m11−m1(d1−βd2)1−β2·X1+m21−m2(d2−βd1)1−β2·X2]<0X1=ξ(ω+r1)−βημ1+χd3ϕ+λπL(ξ−βη+d3ϕ)α+ρ+d3n1n1+n2−1(ψr4+(δr4+λπE)(η+ξ+d3ϕ)d3(α+ρ))X2=η(ω+r2)+χd3ϕ−βμ2ξ+λπR(η+d3ϕ−βξ)α+ρ+d3n2n1+n2−1·(ψr4+(δr4+λπE)(η+ξ+d3ϕ)d3(α+ρ))

Therefore, the cost-sharing mechanism not only increases the pollution control investment of industrial enterprises but also reduces the pollution control expenditure of local governments on both sides and the central government. Due to the complexity of the equilibrium solution, the impact of certain parameters cannot be directly determined. Therefore, it is necessary to conduct numerical analysis to identify the impact patterns of these parameters.

Numerical Simulation

In this section, the water pollution management model is simulated from the perspective of multiple stakeholders using MATLAB software, with exogenous variables assigned specific values. The parameter settings in this study are grounded in real-world scenarios and fundamental assumptions, with the values chosen based on these assumptions and informed by the relevant research, particularly the work of YANG et al. [41]. The initial assumptions for the model’s parameters are as follows: α = 0.01, β = 0.7, ω = 2, r1 = 0.06, r2 = 0.06, r3 = 0.06, r4 = 0.07, φ = 0.18, μ1 = 1, μ2 = 1.2, χ = 3, λ = 1.5, πL = 0.7, πR = 0.75, πE = 0.45, δ = 1.22, n1 = 0.2, n1 = 0.3, m1 = 0.2, m2 = 0.2, k = 2, c = 2, e = 2, ψ = 1.5, η = 0.15, θ = 0.13, ξ = 0.15.

As illustrated in Fig. 1 through Fig. 3, the implementation of the cost-sharing mechanism significantly improves wastewater treatment outcomes and enhances social welfare. Additionally, it leads to substantial increases in the profits of all participating stakeholders. This demonstrates that by distributing the financial burden among the involved parties, the mechanism not only boosts the overall efficiency of environmental management but also fosters the sustainability of the system.

Fig. 4 illustrates that as the competition coefficient increases, the local governments face heightened economic pressures, which may drive them to prioritize economic growth over environmental concerns. This results in a reduced focus on pollution control efforts, further exacerbating the overall decline in pollution management efficiency and profit (Fig. 5). Consequently, the negative impact of increased competition on environmental governance becomes more pronounced.

Fig. 3.

The revenue of each stakeholder under two different modes.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_3_min.jpg
Fig. 1.

Comparison of sewers treatment effects under two different modes.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_1_min.jpg
Fig. 2.

Comparison of social welfare under two different modes.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_2_min.jpg
Fig. 4.

The impact of the economic competition coefficient between the two banks on the governments’ optimal effort.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_4_min.jpg

Fig. 6 illustrates that subsidy policies can significantly increase the pollution control efforts of industrial enterprises and help alleviate the financial pressure on the governments of both sides regarding wastewater treatment.

Fig. 7 indicates that the pollution control subsidy policy for enterprises effectively enhances the overall welfare. However, as the subsidy ratio increases beyond a certain threshold, the revenues of the local governments on both banks begin to decline. This suggests that determining the optimal subsidy coefficient is crucial to balancing the benefits of pollution control with the financial sustainability of government resources.

As shown in Fig. 8, an increase in the cost-sharing ratio effectively incentivizes local governments to increase their investments in pollution control. Moreover, this mechanism reduces the financial burden on the central government, allowing it to allocate more resources to other areas.

Fig. 6.

The impact of the sum of the subsidy coefficients for industrial enterprises by the two banks’ governments on the government’s optimal effort.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_6_min.jpg

Conclusions and Suggestions

This study underscores that cost-sharing mechanisms between central and local governments significantly boost pollution control efforts, improve environmental quality, and enhance economic returns and social welfare. When the central government shares costs, local governments and enterprises face reduced financial burdens, promoting sustainability and more robust governance. By contrast, in the absence of cost-sharing, local governments under strong economic pressure often underinvest in pollution control as competition intensifies, ultimately undermining water quality and long-term benefits.

Fig. 5.

The impact of the economic competition coefficient between the two banks on the governments’ profit.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_5_min.jpg
Fig. 7.

The impact of the sum of the subsidy coefficients for industrial enterprises by the two banks’ governments on the government’s optimal profit.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_7_min.jpg
Fig. 8.

The impact of the central government’s cost-sharing coefficient (m) on the pollution control efforts of each participant by the governments on both banks.

https://www.pjoes.com/f/fulltexts/203358/PolJEnvironStud-2026-35-3-3151-3162_Fig_8_min.jpg

To address these issues, the central government should strengthen oversight to ensure policy consistency across transboundary water bodies. In regions under considerable economic strain, targeted financial and technical support can further reinforce pollution management. Increasing the central government’s cost-sharing ratio and linking funding to pollution control performance can incentivize more effective local strategies. Local governments, for their part, can implement differentiated subsidies based on the pollution control capacities of industrial enterprises, thereby encouraging greater investment and fostering cross-regional collaboration. By sharing technologies and experiences, local governments can balance resource allocation and generate collective benefits. Additionally, implementing environmental tax policies and intensifying the assessment of corporate environmental performance can motivate enterprises to adopt responsible practices and invest more in pollution mitigation.

To further enhance pollution control efforts, it is recommended to introduce transparent monitoring tools, such as real-time water quality monitoring systems and IoT technologies, to improve pollution source tracking and management, thereby enhancing transparency and the effectiveness of governance. Additionally, exploring cutting-edge technologies, such as reverse osmosis membrane technology and AI-driven optimization of wastewater treatment processes, can increase the efficiency of pollution control, reduce energy consumption, and promote innovation in industrial wastewater treatment.

A limitation of this study is its reliance on hypothetical parameters due to limited real-world data, which may affect its generalizability. Moreover, the model does not fully capture the potential nonlinear factors in environmental and economic systems, which is our next research focus.