Code for: Global carbon investment in terrestrial biological nitrogen fixation
Data files
Jul 27, 2026 version files 366.44 KB
Abstract
Code that reads in a variety of data sources and analyzes them for the article with the following abstract. Biological nitrogen (N) fixation (BNF) provides the N needed to produce proteins and other biological building blocks, helping feed humanity and mitigate climate change. Due to its high energetic cost compared to other forms of N acquisition, investment in BNF indicates N limitation. Globally gridded BNF flux data provide an opportunity to determine the energetic investment in BNF across ecosystems, which would help reconcile conflicting indicators of N limitation. Here, we use a new BNF synthesis to quantify the relative importance of BNF in different N-fixing niches, in different biomes, and across the globe by calculating the fraction of net primary productivity (NPP) invested in BNF and the fraction of plant N acquisition provided by BNF. Larger fractions of non-agricultural NPP were invested in BNF in less-productive, higher-latitude biomes. This pattern was driven by biocrusts and mosses. Similarly, non-agricultural symbiotic N-fixing plants invested relatively more of their own NPP in BNF in less-productive, higher-latitude biomes. This symbiotic plant pattern was driven by shrubs and herbs, overriding the opposite pattern in trees. Symbiotic plants also acquired a higher fraction of their N from BNF at higher latitudes and in less productive biomes. Investments in symbiotic BNF were 10x higher in agricultural (2.9% of NPP) vs. natural (0.29%) biomes, providing 26% vs. 3.1% of ecosystem-scale plant N acquisition. These results support the paradigm of strong N limitation at higher latitudes, help understand the rarity of N-fixing trees at higher latitudes, underscore the dominance of human activity, and inform models.
Dataset title: Code for the article "Global carbon investment in terrestrial biological nitrogen fixation," published in Global Biogeochemical Cycles.
Dataset contributors: Duncan N. L. Menge
Version number 1.0.1
Date: July 14, 2026
Dataset DOI: https://doi.org/10.5061/dryad.tqjq2bwbn
Suggested citations:
Dataset citation:
Menge DNL. 2026. Code for:
"Global carbon investment in terrestrial biological nitrogen fixation", Dryad, Dataset, https://doi.org/10.5061/dryad.tqjq2bwbn
Corresponding publication:
Menge DNL, Reis Ely CR, Perakis SS, Kou-Giesbrecht S, Cleveland CC, Reed SC, Taylor BN, Batterman SA, Crews TE, Dynarski KA, Funk JL, Gei M, Griffin KL, Gundale MJ, Herridge DF, Jovan SE, Peoples MB, Piipponen J, Rodríguez-Caballero E, Salmon VG, Soper FM, Staccone AP, Weber B, Wolf AA, Wurzburger N. Global carbon investment in terrestrial biological nitrogen fixation. 2026. Global Biogeochemical Cycles e2026GB009098.
Contact information:
Duncan Menge
Department of Ecology, Evolution, and Environmental Biology, Columbia University
dm2972@columbia.edu
ORCID ID: 0000-0003-4736-9844
Funding sources:
United States Geological Survey Powell Center grant (SSP, DNLM, SCR, CCC);
National Science Foundation DEB-1457650 (DNLM, SSP);
National Science Foundation DEB-1457444 (JLF);
National Science Foundation IOS-2129542 (DNLM)
Dates and locations of data collection:
Dates: This dataset only includes code. The underlying data come from other publications and datasets.
Data and code file overview:
Summary metrics:
File count: 3 files
File breakdown and formats:
3 code files (.R).
Total file size: <0.5 MB
File naming conventions: Code files have a .R suffix.
Code file names:
- Menge_etal_Global_C_investment_in_BNF.R
- Menge_etal_Global_C_investment_in_BNF_SuppTable_Field.R
- Menge_etal_Global_C_investment_in_BNF_SuppTable_Greenhouse.R
NOTE: The data are not included as part of this Dryad submission because they are already published in other sources. The data sources are listed in the paper and in the detailed description below. If any users have trouble finding or processing the exact files from these sources, feel free to contact Duncan Menge.
Data and code file details:
The code files are as follows:
Menge_etal_Global_C_investment_in_BNF.R is the main analysis file. It reads in datafiles, calculates all the metrics
needed, prints values for all of the tables except the supplementary table from field and greenhouse data, and creates all the figures. The data files it reads in are:
- FAO_legume_crops_pct_cover_1degree_adj.tif (from https://www.fao.org/faostat/en/#data)
- FluxCom_GPP_2001_2015_crop_1degree.nc (from https://doi.org/10.5194/bg‐17‐1343‐2020)
- IGBP_2019_mosaic_1degree.tif (from https://doi.org/10.5067/MODIS/MCD12Q1.006)
- IGBP_2019_mosaic_1degree_ha.tif (from https://doi.org/10.5067/MODIS/MCD12Q1.006)
- nfix_1degreeMasked.nc (from https://doi.org/10.1038/s41586‐019‐1128‐0)
- N-fixer_cultivated_pastures_pct_cover_1degree.tif (from https://doi.org/10.1111/jvs.12710)
- NPP_MODIS_crop_1degree.nc (from https://doi.org/10.5067/MODIS/MOD17A3HGF.061)
- Powell_BNF_legume_crops_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_cultivated_pastures_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_herbs_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_biocrusts_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_lichen_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_litter_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_mosses_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_shrubs_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_soil_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_trees_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- Powell_BNF_wood_1_1degree.tif (from https://www.sciencebase.gov/catalog/item/66b53cc6d34eebcf8bb3850e)
- sPLOT_Nfixer_herb_pct_rel_cover_1degree.tif (from https://doi.org/10.1111/jvs.12710)
- sPLOT_Nfixer_shrub_pct_cover_1degree.tif (from https://doi.org/10.1111/jvs.12710)
It also inputs data directly from Reis Ely et al. 2025 Nature and Cleveland et al. 2013 PNAS.
Menge_etal_Global_C_investment_in_BNF_SuppTable_Field.R reads in the data for the field portion of Supplementary Table S3, calculates the metrics needed, and prints the values for Supplementary Table S3. The data files it reads in are:
- NY_FX_SNF.csv (from https://doi.org/10.5061/dryad.fbg79cnxk)
- OR_FX_SNF.csv (from https://doi.org/10.5061/dryad.fbg79cnxk)
- HI_W_FX_SNF.csv (from https://doi.org/10.5061/dryad.fbg79cnxk)
- HI_V_FX_SNF.csv (from https://doi.org/10.5061/dryad.fbg79cnxk)
Menge_etal_Global_C_investment_in_BNF_SuppTable_Greenhouse.R reads in the data for the greenhouse portion of Supplementary Table S3, calculates the metrics needed, and prints the values for Supplementary Table S3. The data file it reads in is:
- GX1_data_master_processed.csv (from https://doi.org/10.5061/dryad.cnp5hqc9s)
