Key crops for processed foods have spatially variable biodiversity impacts not captured in other environmental impact indicators

Key crops for processed foods have spatially variable biodiversity impacts not captured in other environmental impact indicators | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Key crops for processed foods have spatially variable biodiversity impacts not captured in other environmental impact indicators Felicitas Pamatat, Charlotte L. Outhwaite, Amy Molotoks, Carole Dalin, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7380465/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 29 Jan, 2026 Read the published version in Scientific Reports → Version 1 posted 11 You are reading this latest preprint version Abstract We must urgently address unsustainable food system practices to ensure planetary health. Most research quantifying the environmental impacts of food production and consumption has not measured biodiversity impacts directly or spatially, limiting our ability to evaluate how producing and consuming countries are contributing to food system sustainability. Here, we quantify the biodiversity impacts of food crops commonly found in processed foods using a range of accessible metrics. We focus on ingredients found in popular processed foods, such as the chocolate digestive biscuit, and show globally reaching impacts. Land-use and fertiliser impacts were greatest for the food crops considered. We also identified Critically Endangered species at risk from habitat loss associated with crops being supplied to the UK specifically. For instance, cocoa production overlaps with nearly 20% of the range of both the critically endangered Roloway and Miss Waldron’s Red Colobus monkeys. As crop-specific farm and trade data are not publicly available, it is difficult to know the origin of all ingredients, limiting biodiversity impact estimation. If such data are released in future, our approach can be used to further evaluate biodiversity impacts of foods to inform more sustainable decision-making at consumer, business, and government levels. Biological sciences/Ecology Earth and environmental sciences/Ecology Earth and environmental sciences/Environmental sciences sustainable food system biodiversity environmental impact assessment consumption trade Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 INTRODUCTION Food production is a major driver of biodiversity decline 1 due to the shift of land use from natural lands towards agriculture, the intensification of agricultural practices, and the subsequent loss of habitats and species. Of the Earth’s habitable land surface, approximately 10% is used for cropland and 37% for livestock feed and grazing 2 . With the global human population projected to reach 9.7 billion by 2050 3 , there is a trade-off between land required for food production and for biodiversity conservation 4 . This asserts a need for more sustainable food production for planetary health to ensure that both agricultural and conservation-based benefits can be provided in the future 5,6 . Dependence on imports and differing approaches to agriculture across space and among crop groups 7 makes identifying sustainable food-production pathways complex. As such, there is considerable variability in the impacts associated with the production of each food item 7,8 . Several recent studies have assessed the environmental impacts of specific food products, including assessments of water use, carbon production, and land use 8–12 . For example, Poore and Nemecek 8 used a meta-analytical approach to highlight the variability in environmental impacts associated with different foods (e.g., rice, potatoes, nuts, citrus). Similarly, Clark et al. 10 used open data to estimate the impacts of 57,000 food products on greenhouse gas emissions, land use, water stress, and eutrophication potential, which found that more nutritious foods also tend to be more environmentally sustainable. More recently, Jwaideh and Dalin 9 developed a single, multidimensional index of the environmental impacts of crop production, considering water stress and biodiversity loss due to greenhouse gas emissions, water consumption, land occupation, and fertilizer application. Most research focusing on specific foods or diets does not include effects on biodiversity as a measure of sustainability or environmental impact. It is also uncommon for impacts to be represented spatially, showing where crops are produced and where the demand for those crops originates, limiting our understanding of the contributions of different countries to sustainable development (though notable exceptions include Jwaideh and Dalin 9 , Halpern et al. 7 , and Chaudhary and Kastner 13 ). Furthermore, trade-offs and synergies exist among measures of environmental sustainability. For instance, where carbon footprints are low, biodiversity impacts could be high, making it challenging to quantify the overall environmental sustainability of consumption, especially if studies focus on specific impact metrics or group multiple measures into a single index. These are important research gaps in the context of the United Nations Sustainable Development Goals (UN SDGs 14 ), where trade-offs likely exist between goals for ‘Life on Land’, incorporating biodiversity conservation and restoration, and ‘Zero Hunger’, which likely requires the conversion of natural, biodiverse lands to crop- and pasture-lands to meet food demands 15 . Here, we quantify the environmental impacts associated with the key ingredients of processed foods, using digestive biscuits as an example. Focusing on specific ingredients enables us to explore the challenges of understanding sustainable food systems, including representing biodiversity impacts alongside other impact-based metrics, the variation in impact across space, and novel ways to communicate these impacts to consumers. We focus part of our study on the United Kingdom (UK), as biscuits have become both a cultural norm and a staple (e.g., included as an ‘essential’ in Food Bank lists and among Trussell Trust’s most requested items 16 ) and named an ‘affordable treat’ 17 ). An estimated twelve billion biscuits are consumed in the UK per year (estimate based on the weight of a McVities dark chocolate digestive stated publicly as 16.7 grams and the weekly UK household consumption of biscuits from 2006–2020 on statista.com and UK 2019 population estimate from ons.org.uk). Biscuits are generally made using three main ingredients: wheat, sugar, and palm oil, which are also common in many other low-cost, processed foods and foods considered ‘treat’ foods, rather than essential for human health in many cases, though we acknowledge this is country and context dependent. Biscuits vary in composition and can include or exclude chocolate, so there is scope for more sustainable consumer decision-making even within this product group if chocolate (cocoa) has strong environmental impacts. As well as providing a useful entry point for raising consumer awareness of differences in product sustainability and the challenges associated with quantifying this, biscuits and processed goods are commonly excluded from biodiversity-impact studies, where the focus instead tends to be on general food groups such as meat, dairy, fruits and vegetables, and nuts (though see Konstantas et al., 2019 18 , 2018 19 ). Our work addresses this gap. Using crop distribution, trade, and environmental impact data, we: i) summarise and map the global scale environmental impacts of each of the four crops (wheat, sugar, palm oil and cocoa) identifying hotspots of high impact (Q1); ii) compare the impacts in the countries from which the UK sources the majority of these crops (Q2); and iii) investigate the extent to which threatened mammal species could be impacted by the production of these crops (Q3). As well as highlighting the international impacts of processed goods consumed in the UK, our findings can be used as a communication tool to bridge the gap between knowledge of potential impacts and consumer understanding of those impacts to support better-informed sustainable decision-making. METHODS Global-scale environmental impacts of wheat, sugar, palm oil, and cocoa In brief, we first summarised the global-scale environmental impacts of four focal crops (wheat, sugar, oil palm, and cocoa) commonly used in biscuits and other processed foods using global-scale indicators of environmental impact, based on: greenhouse gas emissions (GHGs; Carlson et al. 20 ), land use (Jwaideh and Dalin 9 ), water debt (Tuninetti et al. 21 ), and fertiliser application (nitrogen impacts on marine biodiversity and phosphorus impacts on freshwater fish; Jwaideh et al. 22 ). Using the UK as a case study, we then assessed the environmental impacts of these four crops in the regions from which they are sourced to supply the UK, identified using international food trade data. Finally, to complement the measures of environmental impact already quantified, we extended our analysis to consider especially high-risk regions by overlaying range maps of threatened mammal species 23 with maps of crop production in countries from which the UK sources the most substantial quantities of each crop (Fig. 1 ). Global scale maps of the environmental impact of crop production We used five environmental impact indicators which were developed using spatially explicit environmental models related to climate change, water stress, and terrestrial and aquatic biodiversity loss. Below, we give a brief outline of each metric and of further processing we applied. We recommend referring to Carlson et al. 20 ; Tuninetti et al. 21 ; Chaudhary et al. 24 ; and Jwaideh et al. 22 , respectively, for a more detailed explanation of their development. We mapped the impacts of all metrics resulting from the four key crops production using ‘ggplot2’ 25 and ‘raster’ 26 in R 27 . We make use of available maps for each of the five indicators at a 5 arc-minute resolution and global extent for each of our four crops of interest (though only two indicators could be mapped for cocoa using available data). Each indicator is centered on data from around the year 2000 as this most closely matched the high resolution, spatially explicit, global-scale and crop-specific data (harvested area, yield, fertiliser application, irrigated area) available at the time 28,29 . Dry-weight production values per grid cell, in tonnes, were calculated using dry-matter fractions in Bouwman et al. 30 and multiplying these by production data from Monfreda et al. 28 . Harvested-area data Data on the global harvested area of each of the four crops (wheat, cocoa, sugarcane and sugarbeet, and oil palm) were sourced from Monfreda et al. 28 . These data, which combine cropland data with national and subnational census statistics, are at a 5-arc minute resolution (approximately 9 km by 9 km at the equator) and provide the harvested area, production (tonnes), and yield values for 172 crops and groups of crops (such as berries, citrus, etc. 28 ). The data represent around the year 2000, but were selected for indices of environmental impact because of the comprehensive suite of crops represented at subnational level (e.g., relative to 42 crops available in MapSPAM 31 ) and to ensure temporal coincidence with the fertiliser data, which were only available for phosphorus for the year 2000. To complement environmental impact analyses using these data, we conducted analyses specific to threatened species using more recent global harvested-area data for 2010 for cocoa, sugar, wheat, and oil palm from MapSPAM, the International Food Policy Research Institute (IFPRI) Spatial Production Allocation Model 31,32 . (Note that here we use the term ‘oil palm’ to refer to the crop plant and ‘palm oil’ to refer to the ingredient or product from this plant. In datasets, this will often be ‘oilpalm’.) The threatened species range data represent circa 2012, so 2010 crop data were better temporally aligned for this part of our analysis. Greenhouse gas emissions (GHG) Total greenhouse gas emission intensity metrics (GHG) were quantified in tonnes of carbon dioxide equivalent per year based on Carlson et al. 20 . This metric estimates impact based on the GHG emissions (methane, carbon dioxide, and nitrous oxide) of rice paddy management, peatland draining, and fertiliser application. We multiplied the GHG by dry-weight production per grid cell, in tonnes, to quantify the GHG impact per grid cell, without the impact of the differing water content of different crops. The GHG measure is not a direct measure of biodiversity impact, but will impact biodiversity via climate change, habitat loss, and increasing risk of species extinctions. Occupational Land Use Biodiversity Impacts (LND) Land use, occupation, or degradation impacts on mammal, bird, reptile, amphibian, and vascular plant species were estimated using the occupational land use biodiversity impacts measure (LND) described in Chaudhary et al. 24 . These data use information on the relative change in species richness in different land uses relative to a natural regional reference, which is then used to model species extinctions due to land use. The unit of this measure is Potentially Disappeared Fraction of species per year, per tonne of crop. We produced impact values per grid cell by multiplying these values by dry-weight production per grid cell, in tonnes. Impacts of nitrogen on marine species (Nit) and of phosphorus fertilisers on biodiversity (Pho) Indicators of the impacts of nitrogen fertiliser on marine demersal fish (Nit, abbreviated as such to distinguish the impact indicator from the N for nitrogen more generally) and invertebrate species and of phosphorus fertiliser on freshwater fish species (Pho, abbreviated as such to distinguish the impact indicator from the P for phosphorus more generally) were calculated as Potentially Disappeared Fraction of species per year, per tonne of crop, using the method detailed in Jwaideh et al. 22 . This approach incorporates a spatially-explicit nutrient fate and transport model, with a receptor exposure and effect model to determine the biodiversity impacts of 17 crops on each system linked to the cropland soil emissions via the watershed 22 . Values per grid cell were calculated by multiplying the Nit and Pho data by dry-weight production in tonnes. Water debt (WAT) To measure the impact of the production of each crop on water sustainability, we used the ‘water debt’ measure (unit: years) detailed in Tuninetti et al. 21 . This is not a direct measure of biodiversity impact, but will impact biodiversity via water available to wild plants, freshwater species, and other fauna. The water debt indicator is available at 5 arc minute spatial resolution for the year 2000 (coincident with the Monfreda et al. 28 , crop distribution data used in our analysis). Water debt is calculated using the ratio of the crop-, water-source-, and location-specific annual water footprint in a given cell to the average renewable amount of water available in the cell each year, where the renewability rate is based on a long-term average renewability rate from 1987–2013. When the water debt value is greater than one year, the annual water footprint is unsustainable relative to the local availability of water. For more detailed information on this measure and its calculation, see Tuninetti et al. 21 . For our analyses, we included data values greater than one, to ensure we considered areas where water use is unsustainable and is therefore likely to put pressure on biodiversity. As the values for water debt for all sources can only exceed one where there is groundwater use (overuse) or if the surface water use relies on precipitation falling in other cells upstream, this means that we are representing areas which have been irrigated and do not consider appropriation of soil moisture in this measure. Mapping areas of greatest environmental impact We used maps of environmental impact per grid cell for greenhouse gas emissions (GHG), occupational land use biodiversity impacts (LND), phosphorus (Pho), nitrogen (Nit), and water debt (WAT), to identify areas of greatest impact for each of our study’s four focal crops (wheat, sugar, oil palm, and cocoa). For cocoa, only greenhouse gas emissions (GHG) and occupational land use biodiversity impact (LND) values were available. Sugar data were created by overlaying the sugarcane and sugar beet datasets for both the impact maps and the harvested-area map. Although mostly cultivated in different regions, to account for any overlap of the sugarcane and sugar beet dataset, the values of impact or harvest in occupied cells were summed. Zeros in all impact factor maps were replaced with NA, as the zeros were predominantly in areas of water bodies and were therefore present because terrestrial impacts were being measured rather than representing genuine zero values. We used scaling normalisation (minimum-maximum) to ensure all impact factor values were comparable across indicators. Notably, so we did not lose information on areas of unsustainable water debt while normalising the data, any values below one were removed before normalising, as values below one are deemed sustainable (less water is used than annually available from the source in the cell 21 , although it could be missing unsustainable cases where environmental flow requirements are not satisfied) and we focused on regions with water debt (unsustainable use, with values greater than one). We categorised environmental impact values into low ( 90%) classes and identified areas with the highest impact values as the top 10% of all impacts by cell. By filtering and selecting the highest impact areas, we were able to pinpoint regions with significant environmental concerns. To provide insights into the spatial distribution of each impact, and trade-offs between them, we created individual hotspot maps for GHG, LND, Pho, Nit, and WAT using the high-impact class. In addition to the individual impact maps, two combined maps were created: i) identifying the dominant environmental factor for each cell; and ii) mapping the total environmental impact, calculated by summing all values of the environmental indicators for each cell. We took this approach to determine the regions with the greatest environmental impacts across all available indicators (Q1). As a simpler, consumer-friendly and complementary presentation of potential regions of trade-off, we also mapped areas of overlap between high species richness and harvested area. We use species richness as a measure of biodiversity because richness data are most readily available at the global scale for the greatest number of vertebrates (here, focusing on terrestrial mammals, birds, reptiles, and amphibians). We focus on vertebrates, as these are the best-sampled taxonomic groups on a global scale. We used spatialised species-richness estimates from Etard et al. 33 , computed at 50-km resolution (27 arc minutes) on the global scale using species range-distribution maps from BirdLife International 34 and the IUCN 35 , excluding areas outside of known elevational limits for each species. The species-richness maps, which included separate maps for amphibians, mammals, birds and reptiles, were at a 50-km spatial resolution and given a Behrmann projection. To ensure that all data layers were aligned, the harvested area maps were resampled to match the species richness maps' lower resolution and cropped according to the maximum and minimum spatial extent. We used the 'ggplot2' 25 , ‘dplyr’ 36 , ‘tidyr’ 37 , ‘classInt’ 38 , ‘magrittr’ 39 , 'raster' 26 , and ‘patchwork’ 40 R packages and R 4.0.0 27 for the visualisation and data preparation. Comparing the environmental impacts of countries exporting cocoa, sugar, wheat, and palm oil to the UK We assessed the total impact across all environmental indicators for the focal crops in the countries the UK imports them from (Q2). To do this, we considered the source countries and the quantities of wheat (including wheat products), sugar (sugarbeet and sugarcane), palm oil, and cocoa (including beans, butter, paste, powder and cake) imported to the UK between 2003 and 2013. FAOSTAT trade matrices (FAOSTAT; http://www.fao.org/faostat/en/#data/TM ) were used to determine the trade of each crop and products from each partner country to the UK. These years were selected to best match the crop data (c.2000) and species-richness data (c.2012) available. These data are also representative of recent trade, as opposed to more recent years that are affected by the Covid-19 pandemic and reporting lags. A trade-flow-correction was used to re-assign trade flows to their origin country when they are re-exported, using an origin-tracing algorithm developed in Kastner et al. 41 and updated and applied in Dalin et al. 42 . Using these data, we identified the producers contributing at least 5% of the UK’s imports of the four focal crops (including UK-UK trade, here termed ‘domestic production’). We combined sugarcane and sugarbeet in trade data before consequent analyses, as many items in the source data (e.g., refined sugar) do not specify the crop, so this was the most conservative approach. Using this list of countries (Table 1 ) and the maps of crop-specific environmental impact for each of the five indicators, we assessed and compared the environmental impacts of the crops within each of the source countries. Evaluating the potential pressures being placed on threatened species by cocoa, sugar, wheat, and palm oil Focusing on terrestrial mammals, we quantified the potential for pressure on threatened species according to the IUCN Red List 23 associated with the four focal crops (Q3). We focused on mammals as they are a well-studied and well-known group, so assist in our goal of identifying reproducible ways to communicate sustainability and impact. We also concentrated on the countries providing most of the cocoa, palm oil, sugar and wheat supply to the UK (Table 1 ). We used the International Food Policy Research Institute's (IFPRI) Spatial Production Allocation Model (also known as MapSPAM) 31 , which generates global estimates of crop distributions and provides high-resolution maps (5 arc minute) of 42 crops. We used the physical area where crops are grown in 2010 for this analysis. Although more recent data are now available, we used a version of MapSPAM which is more closely aligned with the time frame covered by species-range data (c. 2012). By summing species ranges (from the IUCN 43 ) for mammals within each pixel, we assume that a species is present across its entire range and compute a spatialised estimate of species richness in 10-km grid cells across the globe. To measure the relative potential pressure placed on threatened mammal species by cocoa and oil palm, we overlaid crop area maps with maps of threatened species ranges using Google Earth Engine 44 . We then computed the proportion of threatened species ranges found within the producing areas of each focal crop to quantify the relative pressure placed on species threatened by habitat loss (potentially due to agricultural expansion). Spatialised results were cross-referenced with information on specific threats to individual Red List species. This analysis was performed for the categories Critically Endangered (CR), Endangered (EN) and Vulnerable (VU). We show results in this paper only for the highest threat category (CR) but all are provided in the Supporting Information . We took a different approach for wheat and sugar because most of these crops are supplied domestically within the UK (as wheat and sugarbeet). There are no threatened mammal species, according to the IUCN Red List, in the UK. Therefore, we used Sciurus vulgaris (the red squirrel) as an example for our case study. Although this species is given the ‘Least Concern’ status by the IUCN Red List globally, red squirrels are considered threatened in the UK. They are a priority species under the UK Biodiversity Action Plan and are the only species listed as endangered by The Mammal Society’s Red List for Britain’s Mammals 45 that have been assessed across England, Scotland and Wales. Spatial distribution data are available for this species from the IUCN 43 , so we overlaid the red squirrel and crop-area data to produce a map showing harvest values within the red squirrels’ range. RESULTS The four key crops in our study showed distinct global patterns in their production and environmental impacts ( Figs. 2 – 5 ) . Wheat was primarily produced in temperate regions such as Europe, Central Asia, North America and parts of China, while sugar production was more geographically widespread, encompassing both temperate (sugarbeet) and subtropical and tropical regions (sugarcane). Cocoa and oil palm were produced mainly in tropical countries, particularly West Africa and Southeast Asia. High levels of species richness and the harvested area of wheat predominantly overlapped in South-East Europe, Southern Russia and India (Fig. 2 c). The environmental impact of wheat was primarily driven by the nitrogen fertiliser runoff (Nit) in Northern Europe and China (Fig. 2 a). Greenhouse gas (GHG) emissions were the most prominent environmental impact factor in areas of Asia and the United States, while land-use pressures (LND) were the most significant issue in India. For sugar (sugarcane and sugar beet combined), the overall distribution of impacts was similar to wheat, with high-impact areas mainly in Europe, plus in the Caribbean and Southeast Asia-Pacific (Fig. 3 a,b). Here, land use (LND) was the main contributing impact (Fig. 3 a). However, the overlap between species richness and sugar harvested area was spatially concentrated and especially strong in India and South America (Fig. 3 c). The environmental impacts of oil palm production were concentrated in tropical regions, especially in South-East Asia & the Pacifics, West Africa and Central America (Fig. 4 a). Therefore, the areas with the greatest species richness and harvested area have the highest environmental impacts (Figs. 4 a,c). Greenhouse gas (GHG) impacts were mainly concentrated in Malaysia, while fertiliser runoff impacts (Nit) extended across Asia (Figs. 4 a and 5 a). High species-richness areas overlapped with regions of extensive oil palm harvesting in the tropics (Fig. 4 c). However, GHG emissions were lower in Côte d’Ivoire than in Indonesia, while nitrogen (Nit) and phosphorus (Pho) fertiliser runoffs remained relatively constant in both areas (Fig. 4 a). The region most impacted by oil palm was the area around Malaysia and Indonesia, with small areas in the Americas and West Africa also having high impact values (Fig. 4 b). However, the land-use (LND) impacts of oil palm were more widespread than other impact factors (Fig. 4 a). Cocoa shared similar spatial patterns to oil palm but extended across a broader geographical range, with the highest impacts in Central America and the north of South America, Madagascar, West Africa, South East Asia and the Pacific islands (Fig. 5 a,b). These areas experienced significant environmental pressures, though only GHG and LND values were available for cocoa (Fig. 5 a). Land use-related biodiversity loss (LND) was the highest impact factor in most areas where cocoa was grown, with only small areas that show a high impact as a result of GHG emissions. Comparing the environmental impacts of biscuit ingredients in countries supplying the majority of UK imports The UK sourced the majority of its wheat, sugar, palm oil, and cocoa from seven countries (Table 1 ). Much was sourced from the UK itself, including large amounts of sugarbeet and wheat (Table 1 ), for which over 7 million and 14 million tonnes are supplied from domestic production, respectively (FAOSTAT; http://www.fao.org/faostat/en/#data/TM ). Cocoa and palm oil were primarily sourced from tropical countries known to have high biodiversity. For cocoa, these countries were Côte d'Ivoire, Ghana, Nigeria and Cameroon and palm oil was primarily sourced from Indonesia and Malaysia. For palm oil, the proportion of the total production of these crops in these countries that was then exported to the UK was very low (~ 2%) yet fulfilled nearly half of UK consumption (Table 1 ). Table 1 Countries that supplied at least 5% of the UK’s consumption of wheat, cocoa, sugar and palm oil. Total supplied is the tonne weight of the total crop imported into the UK, on average, between 2003 and 2013, and the percentage of UK consumption this represents. Note that sugar was assessed by including both sugarcane and sugar beet sources; however, the majority for the UK is domestic sugar beet production. Also presented is the percentage of the exporting country’s total production for that crop. Data sourced from FAO (FAOSTAT; http://www.fao.org/faostat/en/#data/TM ) and adjusted to re-assign trade flows to origin countries rather than representing re-exports using an algorithm created by Kastner et al. 41 . A number of other countries supply less than 5% of the UK's total imports of these crops but we chose to focus on those that supply the majority. The full list of source countries is provided in Supporting Information ). Crop Country Total supplied (tonnes) Percentage of UK consumption Percentage of country’s production exported to the UK Wheat United Kingdom 14,030,964 87% NA Cocoa Côte d'Ivoire 65,916 36% 4.7% Ghana 64,620 35% 9% Nigeria 19,292 10% 4.8% Cameroon 12,280 7% 5.1% Sugar United Kingdom (sugarbeet) 7,276,708 92% NA Oil palm Indonesia 1,944,456 46% 1.8% Malaysia 1,787,703 43% 2.2% Considering sugar, wheat and oil palm production, for which all five indicators could be calculated, the greatest impact was from oil palm production in Malaysia (especially for GHG and LND), followed by oil palm production in Indonesia, then wheat in the UK, for which Nit was the key impact (Fig. 6 a). Sugar in the UK had the lowest total impact, most of which was attributable to Nit (Fig. 6 a), though sugar had a greater per-tonne impact than wheat in the UK (Fig. 6 b). For cocoa, only GHG and land use impact metrics were available. However, cocoa production showed much lower total impacts for these two metrics than the other three crops, with the highest total impact in Cote d’Ivoire (Fig. 6 a). More specifically, land-use impacts relating to Malaysian oil palm production were greatest, with a total standardised impact of 123, followed by Indonesian oil palm production at 74.5 (Fig. 6 a). Per tonne, while other crops followed overall trends, Indonesian oil palm had a substantially greater per-tonne impact than Malaysian oil palm (Fig. 6 b). Nit, LND, and GHG impacts were dominant in Indonesian oil palm impacts per tonne of crop (Fig. 6 b). The other crops showed a relatively low total land-use impact, with values between 0.005 and 1.9 (Ghanaian cocoa and UK wheat production respectively). GHGs were also a major contributor to crop environmental impacts; the greatest impacts again were from palm oil production, with Malaysian oil palm’s standardised total impact value at 74 and Indonesia’s at 42 (Fig. 6 a). GHG impact was a relatively strong contributor to UK wheat and sugar production impact too, with values of 9.7 for wheat and 1.1 for sugar (Fig. 6 a). GHG impacts from cocoa production across the four producing countries included here was relatively low, with values ranging between 0.02 and 0.21 (Fig. 6 a). For the three crops for which impact information from nitrogen fertiliser was available (oil palm, wheat and sugar), the greatest impacts were from UK wheat production (186), followed by Malaysian oil palm (66) (Fig. 6 a). Indonesian oil palm showed a much lower relative impact with a value of 11. For UK sugar production, nitrogen fertiliser impacts make up the greatest proportion of crop-related impact with a value of 14.7. Both phosphorus and water-use impacts for oil palm, wheat and sugar were much lower than the other indicator metrics, with values for phosphorus ranging between 0.0000089 and 0.00138635 and values for water between 0 (for UK sugar production) and 0.00559003 (Fig. 6 a). Quantifying the extent to which threatened species could be impacted by crop production in countries that supply the UK There is considerable overlap of palm oil and cocoa production with threatened species’ ranges (species classified as Critically Endangered - CR, Endangered - EN, or Vulnerable - VU, by the IUCN. For Indonesian palm oil, the ranges of ten Critically Endangered species of mammal overlapped with Indonesian oil palm production (Fig. 7 a). The largest recorded overlap for this highest category of threat (CR) was for Presbytis chrysomelas (6.3%), the Bornean banded langur, closely followed by Manis javanica the Sunda pangolin (Fig. 7 a). The overlap with threatened species (CR, EN and VU) was greatest in western Indonesia, though the distribution of overlap across the region was patchy (Fig. 8 a). For cocoa, the most notable overlap of species ranges with crop production was in southern Côte d’Ivoire (Fig. 8 b). Piliocolobus waldroni - Miss Waldron's red colobus - had the highest overlap with areas of cocoa production here (~ 18%) (Fig. 7 b), closely followed by Cercopithecus roloway - the Roloway Monkey- and both species ranked much higher for overlap than the four other critically endangered (CR) species (Fig. 7 b). Several hotspots of wheat and sugarbeet production overlapped with the known range of Sciurus vulgaris (red squirrel; Fig. 8 ), which is a species already pressured by habitat loss, as well as competition with invasive species and disease and is classed as “near threatened” in the UK (IUCN, 2021; Mathews and Harrower, 2020). The overlap between sugar, wheat, and the red squirrel varied by crop, with hotspots for sugar (Fig. 8 c) and a more widespread distribution for wheat (though East Anglia was a hotspot region for this crop; Fig. 8 d). DISCUSSION Previous studies have shown that the impact of crop production can vary depending on where and how it is produced 8 . However, these impacts are rarely spatially explicit and infrequently contain estimates of how biodiversity will be impacted, often focusing more on indicators such as greenhouse-gas emissions, land-use change or water use 48–50 . Not only does the impact matter for the sake of biodiversity, but because biodiversity is also critical for a sustainable and resilient agricultural system 51,52 . Here, we quantified and mapped the biodiversity impacts associated with the production of palm oil, cocoa, sugar, and wheat. We focus on the four crops in our study: i) given their dominance in ingredient lists for many luxury and ‘treat’ foods; ii) to give a first estimation of biodiversity impacts; and iii) to emphasise the need for biodiversity and broader, complementary environmental impact data in sustainability assessments of food products and their ingredients. Using UK biscuit consumption as a case study, we identified the global-scale biodiversity impacts associated with the crops grown to produce key ingredients. We show that the demand for food including these crops has strong, varied, and wide-reaching environmental impacts that vary spatially and according to country of origin. By assessing the spatial variation in several biodiversity-based impact indicators (i.e. not focusing solely on forest loss 53 ), we were also explore the impact on threatened species and suggest this as a means of communicating biodiversity impacts to consumers. In line with previous work 54,55 , we identify strong variability in impacts according to crop, region and type of indicator. This highlights the importance of assessing multiple metrics for evaluating the biodiversity impacts of different foods. For instance, we show a distinction between hotspots of environmental impact and areas of strong overlap between harvested area and high species richness (Figs. 2 – 5 ). This emphasises a need for complementary information, such as the amounts of a given crop being sourced from an area via international trade (Table 1 ), and its overlap with threatened species’ ranges (Fig. 7 ), to investigate how impacts vary spatially within and among regions. At this time, the coarse resolution of the impact maps we were able to produce at the global scale makes it challenging to identify specific areas at the local scale where impacts may be particularly high or low. This makes providing conservation recommendations challenging, highlighting the need for finer-resolution data. Nevertheless, we found that a number of Critically Endangered species are at risk from habitat loss via agricultural development in countries where the UK sources the focal crops (Fig. 7 ). Indeed, the impacts of the two imported crops - cocoa and oil palm - are strongest in species-rich tropical regions (Figs. 4 & 5 ), particularly Southeast Asia and Central/West Africa. For oil palm, the largest proportion of overlap of harvested area and species range was for Presbytis chrysomelas at over 6% followed by Manis javanica at 5.6%. However, a greater number of threatened species ranges overlapped with crop production areas for this crop in Indonesia ( Supporting Information Table S3 ). While Indonesian orangutans ( Pongo pygmaeus and Pongo abelii ) are often the ‘face’ of oil palm's impact - ‘icons of wild nature’ according to Spehar et al. 56 - our findings reveal that lesser-known species like the Presbytis chrysomelas and the Sunda pangolin are facing significant threats, urging a broader awareness of biodiversity loss. Presbytis primates are relatively understudied, but the chrysomelas species we identify as having strong overlap with oil palm is endemic to Borneo and is losing its forest habitat to logging and conversion to agriculture. This is of particular concern because these monkeys cannot survive in deforested habitats and are unlikely to cross between forest patches 57 . Previous research into the distribution of Manis javanica (the Sunda pangolin) identified forested areas of importance across Sabah in Borneo. Whilst the Sunda pangolin can adapt to secondary forests, cultivated, and urban areas 58,59 , the oil palm mapped in Indonesia (Fig. 7 ) could put pangolins at risk due to habitat fragmentation, which increases the threat of poaching by improving accessibility to humans 58,60 . Whilst our study cannot be used to accurately pinpoint overlap between pangolins and oil palm, due to the data-imposed assumption we must make that pangolins are present across their entire range, we can take from this that any overlap has a potential for direct impact via habitat loss, but also indirect effects via accessibility to humans. Recent studies further support this, highlighting that even palm oil production labelled as ‘sustainable’ has been associated with widespread habitat loss for orangutans and other species (e.g. see Gatti and Velichevskaya 61 ). For cocoa, we identified two primate species for which almost 20% of their range overlapped with crop production - Cercopithecus roloway and Piliocolobus waldroni. Whilst certain species can adapt to changing landscapes and coexist alongside agricultural production (e.g., generalist species and species able to survive in smaller patches 56,62 ), many cannot, and around 90% of Critically Endangered mammals have agricultural land listed as a major driver of their extinction risk 63 . This is the case for Miss Waldron’s red colobus ( Piliocolobus waldroni ) and the Roloway monkey ( Cercopithecus roloway ) 64,65 . Approximately 17% of tree cover has been lost in Cote d’Ivoire in recent decades 66,67 . This has been exacerbated by a large increase in human population and the spread of cocoa farming, including illegal cocoa cultivation within protected areas 65,68 . The red colobus monkey, which had the greatest overlap with cocoa production in Cote d’Ivoire of all critically endangered species analysed (Fig. 7 ), has been noted as potentially extinct in previous years and is endemic only to the two major cocoa producers - Cote d’Ivoire and Ghana 69 . These case studies reveal the ‘hidden’ impacts of global trade in producer countries, often unknown to consumers purchasing the products and difficult for sellers and researchers to quantify due to lengthy and opaque supply chains 70 . Adding to the complexity, our findings suggest UK biscuit consumption has a potentially high contribution to environmental impacts overseas, but not all of the crops we evaluated are used in biscuits. A study investigating the variation of sugar and calorie content in cakes and biscuits sold in the UK in 2016 estimated that the average sugar content is 30/100g 71 . As a case-study example, this is 4.4g per biscuit for dark chocolate digestive biscuits. Based on our estimate of UK biscuit consumption (see Introduction , calculated from public listing by McVities and data from statista.com and ONS.org.uk), this equates to over 52,000 tonnes of sugar, which is equivalent to 0.7% of total sugar produced in the UK. Similarly, for wheat, per biscuit, 96,000 tonnes of wheat is roughly equivalent to 0.7% of total wheat production. Quantities of palm oil in biscuits are generally undeclared. However, for digestive biscuits, roughly 5g of chocolate per biscuit means approximately 60,000 tonnes of chocolate is consumed in biscuits alone per year. This is about a third of the total imported cocoa (assuming ratios of cocoa:chocolate are 1:1). Although generally contributing relatively small amounts to the demand for these crops, with the exception of cocoa, biscuits and other sweet treats within which these crops are staples could provide an opportunity for individuals to make alternative choices to reduce demand for highly impactful crops. For instance, consumers could reduce their consumption of sweet treats or choose options with a potential for lower impact (e.g., an option without a chocolate coating, a different oil, or reduced sugar). Nevertheless, due to the complexities discussed, there is not always a clear ‘optimal’ consumer choice. There remains, for instance, a lack of transparency in global supply chains, especially those used by businesses creating products like biscuits on an industrial scale. Herein lies a major issue with quantifying and communicating biodiversity impacts associated with commonly consumed food products. Without improved supply-chain transparency, the case studies and methods used in our study provide a means to estimate biodiversity impact but have an unavoidable level of uncertainty. For example, we found that the UK sources some of the focal crops from relatively low-impact regions. The proportion of a partner country’s production of the focal crops supplied to the UK are low in most cases (except cocoa from Ghana, which represents 9% of the country’s production) (Table 1 ). This could lead us to believe that the UK’s environmental impact on other countries is low. However, the reality is more complex. For instance, whilst the UK’s cocoa imports represent only ~ 5% of Cote d’Ivoire’s overall cocoa production (Table 1 ), we do not know from where the UK is sourcing the cocoa within Cote d’Ivoire and from which type of farming system. There are attempts to model these supply chains in finer detail (e.g., the TRASE 72 initiative for tracing deforestation and commodity trade pathways) but the data are limited. Additionally, traceability often only extends to distribution warehouses rather than the smallholder farmers that supply them. Thus, the 5% sourced from Cote d’Ivoire could mean UK demand is responsible for a relatively small proportion of Cote d’Ivoire’s potential biodiversity impacts associated with cocoa or a much larger proportion, hidden behind a lack of farm-specific or smaller spatial scale origin-tracing data 73 . This particular issue emphasises the importance of spatialising measures of environmental impact, to better capture the origin and endpoint of impacts. Also impacting our understanding of the relative sustainability of different food products is differing baselines of biodiversity records. Although the UK has no threatened mammal species according to the IUCN Red List, this does not mean that there are no threats to UK biodiversity. Due to the country’s history of industrialization and agricultural intensification, it is likely that threatened mammalian biodiversity was lost before Red List data were collected. This bias is present in many countries with long histories of intense land-use change and influences our understanding of environmental impacts, which appear stronger in the tropics where biodiversity is higher (e.g., as discussed in Ellis et al. 74 ). This highlights the potential for injustice when making sustainable land-use decisions: do we source ingredients from places which have already gained economically via habitat loss from agricultural development rather than potentially less economically developed regions that have higher biodiversity? Furthermore, we assume that a species’ range overlapping with an area of high crop production is placing that species at risk. This is based on the fact that many species on the IUCN Red List are listed as threatened due to cropland expansion and habitat loss for agriculture, and that previous research has shown lower biodiversity in cropland than natural habitats 75,76 . Importantly, there are large differences in the production systems associated with the crops examined in this study. For instance, cocoa can be grown on small plots alongside other vegetation (shaded cocoa) which can enhance tree cover and ecosystem services 77 (though note this depends on shade tree type, land-use history, and other characteristics - e.g., Blaser-Hart et al. 78 and Maney et al. 79 ), whilst oil palm tends to be grown in vast areas of monoculture. These distinctions are lost in global-scale analyses such as ours. Unfortunately, only data for LND and GHG impacts were available for cocoa. As a result, we are underestimating the total environmental impacts of this crop, especially as fertilisers are recommended by Ghana’s cocoa board and associated with yield improvements, for example 80 . According to water footprinting of cocoa grown in Colombia, some cocoa varieties are more water-efficient than others, providing greater yields with the same quantities of rainfall 81 , highlighting a future research need for variety-specific biodiversity-impact measurement, as well as crop-specific studies. We also focus part of our analysis on threatened mammals specifically, since this taxonomic group has the best available data. However, other forms of biodiversity will also be impacted by crop production and we recommend that future work considers multiple groups and various facets of biodiversity when assessing environmental impacts. Whilst beyond the scope of our study, we anticipate that several of the metrics we used could be developed to incorporate more facets of biodiversity (e.g., functional and phylogenetic diversity). As well as impacting specific threatened species, we show that some of the UK’s biscuit ingredients are being sourced from some of the world's most environmentally impacted hotspots, particularly for oil palm and cocoa, where production, impacts, and overlap with species richness are greatest in the tropics (Figs. 2 – 5 ). However, our analysis of products consumed in the UK highlights complexities inherent in evaluating sustainability of snacks like digestive biscuits. Sustainably sourced food is increasingly high on the agenda of consumers in countries including the UK, where there is an increasing trend towards plant-based diets 82 , one of the recommendations for planetary health in the EAT-Lancet Commission report 12 . However, making diets more sustainable comes with several challenges, such as the cost, political economy, cultural norms, justice, and governance issues highlighted by Béné and Abdulai 83 . Furthermore, the complexity inherent in ensuring the overall environmental sustainability of products, particularly when considering biodiversity, makes it difficult to simply and effectively communicate the environmental impacts of food production to consumers. In the future, when more field data on species' responses to specific crops are available and supply chains are more transparent, consumers could be informed of the specific impacts of each dietary decision more easily. For now, by mapping crop-specific environmental impact hotspots and highlighting vulnerable species at risk according to their range overlap with areas of high levels of food crop production, we anticipate that this, and future studies expanding the scope to other commonly consumed food items, will help to reduce the disconnect between consumers and their biodiversity impacts beyond packaging. To date, there has been a consumer focus on reducing plastic waste to protect biodiversity. We propose that further education and industry initiatives are required to build on this established momentum towards more sustainable practices and expand the focus to ingredient origins and how to make more sustainable food choices. DATA AVAILABILITY STATEMENT The data used as inputs for the analyses in this manuscript were openly available for download, as described and cited in the manuscript. The code used for analyses are available from: https://github.com/CharlieOuthwaite/Biscuit_project . Declarations ACKNOWLEDGEMENTS We gratefully acknowledge the participants of workshop and/or ‘hackathon’ sessions for the initial discussion of the topic (Georgina Adams, Paolo Agnolucci, Elizabeth Boakes, Sara Bonetti, Álvaro Calzadilla Rivera, Silvia Ceausu, Henry Ferguson-Gow, Laura Graham, Jonathan Green, Vivienne Groner, Rob Holland, Jacquie McGlade, Monica Ortiz, Richard Pearson, Chrysanthi Rapti, Ralf Seppelt, Pete Smith, Hemant Tripathi, Peter Verburg, and Jessica Williams). FUNDING This work was supported by ‘UN SDGs: Pathways to Achievement’ funding from UCL Grand Challenges and Global Engagement Office. CLO was supported by a UK Natural Environment Research Council grant NE/R010811/1 and funding from Research England. ASAC acknowledges funding from the Wellcome Trust (Our Planet, Our Health programme, 205200, https://doi.org/10.35802/205200 - Sustainable and Healthy Food Systems (SHEFS) programme; and Climate and Health Programme, 227749/Z/23/Z - Sustainable and Healthy Food Systems - Southern Africa (SHEFS-SA)). FP was funded by the EASTBIO DTP. AM was funded by the Global Challenges Research Fund Trade, Development and the Environment Hub project (ES/S008160/1). C.D. also acknowledges funding from the Natural Environment Research Council Fellowship (NERC NE/N01524X/1) and from the European Union (ERC, FLORA, StG 101039402). AUTHOR CONTRIBUTIONS STATEMENT CLO & ASAC applied for the funding. FP, CLO, AM & ASAC conceptualised the initial project idea, contributed equally to the writing of the manuscript, and performed the analysis. ASAC coordinated the project. MJ & CD provided parts of the data and gave extensive comments on the manuscript. All authors have seen the final manuscript and take responsibility for its contents. ADDITIONAL INFORMATION COMPETING INTERESTS STATEMENT The authors declare no competing interests. References Maxwell, S. L., Fuller, R. A., Brooks, T. M. & Watson, J. E. M. Biodiversity: The ravages of guns, nets and bulldozers. Nature 536 , 143-145, doi:10.1038/536143a (2016). Ritchie, H. & Roser, M. Half of the world’s habitable land is used for agriculture (Published online at OurWorldinData.org, 2019). United Nations Department of Economic and Social Affairs. Statistical Papers - United Nations (Ser. A), Population and Vital Statistics Report 376 (2017). Delzeit, R., Zabel, F., Meyer, C. & Václavík, T. Addressing future trade-offs between biodiversity and cropland expansion to improve food security. Regional Environmental Change 17 , 1429-1441, doi:10.1007/s10113-016-0927-1 (2017). Fan, S., Headey, D., Rue, C. & Thomas, T. Food Systems for Human and Planetary Health: Economic Perspectives and Challenges. Annual Review of Resource Economics 13 , 131-156, doi:https://doi.org/10.1146/annurev-resource-101520-081337 (2021). Canavan, C. R., Noor, R. A., Golden, C. D., Juma, C. & Fawzi, W. Sustainable food systems for optimal planetary health. Transactions of The Royal Society of Tropical Medicine and Hygiene 111 , 238-240, doi:10.1093/trstmh/trx038 (2017). Halpern, B. S. et al. The environmental footprint of global food production. Nature Sustainability 5 , 1027-1039, doi:10.1038/s41893-022-00965-x (2022). Poore, J. & Nemecek, T. Reducing food’s environmental impacts through producers and consumers. Science 360 , 987-992, doi:10.1126/science.aaq0216 (2018). Jwaideh, M. A. A. & Dalin, C. The multi-dimensional environmental impact of global crop commodities. Nature Sustainability 8 , 396-410, doi:10.1038/s41893-025-01528-6 (2025). Clark, M. et al. Estimating the environmental impacts of 57,000 food products. Proceedings of the National Academy of Sciences 119 , e2120584119, doi:10.1073/pnas.2120584119 (2022). Harris, F. et al. The Water Footprint of Diets: A Global Systematic Review and Meta-analysis. Adv Nutr 11 , 375-386, doi:10.1093/advances/nmz091 (2020). Willett, W. et al. Food in the Anthropocene: the EAT-Lancet Commission on healthy diets from sustainable food systems. The Lancet 393 , 447-492, doi:10.1016/S0140-6736(18)31788-4 (2019). Chaudhary, A. & Kastner, T. Land use biodiversity impacts embodied in international food trade. Global Environmental Change 38 , 195-204, doi:https://doi.org/10.1016/j.gloenvcha.2016.03.013 (2016). United Nations General Assembly. (2015). Mehrabi, Z., Ellis, E. C. & Ramankutty, N. The challenge of feeding the world while conserving half the planet. Nature Sustainability 1 , 409-412, doi:10.1038/s41893-018-0119-8 (2018). Give Food. Trussell Trust Most Requested Items , (2025). Mintel. Uk Sweet Biscuits Market Report 2024. (2024). Konstantas, A., Stamford, L. & Azapagic, A. Evaluation of environmental sustainability of biscuits at the product and sectoral levels. Journal of Cleaner Production 230 , 1217-1228, doi:https://doi.org/10.1016/j.jclepro.2019.05.095 (2019). Konstantas, A., Jeswani, H. K., Stamford, L. & Azapagic, A. Environmental impacts of chocolate production and consumption in the UK. Food Research International 106 , 1012-1025, doi:https://doi.org/10.1016/j.foodres.2018.02.042 (2018). Carlson, K. M. et al. Greenhouse gas emissions intensity of global croplands. Nature Climate Change 7 , 63-68, doi:10.1038/nclimate3158 (2017). Tuninetti, M., Tamea, S. & Dalin, C. Water Debt Indicator Reveals Where Agricultural Water Use Exceeds Sustainable Levels. Water Resources Research 55 , 2464-2477, doi:https://doi.org/10.1029/2018WR023146 (2019). Jwaideh, M. A. A., Sutanudjaja, E. H. & Dalin, C. Global impacts of nitrogen and phosphorus fertiliser use for major crops on aquatic biodiversity. The International Journal of Life Cycle Assessment 27 , 1058-1080, doi:10.1007/s11367-022-02078-1 (2022). IUCN. (2021). Chaudhary, A., Verones, F., de Baan, L. & Hellweg, S. Quantifying Land Use Impacts on Biodiversity: Combining Species–Area Models and Vulnerability Indicators. Environmental Science & Technology 49 , 9987-9995, doi:10.1021/acs.est.5b02507 (2015). ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics v. 3.5.3 (2025). raster: Geographic Data Analysis and Modeling v. R package version 3.6-23 (2023). R: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, Vienna, Austria, 2023). Monfreda, C., Ramankutty, N. & Foley, J. A. Farming the planet: 2. Geographic distribution of crop areas, yields, physiological types, and net primary production in the year 2000. Global Biogeochemical Cycles 22 , n/a-n/a, doi:10.1029/2007gb002947 (2008). Mueller, N. D. et al. Closing yield gaps through nutrient and water management. Nature 490 , 254-257, doi:10.1038/nature11420 (2012). Bouwman, A. F. et al. Lessons from temporal and spatial patterns in global use of N and P fertilizer on cropland. Scientific Reports 7 , 40366, doi:10.1038/srep40366 (2017). International Food Policy Research, I. MapSPAM (ed Institute International Food Policy Research) (Harvard Dataverse, 2019). Yu, Q. et al. A cultivated planet in 2010 – Part 2: The global gridded agricultural-production maps. Earth Syst. Sci. Data 12 , 3545-3572, doi:10.5194/essd-12-3545-2020 (2020). Etard, A., Morrill, S. & Newbold, T. Global gaps in trait data for terrestrial vertebrates. Global Ecology and Biogeography 29 , 2143-2158, doi:https://doi.org/10.1111/geb.13184 (2020). Birdlife International & NatureServe. (See http://www.birdlife.org/datazone/info/spcdownload., 2012). IUCN. (See http://www.iucnredlist.org/, 2013). dplyr: A Grammar of Data Manipulation (2023). tidyr: Tidy Messy Data (2023). classInt: Choose Univariate Class Intervals (2023). magrittr: A Forward-Pipe Operator for R (2022). patchwork: The Composer of Plots v. 1.3.0 (2024). Kastner, T., Kastner, M. & Nonhebel, S. Tracing distant environmental impacts of agricultural products from a consumer perspective. Ecological Economics 70 , 1032-1040, doi:https://doi.org/10.1016/j.ecolecon.2011.01.012 (2011). Dalin, C., Wada, Y., Kastner, T. & Puma, M. J. Groundwater depletion embedded in international food trade. Nature 543 , 700-704, doi:10.1038/nature21403 (2017). IUCN. (See http://www.iucnredlist.org/, 2012). Gorelick, N. et al. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sensing of Environment 202 , 18-27, doi:https://doi.org/10.1016/j.rse.2017.06.031 (2017). Mathews, F. & Harrower, C. IUCN – compliant Red List for Britain’s Terrestrial Mammals. Assessment by the Mammal Society under contract to Natural England, Natural Resources Wales and Scottish Natural Heritage. . (Natural England, Peterborough, 2020). maps: Draw Geographical Maps v. R package version 3.4.2 (2023). ArcGIS Pro 3.5 (ESRI, Redlands, CA, 2025). Dalin, C. & Rodríguez-Iturbe, I. Environmental impacts of food trade via resource use and greenhouse gas emissions. Environmental Research Letters 11 , 035012, doi:10.1088/1748-9326/11/3/035012 (2016). Molotoks, A., Smith, P. & Dawson, T. P. Impacts of land use, population, and climate change on global food security. Food and Energy Security 10 , e261, doi:https://doi.org/10.1002/fes3.261 (2021). Dinar, A., Tieu, A. & Huynh, H. Water scarcity impacts on global food production. Global Food Security 23 , 212-226, doi:https://doi.org/10.1016/j.gfs.2019.07.007 (2019). Elouafi, I. Why biodiversity matters in agriculture and food systems. Science 386 , eads8197, doi:10.1126/science.ads8197 (2024). Dardonville, M., Bockstaller, C., Villerd, J. & Therond, O. Resilience of agricultural systems: biodiversity-based systems are stable, while intensified ones are resistant and high-yielding. Agricultural Systems 197 , 103365, doi:https://doi.org/10.1016/j.agsy.2022.103365 (2022). Molotoks, A., Green, J., Ribeiro, V., Wang, Y. & West, C. Assessing the value of biodiversity-specific footprinting metrics linked to South American soy trade. People and Nature n/a , doi:https://doi.org/10.1002/pan3.10457 (2023). Chaudhary, A., Pfister, S. & Hellweg, S. Spatially Explicit Analysis of Biodiversity Loss Due to Global Agriculture, Pasture and Forest Land Use from a Producer and Consumer Perspective. Environ Sci Technol 50 , 3928-3936, doi:10.1021/acs.est.5b06153 (2016). Chapman, A. S. A. et al. Quantifying the biodiversity pressures of fruit and vegetable consumption in the United Kingdom, India, and South Africa. Nature Food (2025). Spehar, S. N. et al. Orangutans venture out of the rainforest and into the Anthropocene. Science Advances 4 , e1701422, doi:doi:10.1126/sciadv.1701422 (2018). Ehlers Smith, D. A. The effects of land-use policies on the conservation of Borneo’s endemic Presbytis monkeys. Biodiversity and Conservation 23 , 891-908, doi:10.1007/s10531-014-0639-0 (2014). Nursamsi, I., Amir, Z., Decoeur, H., Moore, J. H. & Luskin, M. S. Sunda pangolins show inconsistent responses to disturbances across multiple scales. Wildlife Letters 1 , 59-70, doi:https://doi.org/10.1002/wll2.12010 (2023). Chong, J. L. et al. in Pangolins (eds Daniel W. S. Challender, Helen C. Nash, & Carly Waterman) 89-108 (Academic Press, 2020). Panjang, E. et al. Mapping the distribution of the Sunda pangolin (Manis javanica) within natural forest in Sabah, Malaysian Borneo. Global Ecology and Conservation 52 , e02962, doi:https://doi.org/10.1016/j.gecco.2024.e02962 (2024). Cazzolla Gatti, R. & Velichevskaya, A. Certified “sustainable” palm oil took the place of endangered Bornean and Sumatran large mammals habitat and tropical forests in the last 30 years. Science of The Total Environment 742 , 140712, doi:https://doi.org/10.1016/j.scitotenv.2020.140712 (2020). Ancrenaz, M. et al. Of Pongo, palms and perceptions: a multidisciplinary assessment of Bornean orang-utans Pongo pygmaeus in an oil palm context. Oryx 49 , 465-472, doi:10.1017/S0030605313001270 (2015). Ripple, W. J. et al. Extinction risk is most acute for the world’s largest and smallest vertebrates. Proceedings of the National Academy of Sciences 114 , 10678-10683, doi:10.1073/pnas.1702078114 (2017). McGraw, W. S. & Oates, J. F. in Primates in Peril: The World's 25 Most Endangered Primates 2014-2016 (eds C. Schwitzer et al. ) 14-16 (IUCN SSC Primate Specialist Group (PSG), International Primatological Society (IPS), Conservation International (CI), and Bristol Zoological Society, 2015). Bitty, E. A., Bi, S. G., Bene, J.-C. K., Kouassi, P. K. & McGraw, W. S. Cocoa Farming and Primate Extirpation Inside Cote D'ivoire's Protected Areas. Tropical Conservation Science 8 , 95-113, doi:10.1177/194008291500800110 (2015). Chatelain, C., Dao, H., Gautier, L. & Spichiger, R. Forest cover changes in Côte d’Ivoire and Upper Guinea. CABI , 15–32, doi:10.1079/9780851997346.0015 (2004). Global Forest Watch. (2019). Barima, Y. S. S. et al. Cocoa crops are destroying the forest reserves of the classified forest of Haut-Sassandra (Ivory Coast). Global Ecology and Conservation 8 , 85-98, doi:https://doi.org/10.1016/j.gecco.2016.08.009 (2016). Linder, J. M. et al. Red colobus (Piliocolobus) conservation action plan 2021–2026. (IUCN, Gland, Switzerland, 2021). Gardner, T. A. et al. Transparency and sustainability in global commodity supply chains. World Development 121 , 163-177, doi:https://doi.org/10.1016/j.worlddev.2018.05.025 (2019). Hashem, K. M., He, F. J., Alderton, S. A. & MacGregor, G. A. Cross-sectional survey of the amount of sugar and energy in cakes and biscuits on sale in the UK for the evaluation of the sugar-reduction programme. BMJ Open 8 , e019075, doi:10.1136/bmjopen-2017-019075 (2018). Trase.earth. Trase.earth website , (2025). Renier, C. et al. Transparency, traceability and deforestation in the Ivorian cocoa supply chain. Environmental Research Letters 18 , 024030, doi:10.1088/1748-9326/acad8e (2023). Ellis, C. J., Yahr, R. & Coppins, B. J. Quantifying the anthropocene loss of bioindicators for an early industrial region: an equitable baseline for biodiversity restoration. Biodiversity and Conservation 27 , 2363-2377, doi:10.1007/s10531-018-1541-y (2018). Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520 , 45-50, doi:10.1038/nature14324 (2015). Millard, J. et al. Global effects of land-use intensity on local pollinator biodiversity. Nature Communications 12 , 2902, doi:10.1038/s41467-021-23228-3 (2021). Jezeer, R. E., Verweij, P. A., Santos, M. J. & Boot, R. G. A. Shaded Coffee and Cocoa – Double Dividend for Biodiversity and Small-scale Farmers. Ecological Economics 140 , 136-145, doi:https://doi.org/10.1016/j.ecolecon.2017.04.019 (2017). Blaser-Hart, W. J. et al. The effectiveness of cocoa agroforests depends on shade-tree canopy height. Agriculture, Ecosystems & Environment 322 , 107676, doi:https://doi.org/10.1016/j.agee.2021.107676 (2021). Maney, C., Sassen, M. & Hill, S. L. L. Modelling biodiversity responses to land use in areas of cocoa cultivation. Agriculture, Ecosystems & Environment 324 , 107712, doi:https://doi.org/10.1016/j.agee.2021.107712 (2022). Arthur, A., Ja, D. & Fk, P. Impact of fertiliser application on cocoa yields in Ghana: An analysis of cocoa bean yields in farmers’ plantations. Australian Journal of Crop Science 18 , 139-144, doi:10.21475/ajcs.24.18.02.PNE3958 (2024). Naranjo-Merino, C. A., Ortíz-Rodriguez, O. O. & Villamizar-G, R. A. Assessing Green and Blue Water Footprints in the Supply Chain of Cocoa Production: A Case Study in the Northeast of Colombia. Sustainability 10 , 38 (2018). Tedstone, A., Duval, D. & Peacock, E. Dietary health and CVD: implications for dietary policy in England. Proceedings of the Nutrition Society 79 , 95-102, doi:10.1017/S0029665119000600 (2020). Béné, C. & Abdulai, A.-R. Navigating the politics and processes of food systems transformation: guidance from a holistic framework. Frontiers in Sustainable Food Systems Volume 8 - 2024 , doi:10.3389/fsufs.2024.1399024 (2024). Additional Declarations No competing interests reported. Supplementary Files SupplementaryInformation.docx Cite Share Download PDF Status: Published Journal Publication published 29 Jan, 2026 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Revision requested 17 Nov, 2025 Reviews received at journal 11 Nov, 2025 Reviewers agreed at journal 02 Nov, 2025 Reviews received at journal 15 Oct, 2025 Reviewers agreed at journal 01 Oct, 2025 Reviewers agreed at journal 08 Sep, 2025 Reviewers invited by journal 02 Sep, 2025 Editor assigned by journal 28 Aug, 2025 Editor invited by journal 21 Aug, 2025 Submission checks completed at journal 20 Aug, 2025 First submitted to journal 20 Aug, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7380465","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":511511373,"identity":"0450958a-64e3-4890-a0bf-5884c299ef04","order_by":0,"name":"Felicitas Pamatat","email":"","orcid":"","institution":"University of Aberdeen","correspondingAuthor":false,"prefix":"","firstName":"Felicitas","middleName":"","lastName":"Pamatat","suffix":""},{"id":511511374,"identity":"1fe0b5b4-5a87-43d5-88db-e5f23bd66085","order_by":1,"name":"Charlotte L. Outhwaite","email":"","orcid":"","institution":"University College London","correspondingAuthor":false,"prefix":"","firstName":"Charlotte","middleName":"L.","lastName":"Outhwaite","suffix":""},{"id":511511375,"identity":"42732844-a258-4869-a679-bd2a2c2fa791","order_by":2,"name":"Amy Molotoks","email":"","orcid":"","institution":"University of York","correspondingAuthor":false,"prefix":"","firstName":"Amy","middleName":"","lastName":"Molotoks","suffix":""},{"id":511511376,"identity":"673acdec-0fab-4868-acc6-d18400c04463","order_by":3,"name":"Carole Dalin","email":"","orcid":"","institution":"University College London","correspondingAuthor":false,"prefix":"","firstName":"Carole","middleName":"","lastName":"Dalin","suffix":""},{"id":511511377,"identity":"7619d5d7-1a36-4aa9-b373-cd72d0da2838","order_by":4,"name":"Mark Jwaideh","email":"","orcid":"","institution":"University of Oxford","correspondingAuthor":false,"prefix":"","firstName":"Mark","middleName":"","lastName":"Jwaideh","suffix":""},{"id":511511378,"identity":"7458abc0-eb88-485f-87c9-308d74536f8c","order_by":5,"name":"Abbie S. A. Chapman","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/ElEQVRIiWNgGAWjYDACCR6GA0AqAcKrYGDgg7ASiNVyhoGBjRgtCBWMbURo4Z/de/DABwa7PP7Z3YmfC+cdlmdjYH74gbEtDbcld84lHJzBkFwscefsZumZ2w4btjGwGUswtuXgdtiNHIPDQO8kNtzI3SDNu+1wAtBhZkAXVuDUIQ/S8geoZf6N3M2/eeeAtLB/w6vFAKSFAahlw43cbdK8DSAtPCBbcDvMEOSXHoPkxI13zm6z5jmWbtjGzFMskXAOt/flbvce/vCjwi5x3u3ezbd5aqzl+dnbN374UJaM2/sQ5wGxBIzDzIA3IpGABGElo2AUjIJRMEIBAEpiVSIb8xR5AAAAAElFTkSuQmCC","orcid":"","institution":"University College London","correspondingAuthor":true,"prefix":"","firstName":"Abbie","middleName":"S. A.","lastName":"Chapman","suffix":""}],"badges":[],"createdAt":"2025-08-15 10:08:12","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7380465/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7380465/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-025-34850-2","type":"published","date":"2026-01-29T15:59:23+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":90813201,"identity":"716903d4-2de1-4c9b-a131-67ac6e1d8d91","added_by":"auto","created_at":"2025-09-08 12:26:33","extension":"jpeg","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":678974,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eUsing global-scale maps representing five environmental impact indicators for four key crops (wheat, sugar, cocoa and oil palm), we assessed three questions (Q1-Q3) relating to sustainable crop production. \u003c/strong\u003eGHG refers to greenhouse gas emissions, LND the occupational land use biodiversity impacts , Nit - nitrogen fertiliser impacts on marine demersal fish and invertebrate species, Pho - phosphorus impacts on freshwater fish species, and WAT the water debt, measuring the sustainability of water use.\u003c/p\u003e","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/6893149a088cbe46d44b03f2.jpeg"},{"id":90813599,"identity":"ea64b34f-20ea-41dd-9c8a-53d322e7e040","added_by":"auto","created_at":"2025-09-08 12:34:33","extension":"jpeg","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1356199,"visible":true,"origin":"","legend":"\u003cp\u003eThe impact of wheat production on: greenhouse gas emissions (GHG), occupational land use biodiversity (LND), marine demersal fish and invertebrate species via nitrogen fertiliser application (Nit), freshwater fish species via phosphorus (Pho) fertiliser application, and the water debt (WAT) (all impacts shown normalised); a) highest impact within each cell; b) top 10% of areas with the highest summed impact values; c) bivariate map of harvested area and species richness. These maps were generated using the ‘raster’\u003csup\u003e26\u003c/sup\u003e and ‘ggplot2’\u003csup\u003e25\u003c/sup\u003e R packages with a ‘maps’ package\u003csup\u003e46\u003c/sup\u003e world outline.\u003c/p\u003e","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/371f88c2ed77bcbac4b6b35b.jpeg"},{"id":90813205,"identity":"4ec95d3a-5eeb-4a1f-a68a-6d67009ba7b5","added_by":"auto","created_at":"2025-09-08 12:26:33","extension":"jpeg","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":1191685,"visible":true,"origin":"","legend":"\u003cp\u003eThe impact of sugar (sugarcane and -beet combined) on: greenhouse gas emissions (GHG), occupational land use biodiversity (LND), marine demersal fish and invertebrate species via nitrogen fertiliser (Nit), freshwater fish species via phosphorus (Pho), and the water debt (WAT) (all impacts shown normalised); a) highest impact within each cell; b) top 10% of areas with the highest summed impact values; c) bivariate map of harvested area and species richness. These maps were generated using the ‘raster’\u003csup\u003e26\u003c/sup\u003e and ‘ggplot2’\u003csup\u003e25\u003c/sup\u003e R packages with a ‘maps’ package\u003csup\u003e46\u003c/sup\u003e world outline.\u003c/p\u003e","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/c08ba02efed80bbb86f7d4c5.jpeg"},{"id":90813597,"identity":"2d02f8d1-33f9-46f3-8e4e-e71f632a0ce4","added_by":"auto","created_at":"2025-09-08 12:34:33","extension":"jpeg","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":1067902,"visible":true,"origin":"","legend":"\u003cp\u003eThe impact of oil palm on: greenhouse gas emissions (GHG), occupational land use biodiversity (LND), marine demersal fish and invertebrate species via nitrogen fertiliser (Nit), freshwater fish species via phosphorus (Pho), and the water debt (WAT) (all impacts shown normalised); a) highest impact within each cell; b) top 10% of areas with the highest summed impact values; c) bivariate map of harvested area and species richness. These maps were generated using the ‘raster’\u003csup\u003e26\u003c/sup\u003e and ‘ggplot2’\u003csup\u003e25\u003c/sup\u003e R packages with a ‘maps’ package\u003csup\u003e46\u003c/sup\u003e world outline.\u003c/p\u003e","description":"","filename":"floatimage4.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/d21bcebcd39cdb092b00c35f.jpeg"},{"id":90814299,"identity":"e0ce333f-494f-4caa-994d-da471b1535f3","added_by":"auto","created_at":"2025-09-08 12:42:33","extension":"jpeg","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1122094,"visible":true,"origin":"","legend":"\u003cp\u003eThe impact of cocoa on greenhouse gas emissions (GHG) and occupational land use biodiversity (LND) (all impacts shown normalised); a) highest impact within each cell; b) top 10% of areas with the highest summed impact values; c) bivariate map of harvested area and species richness. These maps were generated using the ‘raster’\u003csup\u003e26\u003c/sup\u003e and ‘ggplot2’\u003csup\u003e25\u003c/sup\u003e R packages with a ‘maps’ package\u003csup\u003e46\u003c/sup\u003e world outline.\u003c/p\u003e","description":"","filename":"floatimage5.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/a062ec111cc4f6c6cc2bccc1.jpeg"},{"id":90813206,"identity":"ac318e28-1b59-483f-b579-dffc01cb6928","added_by":"auto","created_at":"2025-09-08 12:26:33","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":256451,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eThe environmental impact of growing the crops in countries that supplied at least 5% of the UK’s palm oil, sugar, wheat and cocoa.\u003c/strong\u003e The bar charts show the standardised value of each environmental indicator generated for each crop, in total (\u003cstrong\u003ea\u003c/strong\u003e) and per tonne of crop (\u003cstrong\u003eb\u003c/strong\u003e). GHG - greenhouse gas emissions; LND - occupational land use biodiversity impacts; Nit - nitrogen fertiliser impacts on marine demersal fish and invertebrate species; Pho - phosphorus fertiliser impacts on freshwater fish species; and WAT - water debt. Note that, for cocoa, measures of impact were only available for GHG and LND, these are highlighted with stars.. The x-axis has been square root transformed to better visualise those countries with low impact scores.\u003c/p\u003e","description":"","filename":"floatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/6339f4fb66d730aab9352d9d.png"},{"id":90813600,"identity":"5da159cb-a06b-4a1d-8aa0-916745305e08","added_by":"auto","created_at":"2025-09-08 12:34:34","extension":"jpeg","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":211611,"visible":true,"origin":"","legend":"\u003cp\u003eThe percentage of critically endangered (CR) species’ ranges for a) oil palm in Indonesia and b) cocoa in Cote d’Ivoire which overlapped with areas producing each crop. Species common names are as follows, for reference: \u003cem\u003eRhinoceros sondaicus: \u003c/em\u003eJavan rhinoceros\u003cem\u003e, Zaglossus bruijnii: \u003c/em\u003eWestern long-beaked echidna,\u003cem\u003eSpilocuscus rufoniger: \u003c/em\u003eBlack-spotted cuscus,\u003cem\u003e Pongo abelii:\u003c/em\u003e Sumatran orangutan,\u003cem\u003e Pongo tapanuliensis:\u003c/em\u003e Tapanuli orangutan,\u003cem\u003e Pongo pygmaeus: \u003c/em\u003eBornean orangutan,\u003cem\u003e Nycticebus bancanus: \u003c/em\u003eBangka slow loris,\u003cem\u003e Dicerorhinus sumatrensis: \u003c/em\u003eSumatran rhinoceros,\u003cem\u003e Manis javanica:\u003c/em\u003e Sunda pangolin, \u003cem\u003ePresbytis chrysomelas: \u003c/em\u003eSarawak surili\u003cem\u003e, Hipposideros lamottei:\u003c/em\u003e Lamotte’s roundleaf bat,\u003cem\u003e Diceros bicornis:\u003c/em\u003e black rhinoceros,\u003cem\u003e Crocidura wimmeri: \u003c/em\u003eWimmer’s shrew,\u003cem\u003e Colobus vellerosus: \u003c/em\u003eursine colobus,\u003cem\u003e Cercopithecus roloway: \u003c/em\u003eRoloway monkey,\u003cem\u003e Piliocolobus waldroni: \u003c/em\u003eMiss Waldron’s red colobus.\u003c/p\u003e","description":"","filename":"floatimage7.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/0436b7b5e603afc06ed5c938.jpeg"},{"id":90814300,"identity":"3600b2da-c20d-4c4c-924d-271da8e04cb9","added_by":"auto","created_at":"2025-09-08 12:42:34","extension":"jpeg","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":144154,"visible":true,"origin":"","legend":"\u003cp\u003eMaps of threatened species richness (critically endangered - CR, endangered - EN, and vulnerable - VU) with crop production areas for a) palm oil in Indonesia and b) cocoa in Côte d’Ivoire. For c) and d) the spatial overlap with crop production and the red squirrel’s (\u003cem\u003eSciurus vulgaris\u003c/em\u003e) range in the UK is shown for sugar beet and wheat. Values for a) and b) are the number of critically endangered species whereas for c) and d) they show the area of crop production in hectares. These maps were produced using Google Earth Engine\u003csup\u003e44\u003c/sup\u003e, ArcGIS Pro version 3.5\u003csup\u003e47\u003c/sup\u003e, and basemaps made with Natural Earth version 4.1.0.\u003c/p\u003e","description":"","filename":"floatimage8.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/e0c3001e603a4375f0cd7150.jpeg"},{"id":101691093,"identity":"a8112857-7c88-49bf-80df-056383593297","added_by":"auto","created_at":"2026-02-02 16:12:07","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":7310937,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/4e51739b-f6aa-435e-b7a4-68a51273eba4.pdf"},{"id":90813203,"identity":"32cd8315-122a-40a4-a33b-4e5df1a42d9b","added_by":"auto","created_at":"2025-09-08 12:26:33","extension":"docx","order_by":0,"title":"","display":"","copyAsset":false,"role":"supplement","size":47466,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryInformation.docx","url":"https://assets-eu.researchsquare.com/files/rs-7380465/v1/6091a1b0370469062caeaff7.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Key crops for processed foods have spatially variable biodiversity impacts not captured in other environmental impact indicators","fulltext":[{"header":"INTRODUCTION","content":"\u003cp\u003eFood production is a major driver of biodiversity decline\u003csup\u003e1\u003c/sup\u003e due to the shift of land use from natural lands towards agriculture, the intensification of agricultural practices, and the subsequent loss of habitats and species. Of the Earth\u0026rsquo;s habitable land surface, approximately 10% is used for cropland and 37% for livestock feed and grazing\u003csup\u003e2\u003c/sup\u003e. With the global human population projected to reach 9.7\u0026nbsp;billion by 2050\u003csup\u003e3\u003c/sup\u003e, there is a trade-off between land required for food production and for biodiversity conservation\u003csup\u003e4\u003c/sup\u003e. This asserts a need for more sustainable food production for planetary health to ensure that both agricultural and conservation-based benefits can be provided in the future\u003csup\u003e5,6\u003c/sup\u003e.\u003c/p\u003e\u003cp\u003eDependence on imports and differing approaches to agriculture across space and among crop groups\u003csup\u003e7\u003c/sup\u003emakes identifying sustainable food-production pathways complex. As such, there is considerable variability in the impacts associated with the production of each food item\u003csup\u003e7,8\u003c/sup\u003e. Several recent studies have assessed the environmental impacts of specific food products, including assessments of water use, carbon production, and land use\u003csup\u003e8\u0026ndash;12\u003c/sup\u003e. For example, Poore and Nemecek\u003csup\u003e8\u003c/sup\u003e used a meta-analytical approach to highlight the variability in environmental impacts associated with different foods (e.g., rice, potatoes, nuts, citrus). Similarly, Clark et al.\u003csup\u003e10\u003c/sup\u003e used open data to estimate the impacts of 57,000 food products on greenhouse gas emissions, land use, water stress, and eutrophication potential, which found that more nutritious foods also tend to be more environmentally sustainable. More recently, Jwaideh and Dalin\u003csup\u003e9\u003c/sup\u003e developed a single, multidimensional index of the environmental impacts of crop production, considering water stress and biodiversity loss due to greenhouse gas emissions, water consumption, land occupation, and fertilizer application.\u003c/p\u003e\u003cp\u003eMost research focusing on specific foods or diets does not include effects on biodiversity as a measure of sustainability or environmental impact. It is also uncommon for impacts to be represented spatially, showing where crops are produced and where the demand for those crops originates, limiting our understanding of the contributions of different countries to sustainable development (though notable exceptions include Jwaideh and Dalin\u003csup\u003e9\u003c/sup\u003e, Halpern et al. \u003csup\u003e7\u003c/sup\u003e, and Chaudhary and Kastner\u003csup\u003e13\u003c/sup\u003e). Furthermore, trade-offs and synergies exist among measures of environmental sustainability. For instance, where carbon footprints are low, biodiversity impacts could be high, making it challenging to quantify the overall environmental sustainability of consumption, especially if studies focus on specific impact metrics or group multiple measures into a single index. These are important research gaps in the context of the United Nations Sustainable Development Goals (UN SDGs\u003csup\u003e14\u003c/sup\u003e), where trade-offs likely exist between goals for \u0026lsquo;Life on Land\u0026rsquo;, incorporating biodiversity conservation and restoration, and \u0026lsquo;Zero Hunger\u0026rsquo;, which likely requires the conversion of natural, biodiverse lands to crop- and pasture-lands to meet food demands\u003csup\u003e15\u003c/sup\u003e.\u003c/p\u003e\u003cp\u003eHere, we quantify the environmental impacts associated with the key ingredients of processed foods, using digestive biscuits as an example. Focusing on specific ingredients enables us to explore the challenges of understanding sustainable food systems, including representing biodiversity impacts alongside other impact-based metrics, the variation in impact across space, and novel ways to communicate these impacts to consumers. We focus part of our study on the United Kingdom (UK), as biscuits have become both a cultural norm and a staple (e.g., included as an \u0026lsquo;essential\u0026rsquo; in Food Bank lists and among Trussell Trust\u0026rsquo;s most requested items\u003csup\u003e16\u003c/sup\u003e) and named an \u0026lsquo;affordable treat\u0026rsquo;\u003csup\u003e17\u003c/sup\u003e). An estimated twelve billion biscuits are consumed in the UK per year (estimate based on the weight of a McVities dark chocolate digestive stated publicly as 16.7 grams and the weekly UK household consumption of biscuits from 2006\u0026ndash;2020 on statista.com and UK 2019 population estimate from ons.org.uk). Biscuits are generally made using three main ingredients: wheat, sugar, and palm oil, which are also common in many other low-cost, processed foods and foods considered \u0026lsquo;treat\u0026rsquo; foods, rather than essential for human health in many cases, though we acknowledge this is country and context dependent. Biscuits vary in composition and can include or exclude chocolate, so there is scope for more sustainable consumer decision-making even within this product group if chocolate (cocoa) has strong environmental impacts. As well as providing a useful entry point for raising consumer awareness of differences in product sustainability and the challenges associated with quantifying this, biscuits and processed goods are commonly excluded from biodiversity-impact studies, where the focus instead tends to be on general food groups such as meat, dairy, fruits and vegetables, and nuts (though see Konstantas et al., 2019\u003csup\u003e18\u003c/sup\u003e, 2018\u003csup\u003e19\u003c/sup\u003e). Our work addresses this gap. Using crop distribution, trade, and environmental impact data, we: i) summarise and map the global scale environmental impacts of each of the four crops (wheat, sugar, palm oil and cocoa) identifying hotspots of high impact (Q1); ii) compare the impacts in the countries from which the UK sources the majority of these crops (Q2); and iii) investigate the extent to which threatened mammal species could be impacted by the production of these crops (Q3). As well as highlighting the international impacts of processed goods consumed in the UK, our findings can be used as a communication tool to bridge the gap between knowledge of potential impacts and consumer understanding of those impacts to support better-informed sustainable decision-making.\u003c/p\u003e"},{"header":"METHODS","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\u003ch2\u003eGlobal-scale environmental impacts of wheat, sugar, palm oil, and cocoa\u003c/h2\u003e\u003cp\u003eIn brief, we first summarised the global-scale environmental impacts of four focal crops (wheat, sugar, oil palm, and cocoa) commonly used in biscuits and other processed foods using global-scale indicators of environmental impact, based on: greenhouse gas emissions (GHGs; Carlson et al. \u003csup\u003e20\u003c/sup\u003e), land use (Jwaideh and Dalin\u003csup\u003e9\u003c/sup\u003e), water debt (Tuninetti et al.\u003csup\u003e21\u003c/sup\u003e), and fertiliser application (nitrogen impacts on marine biodiversity and phosphorus impacts on freshwater fish; Jwaideh et al. \u003csup\u003e22\u003c/sup\u003e). Using the UK as a case study, we then assessed the environmental impacts of these four crops in the regions from which they are sourced to supply the UK, identified using international food trade data. Finally, to complement the measures of environmental impact already quantified, we extended our analysis to consider especially high-risk regions by overlaying range maps of threatened mammal species\u003csup\u003e23\u003c/sup\u003e with maps of crop production in countries from which the UK sources the most substantial quantities of each crop (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\n\u003ch3\u003eGlobal scale maps of the environmental impact of crop production\u003c/h3\u003e\n\u003cp\u003eWe used five environmental impact indicators which were developed using spatially explicit environmental models related to climate change, water stress, and terrestrial and aquatic biodiversity loss. Below, we give a brief outline of each metric and of further processing we applied. We recommend referring to Carlson et al.\u003csup\u003e20\u003c/sup\u003e; Tuninetti et al. \u003csup\u003e21\u003c/sup\u003e; Chaudhary et al.\u003csup\u003e24\u003c/sup\u003e; and Jwaideh et al. \u003csup\u003e22\u003c/sup\u003e, respectively, for a more detailed explanation of their development. We mapped the impacts of all metrics resulting from the four key crops production using \u0026lsquo;ggplot2\u0026rsquo;\u003csup\u003e25\u003c/sup\u003e and \u0026lsquo;raster\u0026rsquo;\u003csup\u003e26\u003c/sup\u003e in R\u003csup\u003e27\u003c/sup\u003e. We make use of available maps for each of the five indicators at a 5 arc-minute resolution and global extent for each of our four crops of interest (though only two indicators could be mapped for cocoa using available data). Each indicator is centered on data from around the year 2000 as this most closely matched the high resolution, spatially explicit, global-scale and crop-specific data (harvested area, yield, fertiliser application, irrigated area) available at the time\u003csup\u003e28,29\u003c/sup\u003e. Dry-weight production values per grid cell, in tonnes, were calculated using dry-matter fractions in Bouwman et al.\u003csup\u003e30\u003c/sup\u003e and multiplying these by production data from Monfreda et al.\u003csup\u003e28\u003c/sup\u003e.\u003c/p\u003e\n\u003ch3\u003eHarvested-area data\u003c/h3\u003e\n\u003cp\u003eData on the global harvested area of each of the four crops (wheat, cocoa, sugarcane and sugarbeet, and oil palm) were sourced from Monfreda et al.\u003csup\u003e28\u003c/sup\u003e. These data, which combine cropland data with national and subnational census statistics, are at a 5-arc minute resolution (approximately 9 km by 9 km at the equator) and provide the harvested area, production (tonnes), and yield values for 172 crops and groups of crops (such as berries, citrus, etc.\u003csup\u003e28\u003c/sup\u003e). The data represent around the year 2000, but were selected for indices of environmental impact because of the comprehensive suite of crops represented at subnational level (e.g., relative to 42 crops available in MapSPAM\u003csup\u003e31\u003c/sup\u003e) and to ensure temporal coincidence with the fertiliser data, which were only available for phosphorus for the year 2000.\u003c/p\u003e\u003cp\u003eTo complement environmental impact analyses using these data, we conducted analyses specific to threatened species using more recent global harvested-area data for 2010 for cocoa, sugar, wheat, and oil palm from MapSPAM, the International Food Policy Research Institute (IFPRI) Spatial Production Allocation Model\u003csup\u003e31,32\u003c/sup\u003e. (Note that here we use the term \u0026lsquo;oil palm\u0026rsquo; to refer to the crop plant and \u0026lsquo;palm oil\u0026rsquo; to refer to the ingredient or product from this plant. In datasets, this will often be \u0026lsquo;oilpalm\u0026rsquo;.) The threatened species range data represent circa 2012, so 2010 crop data were better temporally aligned for this part of our analysis.\u003c/p\u003e\n\u003ch3\u003eGreenhouse gas emissions (GHG)\u003c/h3\u003e\n\u003cp\u003eTotal greenhouse gas emission intensity metrics (GHG) were quantified in tonnes of carbon dioxide equivalent per year based on Carlson et al.\u003csup\u003e20\u003c/sup\u003e. This metric estimates impact based on the GHG emissions (methane, carbon dioxide, and nitrous oxide) of rice paddy management, peatland draining, and fertiliser application. We multiplied the GHG by dry-weight production per grid cell, in tonnes, to quantify the GHG impact per grid cell, without the impact of the differing water content of different crops. The GHG measure is not a direct measure of biodiversity impact, but will impact biodiversity via climate change, habitat loss, and increasing risk of species extinctions.\u003c/p\u003e\n\u003ch3\u003eOccupational Land Use Biodiversity Impacts (LND)\u003c/h3\u003e\n\u003cp\u003eLand use, occupation, or degradation impacts on mammal, bird, reptile, amphibian, and vascular plant species were estimated using the occupational land use biodiversity impacts measure (LND) described in Chaudhary et al.\u003csup\u003e24\u003c/sup\u003e. These data use information on the relative change in species richness in different land uses relative to a natural regional reference, which is then used to model species extinctions due to land use. The unit of this measure is Potentially Disappeared Fraction of species per year, per tonne of crop. We produced impact values per grid cell by multiplying these values by dry-weight production per grid cell, in tonnes.\u003c/p\u003e\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\u003ch2\u003eImpacts of nitrogen on marine species (Nit) and of phosphorus fertilisers on biodiversity (Pho)\u003c/h2\u003e\u003cp\u003eIndicators of the impacts of nitrogen fertiliser on marine demersal fish (Nit, abbreviated as such to distinguish the impact indicator from the N for nitrogen more generally) and invertebrate species and of phosphorus fertiliser on freshwater fish species (Pho, abbreviated as such to distinguish the impact indicator from the P for phosphorus more generally) were calculated as Potentially Disappeared Fraction of species per year, per tonne of crop, using the method detailed in Jwaideh et al.\u003csup\u003e22\u003c/sup\u003e. This approach incorporates a spatially-explicit nutrient fate and transport model, with a receptor exposure and effect model to determine the biodiversity impacts of 17 crops on each system linked to the cropland soil emissions via the watershed\u003csup\u003e22\u003c/sup\u003e. Values per grid cell were calculated by multiplying the Nit and Pho data by dry-weight production in tonnes.\u003c/p\u003e\u003c/div\u003e\n\u003ch3\u003eWater debt (WAT)\u003c/h3\u003e\n\u003cp\u003eTo measure the impact of the production of each crop on water sustainability, we used the \u0026lsquo;water debt\u0026rsquo; measure (unit: years) detailed in Tuninetti et al.\u003csup\u003e21\u003c/sup\u003e. This is not a direct measure of biodiversity impact, but will impact biodiversity via water available to wild plants, freshwater species, and other fauna. The water debt indicator is available at 5 arc minute spatial resolution for the year 2000 (coincident with the Monfreda et al.\u003csup\u003e28\u003c/sup\u003e, crop distribution data used in our analysis). Water debt is calculated using the ratio of the crop-, water-source-, and location-specific annual water footprint in a given cell to the average renewable amount of water available in the cell each year, where the renewability rate is based on a long-term average renewability rate from 1987\u0026ndash;2013. When the water debt value is greater than one year, the annual water footprint is unsustainable relative to the local availability of water. For more detailed information on this measure and its calculation, see Tuninetti et al.\u003csup\u003e21\u003c/sup\u003e. For our analyses, we included data values greater than one, to ensure we considered areas where water use is unsustainable and is therefore likely to put pressure on biodiversity. As the values for water debt for all sources can only exceed one where there is groundwater use (overuse) or if the surface water use relies on precipitation falling in other cells upstream, this means that we are representing areas which have been irrigated and do not consider appropriation of soil moisture in this measure.\u003c/p\u003e\n\u003ch3\u003eMapping areas of greatest environmental impact\u003c/h3\u003e\n\u003cp\u003eWe used maps of environmental impact per grid cell for greenhouse gas emissions (GHG), occupational land use biodiversity impacts (LND), phosphorus (Pho), nitrogen (Nit), and water debt (WAT), to identify areas of greatest impact for each of our study\u0026rsquo;s four focal crops (wheat, sugar, oil palm, and cocoa). For cocoa, only greenhouse gas emissions (GHG) and occupational land use biodiversity impact (LND) values were available. Sugar data were created by overlaying the sugarcane and sugar beet datasets for both the impact maps and the harvested-area map. Although mostly cultivated in different regions, to account for any overlap of the sugarcane and sugar beet dataset, the values of impact or harvest in occupied cells were summed. Zeros in all impact factor maps were replaced with NA, as the zeros were predominantly in areas of water bodies and were therefore present because terrestrial impacts were being measured rather than representing genuine zero values. We used scaling normalisation (minimum-maximum) to ensure all impact factor values were comparable across indicators. Notably, so we did not lose information on areas of unsustainable water debt while normalising the data, any values below one were removed before normalising, as values below one are deemed sustainable (less water is used than annually available from the source in the cell\u003csup\u003e21\u003c/sup\u003e, although it could be missing unsustainable cases where environmental flow requirements are not satisfied) and we focused on regions with water debt (unsustainable use, with values greater than one).\u003c/p\u003e\u003cp\u003eWe categorised environmental impact values into low (\u0026lt;\u0026thinsp;10%), medium (11\u0026ndash;89%), and high-impact (\u0026gt;\u0026thinsp;90%) classes and identified areas with the highest impact values as the top 10% of all impacts by cell. By filtering and selecting the highest impact areas, we were able to pinpoint regions with significant environmental concerns. To provide insights into the spatial distribution of each impact, and trade-offs between them, we created individual hotspot maps for GHG, LND, Pho, Nit, and WAT using the high-impact class. In addition to the individual impact maps, two combined maps were created: i) identifying the dominant environmental factor for each cell; and ii) mapping the total environmental impact, calculated by summing all values of the environmental indicators for each cell. We took this approach to determine the regions with the greatest environmental impacts across all available indicators (Q1).\u003c/p\u003e\u003cp\u003eAs a simpler, consumer-friendly and complementary presentation of potential regions of trade-off, we also mapped areas of overlap between high species richness and harvested area. We use species richness as a measure of biodiversity because richness data are most readily available at the global scale for the greatest number of vertebrates (here, focusing on terrestrial mammals, birds, reptiles, and amphibians). We focus on vertebrates, as these are the best-sampled taxonomic groups on a global scale. We used spatialised species-richness estimates from Etard et al.\u003csup\u003e33\u003c/sup\u003e, computed at 50-km resolution (27 arc minutes) on the global scale using species range-distribution maps from BirdLife International\u003csup\u003e34\u003c/sup\u003e and the IUCN\u003csup\u003e35\u003c/sup\u003e, excluding areas outside of known elevational limits for each species. The species-richness maps, which included separate maps for amphibians, mammals, birds and reptiles, were at a 50-km spatial resolution and given a Behrmann projection. To ensure that all data layers were aligned, the harvested area maps were resampled to match the species richness maps' lower resolution and cropped according to the maximum and minimum spatial extent.\u003c/p\u003e\u003cp\u003eWe used the 'ggplot2'\u003csup\u003e25\u003c/sup\u003e, \u0026lsquo;dplyr\u0026rsquo;\u003csup\u003e36\u003c/sup\u003e, \u0026lsquo;tidyr\u0026rsquo;\u003csup\u003e37\u003c/sup\u003e, \u0026lsquo;classInt\u0026rsquo;\u003csup\u003e38\u003c/sup\u003e, \u0026lsquo;magrittr\u0026rsquo;\u003csup\u003e39\u003c/sup\u003e, 'raster'\u003csup\u003e26\u003c/sup\u003e, and \u0026lsquo;patchwork\u0026rsquo;\u003csup\u003e40\u003c/sup\u003e R packages and R 4.0.0\u003csup\u003e27\u003c/sup\u003e for the visualisation and data preparation.\u003c/p\u003e\u003cp\u003e\u003cb\u003eComparing the environmental impacts of countries exporting cocoa, sugar, wheat, and palm oil to the UK\u003c/b\u003e\u003c/p\u003e\u003cp\u003eWe assessed the total impact across all environmental indicators for the focal crops in the countries the UK imports them from (Q2). To do this, we considered the source countries and the quantities of wheat (including wheat products), sugar (sugarbeet and sugarcane), palm oil, and cocoa (including beans, butter, paste, powder and cake) imported to the UK between 2003 and 2013. FAOSTAT trade matrices (FAOSTAT; \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.fao.org/faostat/en/#data/TM\u003c/span\u003e\u003cspan address=\"http://www.fao.org/faostat/en/#data/TM\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e)\u003c/span\u003e were used to determine the trade of each crop and products from each partner country to the UK. These years were selected to best match the crop data (c.2000) and species-richness data (c.2012) available. These data are also representative of recent trade, as opposed to more recent years that are affected by the Covid-19 pandemic and reporting lags. A trade-flow-correction was used to re-assign trade flows to their origin country when they are re-exported, using an origin-tracing algorithm developed in Kastner et al.\u003csup\u003e41\u003c/sup\u003e and updated and applied in Dalin et al.\u003csup\u003e42\u003c/sup\u003e. Using these data, we identified the producers contributing at least 5% of the UK\u0026rsquo;s imports of the four focal crops (including UK-UK trade, here termed \u0026lsquo;domestic production\u0026rsquo;). We combined sugarcane and sugarbeet in trade data before consequent analyses, as many items in the source data (e.g., refined sugar) do not specify the crop, so this was the most conservative approach.\u003c/p\u003e\u003cp\u003eUsing this list of countries (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e) and the maps of crop-specific environmental impact for each of the five indicators, we assessed and compared the environmental impacts of the crops within each of the source countries.\u003c/p\u003e\u003cp\u003e\u003cb\u003eEvaluating the potential pressures being placed on threatened species by cocoa, sugar, wheat, and palm oil\u003c/b\u003e\u003c/p\u003e\u003cp\u003eFocusing on terrestrial mammals, we quantified the potential for pressure on threatened species according to the IUCN Red List\u003csup\u003e23\u003c/sup\u003e associated with the four focal crops (Q3). We focused on mammals as they are a well-studied and well-known group, so assist in our goal of identifying reproducible ways to communicate sustainability and impact. We also concentrated on the countries providing most of the cocoa, palm oil, sugar and wheat supply to the UK (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eWe used the International Food Policy Research Institute's (IFPRI) Spatial Production Allocation Model (also known as MapSPAM)\u003csup\u003e31\u003c/sup\u003e, which generates global estimates of crop distributions and provides high-resolution maps (5 arc minute) of 42 crops. We used the physical area where crops are grown in 2010 for this analysis. Although more recent data are now available, we used a version of MapSPAM which is more closely aligned with the time frame covered by species-range data (c. 2012). By summing species ranges (from the IUCN\u003csup\u003e43\u003c/sup\u003e) for mammals within each pixel, we assume that a species is present across its entire range and compute a spatialised estimate of species richness in 10-km grid cells across the globe.\u003c/p\u003e\u003cp\u003eTo measure the relative potential pressure placed on threatened mammal species by cocoa and oil palm, we overlaid crop area maps with maps of threatened species ranges using Google Earth Engine\u003csup\u003e44\u003c/sup\u003e. We then computed the proportion of threatened species ranges found within the producing areas of each focal crop to quantify the relative pressure placed on species threatened by habitat loss (potentially due to agricultural expansion). Spatialised results were cross-referenced with information on specific threats to individual Red List species. This analysis was performed for the categories Critically Endangered (CR), Endangered (EN) and Vulnerable (VU). We show results in this paper only for the highest threat category (CR) but all are provided in the \u003cb\u003eSupporting Information\u003c/b\u003e.\u003c/p\u003e\u003cp\u003eWe took a different approach for wheat and sugar because most of these crops are supplied domestically within the UK (as wheat and sugarbeet). There are no threatened mammal species, according to the IUCN Red List, in the UK. Therefore, we used \u003cem\u003eSciurus vulgaris\u003c/em\u003e (the red squirrel) as an example for our case study. Although this species is given the \u0026lsquo;Least Concern\u0026rsquo; status by the IUCN Red List globally, red squirrels are considered threatened in the UK. They are a priority species under the UK Biodiversity Action Plan and are the only species listed as endangered by The Mammal Society\u0026rsquo;s Red List for Britain\u0026rsquo;s Mammals\u003csup\u003e45\u003c/sup\u003e that have been assessed across England, Scotland and Wales. Spatial distribution data are available for this species from the IUCN\u003csup\u003e43\u003c/sup\u003e, so we overlaid the red squirrel and crop-area data to produce a map showing harvest values within the red squirrels\u0026rsquo; range.\u003c/p\u003e"},{"header":"RESULTS","content":"\u003cp\u003eThe four key crops in our study showed distinct global patterns in their production and environmental impacts \u003cb\u003e(\u003c/b\u003eFigs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003e\u003cb\u003e)\u003c/b\u003e. Wheat was primarily produced in temperate regions such as Europe, Central Asia, North America and parts of China, while sugar production was more geographically widespread, encompassing both temperate (sugarbeet) and subtropical and tropical regions (sugarcane). Cocoa and oil palm were produced mainly in tropical countries, particularly West Africa and Southeast Asia.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eHigh levels of species richness and the harvested area of wheat predominantly overlapped in South-East Europe, Southern Russia and India (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e2\u003c/span\u003ec). The environmental impact of wheat was primarily driven by the nitrogen fertiliser runoff (Nit) in Northern Europe and China (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). Greenhouse gas (GHG) emissions were the most prominent environmental impact factor in areas of Asia and the United States, while land-use pressures (LND) were the most significant issue in India.\u003c/p\u003e\u003cp\u003eFor sugar (sugarcane and sugar beet combined), the overall distribution of impacts was similar to wheat, with high-impact areas mainly in Europe, plus in the Caribbean and Southeast Asia-Pacific (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e3\u003c/span\u003ea,b). Here, land use (LND) was the main contributing impact (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e3\u003c/span\u003ea). However, the overlap between species richness and sugar harvested area was spatially concentrated and especially strong in India and South America (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e3\u003c/span\u003ec).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eThe environmental impacts of oil palm production were concentrated in tropical regions, especially in South-East Asia \u0026amp; the Pacifics, West Africa and Central America (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). Therefore, the areas with the greatest species richness and harvested area have the highest environmental impacts (Figs.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ea,c). Greenhouse gas (GHG) impacts were mainly concentrated in Malaysia, while fertiliser runoff impacts (Nit) extended across Asia (Figs.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ea and \u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). High species-richness areas overlapped with regions of extensive oil palm harvesting in the tropics (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ec). However, GHG emissions were lower in C\u0026ocirc;te d\u0026rsquo;Ivoire than in Indonesia, while nitrogen (Nit) and phosphorus (Pho) fertiliser runoffs remained relatively constant in both areas (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). The region most impacted by oil palm was the area around Malaysia and Indonesia, with small areas in the Americas and West Africa also having high impact values (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003eb). However, the land-use (LND) impacts of oil palm were more widespread than other impact factors (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003ea).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eCocoa shared similar spatial patterns to oil palm but extended across a broader geographical range, with the highest impacts in Central America and the north of South America, Madagascar, West Africa, South East Asia and the Pacific islands (Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003ea,b). These areas experienced significant environmental pressures, though only GHG and LND values were available for cocoa (Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003ea). Land use-related biodiversity loss (LND) was the highest impact factor in most areas where cocoa was grown, with only small areas that show a high impact as a result of GHG emissions.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\u003ch2\u003eComparing the environmental impacts of biscuit ingredients in countries supplying the majority of UK imports\u003c/h2\u003e\u003cp\u003eThe UK sourced the majority of its wheat, sugar, palm oil, and cocoa from seven countries (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). Much was sourced from the UK itself, including large amounts of sugarbeet and wheat (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), for which over 7\u0026nbsp;million and 14\u0026nbsp;million tonnes are supplied from domestic production, respectively (FAOSTAT; \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.fao.org/faostat/en/#data/TM\u003c/span\u003e\u003cspan address=\"http://www.fao.org/faostat/en/#data/TM\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e).\u003c/span\u003e Cocoa and palm oil were primarily sourced from tropical countries known to have high biodiversity. For cocoa, these countries were C\u0026ocirc;te d'Ivoire, Ghana, Nigeria and Cameroon and palm oil was primarily sourced from Indonesia and Malaysia. For palm oil, the proportion of the total production of these crops in these countries that was then exported to the UK was very low (~\u0026thinsp;2%) yet fulfilled nearly half of UK consumption (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003e\u003cb\u003eCountries that supplied at least 5% of the UK\u0026rsquo;s consumption of wheat, cocoa, sugar and palm oil.\u003c/b\u003e Total supplied is the tonne weight of the total crop imported into the UK, on average, between 2003 and 2013, and the percentage of UK consumption this represents. Note that sugar was assessed by including both sugarcane and sugar beet sources; however, the majority for the UK is domestic sugar beet production. Also presented is the percentage of the exporting country\u0026rsquo;s total production for that crop. Data sourced from FAO (FAOSTAT; \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.fao.org/faostat/en/#data/TM\u003c/span\u003e\u003cspan address=\"http://www.fao.org/faostat/en/#data/TM\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e) and adjusted to re-assign trade flows to origin countries rather than representing re-exports using an algorithm created by Kastner et al.\u003csup\u003e41\u003c/sup\u003e. A number of other countries supply less than 5% of the UK's total imports of these crops but we chose to focus on those that supply the majority. The full list of source countries is provided in \u003cb\u003eSupporting Information\u003c/b\u003e).\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"5\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eCrop\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eCountry\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eTotal supplied (tonnes)\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003ePercentage of UK consumption\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003ePercentage of country\u0026rsquo;s production exported to the UK\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eWheat\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eUnited Kingdom\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e14,030,964\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e87%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003eNA\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"3\" rowspan=\"4\"\u003e\u003cp\u003e\u003cb\u003eCocoa\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eC\u0026ocirc;te d'Ivoire\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e65,916\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e36%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e4.7%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eGhana\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e64,620\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e35%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e9%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eNigeria\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e19,292\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e10%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e4.8%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eCameroon\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e12,280\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e7%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e5.1%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\"\u003e\u003cp\u003e\u003cb\u003eSugar\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eUnited Kingdom (sugarbeet)\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e7,276,708\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e92%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003eNA\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003e\u003cb\u003eOil palm\u003c/b\u003e\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eIndonesia\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e1,944,456\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e46%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e1.8%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eMalaysia\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e\u003cp\u003e1,787,703\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e43%\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e2.2%\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003eConsidering sugar, wheat and oil palm production, for which all five indicators could be calculated, the greatest impact was from oil palm production in Malaysia (especially for GHG and LND), followed by oil palm production in Indonesia, then wheat in the UK, for which Nit was the key impact (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea). Sugar in the UK had the lowest total impact, most of which was attributable to Nit (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea), though sugar had a greater per-tonne impact than wheat in the UK (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003eb). For cocoa, only GHG and land use impact metrics were available. However, cocoa production showed much lower total impacts for these two metrics than the other three crops, with the highest total impact in Cote d\u0026rsquo;Ivoire (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea).\u003c/p\u003e\u003cp\u003eMore specifically, land-use impacts relating to Malaysian oil palm production were greatest, with a total standardised impact of 123, followed by Indonesian oil palm production at 74.5 (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea). Per tonne, while other crops followed overall trends, Indonesian oil palm had a substantially greater per-tonne impact than Malaysian oil palm (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003eb). Nit, LND, and GHG impacts were dominant in Indonesian oil palm impacts per tonne of crop (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003eb).\u003c/p\u003e\u003cp\u003eThe other crops showed a relatively low total land-use impact, with values between 0.005 and 1.9 (Ghanaian cocoa and UK wheat production respectively). GHGs were also a major contributor to crop environmental impacts; the greatest impacts again were from palm oil production, with Malaysian oil palm\u0026rsquo;s standardised total impact value at 74 and Indonesia\u0026rsquo;s at 42 (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea). GHG impact was a relatively strong contributor to UK wheat and sugar production impact too, with values of 9.7 for wheat and 1.1 for sugar (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea). GHG impacts from cocoa production across the four producing countries included here was relatively low, with values ranging between 0.02 and 0.21 (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea).\u003c/p\u003e\u003cp\u003eFor the three crops for which impact information from nitrogen fertiliser was available (oil palm, wheat and sugar), the greatest impacts were from UK wheat production (186), followed by Malaysian oil palm (66) (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea). Indonesian oil palm showed a much lower relative impact with a value of 11. For UK sugar production, nitrogen fertiliser impacts make up the greatest proportion of crop-related impact with a value of 14.7.\u003c/p\u003e\u003cp\u003eBoth phosphorus and water-use impacts for oil palm, wheat and sugar were much lower than the other indicator metrics, with values for phosphorus ranging between 0.0000089 and 0.00138635 and values for water between 0 (for UK sugar production) and 0.00559003 (Fig.\u0026nbsp;\u003cspan refid=\"Fig10\" class=\"InternalRef\"\u003e6\u003c/span\u003ea).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cb\u003eQuantifying the extent to which threatened species could be impacted by crop production in countries that supply the UK\u003c/b\u003e\u003c/p\u003e\u003cp\u003eThere is considerable overlap of palm oil and cocoa production with threatened species\u0026rsquo; ranges (species classified as Critically Endangered - CR, Endangered - EN, or Vulnerable - VU, by the IUCN. For Indonesian palm oil, the ranges of ten Critically Endangered species of mammal overlapped with Indonesian oil palm production (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003ea). The largest recorded overlap for this highest category of threat (CR) was for \u003cem\u003ePresbytis chrysomelas\u003c/em\u003e (6.3%), the Bornean banded langur, closely followed by \u003cem\u003eManis javanica\u003c/em\u003e the Sunda pangolin (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003ea). The overlap with threatened species (CR, EN and VU) was greatest in western Indonesia, though the distribution of overlap across the region was patchy (Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e8\u003c/span\u003ea). For cocoa, the most notable overlap of species ranges with crop production was in southern C\u0026ocirc;te d\u0026rsquo;Ivoire (Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e8\u003c/span\u003eb). \u003cem\u003ePiliocolobus waldroni - Miss Waldron's red colobus -\u003c/em\u003e had the highest overlap with areas of cocoa production here (~\u0026thinsp;18%) (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003eb), closely followed by \u003cem\u003eCercopithecus roloway\u003c/em\u003e - the Roloway Monkey- and both species ranked much higher for overlap than the four other critically endangered (CR) species (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003eb).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003eSeveral hotspots of wheat and sugarbeet production overlapped with the known range of \u003cem\u003eSciurus vulgaris\u003c/em\u003e (red squirrel; Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e8\u003c/span\u003e), which is a species already pressured by habitat loss, as well as competition with invasive species and disease and is classed as \u0026ldquo;near threatened\u0026rdquo; in the UK (IUCN, 2021; Mathews and Harrower, 2020). The overlap between sugar, wheat, and the red squirrel varied by crop, with hotspots for sugar (Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e8\u003c/span\u003ec) and a more widespread distribution for wheat (though East Anglia was a hotspot region for this crop; Fig.\u0026nbsp;\u003cspan refid=\"Fig13\" class=\"InternalRef\"\u003e8\u003c/span\u003ed).\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e"},{"header":"DISCUSSION","content":"\u003cp\u003ePrevious studies have shown that the impact of crop production can vary depending on where and how it is produced\u003csup\u003e8\u003c/sup\u003e. However, these impacts are rarely spatially explicit and infrequently contain estimates of how biodiversity will be impacted, often focusing more on indicators such as greenhouse-gas emissions, land-use change or water use\u003csup\u003e48\u0026ndash;50\u003c/sup\u003e. Not only does the impact matter for the sake of biodiversity, but because biodiversity is also critical for a sustainable and resilient agricultural system\u003csup\u003e51,52\u003c/sup\u003e. Here, we quantified and mapped the biodiversity impacts associated with the production of palm oil, cocoa, sugar, and wheat. We focus on the four crops in our study: i) given their dominance in ingredient lists for many luxury and \u0026lsquo;treat\u0026rsquo; foods; ii) to give a first estimation of biodiversity impacts; and iii) to emphasise the need for biodiversity and broader, complementary environmental impact data in sustainability assessments of food products and their ingredients. Using UK biscuit consumption as a case study, we identified the global-scale biodiversity impacts associated with the crops grown to produce key ingredients. We show that the demand for food including these crops has strong, varied, and wide-reaching environmental impacts that vary spatially and according to country of origin. By assessing the spatial variation in several biodiversity-based impact indicators (i.e. not focusing solely on forest loss\u003csup\u003e53\u003c/sup\u003e), we were also explore the impact on threatened species and suggest this as a means of communicating biodiversity impacts to consumers.\u003c/p\u003e\u003cp\u003eIn line with previous work\u003csup\u003e54,55\u003c/sup\u003e, we identify strong variability in impacts according to crop, region and type of indicator. This highlights the importance of assessing multiple metrics for evaluating the biodiversity impacts of different foods. For instance, we show a distinction between hotspots of environmental impact and areas of strong overlap between harvested area and high species richness (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003e). This emphasises a need for complementary information, such as the amounts of a given crop being sourced from an area via international trade (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), and its overlap with threatened species\u0026rsquo; ranges (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003e), to investigate how impacts vary spatially within and among regions. At this time, the coarse resolution of the impact maps we were able to produce at the global scale makes it challenging to identify specific areas at the local scale where impacts may be particularly high or low. This makes providing conservation recommendations challenging, highlighting the need for finer-resolution data.\u003c/p\u003e\u003cp\u003eNevertheless, we found that a number of Critically Endangered species are at risk from habitat loss via agricultural development in countries where the UK sources the focal crops (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003e). Indeed, the impacts of the two imported crops - cocoa and oil palm - are strongest in species-rich tropical regions (Figs.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e4\u003c/span\u003e \u0026amp; \u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003e), particularly Southeast Asia and Central/West Africa. For oil palm, the largest proportion of overlap of harvested area and species range was for \u003cem\u003ePresbytis chrysomelas\u003c/em\u003e at over 6% followed by \u003cem\u003eManis javanica\u003c/em\u003e at 5.6%. However, a greater number of threatened species ranges overlapped with crop production areas for this crop in Indonesia (\u003cb\u003eSupporting Information Table S3\u003c/b\u003e). While Indonesian orangutans (\u003cem\u003ePongo pygmaeus\u003c/em\u003e and \u003cem\u003ePongo abelii\u003c/em\u003e) are often the \u0026lsquo;face\u0026rsquo; of oil palm's impact - \u0026lsquo;icons of wild nature\u0026rsquo; according to Spehar et al.\u003csup\u003e56\u003c/sup\u003e - our findings reveal that lesser-known species like the \u003cem\u003ePresbytis chrysomelas\u003c/em\u003e and the Sunda pangolin are facing significant threats, urging a broader awareness of biodiversity loss. \u003cem\u003ePresbytis\u003c/em\u003e primates are relatively understudied, but the \u003cem\u003echrysomelas\u003c/em\u003e species we identify as having strong overlap with oil palm is endemic to Borneo and is losing its forest habitat to logging and conversion to agriculture. This is of particular concern because these monkeys cannot survive in deforested habitats and are unlikely to cross between forest patches\u003csup\u003e57\u003c/sup\u003e. Previous research into the distribution of \u003cem\u003eManis javanica\u003c/em\u003e (the Sunda pangolin) identified forested areas of importance across Sabah in Borneo. Whilst the Sunda pangolin can adapt to secondary forests, cultivated, and urban areas\u003csup\u003e58,59\u003c/sup\u003e, the oil palm mapped in Indonesia (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003e) could put pangolins at risk due to habitat fragmentation, which increases the threat of poaching by improving accessibility to humans\u003csup\u003e58,60\u003c/sup\u003e. Whilst our study cannot be used to accurately pinpoint overlap between pangolins and oil palm, due to the data-imposed assumption we must make that pangolins are present across their entire range, we can take from this that any overlap has a potential for direct impact via habitat loss, but also indirect effects via accessibility to humans. Recent studies further support this, highlighting that even palm oil production labelled as \u0026lsquo;sustainable\u0026rsquo; has been associated with widespread habitat loss for orangutans and other species (e.g. see Gatti and Velichevskaya\u003csup\u003e61\u003c/sup\u003e).\u003c/p\u003e\u003cp\u003eFor cocoa, we identified two primate species for which almost 20% of their range overlapped with crop production - \u003cem\u003eCercopithecus roloway\u003c/em\u003e and \u003cem\u003ePiliocolobus waldroni.\u003c/em\u003e Whilst certain species can adapt to changing landscapes and coexist alongside agricultural production (e.g., generalist species and species able to survive in smaller patches\u003csup\u003e56,62\u003c/sup\u003e), many cannot, and around 90% of Critically Endangered mammals have agricultural land listed as a major driver of their extinction risk\u003csup\u003e63\u003c/sup\u003e. This is the case for Miss Waldron\u0026rsquo;s red colobus (\u003cem\u003ePiliocolobus waldroni\u003c/em\u003e) and the Roloway monkey (\u003cem\u003eCercopithecus roloway\u003c/em\u003e)\u003csup\u003e64,65\u003c/sup\u003e. Approximately 17% of tree cover has been lost in Cote d\u0026rsquo;Ivoire in recent decades\u003csup\u003e66,67\u003c/sup\u003e. This has been exacerbated by a large increase in human population and the spread of cocoa farming, including illegal cocoa cultivation within protected areas\u003csup\u003e65,68\u003c/sup\u003e. The red colobus monkey, which had the greatest overlap with cocoa production in Cote d\u0026rsquo;Ivoire of all critically endangered species analysed (Fig.\u0026nbsp;\u003cspan refid=\"Fig12\" class=\"InternalRef\"\u003e7\u003c/span\u003e), has been noted as potentially extinct in previous years and is endemic only to the two major cocoa producers - Cote d\u0026rsquo;Ivoire and Ghana\u003csup\u003e69\u003c/sup\u003e. These case studies reveal the \u0026lsquo;hidden\u0026rsquo; impacts of global trade in producer countries, often unknown to consumers purchasing the products and difficult for sellers and researchers to quantify due to lengthy and opaque supply chains\u003csup\u003e70\u003c/sup\u003e.\u003c/p\u003e\u003cp\u003eAdding to the complexity, our findings suggest UK biscuit consumption has a potentially high contribution to environmental impacts overseas, but not all of the crops we evaluated are used in biscuits. A study investigating the variation of sugar and calorie content in cakes and biscuits sold in the UK in 2016 estimated that the average sugar content is 30/100g\u003csup\u003e71\u003c/sup\u003e. As a case-study example, this is 4.4g per biscuit for dark chocolate digestive biscuits. Based on our estimate of UK biscuit consumption (see \u003cb\u003eIntroduction\u003c/b\u003e, calculated from public listing by McVities and data from statista.com and ONS.org.uk), this equates to over 52,000 tonnes of sugar, which is equivalent to 0.7% of total sugar produced in the UK. Similarly, for wheat, per biscuit, 96,000 tonnes of wheat is roughly equivalent to 0.7% of total wheat production. Quantities of palm oil in biscuits are generally undeclared. However, for digestive biscuits, roughly 5g of chocolate per biscuit means approximately 60,000 tonnes of chocolate is consumed in biscuits alone per year. This is about a third of the total imported cocoa (assuming ratios of cocoa:chocolate are 1:1). Although generally contributing relatively small amounts to the demand for these crops, with the exception of cocoa, biscuits and other sweet treats within which these crops are staples could provide an opportunity for individuals to make alternative choices to reduce demand for highly impactful crops. For instance, consumers could reduce their consumption of sweet treats or choose options with a potential for lower impact (e.g., an option without a chocolate coating, a different oil, or reduced sugar). Nevertheless, due to the complexities discussed, there is not always a clear \u0026lsquo;optimal\u0026rsquo; consumer choice.\u003c/p\u003e\u003cp\u003eThere remains, for instance, a lack of transparency in global supply chains, especially those used by businesses creating products like biscuits on an industrial scale. Herein lies a major issue with quantifying and communicating biodiversity impacts associated with commonly consumed food products. Without improved supply-chain transparency, the case studies and methods used in our study provide a means to estimate biodiversity impact but have an unavoidable level of uncertainty. For example, we found that the UK sources some of the focal crops from relatively low-impact regions. The proportion of a partner country\u0026rsquo;s production of the focal crops supplied to the UK are low in most cases (except cocoa from Ghana, which represents 9% of the country\u0026rsquo;s production) (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). This could lead us to believe that the UK\u0026rsquo;s environmental impact on other countries is low. However, the reality is more complex. For instance, whilst the UK\u0026rsquo;s cocoa imports represent only\u0026thinsp;~\u0026thinsp;5% of Cote d\u0026rsquo;Ivoire\u0026rsquo;s overall cocoa production (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), we do not know from \u003cem\u003ewhere\u003c/em\u003e the UK is sourcing the cocoa \u003cem\u003ewithin\u003c/em\u003e Cote d\u0026rsquo;Ivoire and from which type of farming system. There are attempts to model these supply chains in finer detail (e.g., the TRASE\u003csup\u003e72\u003c/sup\u003e initiative for tracing deforestation and commodity trade pathways) but the data are limited. Additionally, traceability often only extends to distribution warehouses rather than the smallholder farmers that supply them. Thus, the 5% sourced from Cote d\u0026rsquo;Ivoire could mean UK demand is responsible for a relatively small proportion of Cote d\u0026rsquo;Ivoire\u0026rsquo;s potential biodiversity impacts associated with cocoa or a much larger proportion, hidden behind a lack of farm-specific or smaller spatial scale origin-tracing data\u003csup\u003e73\u003c/sup\u003e. This particular issue emphasises the importance of spatialising measures of environmental impact, to better capture the origin and endpoint of impacts.\u003c/p\u003e\u003cp\u003eAlso impacting our understanding of the relative sustainability of different food products is differing baselines of biodiversity records. Although the UK has no threatened mammal species according to the IUCN Red List, this does not mean that there are no threats to UK biodiversity. Due to the country\u0026rsquo;s history of industrialization and agricultural intensification, it is likely that threatened mammalian biodiversity was lost before Red List data were collected. This bias is present in many countries with long histories of intense land-use change and influences our understanding of environmental impacts, which appear stronger in the tropics where biodiversity is higher (e.g., as discussed in Ellis et al.\u003csup\u003e74\u003c/sup\u003e). This highlights the potential for injustice when making sustainable land-use decisions: do we source ingredients from places which have already gained economically via habitat loss from agricultural development rather than potentially less economically developed regions that have higher biodiversity?\u003c/p\u003e\u003cp\u003eFurthermore, we assume that a species\u0026rsquo; range overlapping with an area of high crop production is placing that species at risk. This is based on the fact that many species on the IUCN Red List are listed as threatened due to cropland expansion and habitat loss for agriculture, and that previous research has shown lower biodiversity in cropland than natural habitats\u003csup\u003e75,76\u003c/sup\u003e. Importantly, there are large differences in the production systems associated with the crops examined in this study. For instance, cocoa can be grown on small plots alongside other vegetation (shaded cocoa) which can enhance tree cover and ecosystem services\u003csup\u003e77\u003c/sup\u003e(though note this depends on shade tree type, land-use history, and other characteristics - e.g., Blaser-Hart et al.\u003csup\u003e78\u003c/sup\u003e and Maney et al.\u003csup\u003e79\u003c/sup\u003e), whilst oil palm tends to be grown in vast areas of monoculture. These distinctions are lost in global-scale analyses such as ours. Unfortunately, only data for LND and GHG impacts were available for cocoa. As a result, we are underestimating the total environmental impacts of this crop, especially as fertilisers are recommended by Ghana\u0026rsquo;s cocoa board and associated with yield improvements, for example\u003csup\u003e80\u003c/sup\u003e. According to water footprinting of cocoa grown in Colombia, some cocoa varieties are more water-efficient than others, providing greater yields with the same quantities of rainfall\u003csup\u003e81\u003c/sup\u003e, highlighting a future research need for variety-specific biodiversity-impact measurement, as well as crop-specific studies. We also focus part of our analysis on threatened mammals specifically, since this taxonomic group has the best available data. However, other forms of biodiversity will also be impacted by crop production and we recommend that future work considers multiple groups and various facets of biodiversity when assessing environmental impacts. Whilst beyond the scope of our study, we anticipate that several of the metrics we used could be developed to incorporate more facets of biodiversity (e.g., functional and phylogenetic diversity).\u003c/p\u003e\u003cp\u003eAs well as impacting specific threatened species, we show that some of the UK\u0026rsquo;s biscuit ingredients are being sourced from some of the world's most environmentally impacted hotspots, particularly for oil palm and cocoa, where production, impacts, and overlap with species richness are greatest in the tropics (Figs.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e5\u003c/span\u003e). However, our analysis of products consumed in the UK highlights complexities inherent in evaluating sustainability of snacks like digestive biscuits. Sustainably sourced food is increasingly high on the agenda of consumers in countries including the UK, where there is an increasing trend towards plant-based diets\u003csup\u003e82\u003c/sup\u003e, one of the recommendations for planetary health in the EAT-Lancet Commission report\u003csup\u003e12\u003c/sup\u003e. However, making diets more sustainable comes with several challenges, such as the cost, political economy, cultural norms, justice, and governance issues highlighted by B\u0026eacute;n\u0026eacute; and Abdulai\u003csup\u003e83\u003c/sup\u003e. Furthermore, the complexity inherent in ensuring the overall environmental sustainability of products, particularly when considering biodiversity, makes it difficult to simply and effectively communicate the environmental impacts of food production to consumers. In the future, when more field data on species' responses to specific crops are available and supply chains are more transparent, consumers could be informed of the specific impacts of each dietary decision more easily. For now, by mapping crop-specific environmental impact hotspots and highlighting vulnerable species at risk according to their range overlap with areas of high levels of food crop production, we anticipate that this, and future studies expanding the scope to other commonly consumed food items, will help to reduce the disconnect between consumers and their biodiversity impacts beyond packaging. To date, there has been a consumer focus on reducing plastic waste to protect biodiversity. We propose that further education and industry initiatives are required to build on this established momentum towards more sustainable practices and expand the focus to ingredient origins and how to make more sustainable food choices.\u003c/p\u003e\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e\u003ch2\u003eDATA AVAILABILITY STATEMENT\u003c/h2\u003e\u003cp\u003eThe data used as inputs for the analyses in this manuscript were openly available for download, as described and cited in the manuscript. The code used for analyses are available from: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://github.com/CharlieOuthwaite/Biscuit_project\u003c/span\u003e\u003cspan address=\"https://github.com/CharlieOuthwaite/Biscuit_project\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e.\u003c/p\u003e\u003c/div\u003e"},{"header":"Declarations","content":"\u003ch3\u003e\u003cstrong\u003eACKNOWLEDGEMENTS\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eWe gratefully acknowledge the participants of workshop and/or \u0026lsquo;hackathon\u0026rsquo; sessions for the initial discussion of the topic (Georgina Adams, Paolo Agnolucci, Elizabeth Boakes, Sara Bonetti, \u0026Aacute;lvaro Calzadilla Rivera,\u0026nbsp;Silvia Ceausu, Henry Ferguson-Gow, Laura Graham, Jonathan Green, Vivienne Groner, Rob Holland, Jacquie McGlade, Monica Ortiz, Richard Pearson, Chrysanthi Rapti, Ralf Seppelt, Pete Smith, Hemant Tripathi, Peter Verburg, and Jessica Williams).\u0026nbsp;\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eFUNDING\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eThis work was supported by \u0026lsquo;UN SDGs: Pathways to Achievement\u0026rsquo; funding from UCL Grand Challenges and Global Engagement Office. CLO was supported by a UK Natural Environment Research Council grant NE/R010811/1 and funding from Research England. ASAC acknowledges funding from the Wellcome Trust (Our Planet, Our Health programme, 205200, https://doi.org/10.35802/205200 - Sustainable and Healthy Food Systems (SHEFS) programme; and \u0026nbsp;Climate and Health Programme, 227749/Z/23/Z - Sustainable and Healthy Food Systems - Southern Africa (SHEFS-SA)). FP was funded by the EASTBIO DTP. AM was funded by the Global Challenges Research Fund Trade, Development and the Environment Hub project (ES/S008160/1). C.D. also acknowledges funding from the Natural Environment Research Council Fellowship (NERC NE/N01524X/1) and from the European Union (ERC, FLORA, StG 101039402).\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eAUTHOR CONTRIBUTIONS STATEMENT\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eCLO \u0026amp; ASAC applied for the funding. FP, CLO, AM \u0026amp; ASAC conceptualised the initial project idea, contributed equally to the writing of the manuscript, and performed the analysis. ASAC coordinated the project. MJ \u0026amp; CD provided parts of the data and gave extensive comments on the manuscript. All authors have seen the final manuscript and take responsibility for its contents.\u003c/p\u003e\n\u003ch3\u003e\u003cstrong\u003eADDITIONAL INFORMATION\u003c/strong\u003e\u003c/h3\u003e\n\u003ch3\u003e\u003cstrong\u003eCOMPETING INTERESTS STATEMENT\u003c/strong\u003e\u003c/h3\u003e\n\u003cp\u003eThe authors declare no competing interests.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eMaxwell, S. L., Fuller, R. A., Brooks, T. M. \u0026amp; Watson, J. E. M. Biodiversity: The ravages of guns, nets and bulldozers. \u003cem\u003eNature\u003c/em\u003e \u003cstrong\u003e536\u003c/strong\u003e, 143-145, doi:10.1038/536143a (2016).\u003c/li\u003e\n\u003cli\u003eRitchie, H. \u0026amp; Roser, M. \u003cem\u003eHalf of the world\u0026rsquo;s habitable land is used for agriculture\u003c/em\u003e (Published online at OurWorldinData.org, 2019).\u003c/li\u003e\n\u003cli\u003eUnited Nations Department of Economic and Social Affairs. \u003cem\u003eStatistical Papers - United Nations (Ser. A), Population and Vital Statistics Report\u003c/em\u003e 376 (2017).\u003c/li\u003e\n\u003cli\u003eDelzeit, R., Zabel, F., Meyer, C. \u0026amp; V\u0026aacute;clav\u0026iacute;k, T. Addressing future trade-offs between biodiversity and cropland expansion to improve food security. \u003cem\u003eRegional Environmental Change\u003c/em\u003e \u003cstrong\u003e17\u003c/strong\u003e, 1429-1441, doi:10.1007/s10113-016-0927-1 (2017).\u003c/li\u003e\n\u003cli\u003eFan, S., Headey, D., Rue, C. \u0026amp; Thomas, T. Food Systems for Human and Planetary Health: Economic Perspectives and Challenges. \u003cem\u003eAnnual Review of Resource Economics\u003c/em\u003e \u003cstrong\u003e13\u003c/strong\u003e, 131-156, doi:https://doi.org/10.1146/annurev-resource-101520-081337 (2021).\u003c/li\u003e\n\u003cli\u003eCanavan, C. R., Noor, R. A., Golden, C. D., Juma, C. \u0026amp; Fawzi, W. Sustainable food systems for optimal planetary health. \u003cem\u003eTransactions of The Royal Society of Tropical Medicine and Hygiene\u003c/em\u003e \u003cstrong\u003e111\u003c/strong\u003e, 238-240, doi:10.1093/trstmh/trx038 (2017).\u003c/li\u003e\n\u003cli\u003eHalpern, B. S.\u003cem\u003e et al.\u003c/em\u003e The environmental footprint of global food production. \u003cem\u003eNature Sustainability\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, 1027-1039, doi:10.1038/s41893-022-00965-x (2022).\u003c/li\u003e\n\u003cli\u003ePoore, J. \u0026amp; Nemecek, T. Reducing food\u0026rsquo;s environmental impacts through producers and consumers. \u003cem\u003eScience\u003c/em\u003e \u003cstrong\u003e360\u003c/strong\u003e, 987-992, doi:10.1126/science.aaq0216 (2018).\u003c/li\u003e\n\u003cli\u003eJwaideh, M. A. A. \u0026amp; Dalin, C. The multi-dimensional environmental impact of global crop commodities. \u003cem\u003eNature Sustainability\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 396-410, doi:10.1038/s41893-025-01528-6 (2025).\u003c/li\u003e\n\u003cli\u003eClark, M.\u003cem\u003e et al.\u003c/em\u003e Estimating the environmental impacts of 57,000 food products. \u003cem\u003eProceedings of the National Academy of Sciences\u003c/em\u003e \u003cstrong\u003e119\u003c/strong\u003e, e2120584119, doi:10.1073/pnas.2120584119 (2022).\u003c/li\u003e\n\u003cli\u003eHarris, F.\u003cem\u003e et al.\u003c/em\u003e The Water Footprint of Diets: A Global Systematic Review and Meta-analysis. \u003cem\u003eAdv Nutr\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 375-386, doi:10.1093/advances/nmz091 (2020).\u003c/li\u003e\n\u003cli\u003eWillett, W.\u003cem\u003e et al.\u003c/em\u003e Food in the Anthropocene: the EAT-Lancet Commission on healthy diets from sustainable food systems. \u003cem\u003eThe Lancet\u003c/em\u003e \u003cstrong\u003e393\u003c/strong\u003e, 447-492, doi:10.1016/S0140-6736(18)31788-4 (2019).\u003c/li\u003e\n\u003cli\u003eChaudhary, A. \u0026amp; Kastner, T. Land use biodiversity impacts embodied in international food trade. \u003cem\u003eGlobal Environmental Change\u003c/em\u003e \u003cstrong\u003e38\u003c/strong\u003e, 195-204, doi:https://doi.org/10.1016/j.gloenvcha.2016.03.013 (2016).\u003c/li\u003e\n\u003cli\u003eUnited Nations General Assembly. (2015).\u003c/li\u003e\n\u003cli\u003eMehrabi, Z., Ellis, E. C. \u0026amp; Ramankutty, N. The challenge of feeding the world while conserving half the planet. \u003cem\u003eNature Sustainability\u003c/em\u003e \u003cstrong\u003e1\u003c/strong\u003e, 409-412, doi:10.1038/s41893-018-0119-8 (2018).\u003c/li\u003e\n\u003cli\u003eGive Food. \u003cem\u003eTrussell Trust Most Requested Items\u003c/em\u003e, \u0026lt;https://www.givefood.org.uk/dashboard/trusselltrust/most-requested-items/\u0026gt; (2025).\u003c/li\u003e\n\u003cli\u003eMintel. Uk Sweet Biscuits Market Report 2024. (2024).\u003c/li\u003e\n\u003cli\u003eKonstantas, A., Stamford, L. \u0026amp; Azapagic, A. Evaluation of environmental sustainability of biscuits at the product and sectoral levels. \u003cem\u003eJournal of Cleaner Production\u003c/em\u003e \u003cstrong\u003e230\u003c/strong\u003e, 1217-1228, doi:https://doi.org/10.1016/j.jclepro.2019.05.095 (2019).\u003c/li\u003e\n\u003cli\u003eKonstantas, A., Jeswani, H. K., Stamford, L. \u0026amp; Azapagic, A. Environmental impacts of chocolate production and consumption in the UK. \u003cem\u003eFood Research International\u003c/em\u003e \u003cstrong\u003e106\u003c/strong\u003e, 1012-1025, doi:https://doi.org/10.1016/j.foodres.2018.02.042 (2018).\u003c/li\u003e\n\u003cli\u003eCarlson, K. M.\u003cem\u003e et al.\u003c/em\u003e Greenhouse gas emissions intensity of global croplands. \u003cem\u003eNature Climate Change\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, 63-68, doi:10.1038/nclimate3158 (2017).\u003c/li\u003e\n\u003cli\u003eTuninetti, M., Tamea, S. \u0026amp; Dalin, C. Water Debt Indicator Reveals Where Agricultural Water Use Exceeds Sustainable Levels. \u003cem\u003eWater Resources Research\u003c/em\u003e \u003cstrong\u003e55\u003c/strong\u003e, 2464-2477, doi:https://doi.org/10.1029/2018WR023146 (2019).\u003c/li\u003e\n\u003cli\u003eJwaideh, M. A. A., Sutanudjaja, E. H. \u0026amp; Dalin, C. Global impacts of nitrogen and phosphorus fertiliser use for major crops on aquatic biodiversity. \u003cem\u003eThe International Journal of Life Cycle Assessment\u003c/em\u003e \u003cstrong\u003e27\u003c/strong\u003e, 1058-1080, doi:10.1007/s11367-022-02078-1 (2022).\u003c/li\u003e\n\u003cli\u003eIUCN. (2021).\u003c/li\u003e\n\u003cli\u003eChaudhary, A., Verones, F., de Baan, L. \u0026amp; Hellweg, S. Quantifying Land Use Impacts on Biodiversity: Combining Species\u0026ndash;Area Models and Vulnerability Indicators. \u003cem\u003eEnvironmental Science \u0026amp; Technology\u003c/em\u003e \u003cstrong\u003e49\u003c/strong\u003e, 9987-9995, doi:10.1021/acs.est.5b02507 (2015).\u003c/li\u003e\n\u003cli\u003eggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics v. 3.5.3 (2025).\u003c/li\u003e\n\u003cli\u003eraster: Geographic Data Analysis and Modeling v. R package version 3.6-23 (2023).\u003c/li\u003e\n\u003cli\u003eR: A Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, Vienna, Austria, 2023).\u003c/li\u003e\n\u003cli\u003eMonfreda, C., Ramankutty, N. \u0026amp; Foley, J. A. Farming the planet: 2. Geographic distribution of crop areas, yields, physiological types, and net primary production in the year 2000. \u003cem\u003eGlobal Biogeochemical Cycles\u003c/em\u003e \u003cstrong\u003e22\u003c/strong\u003e, n/a-n/a, doi:10.1029/2007gb002947 (2008).\u003c/li\u003e\n\u003cli\u003eMueller, N. D.\u003cem\u003e et al.\u003c/em\u003e Closing yield gaps through nutrient and water management. \u003cem\u003eNature\u003c/em\u003e \u003cstrong\u003e490\u003c/strong\u003e, 254-257, doi:10.1038/nature11420 (2012).\u003c/li\u003e\n\u003cli\u003eBouwman, A. F.\u003cem\u003e et al.\u003c/em\u003e Lessons from temporal and spatial patterns in global use of N and P fertilizer on cropland. \u003cem\u003eScientific Reports\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, 40366, doi:10.1038/srep40366 (2017).\u003c/li\u003e\n\u003cli\u003eInternational Food Policy Research, I. \u003cem\u003eMapSPAM\u003c/em\u003e (ed Institute International Food Policy Research) (Harvard Dataverse, 2019).\u003c/li\u003e\n\u003cli\u003eYu, Q.\u003cem\u003e et al.\u003c/em\u003e A cultivated planet in 2010 \u0026ndash; Part 2: The global gridded agricultural-production maps. \u003cem\u003eEarth Syst. Sci. Data\u003c/em\u003e \u003cstrong\u003e12\u003c/strong\u003e, 3545-3572, doi:10.5194/essd-12-3545-2020 (2020).\u003c/li\u003e\n\u003cli\u003eEtard, A., Morrill, S. \u0026amp; Newbold, T. Global gaps in trait data for terrestrial vertebrates. \u003cem\u003eGlobal Ecology and Biogeography\u003c/em\u003e \u003cstrong\u003e29\u003c/strong\u003e, 2143-2158, doi:https://doi.org/10.1111/geb.13184 (2020).\u003c/li\u003e\n\u003cli\u003eBirdlife International \u0026amp; NatureServe. (See http://www.birdlife.org/datazone/info/spcdownload., 2012).\u003c/li\u003e\n\u003cli\u003eIUCN. (See http://www.iucnredlist.org/, 2013).\u003c/li\u003e\n\u003cli\u003edplyr: A Grammar of Data Manipulation (2023).\u003c/li\u003e\n\u003cli\u003etidyr: Tidy Messy Data (2023).\u003c/li\u003e\n\u003cli\u003eclassInt: Choose Univariate Class Intervals (2023).\u003c/li\u003e\n\u003cli\u003emagrittr: A Forward-Pipe Operator for R (2022).\u003c/li\u003e\n\u003cli\u003epatchwork: The Composer of Plots v. 1.3.0 (2024).\u003c/li\u003e\n\u003cli\u003eKastner, T., Kastner, M. \u0026amp; Nonhebel, S. Tracing distant environmental impacts of agricultural products from a consumer perspective. \u003cem\u003eEcological Economics\u003c/em\u003e \u003cstrong\u003e70\u003c/strong\u003e, 1032-1040, doi:https://doi.org/10.1016/j.ecolecon.2011.01.012 (2011).\u003c/li\u003e\n\u003cli\u003eDalin, C., Wada, Y., Kastner, T. \u0026amp; Puma, M. J. Groundwater depletion embedded in international food trade. \u003cem\u003eNature\u003c/em\u003e \u003cstrong\u003e543\u003c/strong\u003e, 700-704, doi:10.1038/nature21403 (2017).\u003c/li\u003e\n\u003cli\u003eIUCN. (See http://www.iucnredlist.org/, 2012).\u003c/li\u003e\n\u003cli\u003eGorelick, N.\u003cem\u003e et al.\u003c/em\u003e Google Earth Engine: Planetary-scale geospatial analysis for everyone. \u003cem\u003eRemote Sensing of Environment\u003c/em\u003e \u003cstrong\u003e202\u003c/strong\u003e, 18-27, doi:https://doi.org/10.1016/j.rse.2017.06.031 (2017).\u003c/li\u003e\n\u003cli\u003eMathews, F. \u0026amp; Harrower, C. IUCN \u0026ndash; compliant Red List for Britain\u0026rsquo;s Terrestrial Mammals. Assessment by the Mammal Society under contract to Natural England, Natural Resources Wales and Scottish Natural Heritage. . (Natural England, Peterborough, 2020).\u003c/li\u003e\n\u003cli\u003emaps: Draw Geographical Maps v. R package version 3.4.2 (2023).\u003c/li\u003e\n\u003cli\u003eArcGIS Pro 3.5 (ESRI, Redlands, CA, 2025).\u003c/li\u003e\n\u003cli\u003eDalin, C. \u0026amp; Rodr\u0026iacute;guez-Iturbe, I. Environmental impacts of food trade via resource use and greenhouse gas emissions. \u003cem\u003eEnvironmental Research Letters\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 035012, doi:10.1088/1748-9326/11/3/035012 (2016).\u003c/li\u003e\n\u003cli\u003eMolotoks, A., Smith, P. \u0026amp; Dawson, T. P. Impacts of land use, population, and climate change on global food security. \u003cem\u003eFood and Energy Security\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, e261, doi:https://doi.org/10.1002/fes3.261 (2021).\u003c/li\u003e\n\u003cli\u003eDinar, A., Tieu, A. \u0026amp; Huynh, H. Water scarcity impacts on global food production. \u003cem\u003eGlobal Food Security\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 212-226, doi:https://doi.org/10.1016/j.gfs.2019.07.007 (2019).\u003c/li\u003e\n\u003cli\u003eElouafi, I. Why biodiversity matters in agriculture and food systems. \u003cem\u003eScience\u003c/em\u003e \u003cstrong\u003e386\u003c/strong\u003e, eads8197, doi:10.1126/science.ads8197 (2024).\u003c/li\u003e\n\u003cli\u003eDardonville, M., Bockstaller, C., Villerd, J. \u0026amp; Therond, O. Resilience of agricultural systems: biodiversity-based systems are stable, while intensified ones are resistant and high-yielding. \u003cem\u003eAgricultural Systems\u003c/em\u003e \u003cstrong\u003e197\u003c/strong\u003e, 103365, doi:https://doi.org/10.1016/j.agsy.2022.103365 (2022).\u003c/li\u003e\n\u003cli\u003eMolotoks, A., Green, J., Ribeiro, V., Wang, Y. \u0026amp; West, C. Assessing the value of biodiversity-specific footprinting metrics linked to South American soy trade. \u003cem\u003ePeople and Nature\u003c/em\u003e \u003cstrong\u003en/a\u003c/strong\u003e, doi:https://doi.org/10.1002/pan3.10457 (2023).\u003c/li\u003e\n\u003cli\u003eChaudhary, A., Pfister, S. \u0026amp; Hellweg, S. Spatially Explicit Analysis of Biodiversity Loss Due to Global Agriculture, Pasture and Forest Land Use from a Producer and Consumer Perspective. \u003cem\u003eEnviron Sci Technol\u003c/em\u003e \u003cstrong\u003e50\u003c/strong\u003e, 3928-3936, doi:10.1021/acs.est.5b06153 (2016).\u003c/li\u003e\n\u003cli\u003eChapman, A. S. A.\u003cem\u003e et al.\u003c/em\u003e Quantifying the biodiversity pressures of fruit and vegetable consumption in the United Kingdom, India, and South Africa. \u003cem\u003eNature Food\u003c/em\u003e (2025).\u003c/li\u003e\n\u003cli\u003eSpehar, S. N.\u003cem\u003e et al.\u003c/em\u003e Orangutans venture out of the rainforest and into the Anthropocene. \u003cem\u003eScience Advances\u003c/em\u003e \u003cstrong\u003e4\u003c/strong\u003e, e1701422, doi:doi:10.1126/sciadv.1701422 (2018).\u003c/li\u003e\n\u003cli\u003eEhlers Smith, D. A. The effects of land-use policies on the conservation of Borneo\u0026rsquo;s endemic Presbytis monkeys. \u003cem\u003eBiodiversity and Conservation\u003c/em\u003e \u003cstrong\u003e23\u003c/strong\u003e, 891-908, doi:10.1007/s10531-014-0639-0 (2014).\u003c/li\u003e\n\u003cli\u003eNursamsi, I., Amir, Z., Decoeur, H., Moore, J. H. \u0026amp; Luskin, M. S. Sunda pangolins show inconsistent responses to disturbances across multiple scales. \u003cem\u003eWildlife Letters\u003c/em\u003e \u003cstrong\u003e1\u003c/strong\u003e, 59-70, doi:https://doi.org/10.1002/wll2.12010 (2023).\u003c/li\u003e\n\u003cli\u003eChong, J. L.\u003cem\u003e et al.\u003c/em\u003e in \u003cem\u003ePangolins\u003c/em\u003e (eds Daniel W. S. Challender, Helen C. Nash, \u0026amp; Carly Waterman) 89-108 (Academic Press, 2020).\u003c/li\u003e\n\u003cli\u003ePanjang, E.\u003cem\u003e et al.\u003c/em\u003e Mapping the distribution of the Sunda pangolin (Manis javanica) within natural forest in Sabah, Malaysian Borneo. \u003cem\u003eGlobal Ecology and Conservation\u003c/em\u003e \u003cstrong\u003e52\u003c/strong\u003e, e02962, doi:https://doi.org/10.1016/j.gecco.2024.e02962 (2024).\u003c/li\u003e\n\u003cli\u003eCazzolla Gatti, R. \u0026amp; Velichevskaya, A. Certified \u0026ldquo;sustainable\u0026rdquo; palm oil took the place of endangered Bornean and Sumatran large mammals habitat and tropical forests in the last 30 years. \u003cem\u003eScience of The Total Environment\u003c/em\u003e \u003cstrong\u003e742\u003c/strong\u003e, 140712, doi:https://doi.org/10.1016/j.scitotenv.2020.140712 (2020).\u003c/li\u003e\n\u003cli\u003eAncrenaz, M.\u003cem\u003e et al.\u003c/em\u003e Of Pongo, palms and perceptions: a multidisciplinary assessment of Bornean orang-utans Pongo pygmaeus in an oil palm context. \u003cem\u003eOryx\u003c/em\u003e \u003cstrong\u003e49\u003c/strong\u003e, 465-472, doi:10.1017/S0030605313001270 (2015).\u003c/li\u003e\n\u003cli\u003eRipple, W. J.\u003cem\u003e et al.\u003c/em\u003e Extinction risk is most acute for the world\u0026rsquo;s largest and smallest vertebrates. \u003cem\u003eProceedings of the National Academy of Sciences\u003c/em\u003e \u003cstrong\u003e114\u003c/strong\u003e, 10678-10683, doi:10.1073/pnas.1702078114 (2017).\u003c/li\u003e\n\u003cli\u003eMcGraw, W. S. \u0026amp; Oates, J. F. in \u003cem\u003ePrimates in Peril: The World\u0026apos;s 25 Most Endangered Primates 2014-2016 \u003c/em\u003e(eds C. Schwitzer\u003cem\u003e et al.\u003c/em\u003e) 14-16 (IUCN SSC Primate Specialist Group (PSG), International Primatological Society (IPS), Conservation International (CI), and Bristol Zoological Society, 2015).\u003c/li\u003e\n\u003cli\u003eBitty, E. A., Bi, S. G., Bene, J.-C. K., Kouassi, P. K. \u0026amp; McGraw, W. S. Cocoa Farming and Primate Extirpation Inside Cote D\u0026apos;ivoire\u0026apos;s Protected Areas. \u003cem\u003eTropical Conservation Science\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 95-113, doi:10.1177/194008291500800110 (2015).\u003c/li\u003e\n\u003cli\u003eChatelain, C., Dao, H., Gautier, L. \u0026amp; Spichiger, R. Forest cover changes in C\u0026ocirc;te d\u0026rsquo;Ivoire and Upper Guinea. \u003cem\u003eCABI\u003c/em\u003e, 15\u0026ndash;32, doi:10.1079/9780851997346.0015 (2004).\u003c/li\u003e\n\u003cli\u003eGlobal Forest Watch. \u0026lt;www.globalforestwatch.org\u0026gt; (2019).\u003c/li\u003e\n\u003cli\u003eBarima, Y. S. S.\u003cem\u003e et al.\u003c/em\u003e Cocoa crops are destroying the forest reserves of the classified forest of Haut-Sassandra (Ivory Coast). \u003cem\u003eGlobal Ecology and Conservation\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, 85-98, doi:https://doi.org/10.1016/j.gecco.2016.08.009 (2016).\u003c/li\u003e\n\u003cli\u003eLinder, J. M.\u003cem\u003e et al.\u003c/em\u003e Red colobus (Piliocolobus) conservation action plan 2021\u0026ndash;2026. (IUCN, Gland, Switzerland, 2021).\u003c/li\u003e\n\u003cli\u003eGardner, T. A.\u003cem\u003e et al.\u003c/em\u003e Transparency and sustainability in global commodity supply chains. \u003cem\u003eWorld Development\u003c/em\u003e \u003cstrong\u003e121\u003c/strong\u003e, 163-177, doi:https://doi.org/10.1016/j.worlddev.2018.05.025 (2019).\u003c/li\u003e\n\u003cli\u003eHashem, K. M., He, F. J., Alderton, S. A. \u0026amp; MacGregor, G. A. Cross-sectional survey of the amount of sugar and energy in cakes and biscuits on sale in the UK for the evaluation of the sugar-reduction programme. \u003cem\u003eBMJ Open\u003c/em\u003e \u003cstrong\u003e8\u003c/strong\u003e, e019075, doi:10.1136/bmjopen-2017-019075 (2018).\u003c/li\u003e\n\u003cli\u003eTrase.earth. \u003cem\u003eTrase.earth website\u003c/em\u003e, \u0026lt;https://trase.earth/about\u0026gt; (2025).\u003c/li\u003e\n\u003cli\u003eRenier, C.\u003cem\u003e et al.\u003c/em\u003e Transparency, traceability and deforestation in the Ivorian cocoa supply chain. \u003cem\u003eEnvironmental Research Letters\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, 024030, doi:10.1088/1748-9326/acad8e (2023).\u003c/li\u003e\n\u003cli\u003eEllis, C. J., Yahr, R. \u0026amp; Coppins, B. J. Quantifying the anthropocene loss of bioindicators for an early industrial region: an equitable baseline for biodiversity restoration. \u003cem\u003eBiodiversity and Conservation\u003c/em\u003e \u003cstrong\u003e27\u003c/strong\u003e, 2363-2377, doi:10.1007/s10531-018-1541-y (2018).\u003c/li\u003e\n\u003cli\u003eNewbold, T.\u003cem\u003e et al.\u003c/em\u003e Global effects of land use on local terrestrial biodiversity. \u003cem\u003eNature\u003c/em\u003e \u003cstrong\u003e520\u003c/strong\u003e, 45-50, doi:10.1038/nature14324 (2015).\u003c/li\u003e\n\u003cli\u003eMillard, J.\u003cem\u003e et al.\u003c/em\u003e Global effects of land-use intensity on local pollinator biodiversity. \u003cem\u003eNature Communications\u003c/em\u003e \u003cstrong\u003e12\u003c/strong\u003e, 2902, doi:10.1038/s41467-021-23228-3 (2021).\u003c/li\u003e\n\u003cli\u003eJezeer, R. E., Verweij, P. A., Santos, M. J. \u0026amp; Boot, R. G. A. Shaded Coffee and Cocoa \u0026ndash; Double Dividend for Biodiversity and Small-scale Farmers. \u003cem\u003eEcological Economics\u003c/em\u003e \u003cstrong\u003e140\u003c/strong\u003e, 136-145, doi:https://doi.org/10.1016/j.ecolecon.2017.04.019 (2017).\u003c/li\u003e\n\u003cli\u003eBlaser-Hart, W. J.\u003cem\u003e et al.\u003c/em\u003e The effectiveness of cocoa agroforests depends on shade-tree canopy height. \u003cem\u003eAgriculture, Ecosystems \u0026amp; Environment\u003c/em\u003e \u003cstrong\u003e322\u003c/strong\u003e, 107676, doi:https://doi.org/10.1016/j.agee.2021.107676 (2021).\u003c/li\u003e\n\u003cli\u003eManey, C., Sassen, M. \u0026amp; Hill, S. L. L. Modelling biodiversity responses to land use in areas of cocoa cultivation. \u003cem\u003eAgriculture, Ecosystems \u0026amp; Environment\u003c/em\u003e \u003cstrong\u003e324\u003c/strong\u003e, 107712, doi:https://doi.org/10.1016/j.agee.2021.107712 (2022).\u003c/li\u003e\n\u003cli\u003eArthur, A., Ja, D. \u0026amp; Fk, P. Impact of fertiliser application on cocoa yields in Ghana: An analysis of cocoa bean yields in farmers\u0026rsquo; plantations. \u003cem\u003eAustralian Journal of Crop Science\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, 139-144, doi:10.21475/ajcs.24.18.02.PNE3958 (2024).\u003c/li\u003e\n\u003cli\u003eNaranjo-Merino, C. A., Ort\u0026iacute;z-Rodriguez, O. O. \u0026amp; Villamizar-G, R. A. Assessing Green and Blue Water Footprints in the Supply Chain of Cocoa Production: A Case Study in the Northeast of Colombia. \u003cem\u003eSustainability\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 38 (2018).\u003c/li\u003e\n\u003cli\u003eTedstone, A., Duval, D. \u0026amp; Peacock, E. Dietary health and CVD: implications for dietary policy in England. \u003cem\u003eProceedings of the Nutrition Society\u003c/em\u003e \u003cstrong\u003e79\u003c/strong\u003e, 95-102, doi:10.1017/S0029665119000600 (2020).\u003c/li\u003e\n\u003cli\u003eB\u0026eacute;n\u0026eacute;, C. \u0026amp; Abdulai, A.-R. Navigating the politics and processes of food systems transformation: guidance from a holistic framework. \u003cem\u003eFrontiers in Sustainable Food Systems\u003c/em\u003e \u003cstrong\u003eVolume 8 - 2024\u003c/strong\u003e, doi:10.3389/fsufs.2024.1399024 (2024).\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"sustainable, food system, biodiversity, environmental impact assessment, consumption, trade","lastPublishedDoi":"10.21203/rs.3.rs-7380465/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7380465/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eWe must urgently address unsustainable food system practices to ensure planetary health. Most research quantifying the environmental impacts of food production and consumption has not measured biodiversity impacts directly or spatially, limiting our ability to evaluate how producing and consuming countries are contributing to food system sustainability. Here, we quantify the biodiversity impacts of food crops commonly found in processed foods using a range of accessible metrics. We focus on ingredients found in popular processed foods, such as the chocolate digestive biscuit, and show globally reaching impacts. Land-use and fertiliser impacts were greatest for the food crops considered. We also identified Critically Endangered species at risk from habitat loss associated with crops being supplied to the UK specifically. For instance, cocoa production overlaps with nearly 20% of the range of both the critically endangered Roloway and Miss Waldron\u0026rsquo;s Red Colobus monkeys. As crop-specific farm and trade data are not publicly available, it is difficult to know the origin of all ingredients, limiting biodiversity impact estimation. If such data are released in future, our approach can be used to further evaluate biodiversity impacts of foods to inform more sustainable decision-making at consumer, business, and government levels.\u003c/p\u003e","manuscriptTitle":"Key crops for processed foods have spatially variable biodiversity impacts not captured in other environmental impact indicators","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-09-08 12:26:29","doi":"10.21203/rs.3.rs-7380465/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2025-11-17T10:12:49+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-11-11T16:58:29+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"60405344949612629446954638829358368626","date":"2025-11-02T11:31:21+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2025-10-15T08:24:56+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"283424410495443861306883312943015172449","date":"2025-10-01T11:11:40+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"186587943677245470719652007590959670144","date":"2025-09-08T13:03:35+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-09-02T09:18:12+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-08-28T08:59:23+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2025-08-21T12:53:59+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2025-08-20T13:40:19+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2025-08-20T13:26:39+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"d1134570-de91-48df-ad52-e11948d61270","owner":[],"postedDate":"September 8th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":54324776,"name":"Biological sciences/Ecology"},{"id":54324777,"name":"Earth and environmental sciences/Ecology"},{"id":54324778,"name":"Earth and environmental sciences/Environmental sciences"}],"tags":[],"updatedAt":"2026-02-02T16:07:48+00:00","versionOfRecord":{"articleIdentity":"rs-7380465","link":"https://doi.org/10.1038/s41598-025-34850-2","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2026-01-29 15:59:23","publishedOnDateReadable":"January 29th, 2026"},"versionCreatedAt":"2025-09-08 12:26:29","video":"","vorDoi":"10.1038/s41598-025-34850-2","vorDoiUrl":"https://doi.org/10.1038/s41598-025-34850-2","workflowStages":[]},"version":"v1","identity":"rs-7380465","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7380465","identity":"rs-7380465","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}