Jukuri, open repository of the Natural Resources Institute Finland (Luke) All material supplied via Jukuri is protected by copyright and other intellectual property rights. Duplication or sale, in electronic or print form, of any part of the repository collections is prohibited. Making electronic or print copies of the material is permitted only for your own personal use or for educational purposes. For other purposes, this article may be used in accordance with the publisher’s terms. There may be differences between this version and the publisher’s version. You are advised to cite the publisher’s version. This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Author(s): Albert Muleke, Matthew Tom Harrison, Rowan Eisner et al. Title: Clarifying confusions over carbon conclusions: antecedent soil carbon drives gains realised following intervention Year: 2023 Version: Publisher’s version Copyright: The author(s) 2023 Rights: CC BY 4.0 Rights url: https://creativecommons.org/licenses/by/4.0/ Please cite the original version: Albert Muleke, Matthew Tom Harrison, Rowan Eisner et al. (2023) Clarifying confusions over carbon conclusions: antecedent soil carbon drives gains realised following intervention, Global Environmental Change Advances, DOI: 10.1016/j.gecadv.2023.100001 Global Environmental Change Advances 1 (2023) 100001 Available online 20 August 2023 2950-1385/© 2023 Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Clarifying confusions over carbon conclusions: antecedent soil carbon drives gains realised following intervention Albert Muleke a, Matthew Tom Harrison a,*, Rowan Eisner a, Maria Yanotti p, Peter de Voil b, Shah Fahad c, Wang Fei d, Puyu Feng e, Carla Ferreira f,g, Daniel Forster h, Xionghui Gao i, Ke Liu a, Jianguo Man j, Lixiao Nie k, Jiangwen Nie l, Zhiming Qi m, Narasinha Shurpali h, Weilu Wang n, Rui Yang o, Xiaogang Yin l, Feng Zhang q, Yunbo Zhang r, Jin Zhao s a Tasmanian Institute of Agriculture, University of Tasmania, Newnham Drive, Launceston 7248, Australia b Queensland Alliance for Agriculture and Food Innovation (QAAFI), The University of Queensland, Gatton Campus, Gatton, QLD 4343, Australia c Department of Agronomy, Abdul Wali Khan University Mardan, Khyber Pakhtunkhwa 23200, Pakistan d College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, Hubei, China e College of Land Science and Technology, China Agricultural University, Beijing 100193, China f Department of Physical Geography and Bolin Centre for Climate Research, Stockholm University, Stockholm, Sweden g Research Centre for Natural Resources, Environment and Society (CERNAS), Polytechnic Institute of Coimbra, Coimbra Agrarian Technical School, Coimbra, Portugal h Grasslands and Sustainable Farming, Production systems, Natural Resources Institute Finland (Luke), Halolantie 31 A, 71750 Maaninka, Finland i Chinese Academy of Agriculture Science, Beijing, China j College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070, China k Sanya Nanfan Research Institute of Hainan University, Hainan University, Sanya 572025, China l College of Agronomy and Biotechnology, China Agricultural University, Beijing 100193, China m Department of Bioresource Engineering and Brace Centre for Water Resources Management, MS1–024 Macdonald Campus, McGill University, Sainte-Anne-de- Bellevue, QC H9X 3V9 Canada n The Institutes of Agricultural Science and Technology Development, Yangzhou University, Wenhui East Road, Yangzhou City, China o Hubei Collaborative Innovation Centre for Grain Industry/Agriculture college, Yangtze University, Jingzhou, China p Tasmanian School of Business and Economics, University of Tasmania, Private Bag 98, Hobart, TAS 7001, Australia q State Key Laboratory of Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, China r Hubei Collaborative Innovation Centre for Grain Industry/Agriculture College, Yangtze University, Jingzhou, China s College of Resources and Environmental Sciences, China Agricultural University, 100193 Beijing, China A R T I C L E I N F O Keywords: Carbon Net-zero Emissions intensity Soil Mitigation Adaptation Climate crisis Crop Drought Productivity Offsets Accrual Mineralisation Fertiliser A B S T R A C T Carbon removals associated with incremental gains in soil organic carbon (SOC) at scale have enormous po tential to mitigate global warming, yet confusion over contexts that elicit SOC accrual abound. Here, we examine how bespoke interventions (through irrigation, fertiliser, crop type and rotations), antecedent SOC levels and soil type impact on long-term SOC accrual and greenhouse gas (GHG) emissions. Using a whole farm systems modelling approach informed using participatory research, we discovered an inverse relationship between antecedent SOC stocks and SOC gains realised following intervention, with greater initial SOC levels resulting in lower ex poste change in SOC. We found that SOC accrual was greatest for clays and least for sands, although changes in SOC in sandy loam soils were also low. Diversified whole farm adaptations – implemented through inclusion of grain legumes within wheat/canola crop rotations – were more conducive to improvement in SOC stocks, followed by Intensified systems (implemented through greater rates of irrigation, farm areas under irri gation, nitrogen fertiliser and inclusion of rice and maize in crop rotations). Adaptations that Simplified farm systems by reducing irrigation and fertiliser use resulted in the lowest SOC accrual. In most cases, long-term SOC stocks fell when SOC at the outset was greater than 4–5%, regardless of intervention made, soil or crop type, crop rotation, production system or climate. We contend that (1) management interventions primarily impacted SOC in the soil surface (0–30 cm) and had de minimus impact on deep SOC stocks (30–100 cm), (2) crop rotations including wheat, canola and faba beans were more conducive to improvement in SOC stocks, (3) scenarios with high status quo SOC had little impact on crop productivity, and not necessarily the lowest GHG emissions * Corresponding author. E-mail address: matthew.harrison@utas.edu.au (M.T. Harrison). Contents lists available at ScienceDirect Global Environmental Change Advances journal homepage: www.journals.elsevier.com/global-environmental-change-advances https://doi.org/10.1016/j.gecadv.2023.100001 Received 15 May 2023; Received in revised form 9 July 2023; Accepted 30 July 2023 mailto:matthew.harrison@utas.edu.au www.sciencedirect.com/science/journal/29501385 https://www.journals.elsevier.com/global-environmental-change-advances https://doi.org/10.1016/j.gecadv.2023.100001 https://doi.org/10.1016/j.gecadv.2023.100001 https://doi.org/10.1016/j.gecadv.2023.100001 http://crossmark.crossref.org/dialog/?doi=10.1016/j.gecadv.2023.100001&domain=pdf http://creativecommons.org/licenses/by/4.0/ Global Environmental Change Advances 1 (2023) 100001 2 intensity, (4) productivity and GHG emissions intensity were largely a function of the quantum of nitrogenous fertiliser added, rather than SOC stocks, and (5) aspirations for improving SOC are likely to be futile if antecedent SOC stocks are already high (4–5 %). We conclude that potential for improving SOC stocks exists in contexts where antecedent stocks are low (<1%), which may include regions with land degradation, chronic erosion and/ or other constraints to vegetative ground cover that could be sustainably and consistently alleviated. 1. Introduction The land use sector faces increasing pressure from society to sus tainably intensify food production and improve environmental stew ardship on the one hand, yet also to reduce greenhouse gas (GHG) emissions and engender carbon removals on the other (Ara et al., 2021; Ibrahim et al., 2019; Sándor et al., 2020). Globally, agriculture accounts for around 18.4 % of the GHG inventory, with 5.5 % from croplands and agricultural soils (Balogh, 2022; Ritchie et al., 2020). Ceteris paribus, GHG emissions from the agricultural sector are expected to rise in par allel with expansion of intensified agricultural land use (Hunter et al., 2017; Taylor et al., 2016). While irrigated croplands offer unique op portunities to modulate both agricultural GHGs and soil carbon sinks (McGill et al., 2018; Sapkota et al., 2020), our knowledge of potential GHG mitigation through SOC sequestration associated with irrigated systems is very much in its infancy (McGill et al., 2018). Irrigated sys tems conceivably provide necessary impetus for greater production and presumably SOC sequestration, yet greater production is often driven by higher use of nitrogenous fertilisers (Adnan et al., 2022; Rawnsley et al., 2019), which in turn may increase nitrous oxide (N2O) emissions, causing net GHG emissions to overcompensate for mitigation realised by enhancement of SOC stocks (Bilotto et al., 2021; Christie et al., 2014; Henry et al., 2022). Indeed, recent work has suggested that greater SOC stocks may not ultimately result in lower whole farmnet GHG emissions due to the tight coupling between productivity and GHG emissions (Bilotto et al., 2023; Harrison et al., 2021), since greater soil organic matter and SOC may be conducive to greater crop or pasture production, thereby potentially allowing nitrous oxide and/or methane emissions to rise in parallel with SOC accrual. Contemporary scholarly research of sustainable intensification via irrigation has primarily focused adopted a productivity and/or eco nomic lens (Ash et al., 2017; Monjardino et al., 2022; Muleke et al., 2022a,b; Muleke et al., 2023; Muleke et al., 2022c). For example, Muleke et al. (2023) concluded that productivity and profitability could be sustainably raised through integration of grain legume crops into existing rotations, while Muleke et al. (2022b) show that the number of profitable crop types fell from 35 to 10 under future climates, reflecting an interplay between commodity price, yield, crop water requirements, water costs and variable costs, wherein water use and costs increased, while crop yields decreased under future climates. In their study, effects of climate change on profit were not related to long-term rainfall, with future climates depressing profit by 11–23 % relative to historical cli mates, but rather impacts of future climates were more closely related to crop type and maturity duration, with many crop types that were traditionally profitable under historical climates being no longer prof itable in future. Such studies demonstrate that productivity (e.g., crop yield) and profitability have often been the motivation of previous work, primarily to the exclusion of associated or interacting effects on soils, including implications of systems interventions on soil health, carbon and greenhouse gas (GHG) emissions. Confusion among the scientific and practitioner community abounds as to the contexts with which SOC increases following intervention. Previous research has shown that weather conditions (largely precipi tation), management and soil texture drive SOC sequestration (Emde et al., 2021; Henry et al., 2022; Singhal et al., 2023), primarily through modulating influence on pasture or crop productivity (Langworthy et al., 2018; Yan et al., 2022). Conversely, other field studies have concluded that tillage and residue management practices manifest greater control over SOC sequestration (Chao, 2019; Karstens et al., 2022; Saurabh et al., 2021). Others argue that changes in SOC seques tration are contingent upon type of crop, genotype, and crop rotation (Antón et al., 2022; Sándor et al., 2020; Yang et al., 2023). Such treat ment effects are further confounded by regional differences in global warming, with soil carbon sequestration expected to fall depending on regional specificity and extent of the climate crisis (Bilotto et al., 2023) and frequency of extreme events (Harrison, 2021; Langworthy et al., 2018; Liu et al., 2023). Within such studies, particularised conclusions may well be robust (e.g. Liu et al., 2020), but valid interpretations while bespoke can be voided when extrapolated to new contexts, agroeco logical regions and/or production systems. We posit that much of this confusion can be attributed to initial (antecedent) SOC stocks: experiments beginning from SOC levels close to their floor (e.g., less than 1 %) could be reasonably expected to improve SOC, whereas antecedent SOC at ceiling levels (e.g., >6 %) might be expected to demonstrate less gain, because stocks are already close to their maximum levels. While scientists and practitioners often focus on within-experiment activities to assess the relative promise associated with an intervention – such as treatments and climatic con ditions – they tend to downplay or even completely overlook the influ ence of initial SOC on their results. Experimental results may be further confounded by limited duration of field experiments, which canconfine interpretations drawn to parochial weather conditions realised during the finite window within which the experiment is conducted. Trost et al. (2013) note that it is difficult to estimate the effects of factorial com binations caused by the environment and management on SOC seques tration and, as such, the present study aims to address this knowledge gap. Here, we build on foundational insights of previous work (Muleke et al., 2023) to investigate how initial SOC stocks influence long-term trajectories of SOC accrual using a real irrigated grains farm as a use case. We invoked three whole-farm adaptations and compared direc tional change with the existing system (Baseline) to examine how man agement (crop type, rotation, irrigation, fertilisation) and soil type (sand, clay, sandy loam) impacted on SOC over the long-term. Specif ically, our aims were to (1) examine the impacts of antecedent SOC on sequestration for varying whole-farm adaptations and soil types, and (2), quantify effects of these adaptations and soil types on net farm GHG emissions. 2. Methods 2.1. Study overview An irrigated broadacre farm situated in the Riverina region of New South Wales (NSW), Australia (− 34.8016̊S, 145.8904̊E) was used as a case study. A detailed description of the case study, including climate, soils, farm layout, crop rotation, irrigation methods and source of irri gation water supply is provided in Muleke et al. (2023) and Monjardino et al. (2022). In brief, our study was conducted by initialising our bio physical model (APSIM) according to the cropping production system operating on the case study farm. This step was made to ensure that the simulations were realistic and regionally representative. This included initialisation of farm management (fertiliser and irrigation use, crop rotation etc) and soil physical characteristics. Soil and climatic details were derived from publicly available scientific databases (Sections 2.2 and 2.3). Soil data were used to initialise APSIM, while daily climate A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 3 data drives APSIM (daily computations of crop growth and development are driven by climate inputs, similar to other contemporary cropping systems models). We next developed three generic adaptations scenarios in consulta tion with a regional reference group of experts comprising farmers, farmer groups, consultants and extension staff with the aim of improving SOC stocks. These adaptations comprised Simplification (reducing in puts), Diversification (increasing crop type and genotype diversity in rotation) and Intensification through use of greater inputs and crops requiring irrigation (detailed in Section 2.3.3). Because net greenhouse gas (GHG) emissions depend on the balance of nitrous oxide (N2O) and SOC sequestration, we computed net emissions asociated with each adaptation. Adaptations were compared with the baseline (status quo operation) to assess direction and extent of change. Each adaptation was made in the model by changing farm management (sowing time, irri gation, fertiliser, crop rotation) as shown in Table 1. Thus, our primary focus in this study was on farming systems adaptation and on change relative to the baseline (fractional and percentage), rather than on ab solute values (e.g., actual tonnes SOC/ha). 2.2. Historical climate data We conducted simulations using 110 years of historical climatic data (1 January 1910–31 December 2019). Data for daily maximum and minimum temperature, rainfall and solar insolation were sourced from SILO (Scientific Information for Land Owners) Patched Point Dataset (Jeffrey et al., 2001) (shown in Fig. S1). Historical annual rainfall of the case study location was 403 mm, with 218 mm in winter and 185 mm in summer (Jeffrey et al. (2001). These data have been rigorously screened for anomalies by the data provider using a two-pass interpolation sys tem; erroneous values were excluded from the interpolation, so they did not adversely affect the gridded surface (https://www.longpaddock.qld. gov.au/silo/about/about-data/). 2.3. Simulation of crop growth, yield and soil carbon stocks We adopted the farming systems outlined in Muleke et al. (2023) to examine the impact of initial SOC and whole farm adaptation on SOC accrual and nitrous oxide (S-N₂O) emissions using the Agricultural Production Systems sIMulator (APSIM) v7.10, (Holzworth et al., 2014; Keating et al., 2003). APSIM has been comprehensively validated under a wide range of conditions, with such studies showing that the model performs well in simulating crop growth, soil C dynamics and soil ni trogen (N) cycling (Bilotto et al., 2021; Harrison et al., 2019; Holzworth et al., 2014; Keating et al., 2003; Luo et al., 2011; Mohanty et al., 2020; O’Leary et al., 2016; Wang et al., 2003a). 2.3.1. Surface and soil organic matter and carbon Daily climate data are used in APSIM to compute biophysical values, including plant growth, labile soil carbon inputs, mineralisation, plant senescence and soil water infiltration and extraction. Management in the model further influences growth and carbon inputs, with irrigation improving growth during periods of plant water deficit, and nitrogen improving growth during periods of nitrogen stress. Plant growth, soil inputs and soil carbon in the model are derived by rate-state equations: the rate of any given biophysical process depends on the state of that variable on the prior day. For example, growth is greater when leaf area index on the prior day is relatively high. APSIM simulates surface and soil organic matter, of which soil organic carbon is subset. When surface organic matter is added (e.g., from wheat straw at harvest), that residue will be added to the existing pool of surface organic matter. For all simulations, we assumed an initial surface organic matter pool of 1000 kg DM/ha. Each organic matter pool has separate C:N ratios, decom position rates and specific areas. Leaching of soluble inorganic N (nitrate and ammonium) is determined by rates of rainfall or irrigation. Soil organic matter is partitioned into fresh (FOM), medium (BIOM) and a more stable humus (HUM) decomposition pool, with fractionation be tween these pools dependent on the type of organic matter and crop residues. A subset of the HUM pool is defined as inert. Residue decom position is dynamically determined by soil moisture, temperature, C:N ratio of the substrate, whether the residue is standing or prostrate, and the amount of material contacting the soil surface. Decomposition re sults in some loss of carbon in surface organic matter as CO2 and transfer of C and N to the soil; material with higher C:N ratios creates an immobilisation demand, which is satisfied from mineral N in surface layers; inadequate mineral N can restrict decomposition of residues in surface layers. For detailed accounts of surface and soil organic matter models in APSIM, readers are referred to Dimes and Revanuru (2004) and Probert et al. (2005). Given that soil water holding capacity and growth differs across soil types, we adopted three regionally representative soil types: brown vertosol clay, sand and sandy loam, following soil classification nomenclature of Isbell (2021). We used default hydraulic attributes for saturated water flow in APSoil (Table S2; Dalgliesh et al., 2016). To examine effects of initial SOC stocks, we sourced measured data for bulk density and soil organic matter pools (fresh organic matter pool (FOM), microbial biomass pool (BIOM) and humic pool (HUM); Probert et al., 1998) from the APSoil database (Dalgliesh et al., 2012) (Table S1), then initialised APSIM with six initial SOC concentrations (0.5 %, 2 %, 3 %, 4 %, 5 % and 6 %) within the top 100 cm of soil profile, following Dal gliesh et al. (2016). SOM and SOC in the APSoil database were deter mined by Dalgliesh et al. (2016) using Walkley Black or Leco analytical methods on soil cores taken in the Coleambally region of the case study farm; historical management associated with these data were assumed to be similar to status quo operations on the case study farm (default parameters shown in Table S1). Neither bulk density nor management were altered within any 110-year simulation (except for crop rotation specified in Section 2.3.3). In line with our underpinning aim – to examine the effect of initial SOC on subsequent change in the simulation – and to account for SOC variability due to historical climate and management types, we specified a range of SOC values that were used to initialise APSIM. Ensuing changes in simulated SOC over the 110 years were an emergent property, driven by climate, and modulated by in teractions between antecedent SOC, soil type and management over the course of the simulation. 2.3.2. Management and cropping systems We used various APSIM crop modules to simulate growth and development of wheat (Brown et al., 2014; Wang et al., 2003b), canola (Robertson et al., 1999), faba bean (Turpin et al., 2003), maize (Harrison et al., 2014), rice (Bouman et al., 2001) and cotton (Hearn, 1994). Irrigated crops were sown on fixed dates (mid-May for winter crops, mid-November for summer crops); dryland crops were sown when Table 1 Details of the Baseline and whole-farm adaptations modified after Muleke et al. (2023) and Muleke et al. (2022a). Whole-farm adaptation Crop Sowing date Irrigated area (ha) N applied (kg N/ha/yr) Baseline Canola 17-May 750 100 Wheat 7-Jun 750 150 Diversified Canola 17-May 750 100 Wheat 7-Jun 750 150 Faba bean 29-Mar 750 50 Intensified Canola 17-May 750 150 Maize 29-Dec 188 400 Wheat 7-Jun 750 250 Cotton 1-Oct 188 300 Rice 15-Nov 188 500 Simplified Canola (Dry) 10-May 750 100 Wheat 7-Jun 375 150 A. Muleke et al. https://www.longpaddock.qld.gov.au/silo/about/about-data/ https://www.longpaddock.qld.gov.au/silo/about/about-data/ Global Environmental Change Advances 1 (2023) 100001 4 sufficient autumnal rainfall occurred (25 mm over four days). Sowing dates were set according to information outlined in Muleke et al. (2022a), shown in Table 1. Default soil parameters were adopted from the APSIM soil library following Dalgliesh et al. (2012). Plant available soil water at sowing and application of irrigation and nitrogen (N) fer tiliser for each adaptation was set following Muleke et al. (2023) (details provided in Section 2.3.3, Table 2 and Fig. S2). Maximum potential long-term yields were estimated by identifying optimal flowering pe riods (OFPs), defined as the window that minimises long-term risk of abiotic stress exposure. OFPs for irrigated and dryland crops were computed as the flowering dates corresponding to ⩾95 % of the maximum 15-day running average frost-heat yield following Muleke et al. (2022a), Liu et al. (2020) and Liu et al. (2021). 2.3.3. Whole farm systems adaptations Three whole farm systems adaptations suggested by regional experts were compared with the existing farm system (the Baseline), which comprised 750 ha of irrigated canola− wheat− wheat in serial rotation, with each winter crop followed by a summer fallow (Table 2, Fig. S1). Summer fallows stored water and nitrogen, incurring weed control costs (0–4 herbicide spray events). Thematic adaptations included Intensifi cation (including higher rates of irrigation and fertiliser, greater farm area under irrigation, and inclusion of rice and maize in the rotation); Diversification (including a grain legume to diversify the crop rotation) and Simplification, being the opposite of Intensification. Further details are provided in Tables 1, 2 and Fig. S1 and previous studies (Monjardino et al., 2022; Muleke et al., 2023). Diversified scenarios included a canola–wheat− faba bean rotation, with each winter crop followed by a summer fallow (750 ha). Relative to the Baseline, Simplified scenarios were allocated lower irrigation rates and irrigated farm area (only 50 % of wheat area; 375 ha); dryland canola was sown on 750 ha of the farm area (Table 2 and Fig. S1). Intensified scenarios were irrigated at higher rates (ML/ha), assuming 750 ha for the winter crops (canola, wheat) and 25 % of the farm area (188 ha) for summer crops (maize, cotton, rice), with the remaining 75 % of farm area fallowed in summer. Nitrogen fertiliser rates for each crop type in the Baseline scenario were obtained from the case study region. We then applied these fertiliser rates for corresponding crops in the Diversified and Simplified scenarios. We assumed 50–100 kg N/ha greater application rates for each crop in the Intensified scenario, as our aim for this adaptation was to examine the impact of intensification through greater application of N and irrigation (Table 1). 2.3.4. Greenhouse gas (GHG) emissions and global warming potential (GWP) 2.3.4.1. Soil emissions. To quantify direct net farm soil GHG, we con verted S-N₂O and the change in SOC over the simulation duration (ΔSOC) to carbon dioxide equivalents (CO₂-e). For S-N₂O, we used a 100 year GWP of 298 (Eq.1; IPCC, 2013): Soil N₂O emissions(tonnes CO₂­e) = N₂O(Layer) × 298 1000 (1) Where Soil N₂O emissions represent tonnes of carbon dioxide equivalents (tonnes CO₂-e) from simulated S-N₂O and N₂O(Layer) represents simulated S-N₂O in uppermost (0–30 cm) and deeper (30–100 cm) soil layers. SOC stocks (0–30 and 30–100 cm) were converted into CO2-e by multiplying SOC by 3.67 using Eq. 2 (IPCC, 2013): ΔSOC emissions(tonnes CO₂­e) = SOC(Layer) × 3.67 (2) Where ΔSOC emissions represent tonnes of carbon dioxide equivalents from ΔSOC (difference between SOC at start and end of 110-year period), while SOC(Layer) represents simulated SOC in each of the upper (0–30 cm) and deeper (30–100 cm) soil layers in tonnes C per hectare. Net farm GWP across whole-farm adaptations, soil types and initial SOC levels was calculated as the sum of S-N₂O and ΔSOC using Eq. 3, with positive gains in ΔSOC resulting in negative GHG emissions. Net onfarm GWP(tonnes CO₂e)=Soil N₂O emissions+ΔSOC emissions (3) 2.3.4.2. Emissions intensities of grain yield. Emissions intensity (carbon footprints or CF) were computed as the annual per hectare direct emissions in kilograms of CO₂-e per kilogram of crop yield using equa tion Eq.4 (FAO, 2022). EI(kg CO₂­e / kg crop yield)= Annual GHG emissions(kg CO₂­e/ha) yield(kg/ha) (4) Where EI (kg CO₂-e/kg crop yield) represents emissions intensity, annual GHG emissions represents average direct GHG emissions in kilogram carbon dioxide equivalents (kg CO₂-e) and yield represents average simulated harvested crop yield (kg/ha). 2.3.5. Statistical analysis Regression analyses were performed to test statistical significance of the impact of initial SOC, whole farm adaptation and soil type on rela tive change in long-term SOC storage. We fitted linear models using R 4.2.2 (R Core Team, 2022) to predict change in SOC, using initial SOC, soil type and adaptation as predictor variables. We adopted confidence intervals (CI) of 95 %, a p-value of < 0.05 to denote significance, and adjusted R-squared (R2) of 0.99 (refer to Table S3). We applied omnibus “analysis of variance” (ANOVA) to the linear model to test the statistical significance between initial SOC, soil type, adaptations and their inter action, while Tukey’s post-hoc test was employed to identify disparities between paired level combinations (Tables S4-6). The relative variation contributed by each variable was computed using Eq. 5. Variation contributed = SS associated with the variable Total SS × 100 (5) Where SS represents the sum of squares for each factor (initial SOC, soil type, adaptations and interactions therein) from ANOVA (Table 5); Total SS represents total sum of squares for all dependent factors and their interactions in the regression model. Table 2 Description of the Baseline and whole farm adaptations, after Muleke et al. (2023). Adaptation Description Baseline The current farm system with gravity irrigation. The Baseline represents the current farm situation in terms of agronomy and irrigation. Diversified Diversification of the Baseline system using a grain legume (Faba beans; Vicia faba). This adaptation has similar average water usage and irrigated farm area as the Baseline, but a greater variety of winter crops grown within the year − i.e., wheat, canola and a grain legume – with similar inputs and costs. Intensified This was designed to be a high-input, high-output adaptation. Relative to the Baseline, this scenario applies higher amounts of water per unit area and per year and assumes a larger portion of the farm area is irrigated (i.e., less unirrigated fallow), and has higher inputs per unit area and year. Simplified This adaptation was designed as a low-input, low-output adaptation. Relative to the Baseline, the Simplified scenario was designed to require less irrigation water per unit area and year, has more rainfed crops, and reduced inputs per unit area. A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 5 3. Results 3.1. Antecedent SOC stocks determine SOC gain realised following intervention Across soil types and adaptations, greater initial SOC stocks corre lated inversely with long-term SOC change in the surface 30 cm (Table 3 and Fig. 1). In the deeper soil layers (30–100 cm), rate of change in SOC showed little association with initial SOC stocks (Table 3 and Fig. 3). Across all climate years (110 years), soil types and adaptations, a one unit increase in initial SOC within the upper 30 cm decreased temporal SOC accrual by 0.49 % (10 tonnes C/ha) (Fig. 1) or approximately 20 tonnes C/ha in SOC stocks over the whole soil profile (Table 3 and Fig. 1). Sandy loam and sandy soils had more pronounced reductions in SOC with increasing initial SOC stocks (0.07 to − 1.6 % and 0.09 to − 1.4 %, respectively; Table 3 and Fig. 1), whereas clay soils exhibited the lowest decrement in sequestered SOC (1.8 to − 0.1 %), suggesting higher potential to accrue SOC in clay soils. Simplified scenarios elicited reductions in long-term SOC in the up permost 30 cm in response to increasing SOC levels (from 0.07 % to − 1.3 %; Table 3 and Fig. 1). In contrast, Diversified and Intensified ad aptations showed less decline in SOC accumulation (1.4 to − 0.7 % and 1.2 to − 1 %; Table 3), suggesting superior SOC accrual (Fig. 3) attain able through diversification and intensification. Lower initial SOC (0.5 %) elicited greater gains in SOC accumulation, whereas higher initial SOC concentrations (>2 % for sandy and Sandy loam soils, and >5 % for clay soil; Fig. 2) resulted in considerable SOC loss across adaptations (Table 3 and Fig. 2). This finding suggests that higher antecedent SOC or SOC that exceed ceiling SOC result in lower gains over the duration of a simulation or experimental period. Table 3 Influence of antecedent SOC on long-term change (1910–2019) in soil organic carbon (ΔSOC) in upper (0–30 cm) and lower soil profiles (30–100 cm) as influenced by soil type (grey rows) and whole-farm adaptation (blue rows). A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 6 3.2. Absolute SOC accrual and long term delta SOC Absolute SOC accrual - quantified as the change in SOC stocks over a 110-year timeframe (1910–2019) – decreased with higher initial SOC across soil types and adaptations (Table 4 and Fig. 3). Lower initial SOC (0.5 %) evoked greater SOC accrual in upper layers (0.02 t C/ha /yr), which significantly increased the absolute SOC gains (2.13 t C/ha) resulting in greater SOC accumulation over the 110 year simulation (28 t C/ha) across soil types and adaptations (Tables 4, 5, Figs. 3, 4 and S3). In contrast, initial SOC of 6 % (or above) in the topsoil resulted in SOC loss (− 0.002 tonnes C/ha/yr), considerably diminishing absolute SOC accrual (− 0.2 tonnes C/ha) and long-term SOC storage (− 26 t C/ ha; Table 4, Figs. 3, 4, S3). Deeper soil layers (30–100 cm) exhibited de minimus SOC change per annum (− 0.001 to 0.004 t C/ha/yr; Table 4) and change over time (− 4 to 1 t C/ha; Table 4 and Fig. 2). At lower initial SOC (0.5 %), clay soils had higher rates and absolute SOC change (0.03 t C/ha /yr and 3.4 t C/ha respectively), resulting in larger long-term gains in SOC (50 tonnes C/ha) in the uppermost soil layers. In contrast, sandy soils had lower accrual (0.01 t C/ha/yr and 1.3 t C/ha respectively), as shown in Tables 3, 4, 5 and Fig. 3. Lower SOC accrual with greater antecedent SOC was more evident in sandy soils compared with clay soils (Table 5). Diversified and Intensified adaptations had higher potential for SOC gains at lower initial SOC content (0.5 %) in soil surface (2.8 and 2.4 t C/ ha respectively; Table 4 and S3) relative to the Baseline and Simplified adaptations (1.8 and 1.5 t C/ha respectively; Tables 4, 5 and Fig. 3). The increased SOC accrual resulted in greater gains in long-term SOC under Diversified and Intensified scenarios (39 and 31 t C/ha respectively) compared with the Baseline and Simplified scenarios (23 and 19 t C/ha respectively; Table 3, Fig. 4 and S6). Diversified and Intensified scenarios accrued more SOC on clay soils (72 and 50 t C/ha respectively; Fig. 4), compared with Baseline and Simplified scenarios (44 and 35 t C/ha respectively). The greatest losses in SOC (− 17 to − 44 t C/ha) occurred in Simplified scenarios at high initial SOC (6 %), suggesting that simplifying cropping systems (or de- intensification) may not be an enabler for long-term SOC accrual. High SOC accrual for wheat and canola underpinned larger SOC stock gains in Intensified and Diversified scenarios under clay and sandy soils, and Diversified scenarios under sandy loam soil. Enhanced SOC gains for summer crops (maize, cotton and rice) increased SOC stocks for Intensified scenarios relative to Baseline and Simplified scenarios (Fig. 5). Even so, increments in initial SOC elicited greater SOC losses for summer crops (cotton and rice) cultivated on clay and sandy soils under Inten sified scenarios (Fig. 5). 3.3. Relative impact of initial SOC, adaptation and soil type on SOC stocks Effects of initial SOC, adaptation and soil type on SOC and in teractions therein were significant (Table 5 and S3). Initial SOC accounted for an average of 74 % of the variation in long term SOC stocks, in comparison to 12 % for soil type and 3 % across adaptations (Fig. 6). Clay soils had greater effects on long-term SOC stocks relative to sandy soils (Table S4), whereas Intensified and Diversified scenarios induced greater change in SOC stocks compared with Baseline and Simplified scenarios (Table S5). Low initial SOC (0.5–2 %) were condu cive to considerable gains in SOC storage compared with high initial SOC (3–6 %) (Table S6). 3.4. Global warming potential (GWP) and GHG emissions intensity Despite greater mean yield and significant carbon removal by SOC accrual, Intensified scenarios resulted in higher GWP (2.8–5.8 t CO₂-e) and emissions intensities (0.5–0.9 kg CO₂-e/kg crop yield, relative to the Baseline (− 0.1 to 1.9 t CO₂-e and 0.1–0.5 kg CO₂-e/kg crop yield respectively). Diversified scenarios had the lowest net GWP (− 0.7 to 1.5 t CO₂-e) and carbon footprint (− 0.2 to 0.2 kg CO₂-e/kg crop yield). Conversely, Simplified scenarios evidenced lowest yields and greater emissions intensities (0.2–0.6 kg CO₂-e/kg crop yield) in comparison with the Baseline (Figs. 7 and 8). Summer crop rotations (maize, cotton and rice) accounted for a large portion of farm GWP and emissions in tensities, associated with increased S-N₂O emissions resulting from higher application rates of N-fertilisers and irrigation on sandy soils (Fig. 9, S4 and S5). High N-fertiliser rates elicited greater N leaching in irrigated summer crops (in particular rice) on clay and sandy soils (Fig. S5). These results highlight a tight coupling between crop pro ductivity, S-N₂O emissions and N-leaching, demonstrating that intensi fication with synthetic N fertiliser can increase net GHG emissions. In addition to significant removal of atmospheric CO₂ via SOC accrual, Diversified scenarios were more effective in mitigating GHG emissions compared with Intensification (Figs. 7, 9 and S4), reflecting the potential for Diversification as a prospective climate change adaptation and mitigation intervention through addition of grain legumes (here, faba bean) to the baseline crop rotation, although Diversification had lower yields overall than those modelled for Intensification. Simplification caused relatively little difference to the Baseline (Table 1, Figs. 3, 4 and S5). Sandy and clay soils exhibited greater vulnerability N-leaching compared with sandy loam soils, the former leading to elevated Fig. 1. Relationship between antecedent soil organic carbon (SOC) and long-term change in SOC in the uppermost 30 cm of soil. Colours represent whole farm adaptations and symbols represent soil type. A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 7 emissions intensities (0.4–0.9 kg CO₂-e/kg crop yield) and GWP (1–4.7 t CO₂-e) from increased S-N₂O emissions (Fig. 7 and S5). Only sandy loam soils in Intensified systems generated noteworthy GHG emissions from subsoil layers (0.5–0.7 t CO₂-e) primarily due to leaching and S-N₂O (Fig. S5). Higher initial SOC increased GHG emissions and carbon footprint, but had no discernible impact on crop yield (Figs. 7 and 8). 4. Discussion Past work has shown that sustainable intensification underpinned by higher rates of irrigation, greater farm area under irrigation and higher rates of N application may represent a financially viable alternative in the face of the climate emergency (Muleke et al., 2023). However, this may not always be the case, as (1) irrigation water may not always be physically available or economically affordable (Monjardino et al., Fig. 2. Long-term (1910–2019) change in soil organic carbon (SOC) in the soil surface (0–30 cm; purple and green bars) and subsurface (30–100 cm; red points) at varying antecedent SOC, across soil-types (clay, sandy and sandy loam) and whole-farm adaptations (Baseline, Diversified, Intensified and Simplified). A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 8 2022), (2) intensification per se can increase N leaching and/or elevate N2O emissions, or (3) catalyse other environmental problems, such as salinisation of the ground water and surface water (Bilotto et al., 2021; Christie et al., 2018, 2020; Rawnsley et al., 2019), underscoring a need for development of more comprehensive approaches that holistically account for environmental, economic, biophysical and social trade-offs evoked by purported adaptation or mitigation interventions (Harrison, 2021; Rawnsley et al., 2018; Shahpari et al., 2021). Here, our aim was to Table 4 Absolute soil organic carbon (SOC) accrual in the upper (0–30 cm) and deep (30–100 cm) soil layers for varying antecedent SOC, soil-types (grey rows) and whole-farm adaptations (blue rows). Absolute SOC accrual was estimated as the change in SOC stocks over 110-years (1910–2019) relative to base period (1910). A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 9 disentangle the impacts of whole-farm intensification/simplifica tion/diversification on (1) long-term SOC accrual and (2) global warming potential (GWP) resulting from adaptations across varying antecedent SOC levels and soil types. While we address our aims via use case, our concepts and systems approach underpinning Simplification, Diversification and Intensification could be generically adapted to any farming system or agroecological region; Simplification being aimed at reducing inputs and resource use, Diversification aimed at increasing diversity of enterprises employed on farm (here, crop types), and Intensification designed to increase irrigation and fertilisation per unit area and per farm area. While numerical results obtained for any pro duction system or agroecological region are contextual, the approach we articulate hereto is amenable to generalisation. Although Intensification resulted in superior productivity and improved SOC accrual, this adaptation also rsulted in the highest GWP (5.8 t CO₂-e) and emissions intensity (0.9 kg CO₂-e/kg crop yield; Table 4, Figs. 4 and 7). Diversification manifested greater SOC, GWP mitigation (1.5 t CO₂-e) and carbon footprint (0.2 kg CO₂-e/kg crop yield) consistent with past work (Yang et al., 2023; Zou et al., 2022), but was accompanied by concomitant reductions in yield compared with Intensification (Figs. 7 and 8). Simplification increased emissions in tensities (0.6 kg CO₂-e/kg crop yield), diminished SOC accrual and evoked significant yield penalties relative to the Baseline (Table 4, Figs. 3, 7 and 8). These findings suggest that diversifying irrigated sys tems would seem to be one solution for both adaptation and mitigation (McGill et al., 2018; Meier et al., 2020; Paustian et al., 2019; Sándor et al., 2016), albeit at the expense of yield and profitability (Monjardino et al., 2022; Muleke et al., 2023). On the other hand, intensification can elevate productivity and profitability (Muleke et al., 2023) but would exacerbate climate change through additional GHG (McGill et al., 2018; Phelan et al., 2015) even though Intensification is often the most resilient form of adaptation (Muleke et al., 2023). These findings suggest that the long-term sustainability and endurance of any purported mitigation or adaptation option will be a function of sustainable improvement in multiple dimensions (economic, environmental, social, institutional). It is important to note that the conceptualisation of SOC processes in the model we used (APSIM) quantify change in SOC inputs relative to SOC mineralisation; as APSIM does not account for molecular SOC residence times; we thus refer to “SOC accrual” and “SOC stocks” rather than sequestration per se. SOC residence time, including sequestration and formation of mineral-associated organic carbon, are aspects worthy of future investigation. For example, Begill et al. (2023) assessed a wide range of SOC (5–118 g kg-1) and clay contents (30–770 g kg-1) and found that the proportion of mineral-associated organic carbon (MAOC) was stable across a range of cropland and grassland management types, and across SOC levels. While soil texture influenced the slope of the relationship between SOC and MAOC, no upper limit to MAOC was observed. These insights challenge the notion that MAOC accumulation is limited by the soil fine fraction, and thus may compriseamenable targets for conceptualisation and contrasting within contemporary agro-ecosystemsmodels, such as APSIM. We showed that long-term SOC stocks in the topsoil (0–30 cm) and subsoil (30–100 cm) in irrigated cropping systems was primarily gov erned by initial SOC relative to practice changes, soil type and prevailing climate (74 % effect; Table 5 and Fig. 6). We found evidence suggesting that lower initial SOC (0.5–2 %) are likely elicit greater ex poste SOC accrual compared with contexts with higher initial SOC concentrations (3–6 %). Our findings clarify confusions over previous conclusions and are consistent with cross-site empirical studies demonstrating an inverse relationship between changes in SOC and initial SOC (Arndt et al., 2022; Emde et al., 2021; Slessarev et al., 2023). We showed that clay-rich soils with low antecedent SOC (0.5 %) accrued more SOC (50 t C/ha) in topsoil (30 cm), relative to sandy soils (17 t C/ha; Tables 3–5 and Figs. 2–4). Notwithstanding our discussion above relating to MAOC, this result is consistent with previous research opining that incremental gains in SOC are associated with higher con centrations of silt and clay, wherein fine fractions (<50 µm) tend to protect SOC (Beare et al., 2014; Feng et al., 2013; Matus, 2021), such that C sorbed to mineral surfaces increases with soil specific surface area (Kaiser and Guggenberger, 2003). As such, sand-rich soils are more Fig. 3. Influence of antecedentsoil organic carbon (SOC) on the annual rate and long-term SOC accrual within the upper soil profile layers (0–30 cm) across soil-types and whole-farm adaptations. Green bars represent SOC change (t C/ha) over 110-years (1910–2019). Red points represent mean annual SOC change (t C/ha/yr). A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 10 likely to reach their C equilibrium at lower SOC due to lower SOC protective capacity in comparison to clay soils (Angers et al., 2011; Six et al., 2002; West and Six, 2007). We acknowledge that the presence of a theoretical equilibrium or “saturation” point is contextualised to the soil types in question, and indeed, is the subject of ongoing scholarly debate, with some authors suggesting that “mineral-associated organic carbon formation capacity” would be more appropriate nomenclature than use of the term “saturation” (Emde et al., 2022). Within the conceptual and mathematical boundaries of the present study however, our results suggest that sandy soils reached C saturation levels beginning from lower antecedent SOC levels (2 %) than clayey soils (5 %), evidenced by trends shown in Fig. 4. Incremental levels of initial SOC beyond this upper limit induced greater SOC losses in sandy soils compared with clayey soils. Taken together, our results suggest that management practices aimed at increasing carbon inputs may be futile where SOC is relatively high, as C-retention rate asymptotically diminishes as SOC approaches the equilibrium point (Castellano et al., 2015; Guillaume et al., 2022). Ascertaining SOC equilibrium thresholds across sites, cli mates and soil types would be fertile ground for future research, for knowledge of such thresholds would guide practitioners aspiring to improve SOC as to whether changing management practices would be worth while. Our findings suggest that winter crops (wheat and canola), cultivated on clayey and sandy soils enhanced SOC accrual for Diversified and Intensified systems (72 and 23 t C/ha respectively; Table 3 and Figs. 4–5) whereas, when grown on sandy loam soils, these crops were only likely Fig. 4. Influence of antecedent soil organic carbon (SOC) on long-term SOC stocks in the uppermost 30 cm for three soil types and four whole farm adaptations. A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 11 to enhance SOC accrual under Intensification (25 t C/ha). Summer crops (maize, cotton, and rice) improved SOC accrual in Intensified scenarios relative to Baseline and Simplified. Nevertheless, these crops were prone to N-leaching on clay and sandy soils, and increased S-N₂O emissions on sandy loam soils, resulting in notable yield losses and increased global warming potential (Figs. 4–8 and S4–5). Broadly, this finding underscores the beneficial impact of winter crops on SOC accrual, and the potential integration into crop rotations aimed at reducing GHG emissions through improved C accrual. Deeper soil layers (30–100 cm) exhibited minimal variability in SOC accrual, across initial SOC levels, soil types and whole-farm adaptations (as shown in Table 3 and Fig. 2). Notably, only sandy loam soils in intensified systems showed significant GHG emissions in subsoil layers (0.5–0.7 t CO₂-e), primarily due to elevated S-N₂O emissions (Fig. S5). We found that increasing farm area under irrigation, and/or increasing nitrogen fertilisation rates can and does improve SOC stocks. We note however, that such increase is relative to the status quo farming system in operation and antecedent SOC. In general, we found that over the 110-year simulation, SOC of Simplified systems changed little compared with SOC of the Baseline, suggesting that farm area under irrigation had relatively little impact at the whole farm scale (Simplifi cation irrigated 350 ha of wheat, whereas the Baseline irrigated 750 ha of wheat). When antecedent SOC was low, Diversification and Intensification elicited similarly high SOC increments, however when antecedent SOC was high (e.g., SOC 6 %), Intensified systems lost more SOC than Diver sified systems. Inclusion of legumes (faba bean) into canola-wheat-wheat rotations (i.e. Diversified system) improved soil nitrogen fixation and reduced soil C:N ratios, increasing soil organic matter decomposition rates, improving subsequent temporal gains in SOC. Intensified systems generally improved soil organic matter inputs when initial SOC was low, particularly for the highly productive crops, such as maize, which required an average of 15 ML/ha/year compared with 2.8 ML/ha/year for wheat crops (averaged over the simulation). However, when initial SOC was high (6%), high production and residue crops immobilised more mineral nitrogen, subsequently accelerating SOC decrement over the 110-year simulation. In general, changes in total SOC were primarily attributable to changes in the labile pool (microbial biomass and prod ucts), and to a lesser extent the humic pool, with inert organic carbon pools remaining relatively stable over the simulation. In general, changes in the proportion of fresh and humic pools of soil organic matter relative to the change in total organic matter were relatively consistent across adaptations, soil types and climate years. Fig. 5. Average annual change in soil organic carbon (SOC) storage disaggregated by crop types for three soil types and four whole-farm adaptations. Stacked columns represent cumulative SOC stocks in the surface 30 cm for six initial SOC levels. Table 5 Influence of initial SOC, adaptation and soil type on long-term change in SOC stocks. Source of variation DF Sum Sq Mean Sq F value p-value Initial SOC 5 37.4 7.5 1044.3 < 0.001 Soil type 2 5.9 3.0 412.9 < 0.001 Adaptation 3 1.6 0.5 74.4 < 0.001 Initial SOC × Adaptation 15 1.6 0.1 14.7 < 0.001 Initial SOC × Soil type 10 3.5 0.4 48.9 < 0.001 Adaptation × Soil type 6 0.4 0.07 9.51 < 0.001 Residuals 30 0.2 0.01 DF = Degrees of Freedom, Sum Sq= Sum of Squares, Mean Sq= Mean Square, × = Interaction. Fig. 6. Relative impact of initial SOC, adaptation, soil types and interactions therein on changes in SOC stocks. A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 12 Fig. 7. Net global warming potential (GWP) and emissions intensities as a function of initial SOC, whole-farm adaptations (Baseline, Diversified, Intensified and Simplified) and soil type. Green bars represent emissions intensities estimated as the annual per hectare direct emissions in kg CO₂-e per kilogram of crop yield. Blue points depict net farm carbon dioxide equivalents (t CO₂-e) computed as the sum of direct emissions from soil nitrous oxide (S-N₂O) and net change in soil organic carbon (ΔSOC) stocks in 0–100 cm layers. Fig. 8. Effect of antecedent soil organic carbon (SOC) stocks, whole farm adaptation and soil type on yield. Green violin plots represent distributions of crop yield. Yellow boxplots represent mean long-term crop yield. Black circles represent means of crop yield (averaged across adaptation and soil type). A. Muleke et al. Global Environmental Change Advances 1 (2023) 100001 13 5. Concluding remarks This study has underscored existential challenges landholders face in balancing priorities for agri-food production aginst those of environ mental stewardship and GHG emissions abatement, with Intensification enabling higher productivity and SOC accrual but also evoking greater global warming potential (5.8 t CO₂-e). Conversely, Diversification (via adding grain legumes to crop rotations) was more conducive to SOC accrual and GHG emissions mitigation, albeit with lower yields compared with Intensification. As well, Intensification (via additional irrigation and nitrogen fertiliser) was more likely to catalyse N2O emissions, counterbalancing carbon removals associated with greater SOC stocks. Collectively, these results call for more holistic assessments of GHG emissions (e.g., N2O, CH4, CO2) to account for pollution swap ping associated with intervention. We showed that clay-rich soils have higher asymptotic carbon thresholds (>5 %) and accrue more SOC, while sandy soils have lower SOC ceilings (2 %), truncating SOC pla teaus to lower thresholds, as shown in Fig. 4. We found that antecedent SOC exerted the greatest control over long-term SOC accrual relative to practice change, soil type, crop type and climate. We conclude that (1) contextualised identification and benchmarking of soil C-equilibrium thresholds may elicit more probable improvement in carbon removals and agri-food production, (2) improvement in SOC stocks had little impact on productivity compared with addition of N fertiliser, noting that the latter concurrently increased N2O emissions, and (3) identifi cation of antecedent SOC relative to SOC ceilings would help guide practitioners as to whether changing practices to improve SOC stocks would likely be a fruitful or futile endeavour. Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Data availability Data will be made available on request. Acknowledgements We acknowledge financial support from the Grains Research and Development Corporation (project UOT1906–002RTX), the Tasmanian Institute of Agriculture, the University of Tasmania Graduate Research Scholarship Program, and the Australian Sustainable Agriculture Scholarship (ASAS) Program. Appendix A. Supporting information Supplementary data associated with this article can be found in the online version at doi:10.1016/j.gecadv.2023.100001. References Adnan, M., Fahad, S., Saleem, M.H., Ali, B., Mussart, M., Ullah, R., Amanullah, Arif, M., Ahmad, M., Shah, W.A., Romman, M., Wahid, F., Wang, D., Saud, S., Liu, K., Harrison, M.T., Wu, C., Danish, S., Datta, R., Muresan, C.C., Marc, R.A., 2022. Comparative efficacy of phosphorous supplements with phosphate solubilizing bacteria for optimizing wheat yield in calcareous soils. Sci. Rep. 12 (1), 11997 https://doi.org/10.1038/s41598-022-16035-3. Angers, D.A., Arrouays, D., Saby, N.P.A., Walter, C., 2011. Estimating and mapping the carbon saturation deficit of French agricultural topsoils. Soil Use Manag. 27 (4), 448–452. https://doi.org/10.1111/j.1475-2743.2011.00366.x. Antón, R., Derrien, D., Urmeneta, H., van der Heijden, G., Enrique, A., Virto, I., 2022. Organic carbon storage and dynamics as affected by the adoption of irrigation in a cultivated calcareous mediterranean soil [Original Research]. Front. Soil Sci. 2 https://doi.org/10.3389/fsoil.2022.831775. Ara, I., Turner, L., Harrison, M.T., Monjardino, M., deVoil, P., Rodriguez, D., 2021. Application, adoption and opportunities for improving decision support systems in irrigated agriculture: a review. Agric. Water Manag. 257, 107161 https://doi.org/ 10.1016/j.agwat.2021.107161. Arndt, K.A., Campbell, E.E., Dorich, C.D., Grandy, A.S., Griffin, T.S., Ingraham, P., Perry, A., Varner, R.K., Contosta, A.R., 2022. Initial soil conditions outweigh management in a cool-season dairy farm’s carbon sequestration potential. Sci. Total Environ. 809, 152195 https://doi.org/10.1016/j.scitotenv.2021.152195. Fig. 9. Average emissions intensities disaggregated by crop type for various antecedent soil organic carbon (SOC) content and adaptations. Stacked columns represent cumulative carbon footprint in kg CO₂e per kilogram of yield. A. Muleke et al. https://doi.org/10.1016/j.gecadv.2023.100001 https://doi.org/10.1038/s41598-022-16035-3 https://doi.org/10.1111/j.1475-2743.2011.00366.x https://doi.org/10.3389/fsoil.2022.831775 https://doi.org/10.1016/j.agwat.2021.107161 https://doi.org/10.1016/j.agwat.2021.107161 https://doi.org/10.1016/j.scitotenv.2021.152195 Global Environmental Change Advances 1 (2023) 100001 14 Ash, A., Gleeson, T., Hall, M., Higgins, A., Hopwood, G., MacLeod, N., Paini, D., Poulton, P., Prestwidge, D., Webster, T., Wilson, P., 2017. Irrigated agricultural development in northern Australia: value-chain challenges and opportunities. Agric. Syst. 155, 116–125. https://doi.org/10.1016/j.agsy.2017.04.010. Balogh, J.M., 2022. The impacts of agricultural development and trade on CO2 emissions? Evidence from the Non-European Union countries. Environ. Sci. Policy 137, 99–108. https://doi.org/10.1016/j.envsci.2022.08.012. Beare, M.H., McNeill, S.J., Curtin, D., Parfitt, R.L., Jones, H.S., Dodd, M.B., Sharp, J., 2014. Estimating the organic carbon stabilisation capacity and saturation deficit of soils: a New Zealand case study. Biogeochemistry 120 (1), 71–87. https://doi.org/ 10.1007/s10533-014-9982-1. Begill, N., Don, A., Poeplau, C., 2023. No detectable upper limit of mineral-associated organic carbon in temperate agricultural soils. Glob. Change Biol. (1– 8) https://doi. org/10.1111/gcb.16804. Bilotto, F., Harrison, M.T., Migliorati, M.D.A., Christie, K.M., Rowlings, D.W., Grace, P. R., Smith, A.P., Rawnsley, R.P., Thorburn, P.J., Eckard, R.J., 2021. Can seasonal soil N mineralisation trends be leveraged to enhance pasture growth? Sci. Total Environ. 772, 145031 https://doi.org/10.1016/j.scitotenv.2021.145031. Bilotto, F., Christie-Whitehead, K.M., Malcolm, B., Harrison, M.T., 2023. Carbon, cash, cattle and the climate crisis. Sustain. Sci. https://doi.org/10.1007/s11625-023- 01323-2. Bouman, B.A. M., Kropff, M.J., Wopereis, M.C. S., ten Berge, H.F. M., van Laar, H.H., 2001. ORYZA2000: modeling lowland rice. In. https://doi.org/10.22004/ag. econ.281825. Brown, H.E., Huth, N.I., Holzworth, D.P., Teixeira, E.I., Zyskowski, R.F., Hargreaves, J.N. G., Moot, D.J., 2014. Plant modelling framework: software for building and running crop models on the APSIM platform. Environ. Modell. Softw. 62, 385–398. https:// doi.org/10.1016/j.envsoft.2014.09.005. Castellano, M.J., Mueller, K.E., Olk, D.C., Sawyer, J.E., Six, J., 2015. Integrating plant litter quality, soil organic matter stabilization, and the carbon saturation concept. Global Change Biology 21 (9), 3200–3209. https://doi.org/10.1111/gcb.12982. Chao, P.U., 2019. Residue management induced changes in soil organic carbon and total nitrogen under different tillage practices in the North China Plain. J. Integr. Agric. 18 (6), 1337–1347. https://doi.org/10.1016/S2095-3119(18)62079-9. Christie, K.M., Rawnsley, R.P., Harrison, M.T., Eckard, R.J., 2014. Using a modelling approach to evaluate two options for improving animal nitrogen use efficiency and reducing nitrous oxide emissions on dairy farms in southern Australia. Anim. Prod. Sci. 54 (12), 1960–1970. https://doi.org/10.1071/AN14436. Christie, K.M., Smith, A.P., Rawnsley, R.P., Harrison, M.T., Eckard, R.J., 2018. Simulated seasonal responses of grazed dairy pastures to nitrogen fertilizer in SE Australia: pasture production. Agric. Syst. 166, 36–47. https://doi.org/10.1016/j. agsy.2018.07.010. Christie, K.M., Smith, A.P., Rawnsley, R.P., Harrison, M.T., Eckard, R.J., 2020. Simulated seasonal responses of grazed dairy pastures to nitrogen fertilizer in SE Australia: N loss and recovery. Agric. Syst. 182, 102847 https://doi.org/10.1016/j. agsy.2020.102847. R. Core Team, 2022. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. In https://www.R-project. org/. Dalgliesh, N., Cocks, B., & Horan, H. , 2012. APSoil-providing soils information to consultants, farmers and researchers. Capturing Opportunities and Overcoming Obstacles in Australian Agronomy, Armidale, NSW. Dalgliesh, N., Hochman, Z., Huth, N., & Holzworth, D. , 2016. Field protocol to APSoil characterisations. Canberra, ACT, Australia. CSIRO Agriculture and Food. https:// publications.csiro.au/rpr/download?pid=csiro:EP166550&dsid=DS4. Dimes, J.P., Revanuru, S., 2004. Evaluation of APSIM to simulate plant growth response to applications of organic and inorganic N and P on an Alfisol and Vertisol in India. In: Delve, R.J., Probert, M.E. (Eds.), Modelling Nutrient Management in Tropical Cropping Systems, vol. 114. The Australian Centre for International Agricultural Research (ACIAR), Canberra, pp. 118–125. Emde, D., Hannam, K.D., Most, I., Nelson, L.M., Jones, M.D., 2021. Soil organic carbon in irrigated agricultural systems: A meta-analysis. Glob. Change Biol. 27 (16), 3898–3910. https://doi.org/10.1111/gcb.15680. Emde, D., Hannam, K.D., Midwood, A.J., Jones, M.D., 2022. Estimating mineral- associated organic carbon deficits in soils of the okanagan valley: a regional study with broader implications [Original Research]. Front. Soil Sci. 2 https://doi.org/ 10.3389/fsoil.2022.812249. FAO, 2022. Emissions intensities: Methodological note, release October 2022. Rome, Italy. Food and Agriculture Organization (FAO). https://www.fao.org/faostat/en/ #data/EI. Feng, W., Plante, A.F., Six, J., 2013. Improving estimates of maximal organic carbon stabilization by fine soil particles. Biogeochemistry 112 (1), 81–93. https://doi.org/ 10.1007/s10533-011-9679-7. Guillaume, T., Makowski, D., Libohova, Z., Bragazza, L., Sallaku, F., Sinaj, S., 2022. Soil organic carbon saturation in cropland-grassland systems: storage potential and soil quality. Geoderma 406, 115529. https://doi.org/10.1016/j.geoderma.2021.115529. Harrison, M.T., 2021. Climate change benefits negated by extreme heat. Nature Food. https://doi.org/10.1038/s43016-021-00387-6. Harrison, M.T., Tardieu, F., Dong, Z., Messina, C.D., Hammer, G.L., 2014. Characterizing drought stress and trait influence on maize yield under current and future conditions. Glob. Change Biol. 20 (3), 867–878. https://doi.org/10.1111/ gcb.12381. Harrison, M.T., Roggero, P.P., Zavattaro, L., 2019. Simple, efficient and robust techniques for automatic multi-objective function parameterisation: case studies of local and global optimisation using APSIM. Environ. Modell. Softw. 117, 109–133. https://doi.org/10.1016/j.envsoft.2019.03.010. Harrison, M.T., Cullen, B., Mayberry, D.E., Cowie, A.L., Bilotto, F., Badgery, W.B., Liu, K., Davison, T., Christie, K.M., Muleke, A., Eckard, R.J., 2021. Carbon myopia: the urgent need for integrated social, economic and environmental action in the livestock sector. Glob. Change Biol. https://doi.org/10.1111/gcb.15816. Hearn, A.B., 1994. OZCOT: a simulation model for cotton crop management. Agric. Syst. 44 (3), 257–299. https://doi.org/10.1016/0308-521X(94)90223-3. Henry, B., Dalal, R., Harrison, M.T., Keating, B., 2022. Creating frameworks to foster soil carbon sequestration. In: Rumpel, C. (Ed.), Understanding and Fostering Soil Carbon Sequestration. Burleigh Dodds Science Publishing Limited. https://doi.org/ 10.19103/AS.2022.0106.25. Holzworth, D.P., Huth, N.I., deVoil, P.G., Zurcher, E.J., Herrmann, N.I., McLean, G., Chenu, K., van Oosterom, E.J., Snow, V., Murphy, C., Moore, A.D., Brown, H., Whish, J.P.M., Verrall, S., Fainges, J., Bell, L.W., Peake, A.S., Poulton, P.L., Hochman, Z., Thorburn, P.J., Gaydon, D.S., Dalgliesh, N.P., Rodriguez, D., Cox, H., Chapman, S., Doherty, A., Teixeira, E., Sharp, J., Cichota, R., Vogeler, I., Li, F.Y., Wang, E., Hammer, G.L., Robertson, M.J., Dimes, J.P., Whitbread, A.M., Hunt, J., van Rees, H., McClelland, T., Carberry, P.S., Hargreaves, J.N.G., MacLeod, N., McDonald, C., Harsdorf, J., Wedgwood, S., Keating, B.A., 2014. APSIM – Evolution towards a new generation of agricultural systems simulation. Environ. Modell. Softw. 62, 327–350. https://doi.org/10.1016/j.envsoft.2014.07.009. Hunter, M.C., Smith, R.G., Schipanski, M.E., Atwood, L.W., Mortensen, D.A., 2017. Agriculture in 2050: recalibrating targets for sustainable intensification. BioScience 67 (4), 386–391. https://doi.org/10.1093/biosci/bix010. Ibrahim, A., Harrison, M.T., Meinke, H., Zhou, M., 2019. Examining the yield potential of barley near-isogenic lines using a genotype by environment by management analysis. Eur. J. Agron. 105, 41–51. https://doi.org/10.1016/j.eja.2019.02.003. IPCC, 2013. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (T. F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex, & P. M. Midgley, Eds.). Cambridge University Press. https:// www.ipcc.ch/report/ar5/wg1/. Isbell, R.F. , 2021. The Australian Soil Classification (Third ed.). The National Committee on Soil and Terrain, CSIRO Publishing, Melbourne. https://doi.org/10.1071/ 9781486314782. Jeffrey, S.J., Carter, J.O., Moodie, K.B., Beswick, A.R., 2001. Using spatial interpolation to construct a comprehensive archive of Australian climate data. Environ. Modell. Softw. 16 (4), 309–330. https://doi.org/10.1016/S1364-8152(01)00008-1. Kaiser, K., Guggenberger, G., 2003. Mineral surfaces and soil organic matter. Eur. J. Soil Sci. 54 (2), 219–236. https://doi.org/10.1046/j.1365-2389.2003.00544.x. Karstens, K., Bodirsky, B.L., Dietrich, J.P., Dondini, M., Heinke, J., Kuhnert, M., Müller, C., Rolinski, S., Smith, P., Weindl, I., Lotze-Campen, H., Popp, A., 2022. Management-induced changes in soil organic carbon on global croplands. Biogeosciences 19 (21), 5125–5149. https://doi.org/10.5194/bg-19-5125-2022. Keating, B.A., Carberry, P.S., Hammer, G.L., Probert, M.E., Robertson, M.J., Holzworth, D., Huth, N.I., Hargreaves, J.N.G., Meinke, H., Hochman, Z., McLean, G., Verburg, K., Snow, V., Dimes, J.P., Silburn, M., Wang, E., Brown, S., Bristow, K.L., Asseng, S., Chapman, S., McCown, R.L., Freebairn, D.M., Smith, C.J., 2003. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 18 (3), 267–288. https://doi.org/10.1016/S1161-0301(02)00108-9. Langworthy, A.D., Rawnsley, R.P., Freeman, M.J., Pembleton, K.G., Corkrey, R., Harrison, M.T., Lane, P.A., Henry, D.A., 2018. Potential of summer-active temperate (C3) perennial forages to mitigate the detrimental effects of supraoptimal temperatures on summer home-grown feed production in south-eastern Australian dairying regions. Crop Pasture Sci. 69 (8), 808–820. https://doi.org/10.1071/ CP17291. Liu, K., Harrison, M.T., Hunt, J., Angessa, T.T., Meinke, H., Li, C., Tian, X., Zhou, M., 2020. Identifying optimal sowing and flowering periods for barley in Australia: a modelling approach. Agric. For. Meteorol. 282–283, 107871 https://doi.org/ 10.1016/j.agrformet.2019.107871. Liu, K., Harrison, M.T., Archontoulis, S.V., Huth, N., Yang, R., Liu, D.L., Yan, H., Meinke, H., Huber, I., Feng, P., Ibrahim, A., Zhang, Y., Tian, X., Zhou, M., 2021. Climate change shifts forward flowering and reduces crop waterlogging stress. Environ. Res. Lett. 16 (9), 094017 https://doi.org/10.1088/1748-9326/ac1b5a. Liu, K., Harrison, M.T., Ibrahim, A., Nuruzzaman Manik, S.M., Johnson, P., Tian, X., Meinke, H., Zhou, M., 2020. Genetic factors increasing barley grain yields under soil waterlogging. Food Energy Secur. 9, 1–12. https://doi.org/10.1002/fes3.238 e238. Liu, K., Harrison, M.T., Yan, H., Liu, D.L., Meinke, H., Hoogenboom, G., Wang, B., Peng, B., Guan, K., Jaegermeyr, J., Wang, E., Zhang, F., Yin, X., Archontoulis, S., Nie, L., Badea, A., Man, J., Wallach, D., Zhao, J., Benjumea, A.B., Fahad, S., Tian, X., Wang, W., Tao, F., Zhang, Z., Rötter, R., Yuan, Y., Zhu, M., Dai, P., Nie, J., Yang, Y., Zhang, Y., Zhou, M., 2023. Silver lining to a climate crisis in multiple prospects for alleviating crop waterlogging under future climates. Nat. Commun. 14 (1), 765. https://doi.org/10.1038/s41467-023-36129-4. Luo, Z., Wang, E., Smith, C., & Huth, N. , 2011. Modeling long-term soil carbon dynamics in various agro-ecosystems using APSIM. In: Proceedings of the World Congress on Conservation Agriculture, 26–29 September 2011, Brisbane. Matus, F.J., 2021. Fine silt and clay content is the main factor defining maximal C and N accumulations in soils: a meta-analysis. Sci. Rep. 11 (1), 6438. https://doi.org/ 10.1038/s41598-021-84821-6. McGill, B.M., Hamilton, S.K., Millar, N., Robertson, G.P., 2018. The greenhouse gas cost of agricultural intensification with groundwater irrigation in a Midwest U.S. row cropping system. Glob. Change Biol. 24 (12), 5948–5960. https://doi.org/10.1111/ gcb.14472. Meier, E.A., Thorburn, P.J., Bell, L.W., Harrison, M.T., Biggs, J.S., 2020. Greenhouse gas emissions from cropping and grazed pastures are similar: a simulation analysis in A. Muleke et al. https://doi.org/10.1016/j.agsy.2017.04.010 https://doi.org/10.1016/j.envsci.2022.08.012 https://doi.org/10.1007/s10533-014-9982-1 https://doi.org/10.1007/s10533-014-9982-1 https://doi.org/10.1111/gcb.16804 https://doi.org/10.1111/gcb.16804 https://doi.org/10.1016/j.scitotenv.2021.145031 https://doi.org/10.1007/s11625-023-01323-2 https://doi.org/10.1007/s11625-023-01323-2 https://doi.org/10.1016/j.envsoft.2014.09.005 https://doi.org/10.1016/j.envsoft.2014.09.005 https://doi.org/10.1111/gcb.12982 https://doi.org/10.1016/S2095-3119(18)62079-9 https://doi.org/10.1071/AN14436 https://doi.org/10.1016/j.agsy.2018.07.010 https://doi.org/10.1016/j.agsy.2018.07.010 https://doi.org/10.1016/j.agsy.2020.102847 https://doi.org/10.1016/j.agsy.2020.102847 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref18 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref18 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref18 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref18 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref18 https://doi.org/10.1111/gcb.15680 https://doi.org/10.3389/fsoil.2022.812249 https://doi.org/10.3389/fsoil.2022.812249 https://doi.org/10.1007/s10533-011-9679-7 https://doi.org/10.1007/s10533-011-9679-7 https://doi.org/10.1016/j.geoderma.2021.115529 https://doi.org/10.1038/s43016-021-00387-6 https://doi.org/10.1111/gcb.12381 https://doi.org/10.1111/gcb.12381 https://doi.org/10.1016/j.envsoft.2019.03.010 https://doi.org/10.1111/gcb.15816 https://doi.org/10.1016/0308-521X(94)90223-3 https://doi.org/10.19103/AS.2022.0106.25 https://doi.org/10.19103/AS.2022.0106.25 https://doi.org/10.1016/j.envsoft.2014.07.009 https://doi.org/10.1093/biosci/bix010 https://doi.org/10.1016/j.eja.2019.02.003 https://doi.org/10.1016/S1364-8152(01)00008-1 https://doi.org/10.1046/j.1365-2389.2003.00544.x https://doi.org/10.5194/bg-19-5125-2022 https://doi.org/10.1016/S1161-0301(02)00108-9 https://doi.org/10.1071/CP17291 https://doi.org/10.1071/CP17291 https://doi.org/10.1016/j.agrformet.2019.107871 https://doi.org/10.1016/j.agrformet.2019.107871 https://doi.org/10.1088/1748-9326/ac1b5a https://doi.org/10.1002/fes3.238 https://doi.org/10.1038/s41467-023-36129-4 https://doi.org/10.1038/s41598-021-84821-6 https://doi.org/10.1038/s41598-021-84821-6 https://doi.org/10.1111/gcb.14472 https://doi.org/10.1111/gcb.14472 Global Environmental Change Advances 1 (2023) 100001 15 Australia [Original Research]. Front. Sustain. Food Syst. 3 https://doi.org/10.3389/ fsufs.2019.00121. Mohanty, M., Sinha, N.K., Somasundaram, J., McDermid, S.S., Patra, A.K., Singh, M., Dwivedi, A.K., Reddy, K.S., Rao, C.S., Prabhakar, M., Hati, K.M., Jha, P., Singh, R.K., Chaudhary, R.S., Kumar, S.N., Tripathi, P., Dalal, R.C., Gaydon, D.S., Chaudhari, S. K., 2020. Soil carbon sequestration potential in a Vertisol in central India- results from a 43-year long-term experiment and APSIM modeling. Agric. Syst. 184, 102906 https://doi.org/10.1016/j.agsy.2020.102906. Monjardino, M., Harrison, M.T., DeVoil, P., Rodriguez, D., Sadras, V.O., 2022. Agronomic and on-farm infrastructure adaptations to manage economic risk in Australian irrigated broadacre systems: a case study. Agric. Water Manag. 269, 107740 https://doi.org/10.1016/j.agwat.2022.107740. Muleke, A., Harrison, M.T., de Voil, P., Hunt, I., Liu, K., Yanotti, M., Eisner, R., 2022a. Earlier crop flowering caused by global warming alleviated by irrigation. Environ. Res. Lett. 17 (4), 044032 https://doi.org/10.1088/1748-9326/ac5a66. Muleke, A., Harrison, M.T., Eisner, R., de Voil, P., Yanotti, M., Liu, K., Yin, X., Wang, W., Monjardino, M., Zhao, J., Zhang, F., Fahad, S., Zhang, Y., 2022b. Whole farm planning raises profit despite burgeoning climate crisis. Sci. Rep. 12 (1), 17188 https://doi.org/10.1038/s41598-022-20896-z. Muleke, A., Harrison, M.T., Yanotti, M., Battaglia, M., 2022c. Yield gains of irrigated crops in Australia have stalled: the dire need for adaptation to increasingly volatile weather and market conditions. Curr. Res. Environ. Sustain. 4, 100192 https://doi. org/10.1016/j.crsust.2022.100192. Muleke, A., Harrison, M.T., Eisner, R., Voil, P. d, Yanotti, M., Liu, K., Monjardino, M., Yin, X., Wang, W., Nie, J., Ferreira, C., Zhao, J., Zhang, F., Fahad, S., Shurpali, N., Puyu, F., Zhang, Y., Forster, D., Yang, R., Qi, Z., Fei, W., Gao, X., Man, J., Lixiao, N., 2023. Sustainable intensification with irrigation raises farm profit despite climate emergency. Plants, People, Planet 5 (3), 368–385. https://doi.org/10.1002/ ppp3.10354. O’Leary, G.J., Liu, D.L., Ma, Y., Li, F.Y., McCaskill, M., Conyers, M., Dalal, R., Reeves, S., Page, K., Dang, Y.P., Robertson, F., 2016. Modelling soil organic carbon 1. Performance of APSIM crop and pasture modules against long-term experimental data. Geoderma 264, 227–237. https://doi.org/10.1016/j.geoderma.2015.11.004. Paustian, K., Larson, E., Kent, J., Marx, E., Swan, A., 2019. Soil C sequestration as a biological negative emission strategy [Review]. Front. Clim. 1 https://doi.org/ 10.3389/fclim.2019.00008. Phelan, D.C., Harrison, M.T., Kemmerer, E.P., Parsons, D., 2015. Management opportunities for boosting productivity of cool-temperate dairy farms under climate change. Agric. Syst. 138, 46–54. https://doi.org/10.1016/j.agsy.2015.05.005. Probert, M.E., Dimes, J.P., Keating, B.A., Dalal, R.C., Strong, W.M., 1998. APSIM’s water and nitrogen modules and simulation of the dynamics of water and nitrogen in fallow systems. Agric. Syst. 56 (1), 1–28. https://doi.org/10.1016/S0308-521X(97) 00028-0. Probert, M.E., Delve, R.J., Kimani, S.K., Dimes, J.P., 2005. Modelling nitrogen mineralization from manures: representing quality aspects by varying C:N ratio of sub-pools. Soil Biol. Biochem. 37 (2), 279–287. https://doi.org/10.1016/j. soilbio.2004.07.040. Rawnsley, R., Dynes, R.A., Christie, K.M., Harrison, M.T., Doran-Browne, N.A., Vibart, R., Eckard, R., 2018. A review of whole farm-system analysis in evaluating greenhouse-gas mitigation strategies from livestock production systems. Animal Production Science 58, 980–989. Rawnsley, R.P., Smith, A.P., Christie, K.M., Harrison, M.T., Eckard, R.J., 2019. Current and future direction of nitrogen fertiliser use in Australian grazing systems. Crop and Pasture Sci. 70 (12), 1034–1043. https://doi.org/10.1071/CP18566. Ritchie, H., Roser, M., Rosado, P., 2020. CO₂ and greenhouse gas emissions. Published online at OurWorldInData. org. https://ourworldindata.org/co2-and-greenhouse- gas-emissions. Robertson, M.J., Holland, J.F., Kirkegaard, J.A., & Smith, C.J. , 1999. Simulating growth and development of canola in Australia. In: Proceedings 10th International Rapeseed Congress, 26th -29th September 1999, Canberra, ACT. Sándor, R., Ehrhardt, F., Basso, B., Bellocchi, G., Bhatia, A., Brilli, L., Migliorati, M.D., Doltra, J., Dorich, C., Doro, L., Fitton, N., Giacomini, S.J., Grace, P., Grant, B., Harrison, M.T., Jones, S., Kirschbaum, M.U.F., Klumpp, K., Laville, P., Léonard, J., Liebig, M., Lieffering, M., Martin, R., McAuliffe, R., Meier, E., Merbold, L., Moore, A., Myrgiotis, V., Newton, P., Pattey, E., Recous, S., Rolinski, S., Sharp, J., Massad, R.S., Smith, P., Smith, W., Snow, V., Wu, L., Zhang, Q., Soussana, J.F., 2016. C and N models Intercomparison – benchmark and ensemble model estimates for grassland production. Adv. Anim. Biosci. 7 (3), 245–247. https://doi.org/10.1017/ S2040470016000297. Sándor, R., Ehrhardt, F., Grace, P., Recous, S., Smith, P., Snow, V., Soussana, J.-F., Basso, B., Bhatia, A., Brilli, L., Doltra, J., Dorich, C.D., Doro, L., Fitton, N., Grant, B., Harrison, M.T., Kirschbaum, M.U.F., Klumpp, K., Laville, P., Léonard, J., Martin, R., Massad, R.-S., Moore, A., Myrgiotis, V., Pattey, E., Rolinski, S., Sharp, J., Skiba, U., Smith, W., Wu, L., Zhang, Q., Bellocchi, G., 2020. Ensemble modelling of carbon fluxes in grasslands and croplands. Field Crops Res. 252, 107791 https://doi.org/ 10.1016/j.fcr.2020.107791. Sapkota, A., Haghverdi, A., Avila, C.C.E., Ying, S.C., 2020. Irrigation and greenhouse gas emissions: a review of field-based studies. Soil Syst. 4 (2), 20. 〈https://www.mdpi. com/2571-8789/4/2/20〉. Saurabh, K., Rao, K.K., Mishra, J.S., Kumar, R., Poonia, S.P., Samal, S.K., Roy, H.S., Dubey, A.K., Choubey, A.K., Mondal, S., Bhatt, B.P., Verma, M., Malik, R.K., 2021. Influence of tillage based crop establishment and residue management practices on soil quality indices and yield sustainability in rice-wheat cropping system of Eastern Indo-Gangetic Plains. Soil Tillage Res. 206, 104841 https://doi.org/10.1016/j. still.2020.104841. Shahpari, S., Allison, J., Harrison, M.T., Stanley, R., 2021. An integrated economic, environmental and social approach to agricultural land-use planning. Land 10 (4). https://doi.org/10.3390/land10040364. Singhal, R.K., Fahad, S., Kumar, P., Choyal, P., Javed, T., Jinger, D., Singh, P., Saha, D., Md, P., Bose, B., Akash, H., Gupta, N.K., Sodani, R., Dev, D., Suthar, D.L., Liu, K., Harrison, M.T., Saud, S., Shah, A.N., Nawaz, T., 2023. Beneficial elements: new players in improving nutrient use efficiency and abiotic stress tolerance. Plant Growth Regul. 100 (2), 237–265. https://doi.org/10.1007/s10725-022-00843-8. Six, J., Conant, R.T., Paul, E.A., Paustian, K., 2002. Stabilization mechanisms of soil organic matter: implications for C-saturation of soils. Plant and Soil 241 (2), 155–176. https://doi.org/10.1023/A:1016125726789. Slessarev, E.W., Mayer, A., Kelly, C., Georgiou, K., Pett-Ridge, J., Nuccio, E.E., 2023. Initial soil organic carbon stocks govern changes in soil carbon: reality or artifact? Global Change Biol. 29 (5), 1239–1247. https://doi.org/10.1111/gcb.16491. Taylor, C.A., Harrison, M.T., Telfer, M., Eckard, R., 2016. Modelled greenhouse gas emissions from beef cattle grazing irrigated leucaena in northern Australia. Anim. Prod. Sci. 56 (3), 594–604. https://doi.org/10.1071/AN15575. Trost, B., Prochnow, A., Drastig, K., Meyer-Aurich, A., Ellmer, F., Baumecker, M., 2013. Irrigation, soil organic carbon and N2O emissions. A review. Agron. Sustain. Dev. 33 (4), 733–749. https://doi.org/10.1007/s13593-013-0134-0. Turpin, J., Robertson, M., Haire, C., Bellotti, W., Moore, A., Rose, I., 2003. Simulating fababean development, growth, and yield in Australia. Crop Pasture Sci. 54, 39–52. https://doi.org/10.1071/AR02064. Wang, E., Van Oosterom, E., Meinke, H., Asseng, S., Robertson, M., Huth, N., Keating, B., & Probert, M. , 2003b. The new APSIM-Wheat Model—performance and future improvements. In: Proceedings of the 11th Australian Agronomy ConferenceSolutions for a Better Environment, 2nd - 6th February 2003, Geelong, Victoria, Australia. Wang, E., Van Oosterom, E., Meinke, H., Asseng, S., Robertson, M., Huth, N., Keating, B., Probert, M., 2003a. The new APSIM-Wheat Model—performance and future improvements. In: Proceedings of the 11th Australian Agronomy Conference, West, T.O., Six, J., 2007. Considering the influence of sequestration duration and carbon saturation on estimates of soil carbon capacity. Climatic Change 80 (1), 25–41. https://doi.org/10.1007/s10584-006-9173-8. Yan, H., Harrison, M.T., Liu, K., Wang, B., Feng, P., Fahad, S., Meinke, H., Yang, R., Liu, D.L., Archontoulis, S., Huber, I., Tian, X., Man, J., Zhang, Y., Zhou, M., 2022. Crop traits enabling yield gains under more frequent extreme climatic events. Sci. Total Environ. 808, 152170 https://doi.org/10.1016/j.scitotenv.2021.152170. Yang, Y., Ti, J., Zou, J., Wu, Y., Rees, R.M., Harrison, M.T., Li, W., Huang, W., Hu, S., Liu, K., Wen, X., Chen, F., Yin, X., 2023. Optimizing crop rotation increases soil carbon and reduces GHG emissions without sacrificing yields. Agric. Ecosyst. Environ. 342, 108220 https://doi.org/10.1016/j.agee.2022.108220. Zou, J., Yang, Y., Shi, S., Li, W., Zhao, X., Huang, J., Zhang, H., Liu, K., Harrison, M.T., Chen, F., Yin, X., 2022. Farm-scale practical strategies to reduce carbon footprint and emergy while increasing economic benefits in crop production in the North China plain. J. Clean. Prod. 359, 131996 https://doi.org/10.1016/j.jclepro.2022.131996. A. Muleke et al. https://doi.org/10.3389/fsufs.2019.00121 https://doi.org/10.3389/fsufs.2019.00121 https://doi.org/10.1016/j.agsy.2020.102906 https://doi.org/10.1016/j.agwat.2022.107740 https://doi.org/10.1088/1748-9326/ac5a66 https://doi.org/10.1038/s41598-022-20896-z https://doi.org/10.1016/j.crsust.2022.100192 https://doi.org/10.1016/j.crsust.2022.100192 https://doi.org/10.1002/ppp3.10354 https://doi.org/10.1002/ppp3.10354 https://doi.org/10.1016/j.geoderma.2015.11.004 https://doi.org/10.3389/fclim.2019.00008 https://doi.org/10.3389/fclim.2019.00008 https://doi.org/10.1016/j.agsy.2015.05.005 https://doi.org/10.1016/S0308-521X(97)00028-0 https://doi.org/10.1016/S0308-521X(97)00028-0 https://doi.org/10.1016/j.soilbio.2004.07.040 https://doi.org/10.1016/j.soilbio.2004.07.040 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref55 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref55 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref55 http://refhub.elsevier.com/S2950-1385(23)00001-3/sbref55 https://doi.org/10.1071/CP18566 https://doi.org/10.1017/S2040470016000297 https://doi.org/10.1017/S2040470016000297 https://doi.org/10.1016/j.fcr.2020.107791 https://doi.org/10.1016/j.fcr.2020.107791 https://www.mdpi.com/2571-8789/4/2/20 https://www.mdpi.com/2571-8789/4/2/20 https://doi.org/10.1016/j.still.2020.104841 https://doi.org/10.1016/j.still.2020.104841 https://doi.org/10.3390/land10040364 https://doi.org/10.1007/s10725-022-00843-8 https://doi.org/10.1023/A:1016125726789 https://doi.org/10.1111/gcb.16491 https://doi.org/10.1071/AN15575 https://doi.org/10.1007/s13593-013-0134-0 https://doi.org/10.1071/AR02064 https://doi.org/10.1007/s10584-006-9173-8 https://doi.org/10.1016/j.scitotenv.2021.152170 https://doi.org/10.1016/j.agee.2022.108220 https://doi.org/10.1016/j.jclepro.2022.131996 Kansilehti_Clarifying 1-s2.0-S2950138523000013-main Clarifying confusions over carbon conclusions: antecedent soil carbon drives gains realised following intervention 1 Introduction 2 Methods 2.1 Study overview 2.2 Historical climate data 2.3 Simulation of crop growth, yield and soil carbon stocks 2.3.1 Surface and soil organic matter and carbon 2.3.2 Management and cropping systems 2.3.3 Whole farm systems adaptations 2.3.4 Greenhouse gas (GHG) emissions and global warming potential (GWP) 2.3.4.1 Soil emissions 2.3.4.2 Emissions intensities of grain yield 2.3.5 Statistical analysis 3 Results 3.1 Antecedent SOC stocks determine SOC gain realised following intervention 3.2 Absolute SOC accrual and long term delta SOC 3.3 Relative impact of initial SOC, adaptation and soil type on SOC stocks 3.4 Global warming potential (GWP) and GHG emissions intensity 4 Discussion 5 Concluding remarks Declaration of Competing Interest Data availability Acknowledgements Appendix A Supporting information References