<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mar. Sci.</journal-id>
<journal-title>Frontiers in Marine Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mar. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-7745</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmars.2025.1661382</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Marine Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Port green investment strategy under government subsidization and regulation policy</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Xie</surname>
<given-names>Fengjue</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3126085/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhou</surname>
<given-names>Zhongyin</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3141938/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Intelligent Manufacturing and Control Engineering, Shanghai Polytechnic University</institution>, <addr-line>Shanghai</addr-line>,&#xa0;<country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Gastroenterology, Renmin Hospital of Wuhan University</institution>, <addr-line>Wuhan</addr-line>,&#xa0;<country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Lang Xu, Nanyang Technological University, Singapore</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Le Zhang, Jinan University, China</p>
<p>Mingshuo Cao, Shanghai Maritime University, China</p>
<p>Xu Changyan, Shanghai Maritime University, China</p>
<p>Lingchen Huang, Fudan University, China</p>
<p>Min Zhou, Shanghai University of Electric Power, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Zhongyin Zhou, <email xlink:href="mailto:zjh001012@163.com">zjh001012@163.com</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>06</day>
<month>08</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>12</volume>
<elocation-id>1661382</elocation-id>
<history>
<date date-type="received">
<day>07</day>
<month>07</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>07</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Xie and Zhou.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Xie and Zhou</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The greener port has become an important goal for the maritime industry. This paper developed an economic model to discuss how the investment effectiveness and two commonly used governmental policies (regulation and subsidy policies) affect port&#x2019;s output, green investment, and governmental policy choices. Findings show that compared to no policy, both policies promote green investment, but their impacts on output are opposite: regulation restrains output, while subsidies boost it. Subsidies, as a proactive tool, reduce ports&#x2019; green investment costs, driving both investment and output&#x2014;especially in initial stages with high subsidy rates (0.3-0.5) covering over 50% of costs. Regulation, more passive, forces ports to increase investment but limit production, particularly when investment efficiency is low. Investment effectiveness shapes policy selection: subsidies work better at low or high efficiency, while regulation is preferable at medium efficiency. Both policies outperform no intervention in improving social welfare. These insights inform policy design for sustainable port development.</p>
</abstract>
<kwd-group>
<kwd>port green investment</kwd>
<kwd>regulation</kwd>
<kwd>subsidy</kwd>
<kwd>investment effectiveness</kwd>
<kwd>green policy</kwd>
</kwd-group>
<counts>
<fig-count count="8"/>
<table-count count="4"/>
<equation-count count="10"/>
<ref-count count="60"/>
<page-count count="14"/>
<word-count count="8013"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Marine Affairs and Policy</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>The urgency to accelerate decarbonization in the maritime sector corresponds with the International Maritime Organization&#x2019;s (IMO) updated 2023 strategy, which mandates strengthened emission reduction targets: minimum 20% by 2030 and 70-80% by 2040 relative to 2008 baseline levels (<xref ref-type="bibr" rid="B16">IMO, 2023</xref>). Current data reveals maritime transportation contributes 3.1% of global CO<sub>2</sub> emissions (<xref ref-type="bibr" rid="B46">UTCAD, 2023</xref>), with container vessels disproportionately responsible for 23% of sectoral emissions despite comprising merely 11% of the global fleet. Notably, China&#x2019;s pivotal role in container logistics is evidenced by its ports processing 38.7% of worldwide container throughput (over 310 million TEUs) in 2023 (<xref ref-type="bibr" rid="B9">China Ports &amp; Harbors Association, 2024</xref>). This challenge is further intensified by the sector&#x2019;s capital intensity, requiring over $1.5 trillion investment for energy transition through 2050 (<xref ref-type="bibr" rid="B13">Global Maritime Forum, 2023</xref>), particularly for modernizing port infrastructure.</p>
<p>While green technology adoption remains a conventional corporate approach for CO<sub>2</sub> mitigation, port-related environmental investments face dual financial constraints: elevated costs in both eco-technological innovation and low-carbon infrastructure development. These economic barriers substantially dampen operator motivation to engage in sustainable port initiatives. Therefore, developing effective strategies to incentivize green investments while implementing emissions control mechanisms constitutes a critical priority for transportation policymakers. Governments have developed multi-tiered policy interventions combining regulatory mandates with market incentives to address port-related emissions. For example, the Port of Los Angeles in the United States established a port emission inventory in the early 21st century and developed a real-time database for monitoring ship emissions. Local governments have implemented various regulatory measures such as emission taxes and incentives for green technologies. As a result, the Port of Los Angeles achieved a 98% reduction in sulfur oxide emissions while experiencing a 26% increase in container throughput from 2005 to 2018. Similarly, the Port of Rotterdam, located in Europe, has designated fixed areas around the port and developed an air quality assessment model to monitor emissions from ships and trucks, thereby achieving energy efficiency and emission reduction. In China, regulations and policies such as the &#x201c;Law on Prevention and Control of Air Pollution&#x201d; and the &#x201c;Yangtze River Protection Law&#x201d; have been formulated to address emission issues in the maritime transportation system. Local governments are responsible for enforcing these regulations and are encouraged to implement effective measures, including penalties, to reduce ship emissions. Furthermore, the Ministry of Transport has introduced the &#x201c;Management Measures for Shore Power in Ports and Ships,&#x201d; which stipulates that ships berthed in ports should use shore power. Ports and shipping companies that refuse to comply may face administrative fines imposed by local governments. Ports and shipping companies are also expected to supervise each other to promote energy-saving and emission reduction in the maritime transportation system. The implementation of emission taxes and penalties mentioned above are regulatory measures aimed at controlling carbon emissions. By assigning corresponding environmental costs to emissions, ports are compelled to make efforts to reduce their emissions.</p>
<p>Beyond establishing carbon emission thresholds, governments globally have deployed strategic subsidy programs to accelerate renewable energy adoption. Recent analyses indicate that 78% of OECD nations now integrate renewable subsidies with carbon pricing mechanisms (<xref ref-type="bibr" rid="B30">Organisation for Economic Co-operation and Development (OECD), 2023</xref>). There are a number of companies that have gotten significant help from the government in the form of subsidies (<xref ref-type="bibr" rid="B38">Sendstad et&#xa0;al., 2022</xref>). Prominent corporations exemplify this trend: Tesla secured over $15 billion in cumulative U.S. federal and state incentives through 2023, including revamped EV tax credits under the Inflation Reduction Act, while receiving 8.2 billion CNY in Chinese manufacturing subsidies between 2019-2022. The Dutch government has approved approximately $2.4 billion in subsidies for the world&#x2019;s largest Carbon Capture and Storage (CCS) project, which will be located in the port of Rotterdam. It is anticipated that this project will reduce carbon emissions from the industrial cluster located around Europe&#x2019;s largest seaport by approximately 10% (<xref ref-type="bibr" rid="B17">Janipour et&#xa0;al., 2021</xref>).</p>
<p>Based on the aforementioned government policy cases, we categorize them into subsidy policies and regulatory policies. Most of the existing literature focuses on the impact of a single policy on the greening of ports. However, our study comprehensively considers the important role played by both subsidy and regulatory policies in the port industry and conducts a holistic analysis and comparison. This enriches the literature on green investments in ports. The following are the main research questions: (i) Should the government introduce policies to promote the development of green ports? If so, which is better for social welfare: subsidies or regulations? (ii) How do these policies affect port production and green investment planning compared to the benchmark? (iii) How does the optimal policy change with the improvement of green investment efficiency?</p>
<p>In summary, our research aims to examine the impact of regulatory policies and subsidy policies on port production and green investment planning, and determine the optimal policy considering social welfare. To address these issues specifically, we develop a systematic model to explore investment and output decisions by port operators without government intervention, which serves as a benchmark for comparing the results of government regulations or subsidy policies. The main contributions of this study can be summarized as follows: First, it deals with both port production and green investment planning in an integrated model, as an increase in port throughput will simultaneously increase economic contribution and greenhouse gas emissions. In addition to investment costs, the environmental cost of port emissions is proposed as a standard for government regulation. The daily emissions of the port are directly proportional to its output, and the port needs to bear the corresponding environmental cost for its emissions. When it undertakes green investments, it will reduce the expenditure on this portion of the cost. Furthermore, based on the initial emission level of the port, the environmental cost is divided into high-emission and low-emission parts. The model also considers the effectiveness of green investments. The impact of green investment effectiveness on port decision-making is investigated, and the selection of government policies under different levels of effectiveness is determined.</p>
<p>The rest of this paper is organized as follows. Section 2 reviews the relevant literature, and Section 3 establishes the model and elaborate the game stage. The benchmark case without government policy, and the outcomes under subsidization and regulation are derived and compared in Section 4. In Section 5, numerical simulation is conducted to examine the impacts of several key parameters on port profits. Finally, Section 6 concludes the paper and suggests several future research directions. For clarify, all detailed proofs of the analytical results in this paper are provided in the <xref ref-type="app" rid="app1">
<bold>Appendix 1</bold>
</xref>.</p>
</sec>
<sec id="s2">
<label>2</label>
<title>Literature review</title>
<p>Our work is related to three main research streams of literature: port investment strategy, the greening of port management and government&#x2019;s policy on the green development of ocean transportation. For recent studies about the port investment strategy, some discuss the investment on port capacity, see (<xref ref-type="bibr" rid="B53">Xu et&#xa0;al., 2025a</xref>; <xref ref-type="bibr" rid="B26">Lu et&#xa0;al., 2025</xref>; <xref ref-type="bibr" rid="B4">Balliauw et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B14">Gong et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B52">Xie et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B56">Xu et&#xa0;al., 2021</xref>). Some investigate the port investment in disaster prevention or adaptation (<xref ref-type="bibr" rid="B51">Xia et&#xa0;al., 2024</xref>; <xref ref-type="bibr" rid="B50">Xia and Lindsey, 2021</xref>; <xref ref-type="bibr" rid="B59">Zheng et&#xa0;al., 2021</xref>). Others focus on the green investment in port (<xref ref-type="bibr" rid="B54">Xu et&#xa0;al., 2025b</xref>; <xref ref-type="bibr" rid="B36">Schodler and Saraceni, 2024</xref>; <xref ref-type="bibr" rid="B27">Mahmud et&#xa0;al., 2024</xref>; <xref ref-type="bibr" rid="B11">Finlayson et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B40">Song et&#xa0;al., 2022</xref>);, which is can also be included in the second stream of literature.</p>
<p>The second stream of the literature mainly concentrates on the greening of port management. Port authorities around the world are greening their port management to safeguard their operating licenses and enhance their economic and environmental competitiveness (<xref ref-type="bibr" rid="B29">Notteboom and Lam, 2018</xref>). In past two decades, many scholars have shed light on the green strategies of port activities, which includes reducing energy use, using cleaner electric (electric or LNG) vehicles for port activities, and developing renewable energy sources at the port (<xref ref-type="bibr" rid="B55">Xu et&#xa0;al., 2025c</xref>; <xref ref-type="bibr" rid="B2">Acciaro et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B7">Chen et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B49">Winnes et&#xa0;al., 2015</xref>). The plans for waste management and recycling, including sewage, spills, and ballast water, reducing noise, vibration, and dust from handling activities (<xref ref-type="bibr" rid="B21">Lam and Notteboom, 2014</xref>) and ecosystem protection relating to the marine habitat, wetland, and coastal erosion were also taken into consideration (<xref ref-type="bibr" rid="B8">Chen and Pak, 2017</xref>). Further, Port and shipping companies are also focusing on the application of investment in emission reduction technologies. This includes the use of low-carbon equipment and investment in updated green technologies (<xref ref-type="bibr" rid="B1">Abbas, 2019</xref>; <xref ref-type="bibr" rid="B28">Metzger and SChinas, 2019</xref>; <xref ref-type="bibr" rid="B35">Ren et&#xa0;al., 2019</xref>), considering emissions from a port equipment perspective (<xref ref-type="bibr" rid="B12">Geerlings and van Duin, 2011</xref>; <xref ref-type="bibr" rid="B57">Yang and Lin, 2013</xref>). Some ports invest in shore power systems for ships (<xref ref-type="bibr" rid="B22">Li et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B33">Radwan et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B42">Tan et&#xa0;al., 2023</xref>), while others focus on the use of hydrogen energy in port shipping (<xref ref-type="bibr" rid="B18">Kizielewicz, 2022</xref>). <xref ref-type="bibr" rid="B5">Ballini and Bozzo (2015)</xref> proposed that shore power technology can reduce shipping carbon emissions and improve the port city environment <xref ref-type="bibr" rid="B37">Sciberras et&#xa0;al. (2015)</xref> found that using shore power technology can significantly reduce emissions when ships are docked. (<xref ref-type="bibr" rid="B19">Kotrikla et&#xa0;al., 2017</xref>). found that shore power infrastructure can promote the sustainable development of ports and mitigate negative impacts on the urban environment. Expert scholars have found that investing in port infrastructure to increase ship berthing capacity has positive impacts on reducing transportation costs and CO2 emissions. <xref ref-type="bibr" rid="B24">Liu et&#xa0;al. (2021)</xref> focused on port-led supply chain and analyzed the interest relationships between different members in three investment scenarios. They then analyzed the changes in social welfare under different investment strategies. <xref ref-type="bibr" rid="B39">Sogut and Erdo&#x11f;an (2022)</xref> developed a conceptual framework related to the green concept to support environmental sustainability. <xref ref-type="bibr" rid="B20">Lai et&#xa0;al. (2019)</xref> proposed that port and shipping companies face complex decision-making and multiple stakeholder involvement in sustainable development, influenced by economic, social, and environmental responsibilities. Technological, organizational, market, policy support, and socio-economic aspects also pose challenges. To comply with strict emission control regulations, port and shipping companies should invest in the sustainability of the environment, improve the sustainability of the port supply chain, and enhance the sustainability of the investment environment in the face of market uncertainty (<xref ref-type="bibr" rid="B34">Randrianarisoa and Gillen, 2022</xref>) In this paper, we have taken a unique approach in exploring the greening of port management by focusing on the strategy of green investment and output, rather than specific implementation methods discussed in previous literature. Our goal is to effectively reduce the environmental cost in port daily operation. This fresh perspective offers a new dimension to the existing body of knowledge on port management and sustainability.</p>
<p>Studies on the government&#x2019;s policies on green development of ocean transportation can be further divided into two categories: the emission regulations in ocean shipping (<xref ref-type="bibr" rid="B23">Liu et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B43">Tang et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B45">Tseng and Pilcher, 2019</xref>; <xref ref-type="bibr" rid="B47">Wang et&#xa0;al., 2007</xref>), the policies and regulations of green ports (<xref ref-type="bibr" rid="B3">Ayotamuno and Gobo, 2004</xref>; <xref ref-type="bibr" rid="B6">Birch and Taylor, 2002</xref>; <xref ref-type="bibr" rid="B41">Stojanovic et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B58">Zheng et&#xa0;al., 2017</xref>). Here, we focus on the latter. In the context of policies and regulations on green ports development, <xref ref-type="bibr" rid="B10">Farrell (1986)</xref> has proposed the concept of the subsidization of seaports at an early stage. In recent studies, some scholars have paid attention to the impact of policy instruments on port emission reduction, but most focus on the effect analysis of a single policy. For example, <xref ref-type="bibr" rid="B22">Li et&#xa0;al. (2020)</xref> Study the two-echelon maritime supply chain composed of ports and shipping companies under the government&#x2019;s green subsidy and explore the subsidy mechanism and its impact. Considering the influence of government subsidies, a multi-period dual-objective port shore power deployment optimization model is constructed by <xref ref-type="bibr" rid="B25">Lu and Huang (2021)</xref>, which is aiming at minimizing operating costs and minimizing CO2 emissions. <xref ref-type="bibr" rid="B48">Wang and Jiao (2022)</xref> developed a game model to study the shipping companies&#x2019; choice of low sulfur fuel oil with government subsidy, and the possible impact has been explored. Except for the subsidization, the regulations on port emission have been discussed widely. <xref ref-type="bibr" rid="B60">Zhu and Wang (2021)</xref> evaluated wind pattern and pollutant retention time, the actual effect and spatial distribution of fuel content control in four ports based on the control implementation. In existing studies, some scholars have paid attention to the impact of policy instruments on port emission reduction, but most focus on the effect analysis of a single policy. <xref ref-type="bibr" rid="B31">Park et&#xa0;al. (2018)</xref> constructed an analytical model of government emission standards, compared the emission reduction responses of port operators under multilateral, unilateral and single-country regulatory scenarios, and found that multilateral regulation is more likely to achieve global emission reduction goals, but did not involve a comparison with subsidy policies. <xref ref-type="bibr" rid="B60">Zhu and Wang (2021)</xref> evaluated the actual effects of fuel content regulations in four major ports and found that policy stringency and spatial distribution have a significant impact on emission reduction effectiveness, but did not explore the synergistic effect of economic incentive instruments such as subsidies.</p>
<p>To sum up, the existing literature has adequately analyzed the mechanism of a single policy, but there is a lack of a systematic comparison of the &#x201c;incentive-constraint&#x201d; effects of subsidy and regulatory policies in port green investment, especially the failure to combine dynamic evaluation of investment effectiveness to determine the optimal policy choice. This study addresses this gap by integrating the two policy instruments and constructing a dynamic analysis framework to fill the gap in existing research. By incorporating both subsidization and regulation into our analysis, we can gain a comprehensive understanding of their effects on the green development of seaports. We compare the outcomes under these two different government policies with a scenario in which no government intervention is present. This allows us to identify the optimal policy in light of varying levels of investment effectiveness. Furthermore, we derive insights into the investment and output decisions made by seaports under different policy scenarios.</p>
</sec>
<sec id="s3">
<label>3</label>
<title>Modeling preliminaries</title>
<p>In this part, we establish an analytic model to evaluate port operators&#x2019; decisions to adopt green investments in response to two government policies, namely subsidization and regulation. The decision-making without government interference is used as a benchmark to evaluate the impact of government regulation and subsidies.</p>
<p>In multi-port areas, most ports have cross hinterland, and these ports provide homogeneous services. Therefore, we consider that there are <inline-formula>
<mml:math display="inline" id="im1">
<mml:mi>n</mml:mi>
</mml:math>
</inline-formula> identical ports in the same region, whose demand function has the same form as <xref ref-type="disp-formula" rid="eq1">Equation 1</xref>:</p>
<disp-formula id="eq1">
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the throughputs of Port <inline-formula>
<mml:math display="inline" id="im3">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>; <inline-formula>
<mml:math display="inline" id="im4">
<mml:mi>P</mml:mi>
</mml:math>
</inline-formula> is the generalized price for each port, where <inline-formula>
<mml:math display="inline" id="im5">
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> measures the substitutability of port services, which also reflects the degree of competition between these ports.</p>
<p>As one of the most significant sources of carbon emissions, the daily operation of maritime ports must bear the associated environmental cost. The environmental cost of port emissions is the damage to the environment caused by the emissions in the daily production operations of ports. Using shore power when ships are docked is one of the most effective ways to reduce air pollutant emissions. However, the installation and utilization of onshore power varies from terminal to terminal (<xref ref-type="bibr" rid="B32">Peng et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B33">Radwan et&#xa0;al., 2019</xref>). In addition, the daily operation efficiency of each port is different, and the green concept of the cooperative terminal operators is different, which will lead to the classification of ports into low-emission zones that have implemented certain emission reduction technologies and high-emission zones that maintain the original emission levels. Therefore, the environmental cost of emissions to be borne by port <inline-formula>
<mml:math display="inline" id="im6">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> is defined as <xref ref-type="disp-formula" rid="eq2">Equation 2</xref>:</p>
<disp-formula id="eq2">
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where the proportion of high-emission parts is <inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and the environmental costs that ports have to undertake for each unit of output is <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>; meanwhile, the proportion of low-emission part is <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, and the environmental costs that ports have to undertake for each unit of output is <inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B44">Tian et&#xa0;al., 2019</xref>); <inline-formula>
<mml:math display="inline" id="im11">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> represents the green investment effectiveness. Note that a larger <inline-formula>
<mml:math display="inline" id="im12">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula> implies that carbon emissions cost reduction efficiency increases. It worth noting that the effectiveness of green technology adoption in low-emission part will be reduced to <inline-formula>
<mml:math display="inline" id="im13">
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&lt;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. Further, we assume that <inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im16">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, which means green investment in ports can only reduce the environmental costs, but it won&#x2019;t have any other spill-over effects. Thus, the expected profit of port <inline-formula>
<mml:math display="inline" id="im17">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> has the following form:</p>
<disp-formula id="eq3">
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>S</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>
<inline-formula>
<mml:math display="inline" id="im18">
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the green investment of each port <inline-formula>
<mml:math display="inline" id="im19">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> having a convex form (<xref ref-type="bibr" rid="B15">Hong and Guo, 2019</xref>). Further, the social welfare is considered as the sum of port users&#x2019; surplus and the total ports&#x2019; profits, which is:</p>
<disp-formula id="eq4">
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>a</mml:mi>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>b</mml:mi>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mi>S</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Specifically, when the policy makers choose to regulate port carbon emissions, they will set the minimum requirement, the government regulates the expected environmental cost must be less than a target <inline-formula>
<mml:math display="inline" id="im20">
<mml:mi>R</mml:mi>
</mml:math>
</inline-formula>. We employ superscript <inline-formula>
<mml:math display="inline" id="im21">
<mml:mi>R</mml:mi>
</mml:math>
</inline-formula> to denote the regulation model. The constraint can be specified as <xref ref-type="disp-formula" rid="eq5">Equation 5</xref>:</p>
<disp-formula id="eq5">
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>.</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>.</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x2264;</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>When subsidization is adopted, port will receive subsidy <inline-formula>
<mml:math display="inline" id="im22">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> for every unit of the green investment (denoted by superscript <inline-formula>
<mml:math display="inline" id="im23">
<mml:mi>S</mml:mi>
</mml:math>
</inline-formula>). If the government choose the to subsidize the port&#x2019;s green investment, the cost will be reduced from <inline-formula>
<mml:math display="inline" id="im24">
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> to <inline-formula>
<mml:math display="inline" id="im25">
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>s</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>We model the policy decision and port&#x2019;s action with a multistage game and depict in the <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Multistage game.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g001.tif">
<alt-text content-type="machine-generated">Flowchart illustrating a three-stage process. Stage one depicts a city icon with arrows leading to two options: &#x201c;Subsidization S&#x201d; and &#x201c;Regulation R&#x201d;. Both options direct arrows to multiple shipping port icons labeled &#x201c;Port 1&#x201d;, &#x201c;Port 2&#x201d;, up to &#x201c;Port i&#x201d;. This represents decision-making processes affecting different ports. Stage labels are at the bottom.</alt-text>
</graphic>
</fig>
<p>The government will choose between the two policies in the first stage. Once the subsidy <inline-formula>
<mml:math display="inline" id="im26">
<mml:mrow>
<mml:msup>
<mml:mi>s</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> or the regulation <inline-formula>
<mml:math display="inline" id="im27">
<mml:mrow>
<mml:msup>
<mml:mi>R</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is determined, ports will decide their green investments in stage 2 under either the subsidization or the regulation policy and decide the optimal outputs in stage 3.</p>
</sec>
<sec id="s4">
<label>4</label>
<title>Market equilibriums</title>
<p>To clarify all the notation used in the paper, we describe the terms in advance in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Terms description.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Terms</th>
<th valign="top" align="center">Description</th>
<th valign="bottom" align="center">Unit</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im28">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The throughput of port <inline-formula>
<mml:math display="inline" id="im29">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>.</td>
<td valign="bottom" align="center">TEU</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im30">
<mml:mi>P</mml:mi>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The service charge of port <inline-formula>
<mml:math display="inline" id="im31">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>.</td>
<td valign="bottom" align="center">USD/TEU</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im32">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The environmental cost of port <inline-formula>
<mml:math display="inline" id="im33">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>.</td>
<td valign="bottom" align="center">USD</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im34">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The unit environmental cost of port production in high and low emission areas</td>
<td valign="bottom" align="center">USD</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im35">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">Port i&#x2019;s profits.</td>
<td valign="bottom" align="center">USD</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im36">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">Social welfare.</td>
<td valign="bottom" align="center">USD</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im37">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The effectiveness of port green investment.</td>
<td valign="bottom" align="center">/</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im38">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The level of green investment decided by port <inline-formula>
<mml:math display="inline" id="im39">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>.</td>
<td valign="bottom" align="center">/</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im40">
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The investment cost coefficient.</td>
<td valign="bottom" align="center">/</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im41">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">The proportion of areas with high emissions.</td>
<td valign="bottom" align="center">/</td>
</tr>
</tbody>
</table>
</table-wrap>
<sec id="s4_1">
<label>4.1</label>
<title>Benchmark</title>
<p>The equilibrium outcomes of port&#x2019;s output and green investment in benchmark model are derived and summarized in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>. (The detailed process is presented in the Appendix 1 by using backward induction).</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Optimal port production and green investment decision without government intervention.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Port <italic>i</italic>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im42">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im43">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im44">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im45">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="2" align="left">Benchmark</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im46">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im47">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im48">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im49">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im50">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im51">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im52">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im53">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>After observing the equilibrium outcomes, we have the following lemma.</p>
<p>
<bold>Lemma 1.</bold> <italic>Optimal green investment and output decisions of each port depend on the investment effectiveness</italic> <inline-formula>
<mml:math display="inline" id="im54">
<mml:mrow>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>Lemma 1 demonstrates that port investment and output decisions have a step-like structure in response to a change in investment effectiveness. To better comprehend lemma 1, we assumed the values of a few factors, and <xref ref-type="fig" rid="f2">
<bold>Figures&#xa0;2a, b</bold>
</xref> depicts the ideal port&#x2019;s investment and production. When <inline-formula>
<mml:math display="inline" id="im55">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,we observe that the port&#x2019;s green investment will increase with the investment&#x2019;s effectiveness, but the port&#x2019;s production will not increase correspondingly. When <inline-formula>
<mml:math display="inline" id="im56">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,we find <inline-formula>
<mml:math display="inline" id="im57">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&lt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im58">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, which demonstrates that increasing investment effectiveness will reduce investment while increasing output. This is due to the fact that when the emission reductions from green investment can only cover the environmental cost of the low-emission portion, the cost of the investment exceeds its contribution. When <inline-formula>
<mml:math display="inline" id="im59">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,with the growth of green contribution, the environmental cost of the high-emission portion will be partially covered, the port operator will be more likely to invest in green technology. However, after all the environmental cost is fully covered (<inline-formula>
<mml:math display="inline" id="im60">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>, the port will not pursue additional green investment because it is meaningless to the contribution of green investment.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>The port&#x2019;s optimal output and green investment in benchmark model.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g002.tif">
<alt-text content-type="machine-generated">Graph (a) shows a blue line at a constant level of approximately 1.5, marked by vertical dashed lines at \( \zeta_1 \), \( \zeta_2 \), and \( \zeta_3 \). Graph (b) features a blue line that rises sharply from zero to about 1.5, peaks slightly after \( \zeta_1 \), and then declines beyond \( \zeta_3 \). Both graphs have axes labeled \( q_i^B \) or \( l_i^B \) and \( r \).</alt-text>
</graphic>
</fig>
<p>In addition, we find that the port&#x2019;s investment and output will drop with the growth in the number of rivals in the same region and will increase with the market size, which can be obtained from <inline-formula>
<mml:math display="inline" id="im61">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&lt;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&lt;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>&#xa0;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2217;</mml:mo>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>a</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&gt;</mml:mo>
<mml:mn>0.</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.&#xa0;A less concentrated market, on the other hand, may encourage a greater overall investment and output, both of which fits in with the Cournot model. According to the Cournot model, as the number of competitive enterprises increases, output of each individual participant will decrease. However, the output of the whole market will be improved simultaneously. Additionally, the size of the market drives up production and green investment at each port.</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>Regulation and subsidization</title>
<p>When the government chooses a policy between regulation and subsidization, the policy becomes the premise upon which ports determine the optimal green investment and output. The equilibrium outcomes of the regulation and subsidy are presented in <xref ref-type="table" rid="T2">
<bold>Tables&#xa0;2</bold>
</xref>, <xref ref-type="table" rid="T3">
<bold>3</bold>
</xref>, and the calculation procedure is detailed in Appendix 1.</p>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Optimal port production and green investment decisions under government&#x2019;s regulation.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Port <italic>i</italic>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im62">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im63">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left" rowspan="3">Regulation</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im64">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>a</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im65">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im66">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im67">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im68">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im69">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<sec id="s4_2_1">
<label>4.2.1</label>
<title>Regulation policy</title>
<p>
<bold>Lemma 2.</bold> <italic>Optimal green investment and output decisions under government&#x2019;s regulation of each port depend on the investment effectiveness</italic> <inline-formula>
<mml:math display="inline" id="im70">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula> <italic>and the minimum regulation</italic> <inline-formula>
<mml:math display="inline" id="im71">
<mml:mi>R</mml:mi>
</mml:math>
</inline-formula>.</p>
<p>Lemma 2 and <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref> demonstrate that, based on the regulation model, ports&#x2019; decisions about their green investments and outputs have been split into two portions (when <inline-formula>
<mml:math display="inline" id="im72">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) due to an increase in the effectiveness of the investments. We discover that ports always have the intention of increasing their green investments while simultaneously decreasing their outputs, which can be deduced from <inline-formula>
<mml:math display="inline" id="im73">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im74">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&lt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.This is due to the fact that as the effectiveness of investments increases, green investments will provide a greater contribution to the reduction of the environmental cost, while ports will be required to bear a lower total cost of investments. Moreover, the government has a carbon emission regulation that is less strict in the first part of the region (<inline-formula>
<mml:math display="inline" id="im75">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) and a regulation that is more severe in the second portion of the region (<inline-formula>
<mml:math display="inline" id="im76">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>). Because of the more restrictive standards imposed by the government, ports are being required to cover the environmental cost of both low and high emission parts to the greatest extent possible, which has caused a jump at <inline-formula>
<mml:math display="inline" id="im77">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>To compare the benchmark model with the regulatory model, we have made some simple numerical assumptions and depicted the results in <xref ref-type="fig" rid="f3">
<bold>Figures&#xa0;3a, b</bold>
</xref>. It is obvious that after the government makes the decision to set the maximum expected environmental cost, ports will be willing to produce more while investing less(<inline-formula>
<mml:math display="inline" id="im78">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>&lt;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im79">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>&gt;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>).This may be because increasing investment will hurt ports profits; increasing production and reducing investment will also hurt ports profits if they are operating within regulatory boundaries; however, increasing green investment has a greater potential for loss, which is why ports would rather increase production than expand green investment. This kind of situation will be improved as investment effectiveness increases.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>The port&#x2019;s optimal output and green investment in regulation model.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g003.tif">
<alt-text content-type="machine-generated">Two graphs labeled (a) and (b) display different functions of \( q_i \), \( l_i \), and \( r \). In graph (a), lines \( q_i^R \) and \( q_i^B \) are graphed with red and blue lines, respectively. In graph (b), \( l_i^B \) and \( l_i^R \) are plotted, showing different slopes. Both graphs include vertical dashed lines at points \( \zeta_4 \), \( \zeta_1 \), \( \zeta_2 \), and \( \zeta_3 \), indicating key intervals.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s4_2_2">
<label>4.2.2</label>
<title>Subsidization policy</title>
<p>After observing the equilibrium outcomes, we have the following lemma.</p>
<p>
<bold>Lemma 3.</bold> <italic>Optimal green investment and output decisions under government&#x2019;s subsidization of each port depend on the investment effectiveness</italic> <inline-formula>
<mml:math display="inline" id="im80">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula> <italic>and the subsidy</italic> <inline-formula>
<mml:math display="inline" id="im81">
<mml:mi>s</mml:mi>
</mml:math>
</inline-formula>.</p>
<p>Lemma 3 and <xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref> demonstrate that, according to the subsidization model, the port&#x2019;s investments and outputs similarly exhibit a step structure with rising effectiveness. In addition, it states that the subsidy is always necessary for ports and that it aims to reach as much as feasible in order to fully cover the green investments made by ports because <inline-formula>
<mml:math display="inline" id="im82">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. It is found that <inline-formula>
<mml:math display="inline" id="im83">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>&#x2264;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im84">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>&lt;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, as demonstrated by the outcomes depicted in <xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4a, b</bold>
</xref>.</p>
<table-wrap id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Optimal port production and green investment decision under government&#x2019;s subsidization.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Port <italic>i</italic>
</th>
<th valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im85">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
<th valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im86">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="3" align="left">Subsidization</td>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im87">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im88">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im89">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im90">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im91">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im92">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im93">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im94">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<inline-formula>
<mml:math display="inline" id="im95">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im96">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
<tr>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im97">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im98">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im99">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im100">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td valign="top" align="left">
<inline-formula>
<mml:math display="inline" id="im101">
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im102">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
<inline-formula>
<mml:math display="inline" id="im103">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>S</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>The port&#x2019;s optimal output and green investment in subsidization model.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g004.tif">
<alt-text content-type="machine-generated">Two graphs labeled (a) and (b) compare functions of a variable r. Graph (a) shows two lines labeled \(q_i^S\) and \(q_i^B\) intersecting vertically at points \(\zeta_1\) to \(\zeta_7\). Graph (b) features lines labeled \(l_i^S\) and \(l_i^B\) forming a peak between \(\zeta_5\) and \(\zeta_3\). Orange and black dashed vertical lines highlight specific \(\zeta\) points on both graphs.</alt-text>
</graphic>
</fig>
<p>Comparing the subsidization model to the benchmark reveals that, unlike the regulatory model, government subsidies simultaneously stimulate green port investment and outputs. Therefore, the subsidization policy is a more active strategy to controlling carbon emissions by limiting investment, whereas the regulating policy is a more passive approach for ports.</p>
</sec>
<sec id="s4_2_3">
<label>4.2.3</label>
<title>The optimal policy</title>
<p>In order to find out the government&#x2019;s optimal options in the three models, namely benchmark, regulation and subsidization, we compare the expected social welfare and have the results depicted as follows in <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>The optimal government policy.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g005.tif">
<alt-text content-type="machine-generated">A chart representing policy options along a horizontal axis labeled with increasing numbers from zero to four. Vertical dashed lines at positions one, two, three, and four represent thresholds labeled with symbols &#x3b6;4 to &#x3b6;7. The chart alternates between &#x201c;subsidy&#x201d; and &#x201c;subsidy or regulation&#x201d; at different intervals between these thresholds with SW noted at the top. The lines are color-coded blue, black, and red.</alt-text>
</graphic>
</fig>
<p>Without government policy, port green investments will contribute less to social welfare than they would with policy assistance. The policies stimulate green investment and output to some extent. Even if the government provides subsidies, ports will not increase their outputs to cover the additional environmental costs while the effectiveness of investment is at a relatively low level, that is, <inline-formula>
<mml:math display="inline" id="im104">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. During this time, ports will invest in green technology; but, they will not increase their outputs to bear a greater share of the environmental costs. From <xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>, it can be found that of the two kinds of government policy, subsidy policy still holds the dominant position.</p>
<p>When <inline-formula>
<mml:math display="inline" id="im105">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,the total environmental benefit rises, and the subsidy program not only supports port investment but also boosts port output. And in <inline-formula>
<mml:math display="inline" id="im106">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,when the investment is efficient enough, the port has a strong motivation to fully cover the environmental cost from both high emission and low emission parts, thereby making the minimum regulation unnecessary. Consequently, the subsidies contribute to the enhancement of social welfare.</p>
<p>However, when <inline-formula>
<mml:math display="inline" id="im107">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mtext>U</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>,due to the existence of the minimum regulation on the environmental cost, the port will raise its investment in green technology, but the green benefit will still be less than its cost. In this way, large-scale output will harm the port&#x2019;s profitability, thus the port&#x2019;s output will be lowered. In this case, minimum regulations may be preferable than subsidies. At this point, for port green development, although regulation leads to a decline in port output, the cost savings brought about by increased green investment partially offset the profit loss. Especially when investment effectiveness is at a medium level, the marginal emission reduction cost is lower than the fiscal expenditure on subsidies. Under regulation, the reduction in output lowers consumer surplus, but the improvement in environmental benefits makes up for this loss, resulting in an overall increase in social welfare. In contrast, although subsidy policies boost consumer surplus, high subsidy costs may crowd out public welfare expenditures, leading to a net loss of social welfare in the interval of medium investment effectiveness. In addition, it is indicated that if there are more rivals or a larger market, the players will be less interested in growing their investment; hence, regulation will be more beneficial in such a circumstance.</p>
</sec>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Numerical studies</title>
<p>In Section 4, we explore the impact of investment effectiveness <inline-formula>
<mml:math display="inline" id="im108">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula> on port optimal green investment and output decisions under government&#x2019;s policies. However, it&#x2019;s still missing the exploration of the impact of several key parameters on the ports&#x2019; profits. Therefore, in this section, we conduct the numerical studies to examine the impacts of effectiveness parameter <inline-formula>
<mml:math display="inline" id="im109">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>, the proportion of high-emission part in the port <inline-formula>
<mml:math display="inline" id="im110">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> and the substitutability of port services <inline-formula>
<mml:math display="inline" id="im111">
<mml:mi>b</mml:mi>
</mml:math>
</inline-formula>.We set the following parameter values as: <inline-formula>
<mml:math display="inline" id="im112">
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>10</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>n</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mtext>&#xa0;</mml:mtext>
<mml:mi>k</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>,which can be verified that these values satisfy all the positive assumptions in our model. First, we discuss the ports&#x2019; profits with different levels of <inline-formula>
<mml:math display="inline" id="im113">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>. In this case, we let <inline-formula>
<mml:math display="inline" id="im114">
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0.5</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and depict the results in <xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>The effect of <inline-formula>
<mml:math display="inline" id="im115">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula> on ports&#x2019; profits under three models.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g006.tif">
<alt-text content-type="machine-generated">A graph with three lines representing different datasets. The horizontal axis is labeled &#x201c;r&#x201d; and the vertical axis is labeled &#x201c;&#x3c0;&#x201d;. The brown line (S) increases in steps, the blue line (B) rises gradually, and the red line (R) starts lower and curves upward, intersecting with the blue line around r equals three.</alt-text>
</graphic>
</fig>
<p>
<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref> illustrates that in all three options, the revenues generated by the ports will go up by <inline-formula>
<mml:math display="inline" id="im116">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>. It can be obtained that, from the perspective of port operators, once the government decides to subsidize their green investment, they can achieve the highest profits. However, it can be found that the regulation policy would hurt ports&#x2019; profits to some extent. Only when the investment effectiveness is at a high level, the regulation of government won&#x2019;t be annoying for port operators.</p>
<p>Taking the application of shore power at Shanghai Port as an example, the Shanghai Municipal Transportation Commission began to promote the construction of shore power facilities at the port in 2014 and formulated relevant policies and standards. In 2015, the first batch of shore power facilities at Shanghai Port were put into use, and some ships started to use shore power for power supply. In 2017, Shanghai explicitly required ships of specific types and sizes to install shore power interfaces and use shore power for power supply. In the initial stage of shore power implementation, the utilization rate of shore power was relatively low due to factors such as the need for shore power infrastructure construction and ship adaptability. The construction and maintenance of shore power facilities usually require a certain amount of investment, which also has an impact on the financial status of the port. It was not until 2019 that the Shanghai Municipal Government introduced a subsidy policy, providing subsidies and incentive measures for ships that meet specific conditions to use shore power. The implementation of this subsidy policy further encouraged ships to use shore power, thereby increasing the utilization rate of shore power. After the implementation of the subsidy policy in 2019, the utilization rate of shore power increased from 15% to 40%. According to the 2020 Environmental Report of the Shanghai Port Authority, the proportion of high-emission areas &#x3b8; decreased from 0.6 to 0.3. Combined with model calculations, the effectiveness of green investment r at this time is approximately 1.2, which is between the theoretical thresholds &#x3b6;<sub>4</sub> (about 1.0) and &#x3b6;<sub>5</sub> (about 1.5). This is consistent with the model&#x2019;s prediction that &#x201c;regulatory policies are more optimal under moderate investment effectiveness&#x201d; &#x2014; the mandatory regulations on shore power use (regulatory policies) introduced by Shanghai Port during the same period reduced carbon emissions per unit throughput by 22%, which confirms the theoretical conclusion.</p>
<p>Next, we study the impacts of <inline-formula>
<mml:math display="inline" id="im117">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula>. To distinguish the impact of <inline-formula>
<mml:math display="inline" id="im118">
<mml:mi>r</mml:mi>
</mml:math>
</inline-formula>,we divide it into low <inline-formula>
<mml:math display="inline" id="im119">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> and high values <inline-formula>
<mml:math display="inline" id="im120">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>3</mml:mn>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>, and plot the outcomes in <xref ref-type="fig" rid="f7">
<bold>Figures&#xa0;7a, b</bold>
</xref>.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>The effects of <inline-formula>
<mml:math display="inline" id="im121">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula> on ports&#x2019; profits under three models. <bold>(a)</bold> r=1 <bold>(b)</bold> r=3.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g007.tif">
<alt-text content-type="machine-generated">Two charts labeled (a) and (b) show graphs with lines representing \( B \) (blue), \( S \) (red), and \( R \) (brown) against \( \theta \) on horizontal axes. In chart (a), each line decreases steadily. In chart (b), lines continue to decrease, but the red line has a more pronounced curve, especially as \( \theta \) approaches one.</alt-text>
</graphic>
</fig>
<p>It can be obtained from <xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref> intuitively that ports&#x2019; profit will decrease with the increasing in the proportion of high-emission part of the ports. The subsidy policy always benefits the ports&#x2019; profits. When the investment effectiveness is at a low level, the government&#x2019;s regulation policy will worsen the situation of ports. Further, when the high-emission part accounts for a smaller proportion, and the investment effectiveness is high, the ports situation can be improved, which is in line with Lemma 2. However, once the proportion of high-emission part become larger, it is better for ports&#x2019; profits without regulation policy. The same as <inline-formula>
<mml:math display="inline" id="im122">
<mml:mi>&#x3b8;</mml:mi>
</mml:math>
</inline-formula>, the substitutability of port services <inline-formula>
<mml:math display="inline" id="im123">
<mml:mi>b</mml:mi>
</mml:math>
</inline-formula> also has similar effect on ports&#x2019; profits (see <xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>). Interestingly, when competition between ports is intense, profits tend to be the same in all three models. Moreover, when substitutability between ports is not strong, higher investment effectiveness will be more conducive to port profits.</p>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>The effects of <inline-formula>
<mml:math display="inline" id="im124">
<mml:mi>b</mml:mi>
</mml:math>
</inline-formula> on ports&#x2019; profits under three models.  <bold>(a)</bold> r=1 <bold>(b)</bold> r=3.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmars-12-1661382-g008.tif">
<alt-text content-type="machine-generated">Two graphs labeled (a) and (b) display curves representing B, S, and R plotted against b. In (a), curves are blue (B), brown (S), and red (R), with &#x3c0; on the vertical axis decreasing sharply as b increases from 0 to 1. Graph (b) mirrors this pattern with slight curve adjustments.</alt-text>
</graphic>
</fig>
<p>In recent years, Shanghai Port and Ningbo-Zhoushan Port, both ranking among the top five ports globally in terms of container throughput, have maintained a long-term competitive relationship. Both ports have been committed to the development and improvement of shore power facilities. Ningbo-Zhoushan Port began to promote the construction of shore power facilities in 2017 and implemented them within its port area. Shanghai Port, however, started earlier in advancing shore power facility construction and has a more comprehensive set of shore power facilities. Although both ports have received government support in promoting the development and utilization of shore power facilities, there is a key difference: Shanghai Port explicitly requires ships of certain types and sizes to install shore power interfaces and use shore power supply, while Ningbo-Zhoushan Port has no such regulations. As a result, the shore power utilization rate of Shanghai Port is higher, at approximately 30%, while that of Ningbo-Zhoushan Port is relatively lower, at around 10%. It can be calculated that due to the lack of mandatory regulations (&#x3b8;=0.5, r&#x3b8;=0.5, Ningbo-Zhoushan Port falls into the interval below &#x3b6;<sub>4</sub>. The model indicates that subsidy policies are more effective in this scenario, which is consistent with the actual situation of the port&#x2014;its shore power utilization rate only reaches 10% (lower than that of Shanghai Port) even after relying on subsidies. This reveals the limitations of pure subsidies under low investment effectiveness. Therefore, in practical applications, although subsidies can attract green investment, regulatory policies are a better choice to ensure the sustainable development of green ports.</p>
</sec>
<sec id="s6" sec-type="conclusions">
<label>6</label>
<title>Conclusions</title>
<p>An economic model is developed in this paper to analyze the impact of government policies on green investment in ports. Two commonly used policies, regulatory policy and subsidy policy, are examined and compared. The environmental cost of port emissions is used as a measure of government regulation, which has not been utilized in previous studies. Based on this, the paper analyzes the effects of the two policies on port production and green investments, as well as the optimal choices for the government under different circumstances. Therefore, this paper presents the following theoretical implications and new managerial insights.</p>
<list list-type="roman-lower">
<list-item>
<p>The analysis results show that compared with the situation without government policies, these two policies have the same impact on green investment, but their impacts on port output are opposite. From the perspective of regulation, low investment efficiency is not conducive to green investment in ports. To comply with government regulations, ports will increase investment and restrict production to meet the requirements. The differences in such impacts are particularly evident in the stages of policy implementation: in the initial stage (<inline-formula>
<mml:math display="inline" id="im125">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), to avoid the inhibitory effect of regulation on output, the government will adopt loose regulation, accompanied by high subsidies. For example, Shanghai Port subsidized 0.5 yuan per kilowatt-hour of shore power used in 2019, which not only quickly boosted investment willingness through subsidies but also ensured stable growth of output through loose regulation. This also confirms the inverse relationship between the output-promoting effect of subsidies and the potential output restriction caused by regulation.</p>
</list-item>
<list-item>
<p>Different from government regulation, the existence of subsidies effectively reduces the cost of green investment for ports. Ports always strive to obtain as many subsidies as possible. Therefore, government subsidies can both promote green investment and drive the growth of port output. This effect is most significant in the initial stage (<inline-formula>
<mml:math display="inline" id="im126">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>): a high subsidy rate (s=0.3-0.5) directly covers more than 50% of the green investment cost, greatly lowering the investment threshold for ports. When entering the medium-term stage (<inline-formula>
<mml:math display="inline" id="im127">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&lt;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>): subsidies are gradually reduced (s=0.1-0.2), but the investment foundation laid in the early stage has stabilized the application of green technologies, providing conditions for the subsequent policy transition. This also indicates from the side the key driving role of subsidies in green investment and output during the start-up stage.</p>
</list-item>
<list-item>
<p>When comparing these two mechanisms, subsidies are often a more proactive means, while regulatory policies are more passive. However, both policies can promote social welfare better than no policy intervention. When investment efficiency is low, government subsidies will encourage ports to increase investment. When investment efficiency is high, subsidies serve as rewards for ports. Only when investment efficiency is at a medium level will regulation become a better choice. Although strict regulation will not promote output, ports will increase investment in green technologies. In such cases, subsidy policies may not achieve better results. Specifically, the differences in such policy choices show a clear path in different stages: in the initial stage (<inline-formula>
<mml:math display="inline" id="im128">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) with low investment efficiency, high subsidies (s=0.3-0.5) become the core means, supplemented by loose regulation, which can quickly activate ports&#x2019; investment willingness and ensure stable output growth; in the medium-term stage (<inline-formula>
<mml:math display="inline" id="im129">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>&lt;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> as investment efficiency improves, subsidies are gradually reduced (s=0.1-0.2) and shifted to strict regulation, because the improved investment effectiveness makes regulatory costs lower than subsidy costs; in the mature stage (<inline-formula>
<mml:math display="inline" id="im130">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&gt;</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>7</mml:mn>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), with high investment efficiency, subsidies can be canceled, and market-oriented regulation (such as carbon emission rights trading) becomes the main approach. The Port of Rotterdam, through market-oriented financing for the &#x201c;Carbon Capture and Storage (CCS) Project&#x201d; combined with emission taxes (regulatory policies), achieves the maximization of social welfare.</p>
</list-item>
</list>
<p>Based on these conclusions, we believe that in practical applications, the government should select appropriate policy combinations according to specific circumstances to promote the sustainable development of ports. From the government&#x2019;s perspective, initial subsidy policies can incentivize the steady progress of port greening efforts, while regulatory policies serve as the primary means to consolidate the green development process in subsequent stages. Relying solely on subsidies may undermine social welfare, whereas excessive regulation will impose a burden on port operators. In the early stage, subsidies are necessary to stimulate green investment, but as time goes on and the effects of green measures gradually become apparent, regulatory policies emerge as a better choice. The aforementioned phased policy design (high subsidies + loose regulation in the initial stage, reduced subsidies + strict regulation in the medium term, and no subsidies + market-oriented regulation in the mature stage) is a concrete manifestation of this logic. It not only takes into account the characteristics of investment efficiency in different stages but also achieves a balance between policy costs and social welfare.</p>
<p>Future research can further explore the specific design and implementation of these policies and assess their long-term impacts on the sustainable development of ports. By conducting in-depth studies on the interaction between ports and government policies, we can provide more management insights and decision support for the sustainable development of the port industry. In addition, this paper acknowledges certain unconsidered limitations. Firstly, the government uses public funds to provide subsidies, and these funds could otherwise be allocated to other sectors such as education and healthcare, which might enhance social welfare. The missing element here lies in the shadow price. We can intuitively conclude that once shadow prices are taken into account, the advantages of subsidies over regulation will be significantly diminished. Moreover, when considering regulatory policies, we assume that the regulation of carbon emissions represents the environmental cost; however, such costs are also related to factors such as the social and economic environment, climate, and others. We have also assumed that ports are homogeneous, but in reality, ports vary greatly in terms of scale, throughput, and other aspects, and these factors need to be taken into consideration. These unconsidered issues are important for future research and warrant further investigation.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>FX: Methodology, Writing &#x2013; original draft, Investigation, Conceptualization, Writing &#x2013; review &amp; editing, Formal Analysis. ZZ: Resources, Data curation, Project administration, Methodology, Writing &#x2013; review &amp; editing, Investigation, Conceptualization, Supervision, Funding acquisition.</p>
</sec>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research and/or publication of this article.</p>
</sec>
<sec id="s10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s11" sec-type="ai-statement">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
</sec>
<sec id="s12" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abbas</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Impact of total quality management on corporate green performance through the mediating role of corporate social responsibility</article-title>. <source>J. Cleaner Production</source> <volume>227</volume>, <fpage>118458</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jclepro.2019.118458</pub-id>
</citation></ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Acciaro</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ghiara</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Cusano</surname> <given-names>M. I.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Energy management in seaports: A new role for port authorities</article-title>. <source>Energy Policy</source> <volume>71</volume>, <fpage>4</fpage>&#x2013;<lpage>12</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.enpol.2014.04.013</pub-id>
</citation></ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ayotamuno</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Gobo</surname> <given-names>A. E.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Municipal solid waste management in Port Harcourt, Nigeria: Obstacles and prospects</article-title>. <source>Manage. Environ. Quality: Int. J.</source> <volume>15</volume>, <fpage>389</fpage>&#x2013;<lpage>398</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1108/14777830410540135</pub-id>
</citation></ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balliauw</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Time to build: A real options analysis of port capacity expansion investments under uncertainty</article-title>. <source>Res. Transp. Econ.</source> <volume>83</volume>, <elocation-id>100929</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.retrec.2020.100929</pub-id>
</citation></ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ballini</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Bozzo</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Air pollution from ships in ports: The socio-economic benefit of cold-ironing technology</article-title>. <source>Res. Transportation Business Manage.</source> <volume>17</volume>, <fpage>92</fpage>&#x2013;<lpage>98</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rtbm.2015.10.007</pub-id>
</citation></ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Birch</surname> <given-names>G. F.</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>S. E.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Application of sediment quality guidelines in the assessment and management of contaminated surficial sediments in port jackson (Sydney Harbour), Australia</article-title>. <source>Environ. Manage.</source> <volume>29</volume>, <fpage>860</fpage>&#x2013;<lpage>870</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00267-001-2620-4</pub-id>, PMID: <pub-id pub-id-type="pmid">11992176</pub-id></citation></ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Govindan</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Golias</surname> <given-names>M. M.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Reducing truck emissions at container terminals in a low carbon economy: Proposal of a queueing-based bi-objective model for optimizing truck arrival pattern</article-title>. <source>Transportation Res. Part E: Logistics Transportation Rev.</source> <volume>55</volume>, <fpage>3</fpage>&#x2013;<lpage>22</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tre.2013.03.008</pub-id>
</citation></ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Pak</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>A Delphi analysis on green performance evaluation indices for ports in China</article-title>. <source>Maritime Policy Manage.</source> <volume>44</volume>, <fpage>537</fpage>&#x2013;<lpage>550</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088839.2017.1327726</pub-id>
</citation></ref>
<ref id="B9">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>China Ports &amp; Harbors Association</collab>
</person-group> (<year>2024</year>). <source>2023 Annual report on global container throughput statistics</source>. <publisher-loc>Shanghai</publisher-loc>: <publisher-name>China Ports &amp; Harbors Association</publisher-name>.</citation></ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Farrell</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>The subsidization of seaports: an alternative approach</article-title>. <source>Maritime Policy Manage.</source> <volume>13</volume>, <fpage>177</fpage>&#x2013;<lpage>184</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088838600000061</pub-id>
</citation></ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Finlayson</surname> <given-names>K. A.</given-names>
</name>
<name>
<surname>Limpus</surname> <given-names>C. J.</given-names>
</name>
<name>
<surname>van de Merwe</surname> <given-names>J. P.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Temporal changes in chemical contamination of green turtles (Chelonia mydas) foraging in a heavily industrialised seaport</article-title>. <source>Sci. Total Environ.</source> <volume>817</volume>, <elocation-id>152848</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.scitotenv.2021.152848</pub-id>, PMID: <pub-id pub-id-type="pmid">35007578</pub-id></citation></ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Geerlings</surname> <given-names>H.</given-names>
</name>
<name>
<surname>van Duin</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A new method for assessing CO2-emissions from container terminals: a promising approach applied in Rotterdam</article-title>. <source>J. Cleaner Production</source> <volume>19</volume>, <fpage>657</fpage>&#x2013;<lpage>666</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jclepro.2010.10.012</pub-id>
</citation></ref>
<ref id="B13">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Global Maritime Forum</collab>
</person-group> (<year>2023</year>). <source>Mapping of zero-emission pilots and demonstrations in maritime transport</source> (<publisher-loc>Copenhagen</publisher-loc>: <publisher-name>Global Maritime Forum</publisher-name>).</citation></ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gong</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Seaport investments in capacity and natural disaster prevention</article-title>. <source>Transportation Res. Part D: Transport Environ.</source> <volume>85</volume>, <elocation-id>102367</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2020.102367</pub-id>
</citation></ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hong</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Green product supply chain contracts considering environmental responsibilities</article-title>. <source>Omega</source> <volume>83</volume>, <fpage>155</fpage>&#x2013;<lpage>166</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.omega.2018.02.010</pub-id>
</citation></ref>
<ref id="B16">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>IMO</collab>
</person-group> (<year>2023</year>). Available online at: <uri xlink:href="https://www.imo.org/en/OurWork/Environment/Pages/2023-IMO-Strategy-on-Reduction-of-GHG-Emissions-from-Ships.aspx">https://www.imo.org/en/OurWork/Environment/Pages/2023-IMO-Strategy-on-Reduction-of-GHG-Emissions-from-Ships.aspx</uri> (Accessed <access-date>July 29, 2025</access-date>).</citation></ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Janipour</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Swennenhuis</surname> <given-names>F.</given-names>
</name>
<name>
<surname>de Gooyert</surname> <given-names>V.</given-names>
</name>
<name>
<surname>de Coninck</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Understanding contrasting narratives on carbon dioxide capture and storage for Dutch industry using system dynamics</article-title>. <source>Int. J. Greenhouse Gas Control</source> <volume>105</volume>, <elocation-id>103235</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijggc.2020.103235</pub-id>
</citation></ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kizielewicz</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Monitoring energy efficiency and environmental ship index by cruise seaports in Northern Europe</article-title>. <source>Energies</source> <volume>15</volume>, <elocation-id>4215</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/en15124215</pub-id>
</citation></ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kotrikla</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Lilas</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Nikitakos</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Abatement of air pollution at an aegean island port utilizing shore side electricity and renewable energy</article-title>. <source>Mar. Policy</source> <volume>75</volume>, <fpage>238</fpage>&#x2013;<lpage>248</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.marpol.2016.01.026</pub-id>
</citation></ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lai</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Tao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Sustainability investment in maritime supply chain with risk behavior and information sharing</article-title>. <source>Int. J. production economics</source> <volume>218</volume>, <fpage>16</fpage>&#x2013;<lpage>29</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ijpe.2019.02.021</pub-id>
</citation></ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lam</surname> <given-names>J. S. L.</given-names>
</name>
<name>
<surname>Notteboom</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The greening of ports: A comparison of port management tools used by leading ports in Asia and Europe</article-title>. <source>Transport Rev.</source> <volume>34</volume>, <fpage>169</fpage>&#x2013;<lpage>189</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01441647.2014.891162</pub-id>
</citation></ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Kuang</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Using system dynamics and game model to estimate optimal subsidy in shore power technology</article-title>. <source>IEEE Access</source> <volume>8</volume>, <fpage>116310</fpage>&#x2013;<lpage>116320</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1109/ACCESS.2020.3004183</pub-id>
</citation></ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>T.-K.</given-names>
</name>
<name>
<surname>Sheu</surname> <given-names>H.-Y.</given-names>
</name>
<name>
<surname>Tsai</surname> <given-names>J.-Y.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Sulfur dioxide emission estimates from merchant vessels in a port area and related control strategies</article-title>. <source>Aerosol Air Qual. Res.</source> <volume>14</volume>, <fpage>413</fpage>&#x2013;<lpage>421</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.4209/aaqr.2013.02.0061</pub-id>
</citation></ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Research on investment strategy of maritime supply chain based on green investment efficiency</article-title>. <source>Chin. Manage. Sci.</source> <volume>29</volume>, <fpage>33</fpage>&#x2013;<lpage>44</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.16381/j.cnki.issn1003-207x.2019.1898</pub-id>
</citation></ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Optimization of shore power deployment in green ports considering government subsidies</article-title>. <source>Sustainability</source> <volume>13</volume>, <elocation-id>1640</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/su13041640</pub-id>
</citation></ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Luo</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Optimisation of port investment strategies considering uncertain cargo demand and environmental pollutions</article-title>. <source>Comput. Ind. Eng.</source> <volume>200</volume>, <fpage>110788</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cie.2024.110788</pub-id>
</citation></ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mahmud</surname> <given-names>K. K.</given-names>
</name>
<name>
<surname>Chowdhury</surname> <given-names>M. M. H.</given-names>
</name>
<name>
<surname>Shaheen</surname> <given-names>M. M. A.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Green port management practices for sustainable port operations: a multi method study of Asian ports</article-title>. <source>Maritime Policy Manage.</source> <volume>51</volume>, 1902&#x2013;1937. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088839.2023.2258125</pub-id>
</citation></ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Metzger</surname> <given-names>D.</given-names>
</name>
<name>
<surname>SChinas</surname> <given-names>O.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Fuzzy real options and shared savings: Investment appraisal for green shipping technologies</article-title>. <source>Transportation Res. Part D: Transport Environ.</source> <volume>77</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2019.09.016</pub-id>
</citation></ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Notteboom</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Lam</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The greening of terminal concessions in seaports</article-title>. <source>Sustainability</source> <volume>10</volume>, <elocation-id>3318</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/su10093318</pub-id>
</citation></ref>
<ref id="B30">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>Organisation for Economic Co-operation and Development (OECD)</collab>
</person-group> (<year>2023</year>). <source>Climate Policy Instruments Monitoring Report 2023</source> (<publisher-loc>Paris</publisher-loc>: <publisher-name>OECD Publishing</publisher-name>).</citation></ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Park</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Chang</surname> <given-names>Y.-T.</given-names>
</name>
<name>
<surname>Zou</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Emission control under private port operator duopoly</article-title>. <source>Transportation Res. Part E: Logistics Transportation Rev.</source> <volume>114</volume>, <fpage>40</fpage>&#x2013;<lpage>65</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tre.2018.03.010</pub-id>
</citation></ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Cooperative optimization of shore power allocation and berth allocation: A balance between cost and environmental benefit</article-title>. <source>J. Cleaner Production</source> <volume>279</volume>, <elocation-id>123816</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jclepro.2020.123816</pub-id>
</citation></ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Radwan</surname> <given-names>M. E.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wan</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Hua</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Critical barriers to the introduction of shore power supply for green port development: case of Djibouti container terminals</article-title>. <source>Clean Technol. Environ. Policy</source> <volume>21</volume>, <fpage>1293</fpage>&#x2013;<lpage>1306</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10098-019-01706-z</pub-id>
</citation></ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Randrianarisoa</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Gillen</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Reducing emissions in international transport: A supply chain perspective. Transportation research, Part D</article-title>. <source>Transport Environ.</source> <volume>102</volume>, <elocation-id>103074</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2021.103074</pub-id>
</citation></ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ren</surname> <given-names>S.</given-names>
</name>
<name>
<surname>He</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Symbolic reactions or substantive pro-environmental behaviour? An empirical study of corporate environmental performance under the government&#x2019;s environmental subsidy scheme</article-title>. <source>Business strategy Environ.</source> <volume>28</volume>, <fpage>845</fpage>&#x2013;<lpage>860</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/bse.v28.6</pub-id>
</citation></ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schodler</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Saraceni</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Factors influencing greenhouse gas reduction measures in european ports: implications for sustainable investing</article-title>. <source>J. Risk Financial Manage.</source> <volume>17</volume>, <elocation-id>329</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/jrfm17080329</pub-id>
</citation></ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sciberras</surname> <given-names>E. A.</given-names>
</name>
<name>
<surname>Zahawi</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Atkinson</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Electrical characteristics of cold ironing energy supply for berthed ships</article-title>. <source>Transportation Res. Part D Transport Environ.</source> <volume>39</volume>, <fpage>31</fpage>&#x2013;<lpage>43</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2015.05.007</pub-id>
</citation></ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sendstad</surname> <given-names>L. H.</given-names>
</name>
<name>
<surname>Hagspiel</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Mikkelsen</surname> <given-names>W. J.</given-names>
</name>
<name>
<surname>Ravndal</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Tveitst&#xf8;l</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The impact of subsidy retraction on European renewable energy investments</article-title>. <source>Energy Policy</source> <volume>160</volume>, <elocation-id>112675</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.enpol.2021.112675</pub-id>
</citation></ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sogut</surname> <given-names>M. Z.</given-names>
</name>
<name>
<surname>Erdo&#x11f;an</surname> <given-names>O.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>An investigation on a holistic framework of green port transition based on energy and environmental sustainability</article-title>. <source>Ocean Eng.</source> <volume>266</volume>, <elocation-id>112671</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.oceaneng.2022.112671</pub-id>
</citation></ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>He</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>An</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Green pharmaceutical supply chain coordination considering green investment, green logistics, and government intervention</article-title>. <source>Environ. Sci. pollut. Res. Int.</source> <volume>29</volume>, <fpage>63321</fpage>&#x2013;<lpage>63343</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11356-021-18275-8</pub-id>, PMID: <pub-id pub-id-type="pmid">35451714</pub-id></citation></ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stojanovic</surname> <given-names>T. A.</given-names>
</name>
<name>
<surname>Ormerod Smith</surname> <given-names>H. D.</given-names>
</name>
<name>
<surname>Wooldridge</surname> <given-names>C. F.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>The impact of the Habitats Directive on European port operations and management</article-title>. <source>GeoJournal</source> <volume>65</volume>, <fpage>165</fpage>&#x2013;<lpage>176</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10708-006-0004-2</pub-id>
</citation></ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tan</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Zeng</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Capacity investment of shore power berths for a container port: Environmental incentive and infrastructure subsidy policies</article-title>. <source>Ocean Coast. Manage.</source> <volume>239</volume>, <elocation-id>106582</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ocecoaman.2023.106582</pub-id>
</citation></ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Optimal operation of photovoltaic/battery/diesel/cold-ironing hybrid energy system for maritime application</article-title>. <source>Energy</source> <volume>162</volume>, <fpage>697</fpage>&#x2013;<lpage>714</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.energy.2018.08.048</pub-id>
</citation></ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tian</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Zhen</surname> <given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A new assessment model ofenvironmental cost of carbon and its situation analysis in China</article-title>. <source>J. Cleaner Production</source> <volume>211</volume>, <fpage>1434</fpage>&#x2013;<lpage>1443</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jclepro.2018.11.117</pub-id>
</citation></ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tseng</surname> <given-names>P.-H.</given-names>
</name>
<name>
<surname>Pilcher</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Evaluating the key factors of green port policies in Taiwan through quantitative and qualitative approaches</article-title>. <source>Transport Policy</source> <volume>82</volume>, <fpage>127</fpage>&#x2013;<lpage>137</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tranpol.2018.12.014</pub-id>
</citation></ref>
<ref id="B46">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>UTCAD</collab>
</person-group> (<year>2023</year>). Available online at: <uri xlink:href="https://unctad.org/publication/annual-report-2023">https://unctad.org/publication/annual-report-2023</uri> (Accessed <access-date>July 29, 2025</access-date>).</citation></ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Corbett</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Winebrake</surname> <given-names>J. J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Cost-effectiveness of reducing sulfur emissions from ships</article-title>. <source>Environ. Sci. Technol.</source> <volume>41</volume>, <fpage>8233</fpage>&#x2013;<lpage>8239</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1021/es070812w</pub-id>, PMID: <pub-id pub-id-type="pmid">18200845</pub-id></citation></ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Jiao</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Shipping companies&#x2019; choice of low sulfur fuel oil with government subsidy and different maritime supply chain power structures</article-title>. <source>Maritime Policy Manage.</source> <volume>49</volume>, <fpage>323</fpage>&#x2013;<lpage>346</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088839.2021.1950860</pub-id>
</citation></ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Winnes</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Styhre</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Fridell</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Reducing GHG emissions from ships in port areas</article-title>. <source>Res. Transportation Business Manage.</source> <volume>17</volume>, <fpage>73</fpage>&#x2013;<lpage>82</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.rtbm.2015.10.008</pub-id>
</citation></ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Lindsey</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Port adaptation to climate change and capacity investments under uncertainty</article-title>. <source>Transportation Res. Part B: Methodological</source> <volume>152</volume>, <fpage>180</fpage>&#x2013;<lpage>204</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trb.2021.08.009</pub-id>
</citation></ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Mishra</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Adulyasak</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Seaport adaptation and capacity investments under inter-port competition and climate-change uncertainty</article-title>. <source>Transportation Res. Part D: Transport Environ.</source> <volume>130</volume>, <elocation-id>104183</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2024.104183</pub-id>
</citation></ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Whether to invest in terminal efficiency: A perspective considering customer preference and capital constraint in competitive environment</article-title>? <source>Ocean Coast. Manage.</source> <volume>205</volume>, <elocation-id>105563</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ocecoaman.2021.105563</pub-id>
</citation></ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J</given-names>
</name>
</person-group>. (<year>2025</year>a). <article-title>Evaluation of navigation capacity in the Northeast Arctic passage: evidence from multiple factors</article-title>. <source>Maritime Policy Manage</source> <volume>52</volume>, <fpage>497</fpage>&#x2013;<lpage>513</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088839.2024.2376126</pub-id>
</citation></ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Shen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2025</year>b). <article-title>The impact of the Maritime Silk Road Initiative on the carbon intensity of the participating countries</article-title>. <source>Maritime Economics Logistics</source> <volume>27</volume>, <fpage>349</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1057/s41278-024-00295-z</pub-id>
</citation></ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>C. Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y. Q.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Qiu</surname> <given-names>S. Y.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2025</year>c). <article-title>Research on the coordination of a marine green fuel supply chain considering a cost-sharing contract and a revenue-sharing contract</article-title>. <source>Front. Mar. Sci.</source> <volume>12</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmars.2025.1552136</pub-id>
</citation></ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Passive or proactive capacity sharing? A perspective of cooperation and competition between two regional ports: Article</article-title>. <source>Maritime Policy Manage.</source> <volume>49</volume>, <fpage>492</fpage>&#x2013;<lpage>509</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/03088839.2021.1876938</pub-id>
</citation></ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>Y.-C.</given-names>
</name>
<name>
<surname>Lin</surname> <given-names>C.-L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Performance analysis of cargo-handling equipment from a green container terminal perspective</article-title>. <source>Transportation Res. Part D: Transport Environ.</source> <volume>23</volume>, <fpage>9</fpage>&#x2013;<lpage>11</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trd.2013.03.009</pub-id>
</citation></ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Ge</surname> <given-names>Y.-E.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Nie</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Modeling collusion-proof port emission regulation of cargo-handling activities under incomplete information</article-title>. <source>Transportation Res. Part B: Methodological</source> <volume>104</volume>, <fpage>543</fpage>&#x2013;<lpage>567</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trb.2017.04.015</pub-id>
</citation></ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zheng</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Z.-C.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Chan</surname> <given-names>F. T. S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Subsidy or minimum requirement? Regulation of port adaptation investment under disaster ambiguity</article-title>. <source>Transportation Res. Part B: Methodological</source> <volume>150</volume>, <fpage>457</fpage>&#x2013;<lpage>481</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.trb.2021.04.006</pub-id>
</citation></ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>The effects of fuel content regulation at ports on regional pollution and shipping industry</article-title>. <source>J. Environ. Economics Manage.</source> <volume>106</volume>, <elocation-id>102424</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jeem.2021.102424</pub-id>
</citation></ref>
</ref-list>
<app-group>
<app id="app1">
<title>Appendix 1. Proofs of Market equilibriums in Section 4</title>
<p>4.1 benchmark</p>
<p>To expand the expected profits of ports <inline-formula>
<mml:math display="inline" id="im131">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> in <xref ref-type="disp-formula" rid="eq3">Equation 3</xref>, we have</p>
<disp-formula id="eq6">
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>a</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>For convenience, we assume <inline-formula>
<mml:math display="inline" id="im132">
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>When <inline-formula>
<mml:math display="inline" id="im133">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. Then we maximize the objective function <xref ref-type="disp-formula" rid="eq6">Equation 6</xref> with respect to <inline-formula>
<mml:math display="inline" id="im134">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, we obtain <inline-formula>
<mml:math display="inline" id="im135">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>. Substituting it into <xref ref-type="disp-formula" rid="eq6">Equation 6</xref> and imposing symmetry on <inline-formula>
<mml:math display="inline" id="im136">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im137">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>. Substituting <inline-formula>
<mml:math display="inline" id="im138">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im139">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> into <inline-formula>
<mml:math display="inline" id="im140">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, we obtain <inline-formula>
<mml:math display="inline" id="im141">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>When</p>
<disp-formula id="eq7">
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</disp-formula>
<p>we know that <inline-formula>
<mml:math display="inline" id="im142">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im143">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and the <xref ref-type="disp-formula" rid="eq7">Equation 7</xref> becomes:</p>
<disp-formula id="eq8">
<label>(8)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>a</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Then maximizing the <xref ref-type="disp-formula" rid="eq8">Equation 8</xref>, we have <inline-formula>
<mml:math display="inline" id="im144">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im145">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mi>k</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im146">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>When <inline-formula>
<mml:math display="inline" id="im147">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, the objective function is the same as <xref ref-type="disp-formula" rid="eq8">Equation 8</xref>. Maximizing it with respect to <inline-formula>
<mml:math display="inline" id="im148">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> we obtain <inline-formula>
<mml:math display="inline" id="im149">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and imposing the symmetry, we obtain <inline-formula>
<mml:math display="inline" id="im150">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>. Then substituting <inline-formula>
<mml:math display="inline" id="im151">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im152">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> into <inline-formula>
<mml:math display="inline" id="im153">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im154">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>When <inline-formula>
<mml:math display="inline" id="im155">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, we have</p>
<disp-formula id="eq9">
<label>(9)</label>
<mml:math display="block" id="M9">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>r</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula id="eq10">
<label>(10)</label>
<mml:math display="block" id="M10">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>a</mml:mi>
<mml:mstyle displaystyle="true">
<mml:msubsup>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:msubsup>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>2</mml:mn>
</mml:mfrac>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Substituting <xref ref-type="disp-formula" rid="eq9">Equation 9</xref> into <xref ref-type="disp-formula" rid="eq10">Equation 10</xref>, and maximizing it with respect to <inline-formula>
<mml:math display="inline" id="im156">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im157">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>B</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, and derive <inline-formula>
<mml:math display="inline" id="im158">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo stretchy="false">/</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula> through <inline-formula>
<mml:math display="inline" id="im159">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>4.2.1 Regulation policy</p>
<p>In this section, the specific value of minimum regulation has been explored. Then we present the proofs within the feasible regions of R.</p>
<p>When <inline-formula>
<mml:math display="inline" id="im160">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im161">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mi>t</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im162">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and obtain <inline-formula>
<mml:math display="inline" id="im163">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>. Substituting it into Eq. (4), and maximizing it with respect to <inline-formula>
<mml:math display="inline" id="im164">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im165">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>+</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>a</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mo>;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula>Then we maximize the social welfare with respect to <inline-formula>
<mml:math display="inline" id="im166">
<mml:mi>R</mml:mi>
</mml:math>
</inline-formula>, and we have <inline-formula>
<mml:math display="inline" id="im167">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:msup>
<mml:mi>b</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<mml:math display="inline" id="im179">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
<mml:math display="inline" id="im180">
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">/</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im168">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>a</mml:mi>
<mml:mi>k</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>When <inline-formula>
<mml:math display="inline" id="im169">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>L</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im170">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>&lt;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>k</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>k</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im171">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>r</mml:mi>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and substitute <inline-formula>
<mml:math display="inline" id="im172">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>R</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
<mml:mi>r</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> into <xref ref-type="disp-formula" rid="eq4">Equation 4</xref>. Then maximizing it with respect to <inline-formula>
<mml:math display="inline" id="im173">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, we have <inline-formula>
<mml:math display="inline" id="im174">
<mml:mrow>
<mml:msubsup>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>R</mml:mi>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>R</mml:mi>
</mml:msubsup>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>r</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:mi>&#x3b8;</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im175">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>4</mml:mn>
<mml:msup>
<mml:mi>b</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>4</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>5</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>+</mml:mo>
<mml:mi>t</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>t</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>b</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>+</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mi>t</mml:mi>
<mml:msubsup>
<mml:mi>C</mml:mi>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msubsup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>, where <inline-formula>
<mml:math display="inline" id="im176">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mo>&#x2208;</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3b6;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:mi>&#x221e;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
<p>4.2.2Subsidization policy</p>
<p>Maximizing the expected social welfare of <xref ref-type="disp-formula" rid="eq4">Equation 4</xref> with respect to <inline-formula>
<mml:math display="inline" id="im177">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula>
<mml:math display="inline" id="im178">
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and then substituting them into <xref ref-type="disp-formula" rid="eq4">Equation 4</xref> to obtain the optimal subsidy. The proofs are based on the same logic presented above and omitted.</p>
</app>
</app-group>
</back>
</article>