<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" article-type="research-article" dtd-version="1.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Sustain. Food Syst.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Sustainable Food Systems</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Sustain. Food Syst.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2571-581X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fsufs.2026.1765501</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>New energy development, land system transformation, and grain productivity quality: evidence from China</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Yaqian</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ma</surname>
<given-names>Shenglin</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3050663"/>
<role>reviewer</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhang</surname>
<given-names>Xilin</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/3312254"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="supervision" vocab-term-identifier="https://credit.niso.org/contributor-roles/supervision/">Supervision</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Project administration" vocab-term-identifier="https://credit.niso.org/contributor-roles/project-administration/">Project administration</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Chen</surname>
<given-names>Jiaxian</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Funding acquisition" vocab-term-identifier="https://credit.niso.org/contributor-roles/funding-acquisition/">Funding acquisition</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Project administration" vocab-term-identifier="https://credit.niso.org/contributor-roles/project-administration/">Project administration</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &#x0026; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x0026; editing</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="supervision" vocab-term-identifier="https://credit.niso.org/contributor-roles/supervision/">Supervision</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="resources" vocab-term-identifier="https://credit.niso.org/contributor-roles/resources/">Resources</role>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label><institution>Zhejiang University of Technology</institution>, <city>Hangzhou</city>, <country country="cn">China</country></aff>
<aff id="aff2"><label>2</label><institution>North University of China</institution>, <city>Taiyuan</city>, <country country="cn">China</country></aff>
<aff id="aff3"><label>3</label><institution>University of Gda&#x0144;sk</institution>, <city>Gda&#x0144;sk</city>, <country country="pl">Poland</country></aff>
<aff id="aff4"><label>4</label><institution>South China Agricultural University</institution>, <city>Guangzhou</city>, <country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>&#x002A;</label>Correspondence: Xilin Zhang, <email xlink:href="mailto:18905442878@163.com">18905442878@163.com</email>; Jiaxian Chen, <email xlink:href="mailto:349701061ng@gmail.com">349701061ng@gmail.com</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2026-02-19">
<day>19</day>
<month>02</month>
<year>2026</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2026</year>
</pub-date>
<volume>10</volume>
<elocation-id>1765501</elocation-id>
<history>
<date date-type="received">
<day>11</day>
<month>12</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>28</day>
<month>01</month>
<year>2026</year>
</date>
<date date-type="accepted">
<day>29</day>
<month>01</month>
<year>2026</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2026 Liu, Ma, Zhang and Chen.</copyright-statement>
<copyright-year>2026</copyright-year>
<copyright-holder>Liu, Ma, Zhang and Chen</copyright-holder>
<license>
<ali:license_ref start_date="2026-02-19">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. 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.</license-p>
</license>
</permissions>
<abstract>
<p>Amid global climate change and China&#x2019;s dual&#x2011;carbon goals, the rapid development of new energy is reshaping the resource and technological foundations of grain production. Regional differences in land&#x2011;institution transformation and governance capacity lead to varied effects of energy transition on agriculture.Using panel data from 30 Chinese provinces (2013&#x2013;2022), this study constructs a multidimensional Grain Productivity Quality Index (GPQI) via the entropy&#x2011;weight method, covering technological efficiency, digitalization, green&#x2011;input adoption, production integration, and inclusiveness. We employ a two&#x2011;way fixed&#x2011;effects model and a multi&#x2011;period difference&#x2011;in&#x2011;differences design for causal inference, with robustness checks.Results show that a 1% increase in new energy development raises GPQI by 0.41% (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01). Solar and nuclear energy have significant positive effects, while hydropower and wind power show insignificant impacts. In addition,the impact of new energy development on food neo-productivity exhibits heterogeneity with respect to spatial, temporal, policy, economic, and financial development levels.These differences in land institutions and governance capacity can be viewed as a critical contextual backdrop and explanatory lens for understanding such heterogeneity.This study reveals a synergistic channel between energy and agriculture, indicating that energy transition can support sustainable and resilient grain production. The heterogeneous patterns provide empirical support for regionally differentiated policies coordinating energy and food systems.</p>
</abstract>
<kwd-group>
<kwd>clean energy</kwd>
<kwd>energy&#x2013;land synergy</kwd>
<kwd>entropy-weighting method</kwd>
<kwd>grain productivity quality index</kwd>
<kwd>new energy development</kwd>
</kwd-group>
<funding-group>
<funding-statement>The author(s) declare that financial support was received for the research and/or publication of this article. (1) Key Project of the Ministry of Education under the National Education Science Planning Program (DFA220434): Research on the Mechanism for Cross-Border Transfer of Scientific and Technological Achievements from Universities in Guangdong, Hong Kong, and Macao under the &#x201C;One Country, Two Systems&#x201D; Framework. (2) Pathways and Countermeasures for the High-Quality Development of Modern Agricultural Industrial Clusters, Special Entrusted Project of the Third Batch of Guangdong Province Philosophy and Social Sciences Innovation Project (GD25WTCXGC08).</funding-statement>
</funding-group>
<counts>
<fig-count count="2"/>
<table-count count="10"/>
<equation-count count="2"/>
<ref-count count="21"/>
<page-count count="13"/>
<word-count count="8791"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Land, Livelihoods and Food Security</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec1">
<label>1</label>
<title>Introduction</title>
<p>Frequent extreme weather and geopolitical conflicts have disrupted global food and energy supply chains, making energy and food security urgent global priorities. China, with its large population and limited arable land, has significantly contributed to global security in both fields but faces persistent pressures. By 2023, China&#x2019;s renewable energy installed capacity reached 1.4 billion kilowatts&#x2014;over half of total power capacity&#x2014;yet traditional agricultural energy dependence remains unsustainable. Developing Grain Productivity Quality Index, supported by new energy, has thus become vital for food security. In particular, this study focuses on four major types of new energy&#x2014;solar, wind, hydropower, and nuclear&#x2014;and emphasizes their differentiated applications and integration pathways within agricultural production systems.</p>
<p>Existing research recognizes energy&#x2013;agriculture linkages but suffers from narrow variable characterization and limited mechanism analysis, often focusing on either land competition or emission reductions, while overlooking systemic drivers like electrification, technological spillovers, and institutional coordination. The widespread adoption of agrivoltaic models in China since 2013 highlights the need for a more comprehensive analytical framework.</p>
<p>This study focuses on a balanced panel of 30 Chinese provinces over the period 2013&#x2013;2022. Regions with severe data limitations are excluded to ensure the consistency and comparability of key variables. The selected provinces collectively represent the vast majority of China&#x2019;s grain production, population, and economic activity, while also exhibiting substantial variation in energy structures, land endowments, agricultural practices, and institutional environments. This provincial coverage allows the analysis to capture both national trends and region-specific dynamics in the interaction between new-energy development and grain productivity quality.</p>
<p>Different types of new energy are expected to influence grain productivity quality through distinct mechanisms. Solar energy, which can be deployed in a decentralized manner and integrated with agricultural land use (e.g., agrivoltaic systems), directly supports on-farm electrification and resource-use efficiency. Nuclear energy provides stable, low-carbon baseload power, facilitating large-scale electrification of agricultural production, processing, and logistics. By contrast, hydropower and wind power are typically more centralized and grid-dependent, and their productivity effects may be diluted by transmission constraints, land-use conflicts, and limited direct integration with agricultural activities. These differences suggest heterogeneous impacts of energy types on the Grain Productivity Quality Index (GPQI), which are empirically examined in this study.</p>
<p>Recent scholarship in China has clarified the definition and assessment of &#x201C;new-quality&#x201D; agricultural productivity, emphasizing the need for a multidimensional indicator system that integrates technological, digital, green, integrative, and inclusive elements in line with agricultural modernization. Building on this literature, we measure grain-sector upgrading using a Grain Productivity Quality Index (GPQI) and test whether provincial new-energy development contributes to higher GPQI levels and how such effects vary across contexts. This study uses panel data from 30 Chinese provinces (2013&#x2013;2022) to examine how new energy development impacts Grain Productivity Quality Index. It constructs a multidimensional evaluation system&#x2014;using the share of new energy generation as a key explanatory variable and an entropy-weighted method for the productivity index. A two-way fixed effects model identifies causal effects, and a multi-period DID approach addresses endogeneity, exploiting the clean energy demonstration province policy as an exogenous shock. The analysis further explores heterogeneity across spatial, policy, economic, and financial dimensions.</p>
<p>This study makes three main contributions. First, it constructs a multidimensional Grain Productivity Quality Index (GPQI) to cap ture the technological, digital, green, integrative, and inclusive upgrad ing of grain production. Second, it provides quasi-causal evidence on the impact of new-energy development on GPQI by combining two way fixed effects with a multi-period DID design based on China&#x2019;s clean-energy demonstration policy. Third, it uncovers systematic het erogeneity across energy types and regional contexts, offering differ entiated insights for coordinating energy transition with agricultural modernization. These contributions are connected to a broad stream of related research, including <xref ref-type="bibr" rid="ref22">Zeng et al. (2018)</xref>, <xref ref-type="bibr" rid="ref2">Chi (2022)</xref>, <xref ref-type="bibr" rid="ref23">Zhang et al. (2025)</xref>, <xref ref-type="bibr" rid="ref3">Coomes et al. (2019)</xref>, <xref ref-type="bibr" rid="ref4">Gyamfi et al. (2022)</xref>, <xref ref-type="bibr" rid="ref7">Hu (2023)</xref>, <xref ref-type="bibr" rid="ref15">Power (2010)</xref>, <xref ref-type="bibr" rid="ref16">Shi et al. (2022)</xref>, <xref ref-type="bibr" rid="ref10">Ma et al. (2025)</xref>, and <xref ref-type="bibr" rid="ref21">Ying et al. (2025)</xref>.</p>
</sec>
<sec id="sec2">
<label>2</label>
<title>Literature review</title>
<sec id="sec3">
<label>2.1</label>
<title>The impact of new energy development on agricultural systems</title>
<p>In recent years, the synergy between new energy development and agricultural systems has become a major topic of global academic interest. Mainstream studies indicate that the application of new energy sources, especially solar, wind, and biomass, can not only enhance energy efficiency in agricultural production, but also promote the green transformation of agriculture and reduce carbon emissions (<xref ref-type="bibr" rid="ref6">Haberl et al., 2011</xref>; <xref ref-type="bibr" rid="ref14">Popp et al., 2014</xref>). <xref ref-type="bibr" rid="ref6">Haberl et al. (2011)</xref> systematically reviewed the complex interactions among land use, energy crops, and food production, highlighting the profound impact of new energy utilization models on agricultural ecology and productivity. <xref ref-type="bibr" rid="ref14">Popp et al. (2014)</xref> empirically demonstrated the important role of biomass energy in rural energy transitions, carbon emissions reduction, and enhancing agricultural resilience. Recent studies further show that new energy influences agricultural systems not only through energy substitution, but also via the electrification of irrigation, stabilization of production processes, and integration with land-use systems such as agrivoltaics. These findings suggest that the impact of new energy on agriculture should be evaluated from a systemic productivity-quality perspective rather than solely through output or emission indicators.</p>
</sec>
<sec id="sec4">
<label>2.2</label>
<title>Characterization of grain productivity quality index and new energy variables</title>
<p>In existing studies, new energy development is mostly measured by indicators such as installed capacity, electricity generation, or investment efficiency (<xref ref-type="bibr" rid="ref1">Chang et al., 2020</xref>; <xref ref-type="bibr" rid="ref19">Xu et al., 2024</xref>), while agricultural productivity has gradually evolved from traditional yield or total factor productivity to green total factor productivity (AGTFP), encompassing environmental, resource, and technological innovation dimensions (<xref ref-type="bibr" rid="ref3">Coomes et al., 2019</xref>; <xref ref-type="bibr" rid="ref1">Chang et al., 2020</xref>. In addition, green finance and policy incentives have been proven as important drivers for improving Grain Productivity Quality Index and promoting the application of new energy (<xref ref-type="bibr" rid="ref2">Chi, 2022</xref>). However, energy poverty is believed to significantly reduce agricultural productivity (<xref ref-type="bibr" rid="ref16">Shi et al., 2022</xref>). These multi-dimensional measurement approaches provide a solid foundation for research into the synergistic mechanisms between new energy and the grain system.</p>
</sec>
<sec id="sec5">
<label>2.3</label>
<title>Mechanisms, Chinese evidence, and research gaps</title>
<p>Research on the mechanisms between new energy and agriculture is mainly divided into two categories: one emphasizes the crowding-out effect of new energy development on high-quality arable land, i.e., competition in land resource allocation (<xref ref-type="bibr" rid="ref5">Haberl, 2015</xref>); the other focuses on technological spillovers and green synergies brought by new energy, such as the improvement in land use efficiency through agrivoltaic projects (<xref ref-type="bibr" rid="ref12">Omer et al., 2022</xref>). <xref ref-type="bibr" rid="ref5">Haberl (2015)</xref> discussed the dual pressures of new energy project siting on food security and ecological protection, warning of potential negative externalities. <xref ref-type="bibr" rid="ref12">Omer et al. (2022)</xref> empirically revealed that agrivoltaic projects can increase land economic output, reduce carbon emissions, and promote ecological protection, serving as an innovative agricultural model driven by new energy.</p>
<p>Supporting literature shows an increasing focus on the link between energy structures and agricultural productivity. For example, <xref ref-type="bibr" rid="ref13">Peng et al. (2024)</xref> construct an index of agricultural green productivity using panel data for 30 Chinese provinces, highlighting the spatial&#x2013;temporal characteristics of agricultural green productivity and the influence of energy structure on its performance. <xref ref-type="bibr" rid="ref17">Wang and Qian (2024)</xref> emphasizes the relationship between agricultural productivity and the pursuit of a just energy transition in China, pointing to the need for integrated approaches that align energy transition goals with agricultural development and sustainability. <xref ref-type="bibr" rid="ref9">Li et al. (2022)</xref> find that both renewable energy consumption and demographic transitions have significant long-run impacts on agricultural green total factor productivity across 30 provinces, suggesting that energy use patterns shape productivity outcomes in rural China. Other provincial-level studies relate energy poverty and agricultural technical efficiency, showing that energy constraints can hinder improvements in agricultural productivity under certain conditions. Together, these studies underline that energy structure and energy use dynamics are important determinants of agricultural productivity at the provincial level, and they motivate further investigation into how new energy development interacts with multi-dimensional measures of productivity quality such as GPQI.</p>
<p>Since 2013, new energy projects such as agrivoltaics have rapidly proliferated in China, driving the green transformation of agricultural systems. Relevant studies have focused on policy promotion, regional disparities, and innovative mechanisms. For example, Xinhua News Agency (2024) reported the central government&#x2019;s introduction of the Grain Productivity Quality Index concept, marking a new stage in China&#x2019;s food security strategy. Furthermore, data from the National Energy Administration indicate that by 2023, China&#x2019;s installed capacity of renewable energy surpassed 1.4 billion kilowatts, reflecting the deep integration of new energy and food security under policy-driven initiatives.</p>
<p>In summary, existing domestic and international studies have preliminarily revealed the interactive mechanisms between new energy development and agricultural production. However, further advancements are needed in variable characterization, mechanism identification, and system modeling. Future research should broaden multi-dimensional evaluation indicators, refine heterogeneity analysis, and apply frontier methods such as policy simulation and causal inference to systematically assess the boundary and practical pathways of new energy&#x2019;s impact on Grain Productivity Quality Index. This will not only provide a theoretical basis for the synergy between China&#x2019;s dual carbon and food security strategies, but also contribute valuable Chinese experience to the global transition toward sustainable agriculture.</p>
</sec>
</sec>
<sec id="sec6">
<label>3</label>
<title>Research design</title>
<sec id="sec7">
<label>3.1</label>
<title>Sample selection and data sources</title>
<p>We compile a provincial panel for 30 Chinese provinces, 2013&#x2013;2022. The outcome variable is GPQI. The key explanatory variable measures new-energy development (clean-energy structure and deployment consistent with national statistical practice), and the control vector includes economic structure, market performance, agricultural inputs, and finance-related variables. All monetary variables are deflated to constant prices (base year specified in the data appendix) and log-transformed where appropriate to reduce skewness. All variables are constructed using publicly available data from the China Statistical Yearbook, China Rural Statistical Yearbook, China Energy Statistical Yearbook, and official provincial statistical bulletins. Detailed variable definitions and data sources are summarized in <xref ref-type="table" rid="tab1">Tables 1</xref>, <xref ref-type="table" rid="tab2">2</xref>.</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption>
<p>Evaluation index system of GPQI.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">First-level indicators</th>
<th align="left" valign="top">Secondary indicators</th>
<th align="left" valign="top">Third-level indicators</th>
<th align="left" valign="top">Measuring methods</th>
<th align="left" valign="top">Indicator properties</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle" rowspan="6">Labourer</td>
<td align="left" valign="middle" rowspan="2">Workers&#x2019; skills</td>
<td align="left" valign="middle">Educational level of rural residents</td>
<td align="left" valign="middle">(Number of rural people without schooling&#x202F;&#x00D7;&#x202F;1&#x202F;+&#x202F;Number of primary school graduates&#x202F;&#x00D7;&#x202F;6&#x202F;+&#x202F;Number of junior high school graduates&#x202F;&#x00D7;&#x202F;9&#x202F;+&#x202F;Number of high school and technical secondary school graduates&#x202F;&#x00D7;&#x202F;12&#x202F;+&#x202F;Number of people with college degree or above&#x202F;&#x00D7;&#x202F;16)/The total rural population over the age of 6</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">The proportion of rural residents receiving higher education</td>
<td align="left" valign="middle">Number of rural residents with college degree or above/The total number of people over the age of 6 in rural areas</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">productivity of labour</td>
<td align="left" valign="middle">Per capita output of grain</td>
<td align="left" valign="middle">The ratio of grain production to the labor force employed in growing food crops</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Per capita output value</td>
<td align="left" valign="middle">Per capita GDP</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">Workers&#x2019; awareness</td>
<td align="left" valign="middle">Labour force level</td>
<td align="left" valign="middle">Primary industry employees&#x202F;&#x00D7;&#x202F;(Total value of farm output/Total output value of agriculture, forestry, animal husbandry and fishery)&#x202F;&#x00D7;&#x202F;(Sown area of grain/Crop acreage)</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Labour force structure</td>
<td align="left" valign="middle">The proportion of employees in the tertiary industry in the total number of employees</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="9">Subject of labour</td>
<td align="left" valign="middle">Technical progress</td>
<td align="left" valign="middle">The comprehensive mechanization rate of grain cultivation and harvest</td>
<td align="left" valign="middle">Mechanized farming rate&#x202F;&#x00D7;&#x202F;40%&#x202F;+&#x202F;Mechanized sowing rate&#x202F;&#x00D7;&#x202F;30%&#x202F;+&#x202F;Machine yield&#x202F;&#x00D7;&#x202F;30%</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="4">Food safety</td>
<td align="left" valign="middle">Grain yield per unit area</td>
<td align="left" valign="middle">Total grain output/Sown area of grain</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Multiple-crop index</td>
<td align="left" valign="middle">Total sown area of crops/agricultural acreage</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Volatility of food production</td>
<td align="left" valign="middle">(Annual grain production-Annual average grain yield)/Mean value</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle">Crop disaster rate</td>
<td align="left" valign="middle">Crop disaster area/Crop acreage</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="4">Green development</td>
<td align="left" valign="middle">Fertilizer application intensity</td>
<td align="left" valign="middle">Agricultural fertilizer conversion rate usage/Total sown area of crops</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle">Intensity of pesticide application</td>
<td align="left" valign="middle">Pesticide application/Total sown area of crops</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle">Water use level</td>
<td align="left" valign="middle">Effective irrigation area /Total sown area of crops</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Carbon emissions from food</td>
<td align="left" valign="middle">Carbon emissions per unit sown area</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="5">Means of labor</td>
<td align="left" valign="middle" rowspan="2">Infrastructure</td>
<td align="left" valign="middle">Agricultural machinery and equipment</td>
<td align="left" valign="middle">Total power of agricultural machinery</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Observations from agricultural weather stations</td>
<td align="left" valign="middle">Direct access</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="2">Energy consumption</td>
<td align="left" valign="middle">Rural energy intensity</td>
<td align="left" valign="middle">Electricity consumption in agriculture, forestry, animal husbandry and fishery/Rural population</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
<tr>
<td align="left" valign="middle">Diesel usage</td>
<td align="left" valign="middle">Diesel use for agriculture/Total value of farm output</td>
<td align="left" valign="middle">Negative direction</td>
</tr>
<tr>
<td align="left" valign="middle">Technological innovation</td>
<td align="left" valign="middle">Intensity of agricultural R&#x0026;D investment</td>
<td align="left" valign="middle">Internal expenditure of R&#x0026;D funds&#x202F;&#x00D7;&#x202F;(Total output value of agriculture, forestry, animal husbandry and fishery/Gross regional product)</td>
<td align="left" valign="middle">Forward direction</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption>
<p>Variable definition table.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th>Variable type</th>
<th align="left" valign="top">Variable name</th>
<th align="center" valign="top">Variable abbreviations</th>
<th align="left" valign="top">Variable definition and description</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Explained variable</td>
<td align="left" valign="middle">Grain Productivity Quality Index (GPQI)</td>
<td align="center" valign="middle">GPQI</td>
<td align="left" valign="middle">The evaluation index system of grain GPQI in <xref ref-type="table" rid="tab1">Table 1</xref> was constructed and calculated by entropy method</td>
</tr>
<tr>
<td align="left" valign="middle">Explanatory variable</td>
<td align="left" valign="middle">The development level of new energy</td>
<td align="center" valign="middle">NE</td>
<td align="left" valign="middle">The sum of wind, hydro, solar and nuclear power generation/Total social electricity generation</td>
</tr>
<tr>
<td align="left" valign="middle" rowspan="4">Controlled variable</td>
<td align="left" valign="middle">Level of informatization</td>
<td align="center" valign="middle">Infor</td>
<td align="left" valign="middle">Total volume of postal and telecommunications business/Gross regional product</td>
</tr>
<tr>
<td align="left" valign="middle">Environmental regulation intensity</td>
<td align="center" valign="middle">Envir</td>
<td align="left" valign="middle">Investment in industrial pollution control has been completed/Industrial value added</td>
</tr>
<tr>
<td align="left" valign="middle">Industrialization level</td>
<td align="center" valign="middle">Indus</td>
<td align="left" valign="middle">Industrial value added/Gross regional product</td>
</tr>
<tr>
<td align="left" valign="middle">The level of aging</td>
<td align="center" valign="middle">Old</td>
<td align="left" valign="middle">The elderly population aged 65 and over/Total population</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="sec8">
<label>3.2</label>
<title>Variable description</title>
<sec id="sec9">
<label>3.2.1</label>
<title>Level of new energy development</title>
<p>The analysis considers wind, solar, hydro, and nuclear power as core indicators of new energy. The extent of new energy development in each province is measured by the proportion of total electricity generated from these sources. This metric provides a clear gauge of new energy adoption across provinces, with higher values signifying more advanced development of new energy infrastructure. To ensure the homogeneity of the research subjects, enable a deeper analysis of specific energy challenges, and align with common academic delineations of &#x201C;new energy,&#x201D; this study confines its scope to typical emerging energy sources like solar and wind power, excluding biomass energy from the current research framework (<xref ref-type="bibr" rid="ref4">Gyamfi et al., 2022</xref>).</p>
</sec>
<sec id="sec10">
<label>3.2.2</label>
<title>Measurement of the outcome: GPQI</title>
<p>Recent scholarship in China has clarified the multidimensional assessment of agricultural productivity, emphasizing a multidimensional indicator system that incorporates technological, digital, green, integrative, and inclusive elements in line with the aims and realities of agricultural modernization (<xref ref-type="bibr" rid="ref18">Xiong and He, 2025</xref>). Related empirical work has estimated agricultural green total factor productivity (AGTFP) and examined its determinants, underscoring the roles of technological advancement, energy efficiency, emissions reduction, and market performance while noting persistent regional and factor-based disparities (<xref ref-type="bibr" rid="ref8">Huang et al., 2022</xref>). Building on this literature and the classic three-element framework of productivity (labor, objects of labor, and means of labor), we construct a multidimensional framework for grain production. We operationalize the Grain Productivity Quality Index (GPQI) as a composite measure integrating the technological, digital, green, integrative, and inclusive dimensions, and compute province&#x2013;year GPQI scores using the entropy-weighting method. While the entropy-weighting method assigns weights based on the information variability of each indicator and thus avoids subjective judgment, it may not fully reflect the normative or policy relevance of different dimensions of grain productivity quality. To address this concern, we further examine the sensitivity of the results to alternative weighting schemes. The complete indicator system and variable definitions are reported in <xref ref-type="table" rid="tab1">Table 1</xref>.</p>
</sec>
<sec id="sec11">
<label>3.2.3</label>
<title>Control variables</title>
<p>To account for additional determinants of Grain Productivity Quality Index, several control variables are incorporated into the analysis: (1) Informatization Level, calculated as the ratio of total postal and telecommunications business to regional GDP; (2) Environmental Regulation Intensity, defined by the proportion of investment in industrial pollution management relative to industrial value added; (3) Industrialization Level, expressed as industrial value added divided by regional GDP; (4) Aging Population, measured by the share of individuals aged 65 and above within the total population; and (5) Agricultural Industrial Structure, represented by the ratio of gross agricultural output to the aggregate output of agriculture, forestry, animal husbandry, and fishery. A comprehensive description of each variable is provided in <xref ref-type="table" rid="tab2">Table 2</xref>.</p>
</sec>
</sec>
<sec id="sec12">
<label>3.3</label>
<title>Model specification</title>
<p>To examine the impact of new energy development on the Grain Productivity Quality Index (GPQI), The baseline regression model is specified as <xref ref-type="disp-formula" rid="E1">Equation 1</xref> below.</p>
<disp-formula id="E1">
<label>(1)</label>
<mml:math id="M1">
<mml:mtable columnalign="left" displaystyle="true">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mtext>GPQI</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>NE</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Infor</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Envir</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Indus</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>Old</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mtext>ProvinceFE</mml:mtext>
<mml:mo>+</mml:mo>
<mml:mtext>YearFE</mml:mtext>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>In <xref ref-type="disp-formula" rid="E1">Equation 1</xref> <inline-formula>
<mml:math id="M2">
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:math>
</inline-formula> is the constant term, <inline-formula>
<mml:math id="M3">
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:math>
</inline-formula> is the main regression coefficient, <inline-formula>
<mml:math id="M4">
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>~</mml:mo>
<mml:msub>
<mml:mi>&#x03B2;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
</mml:math>
</inline-formula> represent the coefficients associated with the control variables, Province fixed effects and year fixed effects are included in the specifications, <inline-formula>
<mml:math id="M5">
<mml:mi>&#x03B5;</mml:mi>
</mml:math>
</inline-formula><sub><italic>i</italic>,<italic>t</italic></sub> as the error term.</p>
</sec>
</sec>
<sec id="sec13">
<label>4</label>
<title>Empirical results and analysis</title>
<sec id="sec14">
<label>4.1</label>
<title>Descriptive statistics</title>
<p><xref ref-type="table" rid="tab3">Table 3</xref> presents the descriptive statistics of the variables, providing a detailed analysis of seven key variables. Overall, the means of the variables indicate a moderate level of development within the sample; however, the standard deviations and ranges reveal significant internal variation, indicating substantial regional disparities. Specifically, the mean values of Grain Productivity Quality Index (0.324) and the level of industrialization (0.312) are both in the moderate range, and their relatively low standard deviations (0.087 and 0.077, respectively) suggest that these indicators are fairly stable. In contrast, the standard deviations for new energy generation (mean&#x202F;=&#x202F;0.307, SD&#x202F;=&#x202F;0.388) and informatization level (mean 0.073, SD 0.153) are relatively large, with extremely low minimum values (as low as 0.004 for new energy generation and 0.015 for informatization) and high maximum values (5.447 and 2.520, respectively). This implies regional imbalance, with some sample points exhibiting weak performance in new energy or information technology application, while others perform exceptionally well. The mean for environmental regulation intensity is extremely low (0.003), but the range is wide (from 0.00006 to 0.031), indicating generally weak environmental policy intensity, though a few regions enforce stricter regulation. In addition, the mean proportion of elderly population is 16.33%, reflecting an overall trend of population aging, but its standard deviation (4.396) and range (from 8.75 to 28.77%) highlight significant differences in population structure. Overall, the data reveal development imbalances, especially in new energy, informatization, and demographic structure, suggesting that targeted policies may be necessary to narrow regional gaps.</p>
<table-wrap position="float" id="tab3">
<label>Table 3</label>
<caption>
<p>Descriptive statistics.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Variables</th>
<th align="center" valign="top">Obs</th>
<th align="center" valign="top">Mean</th>
<th align="center" valign="top">Std. dev.</th>
<th align="center" valign="top">Min</th>
<th align="center" valign="top">Max</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">GPQI</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">0.324</td>
<td align="char" valign="middle" char=".">0.087</td>
<td align="char" valign="middle" char=".">0.147</td>
<td align="char" valign="middle" char=".">0.592</td>
</tr>
<tr>
<td align="left" valign="middle">NE</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">0.307</td>
<td align="char" valign="middle" char=".">0.388</td>
<td align="char" valign="middle" char=".">0.004</td>
<td align="char" valign="middle" char=".">5.447</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">0.073</td>
<td align="char" valign="middle" char=".">0.153</td>
<td align="char" valign="middle" char=".">0.015</td>
<td align="char" valign="middle" char=".">2.520</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">0.003</td>
<td align="char" valign="middle" char=".">0.004</td>
<td align="char" valign="middle" char=".">0.000</td>
<td align="char" valign="middle" char=".">0.031</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">0.312</td>
<td align="char" valign="middle" char=".">0.077</td>
<td align="char" valign="middle" char=".">0.101</td>
<td align="char" valign="middle" char=".">0.498</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">300</td>
<td align="char" valign="middle" char=".">16.330</td>
<td align="char" valign="middle" char=".">4.396</td>
<td align="char" valign="middle" char=".">8.750</td>
<td align="char" valign="middle" char=".">28.770</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="sec15">
<label>4.2</label>
<title>Analysis of baseline regression results</title>
<p><xref ref-type="table" rid="tab4">Table 4</xref> displays the results of the baseline regression examining the relationship between new energy development and Grain Productivity Quality Index. Column (1) shows that, in the absence of control variables, the coefficient for the primary explanatory variable (GPQI) is significantly positive at the 5% significance level. When control variables are introduced in column (2), the significance of GPQI increases, reaching the 1% level, and the adjusted <italic>R</italic><sup>2</sup> also rises, reflecting a better model fit. Overall, these results indicate that the advancement of new energy development has a clear positive effect on Grain Productivity Quality Index.</p>
<table-wrap position="float" id="tab4">
<label>Table 4</label>
<caption>
<p>Benchmark regression results.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
</tr>
<tr>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">NE</td>
<td align="center" valign="middle">0.0028&#x002A;&#x002A; (2.1536)</td>
<td align="center" valign="middle">0.0041&#x002A;&#x002A;&#x002A; (3.6320)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td/>
<td align="center" valign="middle">0.0058 (1.3734)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td/>
<td align="center" valign="middle">&#x2212;1.2114&#x002A;&#x002A; (&#x2212;2.4859)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td/>
<td align="center" valign="middle">&#x2212;0.1756&#x002A;&#x002A;&#x002A; (&#x2212;3.4105)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td/>
<td align="center" valign="middle">&#x2212;0.0021&#x002A; (&#x2212;1.7162)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.3226&#x002A;&#x002A;&#x002A; (258.7843)</td>
<td align="center" valign="middle">0.4141&#x002A;&#x002A;&#x002A; (21.2955)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.948</td>
<td align="center" valign="middle">0.952</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec16">
<label>4.3</label>
<title>Robustness tests</title>
<sec id="sec17">
<label>4.3.1</label>
<title>Alternative explanatory variables</title>
<p>To assess the robustness of the findings, this paper substitutes the original new energy development variable (NE) with two alternatives: the aggregate amount of new energy generated (NE1) and the proportion of new energy generation relative to regional GDP (NE2). As displayed in columns (1) and (2) of <xref ref-type="table" rid="tab5">Table 5</xref>, both alternatives yield coefficients that are significantly positive, further reinforcing the conclusion that new energy development substantially enhances Grain Productivity Quality Index.</p>
<table-wrap position="float" id="tab5">
<label>Table 5</label>
<caption>
<p>Results of robustness test.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
<th align="center" valign="top">(3)</th>
<th align="center" valign="top">(4)</th>
</tr>
<tr>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">NE1</td>
<td align="center" valign="middle">0.0372&#x002A;&#x002A; (2.3351)</td>
<td/>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">NE2</td>
<td/>
<td align="center" valign="middle">0.0323&#x002A; (1.8777)</td>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">NE</td>
<td/>
<td/>
<td align="center" valign="middle">0.0037&#x002A;&#x002A;&#x002A; (3.1512)</td>
<td align="center" valign="middle">0.0041&#x002A;&#x002A;&#x002A; (3.7363)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">0.0057 (1.3583)</td>
<td align="center" valign="middle">0.0058 (1.3873)</td>
<td align="center" valign="middle">0.0059 (1.5018)</td>
<td align="center" valign="middle">0.0057 (1.5018)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">&#x2212;1.2176&#x002A;&#x002A; (&#x2212;2.4972)</td>
<td align="center" valign="middle">&#x2212;1.1618&#x002A;&#x002A; (&#x2212;2.3748)</td>
<td align="center" valign="middle">&#x2212;0.9620&#x002A;&#x002A; (&#x2212;1.9714)</td>
<td align="center" valign="middle">&#x2212;1.3386&#x002A;&#x002A;&#x002A; (&#x2212;2.7941)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">&#x2212;0.1742&#x002A;&#x002A;&#x002A; (&#x2212;3.3864)</td>
<td align="center" valign="middle">&#x2212;0.1756&#x002A;&#x002A;&#x002A; (&#x2212;3.4049)</td>
<td align="center" valign="middle">&#x2212;0.1700&#x002A;&#x002A;&#x002A; (&#x2212;2.9514)</td>
<td align="center" valign="middle">&#x2212;0.2054&#x002A;&#x002A;&#x002A; (&#x2212;3.6511)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">&#x2212;0.0020&#x002A; (&#x2212;1.6914)</td>
<td align="center" valign="middle">&#x2212;0.0020&#x002A; (&#x2212;1.7020)</td>
<td align="center" valign="middle">&#x2212;0.0028&#x002A;&#x002A; (&#x2212;2.0395)</td>
<td align="center" valign="middle">&#x2212;0.0015 (&#x2212;1.1958)</td>
</tr>
<tr>
<td align="left" valign="middle">Trade</td>
<td/>
<td/>
<td/>
<td align="center" valign="middle">&#x2212;0.0562&#x002A;&#x002A;&#x002A; (&#x2212;2.8172)</td>
</tr>
<tr>
<td align="left" valign="middle">Income</td>
<td/>
<td/>
<td/>
<td align="center" valign="middle">&#x2212;0.0723&#x002A;&#x002A;&#x002A; (&#x2212;2.9245)</td>
</tr>
<tr>
<td align="left" valign="middle">RD</td>
<td/>
<td/>
<td/>
<td align="center" valign="middle">0.0190&#x002A; (1.6650)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.4135&#x002A;&#x002A;&#x002A; (21.3067)</td>
<td align="center" valign="middle">0.4137&#x002A;&#x002A;&#x002A; (21.2473)</td>
<td align="center" valign="middle">0.4271&#x002A;&#x002A;&#x002A; (20.1093)</td>
<td align="center" valign="middle">0.5771&#x002A;&#x002A;&#x002A; (9.0507)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">260</td>
<td align="center" valign="middle">300</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.956</td>
<td align="center" valign="middle">0.954</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec18">
<label>4.3.2</label>
<title>Excluding municipality samples</title>
<p>Due to the distinctive agricultural structures and policy advantages of centrally administered municipalities&#x2014;which could potentially bias or limit the external validity of the results&#x2014;this analysis removes Beijing, Shanghai, Tianjin, and Chongqing from the dataset and reruns the regression. The findings, presented in column (3) of <xref ref-type="table" rid="tab5">Table 5</xref>, reveal that the coefficient for new energy development remains significantly positive at the 1% threshold, consistent with the primary results of the baseline model.</p>
</sec>
<sec id="sec19">
<label>4.3.3</label>
<title>Including additional control variables</title>
<p>Building on the baseline regression, this study further incorporates additional control variables, including trade dependence (Trade), the urban&#x2013;rural income ratio (Income), and R&#x0026;D investment intensity (RD). As shown in column (4) of <xref ref-type="table" rid="tab5">Table 5</xref>, even after including these additional control variables, new energy development continues to significantly promote the improvement of Grain Productivity Quality Index at the 1% significance level.</p>
</sec>
</sec>
<sec id="sec20">
<label>4.4</label>
<title>Endogeneity test</title>
<p>To further mitigate potential endogeneity issues in the empirical analysis, this study applies a difference-in-differences (DID) framework to reassess the relationship between new energy development and Grain Productivity Quality Index. In December 2016, the National Energy Administration released the &#x201C;13th Five-Year Plan for Energy Development&#x201D; prioritizing the creation of clean energy demonstration provinces and regions as a central initiative for advancing the energy consumption revolution, and identifying several provinces as key pilots. Each selected province subsequently formulated detailed implementation strategies consistent with the national policy agenda. These pilots have focused on building clean, low-carbon, safe, and efficient energy systems, fostering the transformation of growth drivers, and supporting high-quality economic development. Rather than relying on traditional instrumental variables, this study addresses endogeneity concerns by exploiting the staggered implementation of the clean-energy demonstration province policy as a quasi-exogenous shock within a multi-period DID framework.</p>
<p>Within this context, the clean energy demonstration province policy is treated as an external shock to evaluate its influence on Grain Productivity Quality Index. This approach is justified as all types of new energy considered in this study&#x2014;hydropower, wind, solar, and nuclear&#x2014;are classified as clean energy, and related policies are expected to shape their adoption and expansion. The exogeneity of the policy&#x2019;s rollout makes it suitable for addressing endogeneity challenges in the model.</p>
<p>To measure the policy&#x2019;s impact, a multi-period DID strategy is implemented, treating the policy&#x2019;s introduction as a quasi-natural experiment. The DID method reduces selection bias and helps alleviate endogeneity, while the multi-period approach accommodates policy effects that unfold over time. This methodology statistically estimates the average treatment effect by leveraging all possible combinations of periods and treatment groups, allowing for a more nuanced and accurate assessment of the economic impacts associated with staggered policy implementation.</p>
<p>In operationalizing this approach, the actual year of policy implementation in each pilot province is used as the intervention point, dividing the sample into treatment and control groups to compare shifts in Grain Productivity Quality Index before and after the policy. Based on official records, the designated pilot provinces and their respective years of implementation are: Ningxia (2012), Zhejiang (2014), Sichuan (2016), Tibet (2016), Gansu (2016), and Qinghai (2018). These six provinces comprise the treatment group, with the remaining 25 provinces serving as the control. Due to extensive missing data, Tibet is excluded from the analysis, resulting in a final treatment group of five provinces. The model structure is outlined in <xref ref-type="disp-formula" rid="E2">Equation (2).</xref></p>
<disp-formula id="E2">
<label>(2)</label>
<mml:math id="M6">
<mml:mtable columnalign="left" displaystyle="true">
<mml:mtr>
<mml:mtd>
<mml:mi>D</mml:mi>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Policy</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Infor</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>3</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Envir</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>4</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Indus</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>5</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>Old</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03C9;</mml:mi>
<mml:mn>6</mml:mn>
</mml:msub>
<mml:msub>
<mml:mtext>Stru</mml:mtext>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>+</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtext>ProvinceFE</mml:mtext>
<mml:mo>+</mml:mo>
<mml:mtext>YearFE</mml:mtext>
<mml:mo>+</mml:mo>
<mml:msub>
<mml:mi>&#x03B5;</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>In Formula 2 <inline-formula>
<mml:math id="M7">
<mml:msub>
<mml:mtext>Policy</mml:mtext>
<mml:mi mathvariant="italic">it</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> is a dummy variable indicating whether province <italic>i</italic> has implemented the clean energy demonstration province policy in year t; it takes the value of 1 if and only if the pilot province i is in the year of policy implementation or any year thereafter, and 0 otherwise.</p>
<p>In columns (1) and (2) of <xref ref-type="table" rid="tab6">Table 6</xref>, <inline-formula>
<mml:math id="M8">
<mml:msub>
<mml:mtext>Policy</mml:mtext>
<mml:mi mathvariant="italic">it</mml:mi>
</mml:msub>
</mml:math>
</inline-formula> coefficients are all significantly positive, These results demonstrate that the clean energy demonstration province policy exerts a significant positive effect on Grain Productivity Quality Index, regardless of the inclusion of control variables. Since Ningxia implemented the policy in 2012, but the study period starts in 2013, Ningxia is removed from the sample to avoid potential bias, and the regression is conducted again. As presented in column (3) of <xref ref-type="table" rid="tab7">Table 7</xref>, both the magnitude and significance of the coefficients remain consistent. Overall, the findings from the multi-period DID approach align with the primary regression outcomes, further verifying the robustness and reliability of the main results.</p>
<table-wrap position="float" id="tab6">
<label>Table 6</label>
<caption>
<p>Endogeneity test: multi-time DID test.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
<th align="center" valign="top">(3)</th>
</tr>
<tr>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
<th align="center" valign="top">GPQI</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">did</td>
<td align="center" valign="middle">0.0167&#x002A;&#x002A;&#x002A; (2.9150)</td>
<td align="center" valign="middle">0.0145&#x002A;&#x002A; (2.5527)</td>
<td align="center" valign="middle">0.0154&#x002A;&#x002A;&#x002A; (2.6696)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td/>
<td align="center" valign="middle">0.0050 (1.2646)</td>
<td align="center" valign="middle">0.0050 (1.3041)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td/>
<td align="center" valign="middle">&#x2212;1.1382&#x002A;&#x002A; (&#x2212;2.3014)</td>
<td align="center" valign="middle">&#x2212;0.7760 (&#x2212;0.9924)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td/>
<td align="center" valign="middle">&#x2212;0.1720&#x002A;&#x002A;&#x002A; (&#x2212;3.3476)</td>
<td align="center" valign="middle">&#x2212;0.1697&#x002A;&#x002A;&#x002A; (&#x2212;3.0309)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td/>
<td align="center" valign="middle">&#x2212;0.0019 (&#x2212;1.5898)</td>
<td align="center" valign="middle">&#x2212;0.0019 (&#x2212;1.5257)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.3214&#x002A;&#x002A;&#x002A; (243.7443)</td>
<td align="center" valign="middle">0.4095&#x002A;&#x002A;&#x002A; (21.0744)</td>
<td align="center" valign="middle">0.4096&#x002A;&#x002A;&#x002A; (20.3292)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">290</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.949</td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.953</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="tab7">
<label>Table 7</label>
<caption>
<p>Heterogeneity test of development level of different types of new energy.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
<th align="center" valign="top">(3)</th>
<th align="center" valign="top">(4)</th>
</tr>
<tr>
<th align="center" valign="top">Hydroenergy</th>
<th align="center" valign="top">Wind energy</th>
<th align="center" valign="top">Solar energy</th>
<th align="center" valign="top">Nuclear energy</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Water</td>
<td align="center" valign="middle">&#x2212;0.0084 (&#x2212;0.2915)</td>
<td/>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">Wind</td>
<td/>
<td align="center" valign="middle">&#x2212;0.1074 (&#x2212;1.4833)</td>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">Sun</td>
<td/>
<td/>
<td align="center" valign="middle">0.0038&#x002A;&#x002A;&#x002A; (3.4338)</td>
<td/>
</tr>
<tr>
<td align="left" valign="middle">Nuclear</td>
<td/>
<td/>
<td/>
<td align="center" valign="middle">0.0744&#x002A; (1.7664)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">0.0059 (1.3807)</td>
<td align="center" valign="middle">0.0056 (1.2920)</td>
<td align="center" valign="middle">0.0058 (1.3773)</td>
<td align="center" valign="middle">0.0052 (1.2223)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">&#x2212;1.2244&#x002A;&#x002A; (&#x2212;2.4501)</td>
<td align="center" valign="middle">&#x2212;1.3096&#x002A;&#x002A;&#x002A; (&#x2212;2.7321)</td>
<td align="center" valign="middle">&#x2212;1.2193&#x002A;&#x002A; (&#x2212;2.4976)</td>
<td align="center" valign="middle">&#x2212;1.2472&#x002A;&#x002A;&#x002A; (&#x2212;2.6442)</td>
</tr>
<tr>
<td/>
<td/>
<td/>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">&#x2212;0.1736&#x002A;&#x002A;&#x002A; (&#x2212;3.3695)</td>
<td align="center" valign="middle">&#x2212;0.1671&#x002A;&#x002A;&#x002A; (&#x2212;3.1888)</td>
<td align="center" valign="middle">&#x2212;0.1750&#x002A;&#x002A;&#x002A; (&#x2212;3.3966)</td>
<td align="center" valign="middle">&#x2212;0.1778&#x002A;&#x002A;&#x002A; (&#x2212;3.4845)</td>
</tr>
<tr>
<td/>
<td/>
<td/>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">&#x2212;0.0019 (&#x2212;1.5620)</td>
<td align="center" valign="middle">&#x2212;0.0019 (&#x2212;1.5889)</td>
<td align="center" valign="middle">&#x2212;0.0020&#x002A; (&#x2212;1.7029)</td>
<td align="center" valign="middle">&#x2212;0.0018 (&#x2212;1.4726)</td>
</tr>
<tr>
<td/>
<td/>
<td/>
<td/>
<td/>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.4145&#x002A;&#x002A;&#x002A; (21.2812)</td>
<td align="center" valign="middle">0.4155&#x002A;&#x002A;&#x002A; (21.6349)</td>
<td align="center" valign="middle">0.4146&#x002A;&#x002A;&#x002A; (21.2677)</td>
<td align="center" valign="middle">0.4094&#x002A;&#x002A;&#x002A; (20.0982)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
<td align="center" valign="middle">300</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.952</td>
<td align="center" valign="middle">0.952</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
<p>Considering that the validity of the multi-period DID model relies on the parallel trends assumption, <xref ref-type="fig" rid="fig1">Figure 1</xref> presents the coefficient dynamics of the parallel trends test for the impact of new energy development on Grain Productivity Quality Index. Here, &#x201C;current&#x201D; represents the year of policy implementation, while &#x201C;pre 5&#x201D; &#x201C;pre 4&#x201D; &#x201C;pre 3&#x201D; &#x201C;pre 2&#x201D; and &#x201C;pre 1&#x201D; denote five, four, three, two, and 1&#x202F;year(s) before policy implementation, respectively. Similarly, &#x201C;las 1&#x201D; &#x201C;las 2&#x201D; &#x201C;las 3&#x201D; &#x201C;las 4&#x201D; &#x201C;las 5&#x201D; &#x201C;las 6&#x201D; &#x201C;las 7&#x201D; and &#x201C;las 8&#x201D; refer to one to 8&#x202F;years after policy implementation.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption>
<p>Parallel trend test results.</p>
</caption>
<graphic xlink:href="fsufs-10-1765501-g001.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Line graph with error bars displays policy effect versus policy time, showing data points from negative five to eight on the x-axis, a vertical red line at zero, and policy effect values ranging from negative 0.02 to 0.04.</alt-text>
</graphic>
</fig>
<p>The parallel trends test is conducted by selecting &#x201C;pre_1&#x201D; as the baseline period. According to the results in <xref ref-type="fig" rid="fig1">Figure 1</xref>, the treatment effect coefficients for &#x201C;current&#x201D; and subsequent years are all significantly greater than zero, while the 95% confidence intervals for the coefficients prior to &#x201C;current&#x201D; mostly include zero, indicating that the parallel trends assumption holds. In summary, the multi-period DID model satisfies the parallel trends condition, confirming that new energy development has a significant positive effect on Grain Productivity Quality Index.</p>
<p>To further ensure the unbiasedness of the estimated coefficients in the multi-period DID model, a placebo test is conducted in this study. The Bootstrap method is used to simulate random group assignments between the treatment and control groups. As shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, the circles represent the kernel density distribution of the estimated coefficients from the Bootstrap randomization experiments, while the dashed line indicates the estimated coefficient from the multi-period DID model. It can be seen from <xref ref-type="fig" rid="fig2">Figure 2</xref> that the estimated coefficient of the multi-period DID model in this study is unbiased.</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption>
<p>Placebo test results.</p>
</caption>
<graphic xlink:href="fsufs-10-1765501-g002.tif" mimetype="image" mime-subtype="tiff">
<alt-text content-type="machine-generated">Scatter plot and kernel density plot overlay showing blue circles for p_value on the left y-axis and brown line for kdensity beta on the right y-axis against _b[did] on the x-axis, with a red dashed vertical reference line near 0.017.</alt-text>
</graphic>
</fig>
</sec>
<sec id="sec21">
<label>4.5</label>
<title>Heterogeneity analysis by type of new energy</title>
<p>To further analyze the heterogeneous effects of new energy, the explanatory variable is disaggregated, and the development levels of hydropower, wind power, solar power, and nuclear power in each province are measured by the proportion of each type&#x2019;s electricity generation in total energy generation. The impact of each energy type on Grain Productivity Quality Index is then analyzed, with regression results shown in <xref ref-type="table" rid="tab7">Table 7</xref>. First, as indicated in columns (1) and (2), hydropower (Water) and wind power (Wind) do not have significant effects on Grain Productivity Quality Index. For hydropower, this is primarily because resource distribution does not align with major grain-producing areas, and hydropower cannot be directly utilized at agricultural production sites. Hydropower stations are typically located far from major grain-producing regions, with electricity fed into the national grid through centralized, large-scale generation, making it difficult to serve agricultural production directly. In the case of wind power, wind farms are concentrated in coastal and plateau areas rich in wind resources, while major grain-producing regions lack such resources. The high installation and maintenance costs, as well as the technical threshold for wind power, make self-installation unfeasible in most cases. Therefore, both hydropower and wind power exert only indirect influence on Grain Productivity Quality Index, and their direct effects are not significant. The insignificant effects of hydropower and wind power may be attributable to their technological and spatial characteristics. These energy sources are typically deployed in geographically concentrated locations and depend heavily on large-scale grid integration, which may weaken their direct linkage with local agricultural production systems. In addition, hydropower and wind resources are often located in regions with limited arable land or lower grain-production intensity, reducing their capacity to directly support on-farm electrification and productivity-enhancing technologies. By contrast, solar energy can be deployed in a more decentralized manner and integrated with agricultural land use, while nuclear energy provides stable low-carbon baseload power that better supports electrification across production, processing, and logistics (<xref ref-type="bibr" rid="ref7">Hu, 2023</xref>).</p>
<p>In contrast, as shown in columns (3) and (4) of <xref ref-type="table" rid="tab7">Table 7</xref>, solar power (Sun) and nuclear power (Nuclear) exhibit significant positive effects on Grain Productivity Quality Index. Solar power provides stable electricity for agricultural production through photovoltaic generation, reducing energy costs. Particularly in remote areas, photovoltaic power is not dependent on the grid and can be directly applied to agricultural production, supplying energy for irrigation and electric farm machinery, thus realizing &#x201C;zero distance&#x201D; between energy production and consumption. Moreover, integrated agrivoltaic projects improve land use efficiency and extend the agricultural industry chain. Nuclear power, with its high efficiency, stability, and clean electricity generation, offers reliable power for grain processing and cold-chain logistics, reducing production interruptions and losses. Nuclear energy can also provide high-temperature process heat and clean electricity for zero-carbon hydrogen production, significantly reducing carbon emissions in the fertilizer industry. In addition, technologies such as nuclear irradiation breeding and soil improvement accelerate crop variety improvement, enhance soil quality, and improve pest and disease control, further strengthening the sustainability of grain production.</p>
</sec>
<sec id="sec22">
<label>4.6</label>
<title>Heterogeneity analysis</title>
<sec id="sec23">
<label>4.6.1</label>
<title>Spatial heterogeneity</title>
<sec id="sec24">
<label>4.6.1.1</label>
<title>Major grain-producing regions vs. non-major grain-producing regions</title>
<p>As shown in columns (1) and (2) of <xref ref-type="table" rid="tab8">Table 8</xref>, new energy development has a significant positive effect on Grain Productivity Quality Index in major grain-producing regions. This is mainly because resource endowments in these regions are highly aligned with the needs of agricultural production, and new energy technologies can be directly applied to critical processes such as irrigation and greenhouse cultivation, thus reducing production costs and improving efficiency. In addition, the state prioritizes food security in major grain-producing regions, with supporting policies facilitating the development of agricultural new energy projects. The well-developed agricultural infrastructure and industry chains in these regions create synergistic effects with new energy, further boosting Grain Productivity Quality Index. In contrast, in non-major grain-producing regions, arable land is scattered, posing challenges for new energy application, and low grain yields dampen investment incentives for new energy equipment, resulting in an insignificant effect on Grain Productivity Quality Index.</p>
<table-wrap position="float" id="tab8">
<label>Table 8</label>
<caption>
<p>Spatial heterogeneity.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
<th align="center" valign="top">(3)</th>
<th align="center" valign="top">(4)</th>
</tr>
<tr>
<th align="center" valign="top">Major grain producing areas</th>
<th align="center" valign="top">Non-food main producing areas</th>
<th align="center" valign="top">Eastern Region</th>
<th align="center" valign="top">Central and western regions</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">NE</td>
<td align="center" valign="middle">0.0050&#x002A;&#x002A;&#x002A; (3.3438)</td>
<td align="center" valign="middle">0.0471 (1.3298)</td>
<td align="center" valign="middle">0.0028&#x002A; (1.8205)</td>
<td align="center" valign="middle">&#x2212;0.0087 (&#x2212;0.8237)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">0.0027 (0.5565)</td>
<td align="center" valign="middle">0.1331&#x002A;&#x002A; (2.2897)</td>
<td align="center" valign="middle">0.3748&#x002A;&#x002A;&#x002A; (3.8860)</td>
<td align="center" valign="middle">0.0036 (0.9758)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">&#x2212;1.9466 (&#x2212;1.1334)</td>
<td align="center" valign="middle">&#x2212;0.6361 (&#x2212;1.4466)</td>
<td align="center" valign="middle">&#x2212;1.1094 (&#x2212;0.6495)</td>
<td align="center" valign="middle">&#x2212;1.1427&#x002A;&#x002A; (&#x2212;2.3322)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">&#x2212;0.2557&#x002A;&#x002A; (&#x2212;2.2787)</td>
<td align="center" valign="middle">&#x2212;0.1084&#x002A; (&#x2212;1.9459)</td>
<td align="center" valign="middle">&#x2212;0.2243&#x002A; (&#x2212;1.8362)</td>
<td align="center" valign="middle">&#x2212;0.1525&#x002A;&#x002A;&#x002A; (&#x2212;2.6890)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">&#x2212;0.0004 (&#x2212;0.1298)</td>
<td align="center" valign="middle">&#x2212;0.0020 (&#x2212;1.5887)</td>
<td align="center" valign="middle">&#x2212;0.0038&#x002A;&#x002A; (&#x2212;2.5779)</td>
<td align="center" valign="middle">&#x2212;0.0005 (&#x2212;0.2878)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.4620&#x002A;&#x002A;&#x002A; (11.0675)</td>
<td align="center" valign="middle">0.3300&#x002A;&#x002A;&#x002A; (13.2713)</td>
<td align="center" valign="middle">0.4663&#x002A;&#x002A;&#x002A; (10.9947)</td>
<td align="center" valign="middle">0.3712&#x002A;&#x002A;&#x002A; (13.6055)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">130</td>
<td align="center" valign="middle">170</td>
<td align="center" valign="middle">100</td>
<td align="center" valign="middle">200</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.931</td>
<td align="center" valign="middle">0.950</td>
<td align="center" valign="middle">0.963</td>
<td align="center" valign="middle">0.936</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec25">
<label>4.6.1.2</label>
<title>Eastern region vs. central and Western regions</title>
<p>As shown in columns (3) and (4) of <xref ref-type="table" rid="tab8">Table 8</xref>, the regression coefficient for GPQI is significant only in the eastern region, indicating that the positive effect of new energy development on Grain Productivity Quality Index is more pronounced in the east. This is largely attributable to the eastern region&#x2019;s strong economic foundation, advanced technology, and well-developed infrastructure, all of which provide favorable conditions for the application and promotion of new energy technologies. The high level of economic development in the east is accompanied by greater investment intensity in new energy. At the same time, substantial investment in agricultural technological innovation enables better integration of new energy technologies with agricultural production, enhancing the intelligence and efficiency of grain production. Furthermore, the high capacity of the power grid and digital infrastructure supports new energy development. In contrast, the central and western regions are characterized by lower levels of economic development, lagging infrastructure, traditional industrial structures, lower dependence on new energy, and weaker capacity for technological innovation, making it difficult to rapidly adapt to the application of new energy technologies. Additionally, some areas in the central and western regions are non-major grain-producing regions, and the scale effects of agricultural technological innovation and new energy application are difficult to realize, resulting in an insignificant effect of new energy development on Grain Productivity Quality Index.</p>
</sec>
</sec>
<sec id="sec26">
<label>4.6.2</label>
<title>Heterogeneity by level of policy support</title>
<p>As shown in <xref ref-type="table" rid="tab9">Table 9</xref>, under conditions of high fiscal support for agriculture, new energy development has a significant positive impact on Grain Productivity Quality Index. Sufficient financial investment enables the improvement of agricultural infrastructure, promotes technological innovation, and facilitates the widespread application of new energy technologies&#x2014;for example, through subsidies that encourage the use of solar-powered irrigation systems and wind-powered agricultural machinery, thereby enhancing production efficiency and grain yield. In contrast, under low levels of fiscal support, the lack of financial resources leads to weak agricultural infrastructure, making it difficult to introduce new energy technologies on a large scale; traditional production methods remain prevalent, and thus the effect of new energy on Grain Productivity Quality Index is not significant. Therefore, the level of fiscal support for agriculture directly affects the effectiveness of new energy technology adoption in the agricultural sector and thus determines its impact on Grain Productivity Quality Index. Therefore, the level of fiscal support for agriculture directly affects the effectiveness of new energy technology adoption in the agricultural sector and thus determines its impact on Grain Productivity Quality Index (<xref ref-type="bibr" rid="ref22">Zeng et al.,2018</xref>).</p>
<table-wrap position="float" id="tab9">
<label>Table 9</label>
<caption>
<p>Heterogeneity of policy support.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
</tr>
<tr>
<th align="center" valign="top">High fiscal support for agriculture</th>
<th align="center" valign="top">Low fiscal support for agriculture</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">NE</td>
<td align="center" valign="middle">0.0029&#x002A;&#x002A; (2.0623)</td>
<td align="center" valign="middle">&#x2212;0.0084 (&#x2212;0.7539)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">0.0036 (0.9363)</td>
<td align="center" valign="middle">0.2754&#x002A;&#x002A; (2.2124)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">&#x2212;1.3185&#x002A;&#x002A; (&#x2212;2.5566)</td>
<td align="center" valign="middle">1.0659 (0.5562)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">&#x2212;0.2068&#x002A;&#x002A;&#x002A; (&#x2212;2.8186)</td>
<td align="center" valign="middle">&#x2212;0.1430&#x002A; (&#x2212;1.8397)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">0.0008 (0.3959)</td>
<td align="center" valign="middle">&#x2212;0.0040&#x002A;&#x002A;&#x002A; (&#x2212;2.8725)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.3469&#x002A;&#x002A;&#x002A; (13.4153)</td>
<td align="center" valign="middle">0.4648&#x002A;&#x002A;&#x002A; (14.1443)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">170</td>
<td align="center" valign="middle">130</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.931</td>
<td align="center" valign="middle">0.956</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec27">
<label>4.6.3</label>
<title>Heterogeneity by economic and financial development levels</title>
<sec id="sec28">
<label>4.6.3.1</label>
<title>High marketization level vs. low marketization level</title>
<p>According to Columns (1) and (2) in <xref ref-type="table" rid="tab10">Table 10</xref>, the impact of new energy development on grain new quality productivity exhibits distinct regional heterogeneity. In regions with a high degree of marketization, agricultural production operates at a large scale and benefits from mature market mechanisms, which facilitate the effective integration of new energy technologies and enable their substitution for and complementarity with traditional energy sources. Supported by well-established financing and technology diffusion systems, new energy significantly drives the green transformation and efficiency improvement of agricultural production, thereby exerting a significant positive effect on grain new quality productivity.</p>
<table-wrap position="float" id="tab10">
<label>Table 10</label>
<caption>
<p>Heterogeneity of economic and financial development levels.</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Variables</th>
<th align="center" valign="top">(1)</th>
<th align="center" valign="top">(2)</th>
<th align="center" valign="top">(3)</th>
<th align="center" valign="top">(4)</th>
</tr>
<tr>
<th align="center" valign="top">High marketization level</th>
<th align="center" valign="top">Low marketization level</th>
<th align="center" valign="top">High level of financial development</th>
<th align="center" valign="top">Low level of financial development</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">NE</td>
<td align="center" valign="middle">0.0038&#x002A;&#x002A; (2.5092)</td>
<td align="center" valign="middle">0.0524 (1.3262)</td>
<td align="center" valign="middle">0.0502 (1.0475)</td>
<td align="center" valign="middle">0.0031&#x002A;&#x002A; (2.2149)</td>
</tr>
<tr>
<td align="left" valign="middle">Infor</td>
<td align="center" valign="middle">0.0015 (0.3782)</td>
<td align="center" valign="middle">0.0655 (0.7785)</td>
<td align="center" valign="middle">0.2280&#x002A;&#x002A;&#x002A; (3.3249)</td>
<td align="center" valign="middle">0.0036 (0.9764)</td>
</tr>
<tr>
<td align="left" valign="middle">Envir</td>
<td align="center" valign="middle">&#x2212;1.2873&#x002A;&#x002A; (&#x2212;2.6004)</td>
<td align="center" valign="middle">&#x2212;1.2093 (&#x2212;1.0660)</td>
<td align="center" valign="middle">&#x2212;1.0480&#x002A;&#x002A; (&#x2212;2.1374)</td>
<td align="center" valign="middle">0.6971 (0.5829)</td>
</tr>
<tr>
<td align="left" valign="middle">Indus</td>
<td align="center" valign="middle">&#x2212;0.0815&#x002A; (&#x2212;1.9302)</td>
<td align="center" valign="middle">&#x2212;0.3401&#x002A;&#x002A;&#x002A; (&#x2212;2.7378)</td>
<td align="center" valign="middle">&#x2212;0.1210&#x002A;&#x002A; (&#x2212;2.0082)</td>
<td align="center" valign="middle">&#x2212;0.1391&#x002A; (&#x2212;1.7150)</td>
</tr>
<tr>
<td align="left" valign="middle">Old</td>
<td align="center" valign="middle">0.0003 (0.2376)</td>
<td align="center" valign="middle">&#x2212;0.0031 (&#x2212;1.3823)</td>
<td align="center" valign="middle">&#x2212;0.0035&#x002A;&#x002A; (&#x2212;2.5831)</td>
<td align="center" valign="middle">&#x2212;0.0008 (&#x2212;0.4516)</td>
</tr>
<tr>
<td align="left" valign="middle">_cons</td>
<td align="center" valign="middle">0.3646&#x002A;&#x002A;&#x002A; (14.6357)</td>
<td align="center" valign="middle">0.4355&#x002A;&#x002A;&#x002A; (15.7041)</td>
<td align="center" valign="middle">0.3350&#x002A;&#x002A;&#x002A; (12.1546)</td>
<td align="center" valign="middle">0.4098&#x002A;&#x002A;&#x002A; (14.0559)</td>
</tr>
<tr>
<td align="left" valign="middle">Province FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">Year FE</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
<td align="center" valign="middle">YES</td>
</tr>
<tr>
<td align="left" valign="middle">N</td>
<td align="center" valign="middle">160</td>
<td align="center" valign="middle">140</td>
<td align="center" valign="middle">110</td>
<td align="center" valign="middle">190</td>
</tr>
<tr>
<td align="left" valign="middle">adj. <italic>R</italic><sup>2</sup></td>
<td align="center" valign="middle">0.969</td>
<td align="center" valign="middle">0.916</td>
<td align="center" valign="middle">0.898</td>
<td align="center" valign="middle">0.952</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p><italic>t</italic> statistics in parentheses. &#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.1, &#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.05, &#x002A;&#x002A;&#x002A;<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01.</p>
</table-wrap-foot>
</table-wrap>
<p>In contrast, in regions with a low degree of marketization, agricultural production remains dominated by small-scale farmers. Although new energy technologies could theoretically improve production conditions, their application is constrained by underdeveloped market mechanisms, weak infrastructure, and limited funding and technology extension. These barriers prevent new energy from achieving economies of scale, resulting in an insignificant effect on enhancing grain new quality productivity.</p>
</sec>
<sec id="sec29">
<label>4.6.3.2</label>
<title>High financial development level vs. low financial development level</title>
<p>As shown in columns (3) and (4) of <xref ref-type="table" rid="tab10">Table 10</xref>, regions with a high level of financial development have well-developed financial markets and diversified funding channels, making it easier for enterprises and farmers to access support for traditional energy and reducing dependence on new energy technologies. These regions also focus more on the innovation and upgrading of high value-added industries, with relatively less investment in agriculture, thus limiting the application and promotion of new energy technologies in the agricultural sector and resulting in an insignificant effect on Grain Productivity Quality Index.</p>
<p>In contrast, in regions with a low level of financial development, financial markets are underdeveloped and funding for traditional energy is insufficient, making new energy technologies an important alternative. Their application can effectively reduce production costs and improve energy efficiency&#x2014;for example, the promotion of solar-powered irrigation systems and small wind turbines can address energy shortages in remote areas, significantly enhancing Grain Productivity Quality Index.</p>
</sec>
</sec>
<sec id="sec30">
<label>4.6.4</label>
<title>Synthesis of provincial heterogeneity</title>
<p>Taken together, the substantial variation in effects across provinces highlights the highly context-dependent impact of new energy development on grain production quality. This heterogeneity originates from provincial-level structural differences in land institutions, developmental stages, and policy environments, which collectively moderate the influence on the GPQI. In provinces with robust land governance, advanced infrastructure, and stable agricultural systems, including major grain-producing and certain eastern provinces, the institutional framework enables resource reallocation and long-term technological investment. This facilitates channeling energy transition into quality improvements, an effect further amplified by agricultural fiscal support through the adoption of energy-related technologies. Ultimately, land institutional transformation shapes the capabilities and incentives of agricultural operators and the efficiency of policy implementation, providing a critical institutional backdrop for understanding the deep sources of regional divergence.</p>
</sec>
</sec>
</sec>
<sec id="sec31">
<label>5</label>
<title>Research conclusions and policy recommendations</title>
<sec id="sec32">
<label>5.1</label>
<title>Research conclusions</title>
<p>First, the findings demonstrate that new-energy development exerts a significant positive influence on the Grain Productivity Quality Index (GPQI). Quantitatively, a 1% rise in new-energy development corresponds to an average increase of 0.41% in GPQI (<italic>p</italic>&#x202F;&#x003C;&#x202F;0.01).</p>
<p>Second, this positive effect remains robust across multiple verification strategies, including alternative definitions of the explanatory variable, omission of municipality samples, and incorporation of additional control variables. The results also hold when addressing potential endogeneity using the clean-energy demonstration-province policy within a multi-period DID framework.</p>
<p>Third, the heterogeneity analysis by energy type shows that solar and nuclear meaningfully enhance GPQI, whereas hydropower and wind do not display statistically significant effects.</p>
<p>Fourth, additional heterogeneity analyses indicate that the influence of new-energy development varies across spatial, policy, economic, and financial contexts. In particular, a pronounced positive effect is found in major grain-producing areas, eastern regions, the latter half of the sample period (approximately 2018&#x2013;2022), provinces with strong fiscal support for agriculture, and regions characterized by higher levels of marketization and lower levels of financial development. These findings deepen the understanding of energy&#x2013;agriculture linkages by demonstrating how energy transition can enhance the quality dimension of grain productivity, and they provide empirical support for aligning China&#x2019;s dual-carbon goals with food security and agricultural modernization strategies.</p>
</sec>
<sec id="sec33">
<label>5.2</label>
<title>Policy recommendations</title>
<p>1. Coordinate region-specific new-energy deployment to foster energy&#x2013;agriculture synergy and raise GPQI.</p>
<p>Based on the empirical findings and heterogeneity analyses, the following policy recommendations are proposed to promote coordinated development between new-energy transition and grain productivity quality across regions. Regions should formulate scientific new-energy strategies aligned with local resource endowments, industrial bases, and ecological carrying capacities. In economically developed eastern provinces, it is recommended to accelerate the deployment of efficient clean energy sources, particularly solar and nuclear energy, thereby enhancing agricultural performance through energy-structure optimization. In central and western regions as well as resource based areas, wind and hydropower should be developed in accordance with local conditions, while the expansion of agrivoltaics and agroforestry and solar integration should be promoted to enhance land use efficiency and grain production capacity. In addition, local governments are encouraged to establish demonstration zones and innovation platforms for deep integration between the energy and agricultural sectors, thereby facilitating coordinated regional upgrading of both systems.</p>
<p>2. Strengthen green technological innovation in agriculture and promote adoption of clean-energy equipment to enhance GPQI&#x2019;s multidimensional content.</p>
<p>Increase public investment in green-agriculture R&#x0026;D and agricultural clean-energy equipment, with targeted support for breakthroughs in the electrification and digitalization of production. Encourage industry&#x2013;academia&#x2013;research collaboration to accelerate demonstration and diffusion of technologies such as solar-powered irrigation and biogas utilization. Improve technology-transfer and extension mechanisms to involve innovative firms and cooperatives in green-equipment development, thereby raising the technological content, sustainability, and resilience captured by GPQI.</p>
<p>3. Tighten land and ecological governance to ensure coordinated development of agriculture and new energy without compromising grain security.</p>
<p>While vigorously advancing new-energy projects, strictly protect arable land and the ecological environment to safeguard the baseline of grain security. Establish clear land-use standards and ecological-compensation rules for new-energy siting; promote eco-friendly compound models (e.g., agrivoltaics, agroforestry&#x2013;solar) and encourage multiple-use land arrangements. Strengthen dynamic supervision of facility layout and operations, and enhance ecological restoration and green-compensation efforts to prevent disorderly land competition and ecological degradation stemming from rapid new-energy expansion. This approach supports a virtuous interaction among energy transition, grain security, and ecological protection <xref ref-type="bibr" rid="ref15">Power (2010)</xref>.</p>
<p>4. Policy Recommendations with Provincial Differentiation.</p>
<p>The policy implications of this study vary across provinces with different development conditions. In major grain-producing regions and provinces with relatively strong land governance capacity, policies should prioritize the deep integration of new energy with agricultural production, such as promoting agrivoltaic systems, electrified irrigation, and energy-efficient processing facilities. In economically developed eastern provinces, where energy infrastructure is more advanced, efforts should focus on leveraging clean and stable energy sources to support agricultural upgrading along the entire value chain. By contrast, in less developed regions and provinces with weaker fiscal capacity, policy support should emphasize basic energy access, targeted subsidies, and institutional coordination to ensure that new-energy deployment effectively translates into improvements in grain productivity quality rather than exacerbating land-use conflicts.</p>
</sec>
<sec id="sec34">
<label>5.3</label>
<title>Future research direction</title>
<p>Despite its contributions, this study has several limitations that warrant further investigation. First, the analysis is conducted at the provincial level, which effectively captures regional aggregate characteristics but may mask intra-provincial heterogeneity in land-use patterns, agricultural practices, and new-energy deployment. Second, the sample period from 2013 to 2022 primarily reflects a specific stage of new-energy development and does not explicitly distinguish dynamic differences across development phases. Third, as the empirical analysis is based on macro-level data, the micro-level mechanisms through which new-energy development affects grain productivity quality remain to be further explored using finer-scale data. In addition, this study does not explicitly model potential spatial spillover effects or policy interaction mechanisms between energy transition and agricultural systems. Future research could address these limitations by employing more disaggregated data, extending the time horizon, and incorporating spatial or multi-policy analytical frameworks to deepen understanding of the coordinated evolution between new-energy development and grain production systems.</p>
</sec>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="sec35">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec sec-type="author-contributions" id="sec36">
<title>Author contributions</title>
<p>YL: Writing &#x2013; original draft, Software, Conceptualization, Methodology, Data curation. SM: Writing &#x2013; review &#x0026; editing. XZ: Supervision, Writing &#x2013; review &#x0026; editing, Project administration, Validation, Methodology. JC: Funding acquisition, Project administration, Validation, Writing &#x2013; review &#x0026; editing, Supervision, Resources.</p>
</sec>
<sec sec-type="COI-statement" id="sec37">
<title>Conflict of interest</title>
<p>The author(s) declared that this work 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="ai-statement" id="sec38">
<title>Generative AI statement</title>
<p>The author(s) declared that Generative AI was not used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="sec39">
<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="ref1"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chang</surname><given-names>K.</given-names></name> <name><surname>Wan</surname><given-names>Q.</given-names></name> <name><surname>Lou</surname><given-names>Q.</given-names></name> <name><surname>Chen</surname><given-names>Y.</given-names></name> <name><surname>Wang</surname><given-names>W.</given-names></name></person-group> (<year>2020</year>). <article-title>Green fiscal policy and firms&#x2019; investment efficiency: new insights into firm-level panel data from the renewable energy industry in China</article-title>. <source>Renew. Energy</source> <volume>151</volume>, <fpage>589</fpage>&#x2013;<lpage>597</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.renene.2019.11.064</pub-id></mixed-citation></ref>
<ref id="ref2"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chi</surname><given-names>N. T. K.</given-names></name></person-group> (<year>2022</year>). <article-title>Driving factors for green innovation in agricultural production: an empirical study in an emerging economy</article-title>. <source>J. Clean. Prod.</source> <volume>368</volume>:<fpage>132965</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.jclepro.2022.132965</pub-id></mixed-citation></ref>
<ref id="ref3"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coomes</surname><given-names>O. T.</given-names></name> <name><surname>Barham</surname><given-names>B. L.</given-names></name> <name><surname>MacDonald</surname><given-names>G. K.</given-names></name> <name><surname>Ramankutty</surname><given-names>N.</given-names></name> <name><surname>Chavas</surname><given-names>J. P.</given-names></name></person-group> (<year>2019</year>). <article-title>Leveraging total factor productivity growth for sustainable and resilient farming</article-title>. <source>Nat. Sustain.</source> <volume>2</volume>, <fpage>22</fpage>&#x2013;<lpage>28</lpage>. doi: <pub-id pub-id-type="doi">10.1038/s41893-018-0200-3</pub-id></mixed-citation></ref>
<ref id="ref4"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gyamfi</surname><given-names>B. A.</given-names></name> <name><surname>Agozie</surname><given-names>Q.</given-names></name> <name><surname>Bein</surname><given-names>M. A.</given-names></name> <name><surname>Bekun</surname><given-names>F. V.</given-names></name> <name><surname>Fatai Adedoyin</surname><given-names>F.</given-names></name></person-group> (<year>2022</year>). <article-title>Unlocking the investment impact of biomass energy utilization on environmental degradation for an isolated island</article-title>. <source>Int. J. Energy Sect. Manage.</source> <volume>16</volume>, <fpage>585</fpage>&#x2013;<lpage>604</lpage>. doi: <pub-id pub-id-type="doi">10.1108/IJESM-11-2020-0007</pub-id></mixed-citation></ref>
<ref id="ref5"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haberl</surname><given-names>H.</given-names></name></person-group> (<year>2015</year>). <article-title>Competition for land: a sociometabolic perspective</article-title>. <source>Ecol. Econ.</source> <volume>119</volume>, <fpage>424</fpage>&#x2013;<lpage>431</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ecolecon.2014.10.002</pub-id></mixed-citation></ref>
<ref id="ref6"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haberl</surname><given-names>H.</given-names></name> <name><surname>Erb</surname><given-names>K. H.</given-names></name> <name><surname>Krausmann</surname><given-names>F.</given-names></name> <name><surname>Bondeau</surname><given-names>A.</given-names></name> <name><surname>Lauk</surname><given-names>C.</given-names></name> <name><surname>M&#x00FC;ller</surname><given-names>C.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Global bioenergy potentials from agricultural land in 2050: sensitivity to climate change, diets and yields</article-title>. <source>Biomass Bioenergy</source> <volume>35</volume>, <fpage>4753</fpage>&#x2013;<lpage>4769</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.biombioe.2011.04.035</pub-id>, <pub-id pub-id-type="pmid">22211004</pub-id></mixed-citation></ref>
<ref id="ref7"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hu</surname><given-names>Z.</given-names></name></person-group> (<year>2023</year>). <article-title>Towards solar extractivism? A political ecology understanding of the solar energy and agriculture boom in rural China</article-title>. <source>Energy Res. Soc. Sci.</source> <volume>98</volume>:<fpage>102988</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.erss.2023.102988</pub-id></mixed-citation></ref>
<ref id="ref8"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname><given-names>X.</given-names></name> <name><surname>Feng</surname><given-names>C.</given-names></name> <name><surname>Qin</surname><given-names>J.</given-names></name> <name><surname>Wang</surname><given-names>X.</given-names></name> <name><surname>Zhang</surname><given-names>T.</given-names></name></person-group> (<year>2022</year>). <article-title>Measuring China's agricultural green total factor productivity and its drivers during 1998&#x2013;2019</article-title>. <source>Sci. Total Environ.</source> <volume>829</volume>:<fpage>154477</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.scitotenv.2022.154477</pub-id>, <pub-id pub-id-type="pmid">35304138</pub-id></mixed-citation></ref>
<ref id="ref9"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H.</given-names></name> <name><surname>Zhou</surname><given-names>X.</given-names></name> <name><surname>Tang</surname><given-names>M.</given-names></name> <name><surname>Guo</surname><given-names>L.</given-names></name></person-group> (<year>2022</year>). <article-title>Impact of population aging and renewable energy consumption on agricultural green total factor productivity in rural China: evidence from panel var approach</article-title>. <source>Agriculture</source> <volume>12</volume>:<fpage>715</fpage>. doi: <pub-id pub-id-type="doi">10.3390/agriculture12050715</pub-id></mixed-citation></ref>
<ref id="ref10"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ma</surname><given-names>S.</given-names></name> <name><surname>Yan</surname><given-names>H.</given-names></name> <name><surname>Li</surname><given-names>D.</given-names></name> <name><surname>Liu</surname><given-names>H.</given-names></name> <name><surname>Zeng</surname><given-names>H.</given-names></name></person-group> (<year>2025</year>). <article-title>The impact of agricultural mechanisation on agricultural carbon emission intensity: evidence from China</article-title>. <source>Pak. J. Agric. Sci.</source> <volume>62</volume>, <fpage>99</fpage>&#x2013;<lpage>110</lpage>. doi: <pub-id pub-id-type="doi">10.21162/PAKJAS/25.385</pub-id></mixed-citation></ref>
<ref id="ref12"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Omer</surname><given-names>A. A. A.</given-names></name> <name><surname>Liu</surname><given-names>W.</given-names></name> <name><surname>Li</surname><given-names>M.</given-names></name> <name><surname>Zheng</surname><given-names>J.</given-names></name> <name><surname>Zhang</surname><given-names>F.</given-names></name> <name><surname>Zhang</surname><given-names>X.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Water evaporation reduction by the agrivoltaic systems development</article-title>. <source>Sol. Energy</source> <volume>247</volume>, <fpage>13</fpage>&#x2013;<lpage>23</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.solener.2022.10.022</pub-id></mixed-citation></ref>
<ref id="ref13"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peng</surname><given-names>L.</given-names></name> <name><surname>Chen</surname><given-names>L.</given-names></name> <name><surname>Dai</surname><given-names>H.</given-names></name></person-group> (<year>2024</year>). <article-title>The impact of energy structure on agricultural green productivity in China</article-title>. <source>Sci. Rep.</source> <volume>14</volume>:<fpage>27938</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-024-78876-4</pub-id>, <pub-id pub-id-type="pmid">39537689</pub-id></mixed-citation></ref>
<ref id="ref14"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Popp</surname><given-names>J.</given-names></name> <name><surname>Lakner</surname><given-names>Z.</given-names></name> <name><surname>Harangi-R&#x00E1;kos</surname><given-names>M.</given-names></name> <name><surname>Fari</surname><given-names>M.</given-names></name></person-group> (<year>2014</year>). <article-title>The effect of bioenergy expansion: food, energy, and environment</article-title>. <source>Renew. Sust. Energ. Rev.</source> <volume>32</volume>, <fpage>559</fpage>&#x2013;<lpage>578</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.rser.2014.01.056</pub-id></mixed-citation></ref>
<ref id="ref15"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Power</surname><given-names>A. G.</given-names></name></person-group> (<year>2010</year>). <article-title>Ecosystem services and agriculture: tradeoffs and synergies</article-title>. <source>Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci.</source> <volume>365</volume>, <fpage>2959</fpage>&#x2013;<lpage>2971</lpage>.</mixed-citation></ref>
<ref id="ref16"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Shi</surname><given-names>H.</given-names></name> <name><surname>Xu</surname><given-names>H.</given-names></name> <name><surname>Gao</surname><given-names>W.</given-names></name> <name><surname>Zhang</surname><given-names>J.</given-names></name> <name><surname>Chang</surname><given-names>M.</given-names></name></person-group> (<year>2022</year>). <article-title>The impact of energy poverty on agricultural productivity: the case of China</article-title>. <source>Energy Policy</source> <volume>167</volume>:<fpage>113020</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.enpol.2022.113020</pub-id></mixed-citation></ref>
<ref id="ref17"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>Y.</given-names></name> <name><surname>Qian</surname><given-names>Y.</given-names></name></person-group> (<year>2024</year>). <article-title>Driving factors to agriculture total factor productivity and its contribution to just energy transition</article-title>. <source>Environ. Impact Assess. Rev.</source> <volume>105</volume>:<fpage>107369</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.eiar.2023.107369</pub-id></mixed-citation></ref>
<ref id="ref18"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xiong</surname><given-names>C.</given-names></name> <name><surname>He</surname><given-names>X.</given-names></name></person-group> (<year>2025</year>). <article-title>Construction of the evaluation index system, level measurement, and regional differences of China&#x2019;s new quality agricultural productivity</article-title>. <source>Sci. Tech. Manag. Res.</source> <volume>6</volume>, <fpage>53</fpage>&#x2013;<lpage>62</lpage>.</mixed-citation></ref>
<ref id="ref19"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xu</surname><given-names>G.</given-names></name> <name><surname>Yang</surname><given-names>M.</given-names></name> <name><surname>Li</surname><given-names>S.</given-names></name> <name><surname>Jiang</surname><given-names>M.</given-names></name> <name><surname>Rehman</surname><given-names>H.</given-names></name></person-group> (<year>2024</year>). <article-title>Evaluating the effect of renewable energy investment on renewable energy development in China with panel threshold model</article-title>. <source>Energy Policy</source> <volume>187</volume>:<fpage>114029</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.enpol.2024.114029</pub-id></mixed-citation></ref>
<ref id="ref21"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ying</surname><given-names>J.</given-names></name> <name><surname>Su</surname><given-names>H.</given-names></name> <name><surname>He</surname><given-names>S.</given-names></name> <name><surname>Qiu</surname><given-names>G.</given-names></name> <name><surname>Chen</surname><given-names>X.</given-names></name></person-group> (<year>2025</year>). <article-title>Belief dispersion in credit markets: evidence from CDS-bond basis</article-title>. <source>Financ. Res. Lett.</source> <volume>86</volume>:<fpage>108076</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.frl.2025.108076</pub-id></mixed-citation></ref>
<ref id="ref22"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zeng</surname><given-names>S.</given-names></name> <name><surname>Jiang</surname><given-names>C.</given-names></name> <name><surname>Ma</surname><given-names>C.</given-names></name> <name><surname>Su</surname><given-names>B.</given-names></name></person-group> (<year>2018</year>). <article-title>Investment efficiency of the new energy industry in China</article-title>. <source>Energy Econ.</source> <volume>70</volume>, <fpage>536</fpage>&#x2013;<lpage>544</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.eneco.2017.12.023</pub-id></mixed-citation></ref>
<ref id="ref23"><mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>X.</given-names></name> <name><surname>Li</surname><given-names>G.</given-names></name> <name><surname>Wu</surname><given-names>R.</given-names></name> <name><surname>Zeng</surname><given-names>H.</given-names></name> <name><surname>Ma</surname><given-names>S.</given-names></name></person-group> (<year>2025</year>). <article-title>Impact of carbon emissions, green energy, artificial intelligence and high-tech policy uncertainty on China&#x2019;s financial market</article-title>. <source>Financ. Res. Lett.</source> <volume>82</volume>:<fpage>107599</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.frl.2025.107599</pub-id></mixed-citation></ref>
</ref-list>
<fn-group>
<fn fn-type="custom" custom-type="edited-by" id="fn0001">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1843811/overview">Taiyi He</ext-link>, Southwestern University of Finance and Economics, China</p>
</fn>
<fn fn-type="custom" custom-type="reviewed-by" id="fn0002">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2974169/overview">Shenghao Bi</ext-link>, Beijing Normal University, China</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3315348/overview">Zubaria Andlib</ext-link>, Federal Urdu University of Arts, Sciences and Technology Islamabad, Pakistan</p>
</fn>
</fn-group>
</back>
</article>