<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Environ. Sci.</journal-id>
<journal-title>Frontiers in Environmental Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Environ. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-665X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1095394</article-id>
<article-id pub-id-type="doi">10.3389/fenvs.2022.1095394</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Environmental Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Study on the coupling coordination and pattern evolution of green investment and ecological development: Based on spatial econometric model and China&#x2019;s provincial panel data</article-title>
<alt-title alt-title-type="left-running-head">Peng et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenvs.2022.1095394">10.3389/fenvs.2022.1095394</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Peng</surname>
<given-names>Ruikun</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2094128/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Huichun</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ge</surname>
<given-names>Jiaojiao</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yan</surname>
<given-names>Weimin</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Finance</institution>, <institution>Nanjing Agricultural University</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Management Science and Engineering</institution>, <institution>Nanjing University of Information Science and Technology</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1532718/overview">Yuantao Xie</ext-link>, University of International Business and Economics, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1951026/overview">Danish Ahmed</ext-link>, Shanghai Jianqiao University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1846673/overview">Maham Furqan</ext-link>, Oregon State University, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Ruikun Peng, <email>tianya_1978@126.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Environmental Economics and Management, a section of the journal Frontiers in Environmental Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>03</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>1095394</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>30</day>
<month>11</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Peng, Huang, Ge and Yan.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Peng, Huang, Ge and Yan</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 study determines the coupling degree of green investment and green ecology in China using kernel density estimation, spatial autocorrelation analysis, and standard deviation ellipse model to empirically evaluate the data of 30 Chinese provinces from 2005 to 2019. Moreover, the study investigates the temporal evolution trend, spatial clustering characteristics, and spatial evolution trend of coupling degree. Although the coupling coordination value of green investment and green ecology gradually increased, it is generally at a low coordination stage. At the same time, regional disparities narrowed with the most significant variability in the eastern region. Moreover, results found that the aggregation effect of the coupling and coordination of green investment and green ecology is more significant, and the high-value aggregation area extended from the lower reaches of the Yangtze River to the midstream region, while the western region is dominated by low-value aggregation. Similarly, the degree of synergy between green investment and green ecology is increased; however, the fragmentation trend is inevitable. At the same time, the center of gravity of coupling coordination shifted to the western regions, and the spatial pattern gradually weakened in the &#x201c;northeast-southwest&#x201d; direction. The findings of the study stress that local governments of China should improve the green investment system for green ecological development in the surrounding areas of the Yangtze River Economic Belt. Moreover, it is required to focus on the Northeast revitalization and Western development strategies to promote the synergistic development of green investment and green ecology.</p>
</abstract>
<kwd-group>
<kwd>green investment</kwd>
<kwd>green ecology</kwd>
<kwd>coupling coordination</kwd>
<kwd>spatial autocorrelation</kwd>
<kwd>China</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<sec id="s1-1">
<title>1.1 Background</title>
<p>In parallel with global economic growth, environmental issues have increased (<xref ref-type="bibr" rid="B23">Khan et al., 2020</xref>). About 8.3 million people die each year due to environmental pollution (<xref ref-type="bibr" rid="B22">Khan &#x26; Ozturk, 2020</xref>) and most cases were found in China as it has ranked second in the world (<xref ref-type="bibr" rid="B28">Landrigan et al., 2018</xref>; <xref ref-type="bibr" rid="B17">Gu et al., 2019</xref>; <xref ref-type="bibr" rid="B18">Gu et al., 2020</xref>). Therefore, it is indispensable to advocate green ecological development in China, as an important for social welfare and human health (<xref ref-type="bibr" rid="B43">Wang et al., 2020</xref>). Green investment has involved in the field of renewable energy (<xref ref-type="bibr" rid="B45">Yao et al., 2016</xref>; <xref ref-type="bibr" rid="B10">Elahi et al., 2019</xref>; <xref ref-type="bibr" rid="B9">Elahi et al., 2022</xref>; <xref ref-type="bibr" rid="B11">Elahi et al., 2022</xref>) and promotes green economy (<xref ref-type="bibr" rid="B8">Ekeh et al., 2007</xref>), which also caters to the strong public demand (<xref ref-type="bibr" rid="B31">Liao and Shi, 2018</xref>). May studies have found a positive nexus between green investment and green ecology. <xref ref-type="bibr" rid="B39">Shen et al. (2021)</xref> studied the role of natural resource rents, green investment, financial development, and energy consumption in carbon emission reduction using panel data of Chinese provinces in recent years, and concluded that green investment is negatively related to carbon emissions, suggesting that green investment should be promoted to control carbon emissions. <xref ref-type="bibr" rid="B52">Zhu et al. (2014)</xref> suggested that environmental regulations as well as green investment can motivate pollution-intensive enterprises in China to maintain green ecology through technological innovation, industrial structure upgrading, and plant relocation. However, some scholars argue that there is a non-linear correlation between the two, that the effect of carbon neutrality and green transformation has a lag, that short-term green investment initiatives will increase the proportion of national eco-industrial R&#x26;D investment and weaken green ecological protection, and that the impact of green investment efforts on China&#x2019;s green ecological quality shows a &#x201c;U&#x201d; shaped change (<xref ref-type="bibr" rid="B14">Fuyong and Xiang, 2022</xref>). There are different views on the study of the coupling mechanism between green investment and green ecology. Many gaps in the quantitative studies on coupling coordination, such as coupling trends and the evolution of the dynamics, have been found. To the best of our knowledge, this is the first attempt to explore the evolutionary trend of green investment and green ecology coupling from temporal and spatial dimensions. Moreover, it provides ideas for green investment and green ecology synergistic development policies in China.</p>
</sec>
<sec id="s1-2">
<title>1.2 Literature review and contribution of the study</title>
<sec id="s1-2-1">
<title>1.2.1 Green investment promotes green ecology</title>
<p>Some scholars interpret green investment as &#x201c;socially responsible investment&#x201d;, i.e., investment based on a combination of environmental, social, and economic criteria, which is part of &#x201c;green finance.&#x201d; Others also innovatively interpreted green investment as the &#x201c;environmental protection investment.&#x201d; <xref ref-type="bibr" rid="B48">Zhang et al. (2022)</xref> used the entropy method and undesirable-SE-SBM model to measure provincial green finance and green development efficiency respectively from 2008 to 2018, and the results showed that green finance can significantly promote green development and improve green ecology after a certain percentage of R&#x26;D investment. <xref ref-type="bibr" rid="B34">Mesjasz-Lech (2017)</xref> analyzed the relationship between environmental protection expenditures and environmental governance effectiveness for sustainable development using principal component analysis for EU industrial sector related firms and observed that the increase in environmental protection expenditures in most countries was accompanied by improved environmental governance effectiveness. Green investment plays its capital allocation role to provide financial support to green industries for technological innovation, which leads to energy consumption reduction and pollution reduction, and improvement of green ecosystem.</p>
</sec>
<sec id="s1-2-2">
<title>1.2.2 Impact of green ecology on the efficiency of green investment</title>
<p>Green ecology is an important development approach for social progress and green investment efficiency. Most of the existing literature focused on the impact of green ecology on economic growth (<xref ref-type="bibr" rid="B35">Peng et al., 2021</xref>; <xref ref-type="bibr" rid="B37">Peng et al., 2021</xref>; <xref ref-type="bibr" rid="B49">Zhao et al. 2020</xref>; <xref ref-type="bibr" rid="B5">Aldieri et al., 2020</xref>) reported that environmental regulations impacted the distribution of financial resources. <xref ref-type="bibr" rid="B42">Vogel (2000)</xref> justified that economic integration and green ecological regulation are compatible and that economic integration can provide a certain degree of an effective mechanism for raising regulatory standards. From the perspective of enterprises, <xref ref-type="bibr" rid="B29">Leiter et al. (2011)</xref>, <xref ref-type="bibr" rid="B41">Testa et al. (2011)</xref>, and <xref ref-type="bibr" rid="B52">Zhu et al. (2014)</xref> argued that a series of government management tools such as environmental controls, financial subsidies, and collection of emission fees can significantly optimize the production methods of enterprises and strengthen their level of green investment. Improving the level of green ecology can effectively enhance corporate investment in energy conservation and emission reduction, thus enhancing environmental protection investment, while green ecology can also enhance corporate environmental protection revenue by increasing corporate investment in energy conservation and emission reduction (<xref ref-type="bibr" rid="B7">Chesney et al. 2011</xref>; <xref ref-type="bibr" rid="B20">Jaraite et al. 2014</xref>).</p>
</sec>
<sec id="s1-2-3">
<title>1.2.3 Coordination between green investment and green ecology</title>
<p>Despite the green investment is a subdivision of the economy in the field of capital financing, a few studies focused on the coordination relationship between green investment and green ecology. Most of the existing literature focuses on the harmonious relationship between economic growth and green ecology. <xref ref-type="bibr" rid="B26">Kline (2000)</xref> believed that high-level economic growth and a high-quality environment can coexist in ecological cities. <xref ref-type="bibr" rid="B13">Feiock and Stream (2001)</xref> found that economic growth and the environment should be integrated and that they should have a harmonious development relationship to promote each other. <xref ref-type="bibr" rid="B15">Gan and Bu (2020)</xref> revealed that the coupling coordination between green finance and the eco-environment is developing from an uncoordinated level to a well-coordinated level, but failed to analyze their coupled coordination relationship in depth. On the one hand, <xref ref-type="bibr" rid="B32">Liu et al. (2018)</xref> argued that the social welfare brought by economic development can significantly improve the quality of the environment and gradually lead to the benign development of the overall environment. On the other hand, the gradual improvement of the environment can optimize the economic structure, stimulate economic vitality and adjust the economic scale.</p>
<p>By comparison of the above studies, it is found that most of the existing studies on green investment and green ecology focus on the improvement effect of green investment on green ecology and the promotion effect of green ecology on green investment efficiency, and there are still research gaps on issues of the coupling and coordination of green investment and green ecology, which need to be further explored. The existing studies on the spatial distribution of coupling and coordination mostly focused on the field of spatial clustering, and lack of in-depth exploration of the evolution of spatial and temporal patterns. Therefore, the current study selects general indicators of green investment and green ecology in 30 Chinese provinces and explores the coupling characteristics of green investment and green ecology in China from 2005 to 2019. While studying the temporal evolution trend of coupling coordination, it also delves into the spatial evolution characteristics of coupling coordination and examines whether there is a coupling effect between green investment and green ecology. Furthermore, the study provides ideas for synergistic development of green investment and green ecology in China.</p>
</sec>
</sec>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Research methods</title>
<p>Common methods of variable modeling studies include panel data (<xref ref-type="bibr" rid="B24">Khelfaoui et al., 2022</xref>; <xref ref-type="bibr" rid="B25">Khelfaoui et al., 2022</xref>), bias correction (<xref ref-type="bibr" rid="B3">Ahmed et al, 2021</xref>), structural equation modeling (<xref ref-type="bibr" rid="B2">Ahmed et al., 2021</xref>), and sample analysis (<xref ref-type="bibr" rid="B1">Ahmed et al., 2019</xref>), among others. The coordination model is used to estimate the degree of interaction between green investment and green ecology. Using the kernel density function, another estimation method is constructed to study the time-series evolution of the interaction between green investment and green ecological development. In addition, we used the spatial autocorrelation analysis method (<italic>Moran&#x2019;s I index, Getis-Ord Gi&#x2a; index</italic>) for exploratory spatial data analysis (ESDA) and the spatial standard deviational ellipse model to analyze the spatial global characteristics and local spatial clustering differences in the coordination of green investment and green ecology coupling.</p>
<sec id="s2-1-1">
<title>2.1.1 Pre-processing of data</title>
<p>As the indicators of green investment and green ecology differ in terms of units and areas of expertise. Therefore, to eliminate the influence of the innate defects of the system indicators on the empirical results, the values of the indicators should be carried out using dimensionless processing. Based on this theory, the indicators were standardized into positive and negative indicators using maximum difference denormalization (<xref ref-type="bibr" rid="B44">Wang et al., 2019</xref>).<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">Min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">Max</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">Min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">h</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">p</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">v</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="italic">Max</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">Max</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="italic">Min</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">h</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">v</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the value of indicator <italic>i</italic> in the <italic>j-</italic>th year. <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:mi>M</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">x</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> represent the minimum and maximum values of the indicator among all indicators in the <italic>j-</italic>th year, respectively. <inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the final standard value. To eliminate the interference of human factors on the index values, the standardized index <inline-formula id="inf5">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>Z</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is processed by the entropy method to obtain the weights <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B6">Cao et al., 2020</xref>). The results of the weights are given in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Construction of a systematic index of green ecology and green investment.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left"/>
<th align="left">System layer</th>
<th align="left">Element layer</th>
<th align="left">Indicator layer (indicator code)</th>
<th align="left">Unit</th>
<th align="left">Indicator attribute</th>
<th align="left">Weight</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="22" align="left">Green ecology</td>
<td rowspan="12" align="left">Green growth</td>
<td rowspan="2" align="left">Economic growth</td>
<td align="left">GDP <italic>per capita</italic> (X1)</td>
<td align="left">Yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0367</td>
</tr>
<tr>
<td align="left">Local fiscal revenue <italic>per capita</italic> (X2)</td>
<td align="left">Yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0336</td>
</tr>
<tr>
<td rowspan="3" align="left">Ecological economy</td>
<td align="left">Energy consumption per tenthousand yuan of GDP (X3)</td>
<td align="left">Tons of standard coal/million yuan</td>
<td align="left">&#x2212;</td>
<td align="left">0.0401</td>
</tr>
<tr>
<td align="left">Industrial SO2 emissions per ten thousand yuan of GDP (X4)</td>
<td align="left">kg/million yuan</td>
<td align="left">&#x2212;</td>
<td align="left">0.0403</td>
</tr>
<tr>
<td align="left">Industrial wastewater discharge per ten thousand yuan of GDP (X5)</td>
<td align="left">kg/million yuan</td>
<td align="left">&#x2212;</td>
<td align="left">0.0402</td>
</tr>
<tr>
<td rowspan="4" align="left">Output level</td>
<td align="left">Labor productivity in primary industry (X6)</td>
<td align="left">Yuan/person</td>
<td align="left">&#x2b;</td>
<td align="left">0.0377</td>
</tr>
<tr>
<td align="left">Labor productivity in secondary industry (X7)</td>
<td align="left">Yuan/person</td>
<td align="left">&#x2b;</td>
<td align="left">0.0373</td>
</tr>
<tr>
<td align="left">Labor productivity of tertiary industry (X8)</td>
<td align="left">Yuan/person</td>
<td align="left">&#x2b;</td>
<td align="left">0.0358</td>
</tr>
<tr>
<td align="left">Value-added ratio of tertiary industry (X9)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0378</td>
</tr>
<tr>
<td rowspan="3" align="left">Standard of living</td>
<td align="left">Per capita disposable income of urban residents (X10)</td>
<td align="left">Yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0364</td>
</tr>
<tr>
<td align="left">Net income of rural residents (X11)</td>
<td align="left">Yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0365</td>
</tr>
<tr>
<td align="left">Per capita housing area (X12)</td>
<td align="left">m<sup>2</sup>
</td>
<td align="left">&#x2b;</td>
<td align="left">0.0376</td>
</tr>
<tr>
<td rowspan="6" align="left">Green welfare</td>
<td rowspan="4" align="left">Public service</td>
<td align="left">Per capita park green space (X13)</td>
<td align="left">m<sup>2</sup>/person</td>
<td align="left">&#x2b;</td>
<td align="left">0.0388</td>
</tr>
<tr>
<td align="left">Green coverage rate of built-up areas (X14)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0400</td>
</tr>
<tr>
<td align="left">Decontamination rate of urban refuse (X15)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0381</td>
</tr>
<tr>
<td align="left">Comprehensive utilization rate of industrialsolid waste (X16)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0396</td>
</tr>
<tr>
<td rowspan="2" align="left">Resource consumption</td>
<td align="left">Per capita domestic water consumption (X17)</td>
<td align="left">L/person</td>
<td align="left">&#x2212;</td>
<td align="left">0.0404</td>
</tr>
<tr>
<td align="left">Per capita energy consumption (X18)</td>
<td align="left">kg standard coal/person</td>
<td align="left">&#x2212;</td>
<td align="left">0.0404</td>
</tr>
<tr>
<td rowspan="4" align="left">Resources and environment</td>
<td rowspan="2" align="left">Environmental pressure</td>
<td align="left">Industrial SO2 emission per unit land area (X19)</td>
<td align="left">t/km<sup>2</sup>
</td>
<td align="left">&#x2212;</td>
<td align="left">0.0404</td>
</tr>
<tr>
<td align="left">Industrial wastewater emission per unit land area (X20)</td>
<td align="left">t/km<sup>2</sup>
</td>
<td align="left">&#x2212;</td>
<td align="left">0.0401</td>
</tr>
<tr>
<td rowspan="2" align="left">Ecological quality</td>
<td align="left">Forest coverage rate (X21)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0371</td>
</tr>
<tr>
<td align="left">The ratio of the area of nature reservesto the area under the jurisdiction (X22)</td>
<td align="left">%</td>
<td align="left">&#x2b;</td>
<td align="left">0.0362</td>
</tr>
<tr>
<td rowspan="5" align="left">Green investment</td>
<td rowspan="3" align="left">Environmental pollutioncontrol investment</td>
<td colspan="2" align="left">Investment in industrial pollution control (Y1)</td>
<td align="left">Million yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0351</td>
</tr>
<tr>
<td colspan="2" align="left">Investment in urban environmental infrastructureconstruction (Y2)</td>
<td align="left">Million yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0337</td>
</tr>
<tr>
<td colspan="2" align="left">&#x201c;Three simultaneous&#x201d; project investment inenvironmental protection engineering (Y3)</td>
<td align="left">Million yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0258</td>
</tr>
<tr>
<td align="left">Water conservancyconstruction investment</td>
<td colspan="2" align="left">Water construction investment (Y4)</td>
<td align="left">Million yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0303</td>
</tr>
<tr>
<td align="left">Forestry investment</td>
<td colspan="2" align="left">Forestry investment (Y5)</td>
<td align="left">Million yuan</td>
<td align="left">&#x2b;</td>
<td align="left">0.0338</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-1-2">
<title>2.1.2 Coupling coordination model</title>
<p>The concept of coupling is taken from physics and it refers to the process in which two or more systems interact and thus influence each other. Many scholars used this concept in the fields of economics, ecology, and sociology. Similarly, in this study, we considered the &#x201c;green investment-green ecology&#x201d; system as the research object. Following <xref ref-type="bibr" rid="B40">Tang (2015)</xref>, a coupling coordination model consisting of a coupling degree model and a coordination degree model was constructed.<disp-formula id="e3">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x2003;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m10">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mfrac bevelled="true">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
</mml:msqrt>
<mml:mo>&#x2003;</mml:mo>
<mml:mi>T</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2202;</mml:mo>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>f</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mi>U</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m11">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf7">
<mml:math id="m12">
<mml:mrow>
<mml:mi mathvariant="normal">C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the coupling degree, <inline-formula id="inf8">
<mml:math id="m13">
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the integrated coordination index of green investment and green ecosystem. <inline-formula id="inf9">
<mml:math id="m14">
<mml:mrow>
<mml:mo>&#x2202;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf10">
<mml:math id="m15">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> are the undetermined coefficients of green investment and green ecosystem, respectively, which are listed at the same level in this paper. Therefore, both are brought in with a value of .5 (<xref ref-type="bibr" rid="B4">Ai et al., 2016</xref>). <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the integrated coordination degree of green investment and green ecosystem.</p>
</sec>
<sec id="s2-1-3">
<title>2.1.3 Kernel density estimation</title>
<p>Taking the coupling coordination value as the research object, the kernel density estimation method is applied to explore the time evolution pattern of the coupling relationship between green investment and green ecology. The kernel density estimation is based on non-parametric tests, which evaluate and analyzes the probability densities of variables, and ultimately investigate the evaluation method of data distribution dynamics. While in the fields of management, geography, and economics, the probability density distribution characteristics of the specified years within the study period of the sample can be used to compare the temporal evolution pattern of the study object using a given function.<disp-formula id="e6">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>n</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mi>h</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mi>h</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>h</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where <inline-formula id="inf12">
<mml:math id="m18">
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi mathvariant="normal">h</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the kernel function with bandwidth. In this study, the default bandwidth of Stata software is used, while the Gaussian kernel function is used. For the implementation of kernel density, we have followed <xref ref-type="bibr" rid="B21">Jiang et al. (2019)</xref>.<disp-formula id="e7">
<mml:math id="m19">
<mml:mrow>
<mml:mi>G</mml:mi>
<mml:mi mathvariant="normal">ausian</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:msup>
<mml:mi>t</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-1-4">
<title>2.1.4 Convergence function</title>
<p>To explore the variability characteristics of a certain research object in the research region, the convergence function (<italic>&#x3c3;</italic>) can measure the variability of a certain index. In this study, the coupling coordination value is brought into the convergence function for a preliminary study on the variability of the coupling coordination between green investment and green ecology in each region of China. Following the studies of <xref ref-type="bibr" rid="B46">Wu et al. (2020)</xref>, the specific equation can be written as:<disp-formula id="e8">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>log</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where <inline-formula id="inf13">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> shows the coupling coordination value of province <inline-formula id="inf14">
<mml:math id="m22">
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> in year <inline-formula id="inf15">
<mml:math id="m23">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf16">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>t</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represents the log standard deviation of <inline-formula id="inf17">
<mml:math id="m25">
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> provinces in year <inline-formula id="inf18">
<mml:math id="m26">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>.</p>
</sec>
<sec id="s2-1-5">
<title>2.1.5 Spatial clustering analysis</title>
<p>The spatial autocorrelation method is a type of exploratory spatial data analysis to explore the interconnections and dependencies of variables among regions within the research object that consisted of two theoretical models i.e., global spatial autocorrelation and local spatial autocorrelation (<xref ref-type="bibr" rid="B30">Li et al., 2020</xref>). In this paper, we used <italic>Moran&#x2019;s I</italic> index and <italic>Getis-Ord Gi&#x2a;</italic> index, which are commonly used in the application of indices to focus on the spatial local clustering characteristics among provinces, autonomous regions, and municipalities that are directly under the central government.</p>
<p>Global spatial autocorrelation (<italic>Global Moran&#x2019;s I</italic>):<disp-formula id="e9">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>g</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Local spatial autocorrelation (<italic>Local Moran&#x2019;s I</italic>):<disp-formula id="e10">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>X</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>
<italic>Getis-Ord Gi&#x2a;</italic> index:<disp-formula id="e11">
<mml:math id="m29">
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi>G</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2211;</mml:mo>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mo>&#x2211;</mml:mo>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <inline-formula id="inf19">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi mathvariant="normal">g</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the global Moran index, <inline-formula id="inf20">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the local Moran index, <inline-formula id="inf21">
<mml:math id="m32">
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is the sample size, <inline-formula id="inf22">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf23">
<mml:math id="m34">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denote the coupling coordination values of the <inline-formula id="inf24">
<mml:math id="m35">
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf25">
<mml:math id="m36">
<mml:mrow>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> provinces and municipalities respectively. <inline-formula id="inf26">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">j</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the spatial weight that indicates the spatial relationship between regions within the research object. In this paper, we used the adjacent spatial weight as the operational value, i.e., adjacent is recorded as 1 and non-adjacent is recorded as 0.<disp-formula id="e12">
<mml:math id="m38">
<mml:mrow>
<mml:mi>Z</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>I</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>E</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:msqrt>
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:mi>A</mml:mi>
<mml:mi>R</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
</p>
<p>If <inline-formula id="inf27">
<mml:math id="m39">
<mml:mrow>
<mml:mi>Z</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and passes the <inline-formula id="inf28">
<mml:math id="m40">
<mml:mrow>
<mml:mi>Z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-value significance test, it indicates that the coupling coordination value is significant and positively correlated, while if <inline-formula id="inf29">
<mml:math id="m41">
<mml:mrow>
<mml:mi>Z</mml:mi>
<mml:mo>&#x3c;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula> and passes the <inline-formula id="inf30">
<mml:math id="m42">
<mml:mrow>
<mml:mi>Z</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-value significance test, it indicates that the coupling coordination value is significant and negatively correlated; otherwise, the correlation is not significant.</p>
</sec>
<sec id="s2-1-6">
<title>2.1.6 Spatial standard deviational ellipse</title>
<p>This paper carries out an empirical study using a spatial standard deviational ellipse to investigate the spatial global characteristics for the coupling and coordination characteristics of green investment and green ecology. The standard deviational ellipse focuses on revealing the global characteristics of geographical elements&#x2019; spatial distribution and belongs to the statistical analysis method of spatial pattern. The model mainly describes the spatial distribution dynamics of geographical elements quantitatively with parameters such as the center of gravity, rotation angle <italic>&#x3b8;</italic>, <italic>x</italic>-axis standard deviation, and <italic>y</italic>-axis standard deviation (<xref ref-type="bibr" rid="B16">Gong, 2002</xref>).<disp-formula id="e13">
<mml:math id="m43">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mo>;</mml:mo>
<mml:mi>Y</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>
<disp-formula id="e14">
<mml:math id="m44">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="italic">X</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="italic">n</mml:mi>
</mml:munderover>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
<mml:mo>;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mfrac>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="italic">X</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:msup>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="italic">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="italic">W</mml:mi>
<mml:mi mathvariant="italic">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
<disp-formula id="e15">
<mml:math id="m45">
<mml:mrow>
<mml:mi mathvariant="italic">tan</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b8;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">Y</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">Y</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>4</mml:mn>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">Y</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mo>&#x2a;</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="normal">n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:msup>
<mml:msub>
<mml:mi mathvariant="normal">Y</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>where <inline-formula id="inf31">
<mml:math id="m46">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>X</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>Y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is the barycentric coordinate of the coupling coordination value of the green investment and green ecology. <inline-formula id="inf32">
<mml:math id="m47">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is the spatial coordinate of the study area, and <inline-formula id="inf33">
<mml:math id="m48">
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msup>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>Y</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2a;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
</inline-formula> is the relative coordinate of each point in the study area from the regional coupling coordination center. <inline-formula id="inf34">
<mml:math id="m49">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the value of coupling coordination between green investment and green ecology in each region. <inline-formula id="inf35">
<mml:math id="m50">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf36">
<mml:math id="m51">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mi mathvariant="normal">y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the standard deviation along the <inline-formula id="inf37">
<mml:math id="m52">
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-axis and <inline-formula id="inf38">
<mml:math id="m53">
<mml:mrow>
<mml:mi>Y</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>-axis, respectively. Similarly, <inline-formula id="inf39">
<mml:math id="m54">
<mml:mrow>
<mml:mi>&#x3b8;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> shows the deflection angle.</p>
</sec>
</sec>
<sec id="s2-2">
<title>2.2 Sources of indicator and classification</title>
<sec id="s2-2-1">
<title>2.2.1 Construction of indicator system and sources</title>
<p>Green investment indicators and green ecological indicators are complex systems, including social, economic, environmental, and other factors (<xref ref-type="bibr" rid="B36">Peng et al, 2018</xref>; <xref ref-type="bibr" rid="B38">Peng et al, 2020</xref>; <xref ref-type="bibr" rid="B50">Zhao et al, 2021</xref>; <xref ref-type="bibr" rid="B51">Zhong et al, 2021</xref>). This paper upholds the principles of indicator availability, comprehensiveness, and scientific nature. Moreover, following <xref ref-type="bibr" rid="B12">Fei et al. (2020)</xref>, <xref ref-type="bibr" rid="B19">He et al. (2017)</xref>, and <xref ref-type="bibr" rid="B33">Lu et al. (2017)</xref> we selected a total of three first-level indicators for green ecosystems such as green growth, green welfare, and resource environment. Similarly, a total of three first-level indicators are selected from green investment such as investment in environmental pollution control, investment in water conservancy construction, and investment in forestry. The detail of the selection of indicators is given in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<p>The data were obtained from the China Statistical Yearbook on Environment, the China Water Statistical Yearbook, the China Forestry Yearbook, and the statistical yearbooks of each province from 2005 to 2019. The vector map of China was used to visualize the conclusions which it is derived by scanning the administrative boundaries of the region, and the individual missing values are filled in by linear regression interpolation.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Coupling feature classification</title>
<p>Coordination degree refers to studying the degree of coupling while considering the coordination degree between systems, which is used to explore the good or bad coupling degree between green investment and green ecosystem, usually, the value of D is in the range of 0 and 1. This paper draws on the research results of <xref ref-type="bibr" rid="B15">Gan and Bu (2020)</xref> to make the following classification of coordination degree level (<xref ref-type="table" rid="T2">Table 2</xref>).</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Classification of coupling coordination level.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Range of D value</th>
<th align="left">Types of classification</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">0.81 &#x3c; D &#x2264; 1.00</td>
<td align="left">Complete coordination stage (V)</td>
</tr>
<tr>
<td align="left">0.61 &#x3c; D &#x2264; 0.80</td>
<td align="left">Good coordination stage (IV)</td>
</tr>
<tr>
<td align="left">0.41 &#x3c; D &#x2264; 0.60</td>
<td align="left">Basic coordination stage (III)</td>
</tr>
<tr>
<td align="left">0.21 &#x3c; D &#x2264; 0.40</td>
<td align="left">Low coordination stage (II)</td>
</tr>
<tr>
<td align="left">0.00 &#x3c; D &#x2264; 0.20</td>
<td align="left">No coordination stage (I)</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Analysis of coupling coordination</title>
<p>The coupling coordination values of green investment and green ecology were estimated based on the coupling coordination model. In particular, we selected 2005, 2010, 2014, and 2019 as time points to study the temporal characteristics of the coupling coordination in various provinces of China (<xref ref-type="fig" rid="F1">Figure 1</xref>). Furthermore, the provinces were divided into eastern, central, and western regions to determine the spatial distribution of the coupling coordination between green investment and green ecology.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Evolution values of coupling coordination.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g001.tif"/>
</fig>
<p>The results depict that from 2005 to 2019, the coupling coordination values of provinces showed a continuously increasing trend, but all of them were at the stage of serious disorder, and the coupling coordination values also vary widely between regions. It shows that the combination of green investment and green ecology in China has become deeper during the study period, but the interaction between the two is relatively chaotic and the development trends are different between regions. For instance, from 2005 to 2010, the coupling coordination values were generally larger in the eastern provinces of China particularly, Guangdong, Liaoning, and Shandong. It indicates that the superior economic advantages and infrastructure in the eastern regions provide a solid foundation for the combination of green investment and green ecology. Similarly, from 2010 to 2019, the coupling coordination values of developed eastern coastal provinces such as Jiangsu, Shanghai, Zhejiang, and Guangdong still ranked in the first echelon, but the mean value of coupling coordination values in the central region is higher than that in the eastern region during the period. It means that the coupling coordination values of eastern provinces particularly, Liaoning and Tianjin are too low and there is a certain degree of polarization.</p>
<p>Regarding types of coordination (<xref ref-type="table" rid="T3">Table 3</xref>), the types for coordinating green investment and green ecology in 2005&#x2013;2019 are only the no coordination type and the low coordination type. Specifically, the type of coordination between green investment and green ecology in China from 2005 to 2010 is dominated by the type of no coordination, and the coordination relationship between green investment and green ecology is generally in a chaotic state during the period. From 2014 to 2019, the number of provinces in the low coordination type gradually increased by 66.7% and 76.7% of the total in 2014 and 2019, respectively. This implies that during this period China&#x2019;s green investment and green ecology started the coordination phenomenon; however, the overall coordination level is still not high.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Evolution of types of coupling coordination.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Provinces</th>
<th align="center">2005</th>
<th align="center">2010</th>
<th align="center">2014</th>
<th align="center">2019</th>
<th align="left">Provinces</th>
<th align="center">2005</th>
<th align="center">2010</th>
<th align="center">2014</th>
<th align="center">2019</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Heilongjiang1</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="left">Guangxi</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Xinjiang</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Guangdong</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Shanxi</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Hainan</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Ningxia</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="left">Jilin</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Shandong</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Liaoning</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Henan</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Tianjin</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Jiangsu</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Qinghai</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Anhui</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Gansu</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
</tr>
<tr>
<td align="left">Hubei</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Shaanxi</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Zhejiang</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Inner Mongolia</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Jiangxi</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Chongqing</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Hunan</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Hebei</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Yunnan</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Shanghai</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Guizhou</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Beijing</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
</tr>
<tr>
<td align="left">Fujian</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
<td align="left">Sichuan</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2160;</td>
<td align="center">&#x2161;</td>
<td align="center">&#x2161;</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>
<inline-graphic xlink:href="FENVS_fenvs-2022-1095394_wc_tfx1.tif"/>Shows no coordination stage (I). <inline-graphic xlink:href="FENVS_fenvs-2022-1095394_wc_tfx2.tif"/>Shows low coordination stage (II).</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3-2">
<title>3.2 Time evolution</title>
<p>We used the kernel density estimation method and determined the evolution characteristics of the coupling coordination between green investment and green ecology. <xref ref-type="fig" rid="F2">Figure 2</xref> illustrates that from 2005 to 2019, there is an overall shift to the right in the center of the kernel density curve of the coupling coordination between green investment and green ecology, and the magnitude of the shift gradually increases. However, the change interval has decreased year by year. It suggests that the coupling coordination value of Chinese provinces has increased continuously over 15&#xa0;years and the regional gap has narrowed.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Kernel density distribution of coupling coordination.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g002.tif"/>
</fig>
<p>From 2005 to 2019, the kernel density curve shows a shift from multiple peaks to a single peak. In 2005, there are two peaks, of which the kernel density value corresponding to the first peak is smaller than the second one, indicating a clear polarization of the coupling coordination between green investment and green ecology in each province. In 2010, the kernel density curve consists of the main peak and two side peaks, suggesting that the polarization characteristics of the coupling coordination of green investment and green ecology have been expanded gradually. The shape of the kernel density curve in 2014 is similar to that in 2005; however, the number of peaks decreased as compared to 2010 until the kernel density curve developed into a single peak in 2019. It means that the regional differences in the coupling coordination of green investment and green ecology in each province decreased gradually. This indicates that a higher level of agglomeration phenomenon and polarization characteristics disappeared.</p>
<p>From 2005 to 2019, the kurtosis of the kernel density curve of green investment and green ecological coupling coordination shifts to the right. It shows a shift from a &#x201c;sharp peak&#x201d; to a &#x201c;broad peak.&#x201d; From 2005 to 2014, the left peak of the kernel density curve gradually increased, and the spike feature of the crest become more and more significant, which means that the coupling coordination value of the green investment and green ecology increased in each province of China. Furthermore, <xref ref-type="fig" rid="F2">Figure 2</xref> illustrates that the kurtosis of the kernel density curve decreased year by year from 2014 to 2019, while the crest pattern gradually flattened out and evolved toward a broad peak feature. This implies that the number of regions with the coupling coordination value of the green investment and green ecology corresponding to each crest was expanded.</p>
</sec>
<sec id="s3-3">
<title>3.3 Convergence analysis</title>
<p>The convergence function (<italic>&#x3c3;</italic>) was used to estimate the differences in the coupling coordination of green investment-green ecology among various regions of China (<xref ref-type="fig" rid="F3">Figure 3</xref>). At the national level, from 2005 to 2014, the <italic>&#x3c3;</italic> values continue to decline, indicating a gradual decrease in the variability of the coupling coordination between green investment and green ecology in China. <xref ref-type="fig" rid="F3">Figure 3</xref> illustrates that the value of <italic>&#x3c3;</italic> increased slightly from 2014 to 2019, but it is still smaller than the <italic>&#x3c3;</italic> values in 2005. It reveals an overall decrease in the variability of coupling coordination over the study period. At the regional level particularly in eastern and western regions, the values of <italic>&#x3c3;</italic> experienced an evolutionary trend of slowly decreasing and then slowly increasing. The overall coupling coordination showed a convergence trend. From 2005 to 2014, the values of <italic>&#x3c3;</italic> in the central region declined rapidly and then increased slowly from 2014 to 2019, which indicates a significant overall decrease in the variability of coupling coordination in the region. The order of <italic>&#x3c3;</italic> values in the three regions is the highest in the eastern region, the second-highest in the western region, and the smallest in the central region, which shows that the central region has the least variability in the coupling coordination characteristics of green investment and green ecology.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Convergence analysis of the coupling coordination.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g003.tif"/>
</fig>
</sec>
<sec id="s3-4">
<title>3.4 Spatial clustering analysis</title>
<p>The coupling coordination value of the green investment and green ecology were taken as the research variable and calculated Moran&#x2019;s I index for 2005, 2010, 2014, and 2019. We performed 9999 randomizations to test the significance level (<xref ref-type="table" rid="T4">Table 4</xref>). In 2005, 2014, and 2019, the <italic>Moran&#x27;I</italic> index was greater than zero, and it showed an evolutionary trend of first decreasing and then increasing. It demonstrates that the coupling coordination of green investment and green ecology has spatial clustering characteristics at this time, and the clustering characteristics were continuously strengthened in 2014 and 2019. In 2010, the <italic>Moran&#x27;I</italic> index is less than zero, which means that the coupling coordination at this time does not have the characteristics of spatial agglomeration, but presents a certain degree of dispersion. Except in 2005, when the <italic>Z</italic>-value of coupling coordination passes the 5% significance test, the rest of the time points are not significant at the 5% or 10% statistical level. However, it is not proved that there is a lack of autocorrelation of coupling coordination between green investment and green ecology in China because the selection of spatial weights and the differences between regions may lead to insignificant coupling coordination values. Therefore, the spatial correlation analysis using the local <italic>Moran&#x27;I</italic> index is required.</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Global <italic>Moran&#x27;I</italic> indexes.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Year</th>
<th align="left">
<italic>Moran&#x27;I</italic>
</th>
<th align="left">
<italic>E(I)</italic>
</th>
<th align="left">
<italic>Z(I)</italic>
</th>
<th align="left">
<italic>P</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">2005</td>
<td align="left">0.1228</td>
<td align="left">&#x2212;0.0345</td>
<td align="left">2.0591</td>
<td align="left">0.0395</td>
</tr>
<tr>
<td align="left">2010</td>
<td align="left">&#x2212;0.0153</td>
<td align="left">&#x2212;0.0345</td>
<td align="left">0.2571</td>
<td align="left">0.7971</td>
</tr>
<tr>
<td align="left">2014</td>
<td align="left">0.0104</td>
<td align="left">&#x2212;0.0345</td>
<td align="left">0.5941</td>
<td align="left">0.5524</td>
</tr>
<tr>
<td align="left">2019</td>
<td align="left">0.0901</td>
<td align="left">&#x2212;0.0345</td>
<td align="left">1.6215</td>
<td align="left">0.1049</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The local Moran index can reasonably analyze the spatial autocorrelation of the study area. The local indicators of spatial association (<italic>LISA</italic>) clustering maps are plotted at a 5% significance level. <xref ref-type="fig" rid="F4">Figure 4</xref> illustrates that the coupling coordination of green investment and green ecology in China from 2005 to 2019 possesses four types of spatial agglomeration: 1) Significantly high-high type, i.e., the coupling coordination values of both itself and the surrounding areas are high, and the differences between regions are small. It reveals the phenomenon of high-value aggregation. In terms of its spatial distribution, it was mainly concentrated in Shandong, Jiangsu, Shanghai, Zhejiang, and Henan in 2005, while the number of high-value aggregation areas gradually disappeared in 2010. The number of significantly high-high type areas continued to increase from 2010 to 2019, extending to the middle and lower reaches of the Yangtze River provinces and their surrounding areas. 2) Significantly high-low type, i.e., the self-coupling coordination value is high, but the surrounding areas are small, and the spatial differences are more obvious. The regional distribution of this spatial agglomeration type is limited and existed only in Sichuan Province in 2005. 3) Significantly low-high type, i.e., the self-coupling coordination value is low, but the coupling coordination values of surrounding areas are high and the spatial differences are more significant. From 2005 to 2014, the spatial distribution is more regionally concentrated, mostly located in Anhui and Jiangxi provinces, and occasionally in Hainan and Shanghai in individual years. 4) Significantly low-low type i.e., the coupling coordination value between itself and the surrounding areas is low, and the spatial difference is small, showing a low-value aggregation. In 2005, a significantly low-low type is distributed only in Xinjiang and Gansu. After 2005, there are few low-value aggregation regions in Gansu Province.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Local indicators of spatial association of the coupling coordination. <bold>(A)</bold> 2005, <bold>(B)</bold> 2010, <bold>(C)</bold> 2014, <bold>(D)</bold> 2019.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g004.tif"/>
</fig>
<p>The <italic>Getis-Ord Gi&#x2a;</italic> index can be calculated for each element of the research data, and the statistical method is different from the local Moran index, which is an important complement to local spatial autocorrelation. Therefore, this paper further explores the clustering characteristics of the coupling coordination of green investment and green ecology in China and divides them into cold spot areas, sub-cold spot areas, sub-hot spot areas, and hot spot areas using the natural break method (<xref ref-type="fig" rid="F5">Figure 5</xref>). The results can be reported as: 1) The number of hot spot areas increased year by year, gradually developing from the coastal provinces in the lower reaches of the Yangtze River to the provinces in the middle and lower reaches of the Yangtze River and its surrounding areas, demonstrating that these provinces have higher coupling coordination values, have advantages in the combination of green investment and green ecology, and presented a high-value aggregation trend. 2) During 2005 and 2014, the number of sub-hot spot areas decreased gradually, initially concentrated in the central provinces of East China and some central provinces, and gradually shrinking to the coastal areas of East China. In 2014 and 2019, the number of sub-hot spot areas grew rapidly to 18 and was distributed in the eastern, central, and western regions. There are large differences in coupling coordination values between regions. 3) The development trend of the number of sub-cold spot areas is opposite to that of sub-hot spot areas, experiencing a development trend of gradual expansion followed by rapid shrinkage, and in 2019 there is Heilongjiang Province belonged to sub-cold spot areas. The regional differences are large. 4) The number of cold spot areas gradually decreased and is concentrated in the western region. After 2014, only Qinghai and Gansu existed. It shows that the low-value aggregation area of the coupling coordination value is gradually shrinking in the western region, and the regional differences are small.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Cold and hot spot areas of the coupling coordination degree. <bold>(A)</bold> 2010, <bold>(B)</bold> 2005, <bold>(C)</bold> 2014, <bold>(D)</bold> 2019.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g005.tif"/>
</fig>
</sec>
<sec id="s3-5">
<title>3.5 Analysis of spatial evolution</title>
<p>After proving the existence of spatial clustering characteristics of the coupling coordination of green investment and green ecology in China, this study further applies the center of gravity-standard deviation ellipse method to investigate the spatial evolutionary dynamics of the coupling coordination values. And, 2005, 2010, 2014, and 2019 are selected to study the spatial distribution characteristics of the coupling coordination of green investment and green ecology (<xref ref-type="table" rid="T5">Table 5</xref>; <xref ref-type="fig" rid="F6">Figure 6</xref>).</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Shift of the center of gravity for coupling coordination values and standard deviation ellipse parameters.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">Year</th>
<th colspan="2" align="left">Barycentric coordinates</th>
<th colspan="3" align="left">Center of gravity</th>
<th colspan="2" align="left">Standard deviation (km)</th>
<th rowspan="2" align="left">Angle of rotation</th>
<th align="left">Oblateness</th>
</tr>
<tr>
<th align="left">Longitude (<sup>&#xba;</sup>E)</th>
<th align="left">Latitude (<sup>&#xba;</sup>N)</th>
<th align="left">Direction</th>
<th align="left">Distance (<italic>km</italic>)</th>
<th align="left">Velocity (<italic>km/a</italic>)</th>
<th align="left">Minor axis</th>
<th align="left">Major axis</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">2005</td>
<td align="left">112.59</td>
<td align="left">33.49</td>
<td align="left">&#x2014;</td>
<td align="left">&#x2014;</td>
<td align="left">&#x2014;</td>
<td align="left">964.93</td>
<td align="left">1,098.26</td>
<td align="left">28.99</td>
<td align="left">1.138</td>
</tr>
<tr>
<td align="left">2010</td>
<td align="left">112.34</td>
<td align="left">33.43</td>
<td align="left">Southwest</td>
<td align="left">40.02</td>
<td align="left">8.04</td>
<td align="left">959.13</td>
<td align="left">1,130.56</td>
<td align="left">28.19</td>
<td align="left">1.179</td>
</tr>
<tr>
<td align="left">2014</td>
<td align="left">112.29</td>
<td align="left">33.52</td>
<td align="left">Northwest</td>
<td align="left">18.18</td>
<td align="left">4.55</td>
<td align="left">987.36</td>
<td align="left">1,119.94</td>
<td align="left">31.01</td>
<td align="left">1.134</td>
</tr>
<tr>
<td align="left">2019</td>
<td align="left">112.27</td>
<td align="left">33.38</td>
<td align="left">Southwest</td>
<td align="left">26.10</td>
<td align="left">5.22</td>
<td align="left">981.72</td>
<td align="left">1,108.10</td>
<td align="left">28.35</td>
<td align="left">1.129</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Spatial distribution of evolutionary trend of the coupling coordination degree.</p>
</caption>
<graphic xlink:href="fenvs-10-1095394-g006.tif"/>
</fig>
<p>Examining the distribution, distance, and rate of movement of the center of gravity in the center of gravity model, the centers of gravity of the coupling coordination between green investment and green ecology in China are all distributed in Henan Province. It indicates that in the generalized east-west direction, the combination of green investment and green ecology in eastern China is generally better than that in the western region. It is found that the center of gravity of the coupling coordination value in 2005 is located in Pingdingshan City, Henan Province. From 2005 to 2010, the center of gravity of the coupling coordination value shifted to the southwest from Baofeng County (Pingdingshan City) to Lushan County (Pingdingshan City), moving 40.02&#xa0;km and 8.04&#xa0;km/a, respectively. It is the farthest distance and the fastest moving rate. Similarly, the center of gravity shifted to the northwest from 2010 to 2014, moving from Baofeng County to Song County in Luoyang City, but the distance and rate of movement decreased to 18.18&#xa0;km and 4.545&#xa0;km/a, respectively. Moreover, the center of gravity of the coupling coordination value shifted to the southwest again in 2014 and 2019 and approached the edge of Nanyang City, at this time the moving distance increased slightly and the moving rate increased slowly. The above analysis shows that the center of gravity of China&#x2019;s green investment and green ecological coupling coordination generally showed a westward shift in the east-west direction, and shifts northward and then southward in the north-south direction. It is gradually developed towards the middle reaches of the Yangtze River at the geographical level. The moving speed of the center of gravity generally showed a trend of &#x201c;accelerating to gradually decreasing to slowly increasing.&#x201d;</p>
<p>From the evolution of the rotation angle, the spatial pattern of the coupling coordination between green investment and green ecology in China from 2005 to 2019 showed a &#x201c;northeast-southwest&#x201d; phenomenon and a shift toward the &#x201c;north-south&#x201d; direction, while the rotation angle <italic>&#x3b8;</italic> presents an overall trend of fluctuating declined. This indicated that the &#x201c;northeast-southwest&#x201d; spatial distribution pattern has a certain degree of tendency to weaken. Specifically, the rotation angle <italic>&#x3b8;</italic> changed by 0.8&#xb0; from 2005 to 2010, which showed that the spatial pattern of coupling coordination between green investment and green ecology has shifted from &#x201c;northeast-southwest&#x201d; to &#x201c;due north-south&#x201d; by 0.8&#xb0;. From 2010 to 2014, the rotation angle <italic>&#x3b8;</italic> increased to 31.01&#xb0;, and the spatial distribution of the coupling coordination type shifted to the &#x201c;due north-south&#x201d; direction again by 2.66&#xb0; from 2014 to 2019.</p>
<p>From the analysis of the length of the main axis, the standard deviation of both the major axis and the minor axis of the standard deviation ellipse from 2005 to 2019 showed a fluctuation with an increasing trend. It indicated that the coupling coordination between green investment and green ecology in China is gradually dispersed in the &#x201c;northeast-southwest&#x201d; and &#x201c;southeast-northwest&#x201d; directions. From the change of the major axis, the standard deviation of the major axis in 2005 gradually expanded from 1,098.26 to 1,130.56&#xa0;km in 2010, suggesting that the coupling coordination is dispersed in the spatial direction of &#x201c;northeast-southwest&#x201d;. The standard deviation of the major axis decreased from 1,130.56 to 1,108.10&#xa0;km from 2010 to 2019. It showed that the coupling coordination gradually polarized in the &#x201c;northeast-southwest&#x201d; direction during this period. From the change of the minor axis, the standard deviation of the minor axis decreased from 2005 to 2010 and from 2014 to 2019, which means that the coupling coordination is polarized in the &#x201c;southeastward-northwestward&#x201d; direction. The standard deviation gradually increased from 2010 to 2014. It indicated that the coupling coordination is dispersed in the direction of &#x201c;southeastward - northwestward.&#x201d; The change of the flatness rate fluctuates from 1.138 to 1.129 during the study period. This implies that the coupling coordination between green investment and green ecology in China showed a local aggregation trend, and the overall aggregation degree gradually decreased.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>In this study, the coupling coordination characteristics of green investment and green ecology are studied in spatial and temporal dimensions. In particular, the coupling coordination is explored in three aspects: temporal evolution characteristics, spatial clustering analysis, and spatial pattern evolution.</p>
<p>Regarding the temporal evolution of coupling coordination, the national coupling coordination value developed from 0.1432 to 0.2314 from 2005 to 2019, and the coupling coordination value of the green investment and green ecology gradually increased. However, it is still at the stage of serious unbalanced development. Similar to the study of <xref ref-type="bibr" rid="B15">Gan and Bu (2020)</xref>, who concluded that the coupling coordination of the green financial system and ecology was changed from the uncoordinated stage to the well-coordinated stage. The center of the kernel density curve shifted to the right, but the changed interval decreased year by year. At the same time, the number of peaks gradually decreased, showing a &#x201c;spike-broad peak&#x201d; transformation, indicating that the coupling coordination values of Chinese provinces have narrowed during the study period. Therefore, the clustering phenomenon is more significant, and the polarization characteristics have disappeared significantly. The results of the convergence analysis revealed that the differences in coupling coordination decreased across Chinese provinces, but the most significant difference was found in the eastern region. Probably because China first proposed the concept of ecological civilization in 2007, thus the coupling coordination value of green investment and green ecology gradually increased. Until 2016 when China announced the industrial green development and the green financial system planning, the green investment was formally involved in the green ecological construction, while the coupling coordination was still in the primary stage. However, green development policies vary from region to region, and regional characteristics were formed with the evolution over time. Thus the gap in regional coupling coordination was decreased. Among them, the economic base and geopolitical advantages of the eastern region are better than those of the central and western regions. The polarization phenomenon among the provinces in the eastern region is also the most significant, and therefore there is a greater variability of coupling coordination is found.</p>
<p>Regarding the analysis of spatial cluster, whether it is the local <italic>Moran&#x27;I</italic> index or the <italic>Getis-Ord Gi&#x2a;</italic> index, the spatial distribution characteristics showed that the high-value agglomeration area developed gradually from the Yangtze River inlet and its vicinity to the middle reaches of the Yangtze River, making it a highland for the synergistic development of green investment and green ecology, and the radiation effect on the neighboring provinces has always existed, which is similar to the conclusion of other studies, for instance, <xref ref-type="bibr" rid="B47">Yao et al. (2014)</xref>. In these studies, the significant economic and ecological agglomeration characteristics in the eastern part of China were found because China proposed the strategy of Yangtze River Golden Waterway, the Yangtze River estuary area represented by Shanghai, Jiangsu, Anhui, and Zhejiang which has a strong economic advantage. The amount of related environmental pollution control investment, water conservancy construction investment, and forestry investment is also more than sufficient. The geographical situation is mostly plain, coupled with a variety of policy resources, which can ensure that green investment is efficiently and sufficiently invested in the field of green ecology. In China, the combination of green ecology and green ecology is in the first echelon. However, low-value aggregation regions are only distributed in and around Gansu with the least obvious coupling coordination characteristics. Perhaps because the economic base of the western region represented by Gansu is relatively weak, and the ecological environment is relatively poor, the combination of the two has a long way to go.</p>
<p>Regarding spatial pattern evolution, the center of gravity of the coupling coordination of green investment and green ecology generally shifted 53.84&#xa0;km to the southwest. The results are in line with the studies of <xref ref-type="bibr" rid="B27">Lai et al. (2020)</xref>, who found the coupling coordination of ecology and economy in the central and western provinces gradually increased. Probably because the implementation of the Belt and Road and the Western Development Strategy proposed by China, led to the strengthening of the state&#x2019;s support for the western and southwestern regions, providing an important policy backing for the economic development of the central and western regions and laying the foundation for the implementation of green investments. This also confirmed the conclusion that the coupled and coordinated high-value agglomeration region evolved towards the middle reaches of the Yangtze River. The deflection angle of the standard deviation ellipse fluctuated from 28.99&#xb0; to 28.35&#xb0;, with an overall shift in the &#x201c;due north-south&#x201d; direction, similar to the conclusion that the center of gravity shifted southward in the north-south direction, while the principal axis length of the standard deviation ellipse increased during the study period and the flattening of the ellipsoid decreased from 1.138 to 1.129, indicating that the distribution of coupling coordination expanded during 2005 and 2019, and the degree of synergy between green investment and green ecology increased, but its fragmentation was inevitable. This may be due to the gradual improvement of the green financial system, and the level of green investment is also tilted by policy resources. The provinces can put the amount of green investment into relevant fields, which makes the degree of synergistic development of green investment and green ecology gradually rise. However, the convergent policy requirements also restrict the development of each province, making it difficult to form regional characteristics, and thus creating a fragmented situation.</p>
<p>The findings of the study show a more significant coupling effect between green investment and green ecology in China, at the level of temporal evolution analysis, spatial clustering analysis, and spatial pattern evolution analysis.</p>
</sec>
<sec id="s5">
<title>5 Conclusion and policy implications</title>
<p>The study focuses on the coupling coordination of green investment and green ecology in 30 Chinese provinces and cities from spatial and temporal dimensions using the coupling coordination model, kernel density estimation, spatial autocorrelation, and spatial standard deviation ellipse methods. The main findings of the study can be summarized as:</p>
<p>Overall, the study found that the coupling effect of green investment and green ecology in China is generally at a preliminary stage. In particular, the coupling coordination value of the green investment and green ecology is gradually increased, but the overall state is in a serious imbalance. Meanwhile, the regional gap of coupling coordination narrowed with the most significant variability in the eastern region. Specifically, the coupling coordination value of the green investment and green ecology in China developed from 0.143 to 0.231, but the coordination type is dominated by low-level coordination, while among all regions the difference coefficient in the eastern region is the highest.</p>
<p>The aggregation effect of the coupling coordination of green investment and green ecology is significant. The high-value aggregation area extends from the lower reaches of the Yangtze River to the middle reaches of the Yangtze River, while the western region is dominated by low-value aggregation. The specific performance of the significant high-high type and hot spot areas are from the Yangtze River estuary to the middle reaches of the Yangtze River transfer trend, but the western region such as Gansu has shown low-value aggregation characteristics. Therefore, it is suggested to highlight the dominant position of the Yangtze River Economic Belt, strengthen cooperation with neighboring regions, share green investment patterns and green ecological governance experience, and achieve the goal of coordinated development of green investment and green ecology in internationalized city clusters.</p>
<p>The degree of synergistic development between green investment and green ecology is increased, but fragmentation is inevitable. At the same time, the center of gravity of coupling coordination moved to the southwest, and the spatial pattern gradually weakened in the direction of &#x201c;northeast-southwest.&#x201d; Moreover, the distribution range of the standard deviation ellipse is expanded and the flattened ellipsoid is decreased. The center of gravity of the coupling coordination value is shifted to the southwest in Henan, and the deflection angle is also fluctuated down.</p>
<p>It is suggested that China should continue to improve the green financial system, refine the investment scope of green finance, provide important policy support for the in-depth participation of green investment in the fields of green welfare, green economy, and natural resources, and actively guide the flow of social capital to green ecology. Therefore, it is recommended that Chinese provinces give full play to their geographical advantages, actively carry out the Northeast Revitalization Strategy and Belt and Road Development Strategy, and develop green ecological development strategies according to local conditions.</p>
<p>This paper finds a coupling effect between green investment and green ecology through a quantitative study of coupling characteristics, indicating that there is a degree of positive effect between green investment and green ecology in China. On the other hand, this paper finds that the center of gravity of the combination of green investment and green ecology is close to the middle reaches of the Yangtze River, which shows the excellence of the Yangtze River Economic Belt strategy and provides a regional model for realizing the synergistic development of green investment and green ecology. Due to the complexity of the facilitation between green investment and green ecosystem, this study conducted a macro empirical analysis on the coupling coordination between green investment and green ecology at the provincial level, and ignored at the city level. Moreover, the study did not find the influencing factors of coupling and coordination. The above deficiencies and weaknesses can be overcome in future studies to make more valuable research on green investment and green ecology.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<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 sec-type="disclaimer" id="s9">
<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>Ahmed</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Shahab</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ullah</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Investor sentiment and insurers&#x27; financial stability: Do sovereign rating matter? Geneva papers on risk and insurance - issues and practice</article-title>. <source>Geneva Pap. Risk insur. Issues Pract.</source> <volume>45</volume>, <fpage>281</fpage>&#x2013;<lpage>312</lpage>. <pub-id pub-id-type="doi">10.1057/s41288-020-00160-z</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahmed</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bhutta</surname>
<given-names>U.</given-names>
</name>
</person-group> (<year>2021a</year>). <article-title>Insurers&#x27; risk management as a business process: A prospective competitive advantage or not?</article-title> <source>Eur. J. Manag. Bus. Econ.</source> <volume>31</volume>, <fpage>345</fpage>&#x2013;<lpage>366</lpage>. <comment>ahead-of-print</comment>. <pub-id pub-id-type="doi">10.1108/EJMBE-08-2021-0221</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ahmed</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Issam</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2021b</year>). <article-title>Investor confidence and life insurance demand: Can economic condition limit life insurance business?</article-title> <source>Int. J. Emerg. Mark.</source> <volume>6</volume>, <fpage>0650</fpage>. <comment>ahead-of-print(ahead-of-print)</comment>. <pub-id pub-id-type="doi">10.1108/ijoem-06-2020-0650</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ai</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Exploring coupling coordination between urbanization and ecosystem quality (1985&#x2013;2010): A case study from lianyungang city, China</article-title>. <source>Front. Earth Sci.</source> <volume>10</volume> (<issue>3</issue>), <fpage>527</fpage>&#x2013;<lpage>545</lpage>. <pub-id pub-id-type="doi">10.1007/s11707-015-0531-6</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aldieri</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Makkonen</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Vinci</surname>
<given-names>C. P.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Environmental knowledge spillovers and productivity: A patent analysis for large international firms in the energy, water and land resources fields</article-title>. <source>Resour. Policy</source> <volume>69</volume>, <fpage>101877</fpage>. <pub-id pub-id-type="doi">10.1016/j.resourpol.2020.101877</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Multi-objective optimization of injection-molded plastic parts using entropy weight, random forest, and genetic algorithm methods</article-title>. <source>J J. Polym. Eng.</source> <volume>40</volume> (<issue>4</issue>), <fpage>360</fpage>&#x2013;<lpage>371</lpage>. <pub-id pub-id-type="doi">10.1515/polyeng-2019-0326</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chesney</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Taschini</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Experimental comparison between markets on dynamic permit trading and investment in irreversible abatement with and without non-regulated companies</article-title>. <source>SSRN J.</source> <volume>46</volume>, <fpage>1719026</fpage>. <pub-id pub-id-type="doi">10.2139/ssrn.1719026</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ekeh</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Wara</surname>
<given-names>S. T.</given-names>
</name>
<name>
<surname>Orovwode</surname>
<given-names>H. E.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Management of existing capacity of electric power with energy saving devices</article-title>. <source>Adv. Mat. Res.</source> <volume>19</volume>, <fpage>117</fpage>&#x2013;<lpage>124</lpage>. <pub-id pub-id-type="doi">10.4028/www.scientific.net/amr.18-19.117scientific.net/AMR.18-19.117</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Khalid</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2022a</year>). <article-title>Understanding farmers&#x2019; intention and willingness to install renewable energy technology: A solution to reduce the environmental emissions of agriculture</article-title>. <source>Appl. Energy</source> <volume>309</volume>, <fpage>118459</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2021.118459</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Weijun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Jha</surname>
<given-names>S. K.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Estimation of realistic renewable and non-renewable energy use targets for livestock production systems utilising an artificial neural network method: A step towards livestock sustainability</article-title>. <source>Energy</source> <volume>183</volume>, <fpage>191</fpage>&#x2013;<lpage>204</lpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2019.06.084</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Zhixin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Khalid</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022b</year>). <article-title>Application of an artificial neural network to optimise energy inputs: An energy-and cost-saving strategy for commercial poultry farms</article-title>. <source>Energy</source> <volume>244</volume>, <fpage>123169</fpage>. <pub-id pub-id-type="doi">10.1016/j.energy.2022.123169</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fei</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Can technology R&#x26;D continuously improve green development level in the open economy? Empirical evidence from China&#x2019;s industrial sector</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>27</volume> (<issue>27</issue>), <fpage>34052</fpage>&#x2013;<lpage>34066</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-020-09357-0</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feiock</surname>
<given-names>R. C.</given-names>
</name>
<name>
<surname>Stream</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Environmental protection versus economic development: A false trade-off?</article-title> <source>Public Adm. Rev.</source> <volume>61</volume> (<issue>3</issue>), <fpage>313</fpage>&#x2013;<lpage>321</lpage>. <pub-id pub-id-type="doi">10.1111/0033-3352.00032</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fuyong</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xiang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Regional environmental pressures and corporate green innovation: An examination of low carbon city pilot policies</article-title>. <source>Commun. Finance Account.</source> <volume>21</volume>, <fpage>49</fpage>&#x2013;<lpage>53&#x2b;59</lpage>. <pub-id pub-id-type="doi">10.16144/j.cnki.issn1002-8072.2022.21.027</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Bu</surname>
<given-names>Y. J. E. S. W. C.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A coordinated assessment between green finance and ecological environment along the Yangtze River basin</article-title>. <source>E3S Web Conf.</source> <volume>194</volume>, <fpage>05038</fpage>. <pub-id pub-id-type="doi">10.1051/e3sconf/202019405038</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gong</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Clarifying the standard deviational ellipse</article-title>. <source>Geogr. Anal.</source> <volume>34</volume> (<issue>2</issue>), <fpage>155</fpage>&#x2013;<lpage>167</lpage>. <pub-id pub-id-type="doi">10.1111/j.1538-4632.2002.tb01082.x</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Jha</surname>
<given-names>S. K.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Human health damages related to air pollution in China</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>26</volume> (<issue>13</issue>), <fpage>13115</fpage>&#x2013;<lpage>13125</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-019-04708-y</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Air pollution risks human mental health: An implication of two-stages least squares estimation of interaction effects</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>27</volume> (<issue>2</issue>), <fpage>2036</fpage>&#x2013;<lpage>2043</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-019-06612-x</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Examining the relationship between urbanization and the eco-environment using a coupling analysis: Case study of Shanghai, China</article-title>. <source>Ecol. Indic.</source> <volume>77</volume>, <fpage>185</fpage>&#x2013;<lpage>193</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2017.01.017</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jaraite</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kazukauskas</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lundgren</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The effects of climate policy on environmental expenditure and investment: Evidence from Sweden</article-title>. <source>J. Environ. Econ. Policy</source> <volume>3</volume> (<issue>2</issue>), <fpage>148</fpage>&#x2013;<lpage>166</lpage>. <pub-id pub-id-type="doi">10.1080/21606544.2013.875948</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A novel probabilistic wind speed prediction approach using real time refined variational model decomposition and conditional kernel density estimation</article-title>. <source>Energy Convers. Manag.</source> <volume>185</volume>, <fpage>758</fpage>&#x2013;<lpage>773</lpage>. <pub-id pub-id-type="doi">10.1016/j.enconman.2019.02.028</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khan</surname>
<given-names>M. A.</given-names>
</name>
<name>
<surname>Ozturk</surname>
<given-names>I.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Examining foreign direct investment and environmental pollution linkage in Asia</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>27</volume> (<issue>7</issue>), <fpage>7244</fpage>&#x2013;<lpage>7255</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-019-07387-x</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Umar</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kirikkaleli</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Jiao</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Consumption-based carbon emissions and International trade in G7 countries: The role of Environmental innovation and Renewable energy</article-title>. <source>Sci. Total Environ.</source> <volume>730</volume>, <fpage>138945</fpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2020.138945</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khelfaoui</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hafeez</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ahmed</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Degha</surname>
<given-names>H. E.</given-names>
</name>
<name>
<surname>Meskher</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022a</year>). <article-title>Information communication technology and infant mortality in low-income countries: Empirical study using panel data models</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>19</volume> (<issue>12</issue>), <fpage>7338</fpage>. <pub-id pub-id-type="doi">10.3390/ijerph19127338</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khelfaoui</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hafeez</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ahmed</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Houssem Eddine</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Meskher</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2022b</year>). <article-title>Effects of health shocks, insurance, and education on Income: Fresh analysis using CHNS panel data</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>19</volume>, <fpage>8298</fpage>. <pub-id pub-id-type="doi">10.3390/ijerph19148298</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kline</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Planning and Creating Eco-cities: Indicators as a tool for shaping development and measuring progress</article-title>. <source>Local Environ.</source> <volume>5</volume> (<issue>3</issue>), <fpage>343</fpage>&#x2013;<lpage>350</lpage>. <pub-id pub-id-type="doi">10.1080/13549830050134275</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lai</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yue</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Coupling coordination between environment, economy and tourism: A case study of China</article-title>. <source>PLoS One</source> <volume>15</volume> (<issue>2</issue>), <fpage>e0228426</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0228426</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Landrigan</surname>
<given-names>P. J.</given-names>
</name>
<name>
<surname>Fuller</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Acosta</surname>
<given-names>N. J. R.</given-names>
</name>
<name>
<surname>Adeyi</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Arnold</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Basu</surname>
<given-names>N. N.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>The Lancet Commission on pollution and health</article-title>. <source>Lancet</source> <volume>391</volume> (<issue>10119</issue>), <fpage>462</fpage>&#x2013;<lpage>512</lpage>. <pub-id pub-id-type="doi">10.1016/s0140-6736(17)32345-0</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Leiter</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Parolini</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Winner</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Environmental regulation and investment: Evidence from European industry data</article-title>. <source>Ecol. Econ.</source> <volume>70</volume> (<issue>4</issue>), <fpage>759</fpage>&#x2013;<lpage>770</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecolecon.2010.11.013</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Spatial characteristics of CO2 emissions and PM2.5 concentrations in China based on gridded data</article-title>. <source>Appl. Energy</source> <volume>266</volume>, <fpage>114852</fpage>. <pub-id pub-id-type="doi">10.1016/j.apenergy.2020.114852</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Public appeal, environmental regulation and green investment: Evidence from China</article-title>. <source>Energy Policy</source> <volume>119</volume>, <fpage>554</fpage>&#x2013;<lpage>562</lpage>. <pub-id pub-id-type="doi">10.1016/j.enpol.2018.05.020</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Liyuan</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>The differentiation mechanism of coupling degree among space exploitation, economy development and environment evolution in harbin-dalian giant urban Belt</article-title>. <source>Sci. Geogr. Sin.</source> <volume>38</volume> (<issue>5</issue>), <fpage>662</fpage>&#x2013;<lpage>671</lpage>. <pub-id pub-id-type="doi">10.13249/j.cnki.sgs.2018.05.003</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>An</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Degree of coupling and coordination of eco-economic system and the influencing factors: A case study in yanchi county, ningxia hui autonomous region, China</article-title>. <source>J. Arid. Land</source> <volume>9</volume> (<issue>3</issue>), <fpage>446</fpage>&#x2013;<lpage>457</lpage>. <pub-id pub-id-type="doi">10.1007/s40333-017-0098-z</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mesjasz-Lech</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Environmental protection expenditures and effects of environmental governance of sustainable development in manufacture enterprise</article-title>. <source>Environmental Economics and Managemen</source> <volume>05</volume>, <fpage>96647</fpage>. <pub-id pub-id-type="doi">10.3389/fenvs.2022.96647</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021a</year>). <article-title>Can corporate environmental responsibility improve environmental performance? An inter-temporal analysis of Chinese chemical companies</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>28</volume>, <fpage>12190</fpage>&#x2013;<lpage>12201</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-020-11636-9</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Tu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Extended Producer Responsibility and corporate performance: Effects of environmental regulation and environmental strategy</article-title>. <source>J. Environ. Manag.</source> <volume>218</volume>, <fpage>181</fpage>&#x2013;<lpage>189</lpage>. <pub-id pub-id-type="doi">10.1016/j.jenvman.2018.04.068</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021b</year>). <article-title>Evolution of spatial&#x2013;temporal characteristics and financial development as an influencing factor of green ecology</article-title>. <source>Environ. Dev. Sustain.</source> <volume>24</volume>, <fpage>789</fpage>&#x2013;<lpage>809</lpage>. <pub-id pub-id-type="doi">10.1007/s10668-021-01469-y</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The cultivation mechanism of green technology innovation in manufacturing industry: From the perspective of ecological niche</article-title>. <source>J. Clean. Prod.</source> <volume>252</volume>, <fpage>119711</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2019.119711</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>Z.-W.</given-names>
</name>
<name>
<surname>Malik</surname>
<given-names>M. Y.</given-names>
</name>
<name>
<surname>Umar</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Khan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Khan</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Does green investment, financial development and natural resources rent limit carbon emissions? A provincial panel analysis of China</article-title>. <source>Sci. Total Environ.</source> <volume>755</volume>, <fpage>142538</fpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2020.142538</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>An integrated approach to evaluating the coupling coordination between tourism and the environment</article-title>. <source>Tour. Manag.</source> <volume>46</volume>, <fpage>11</fpage>&#x2013;<lpage>19</lpage>. <pub-id pub-id-type="doi">10.1016/j.tourman.2014.06.001</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Testa</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Iraldo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Frey</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>The effect of environmental regulation on firms&#x27; competitive performance: The case of the building &#x26; construction sector in some EU regions</article-title>. <source>J. Environ. Manage.</source> <volume>92</volume> (<issue>9</issue>), <fpage>2136</fpage>&#x2013;<lpage>2144</lpage>. <pub-id pub-id-type="doi">10.1016/j.jenvman.2011.03.039</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vogel</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Environmental regulation and economic integration</article-title>. <source>J. Int. Econ. LawJournal of International Economic Law</source> <volume>3</volume> (<issue>2</issue>), <fpage>265</fpage>&#x2013;<lpage>279</lpage>. <pub-id pub-id-type="doi">10.1093/jiel/3.2.265</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>C. W.</given-names>
</name>
<name>
<surname>Ali</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>H. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>How China is fostering sustainable growth: The interplay of green investment and production-based emission</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>27</volume> (<issue>31</issue>), <fpage>39607</fpage>&#x2013;<lpage>39618</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-020-09933-4</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Comprehensive evaluation and spatial difference analysis of regional ecological carrying capacity: A case study of the Yangtze River urban agglomeration</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>16</volume> (<issue>18</issue>), <fpage>3499</fpage>. <pub-id pub-id-type="doi">10.3390/ijerph16183499</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhi</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The role of green finance in environmental protection: Two aspects of market mechanism and policies</article-title>. <source>Energy Procedia</source> <volume>104</volume>, <fpage>311</fpage>&#x2013;<lpage>316</lpage>. <pub-id pub-id-type="doi">10.1016/j.egypro.2016.12.053</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>w.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ting</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The distortion degree of relative price of factors and the convergence of regional differences</article-title>. <source>Soc. Sci. Res.</source> <volume>2</volume>, <fpage>75</fpage>&#x2013;<lpage>83</lpage>.</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Jiankang</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Spatial pattern and regional disparity of regional economic integration and ecology-environmental stress in China</article-title>. <source>World Reg. Stud.</source> <volume>23</volume> (<issue>1</issue>), <fpage>111</fpage>&#x2013;<lpage>121</lpage>.</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Ding</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Will green finance promote green development: Based on the threshold effect of R&#x26;D investment</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>29</volume> (<issue>40</issue>), <fpage>60232</fpage>&#x2013;<lpage>60243</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-022-20161-w</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Optimization of Chinese coal-fired power plants for cleaner production using Bayesian network</article-title>. <source>J. Clean. Prod.</source> <volume>273</volume>, <fpage>122837</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2020.122837</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Does the extended producer responsibility system promote the green technological innovation of enterprises? An empirical study based on the difference-in-differences model</article-title>. <source>J. Clean. Prod.</source> <volume>319</volume>, <fpage>128631</fpage>. <pub-id pub-id-type="doi">10.1016/j.jclepro.2021.128631</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhong</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Elahi</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Spatial and temporal pattern evolution and influencing factors of energy&#x2013;environmental efficiency: A case study of Yangtze River urban agglomeration in China</article-title>. <source>Energy &#x26; Environ.</source> <volume>32</volume> (<issue>2</issue>), <fpage>242</fpage>&#x2013;<lpage>261</lpage>. <pub-id pub-id-type="doi">10.1177/0958305x20923114</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Going green or going away: Environmental regulation, economic geography and firms&#x2019; strategies in China&#x2019;s pollution-intensive industries</article-title>. <source>Geoforum</source> <volume>55</volume>, <fpage>53</fpage>&#x2013;<lpage>65</lpage>. <pub-id pub-id-type="doi">10.1016/j.geoforum.2014.05.004</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>