DOI: 10.5281/zenodo.17586378
title: Data for: A fundamental trade-off among resilience, resistance, efficiency, and redundancy in tidal wetlands

description: We filtered the raw NASA-MODIS (MOD13Q1) Enhanced Vegetation Index (EVI) dataset to only inlcude pixels with high tidal wetland class purity and Quality Assurance (QA) reliability scores. We filtered 782,693 tidal wetland pixels with coverage spanning the entire contiguous United States to only include those with greater than 90% tidal wetland class purity. We then further filtered these pixels to only include those where data dropouts in the EVI or QA layer occured fewer than 10% of the time. In the end, we used 145,871 pixels in our analysis. Tidal wetland gross primary production (GPP) was calculated by pixel for the dates 3/5/2000 to 12/2/2020 at 16-day intervals using the algorithms published in: 

R. A. Feagin, I. Forbrich, T.P. Huff, J.G. Barr, J. Ruiz-plancarte, J.D Fuentes, R.G. Najjar, R. Vargas, A. Vazquez-lule, L. Windham-Myers, K. Kroeger, E.J. Ward, G.W. Moore, M. Leclerc, K.W. Krauss, C.L. Stagg, M. Alber, S.H. Knox, K.V.R. Schafer, T.S., Bianchi, J.A. Hutchings, H.B. Nahrawi, A. Noormets, B. Mitra, A. Jaimes, A.L. Hinson, B. Bergamaschi, J. King, and G. Miao., Tidal wetland gross primary production across the continental United States, 2000–2019. Global Biogeochemical Cycles 34, e2019GB006349 (2020).

The file named "rawEVI.csv" contains the raw EVI database prior to filtering and gap-filling. File named "rawSWR.csv" contains the raw SWR database prior to filtering. File named "rawAirTemp.csv" contains the air temperature database prior to filtering. File named "rawQA.csv" contains the QA layer of the MOD13 satellite product prior to filtering. These files contain the raw data for all 782,693 tidal wetland pixel locations. These raw files can also be accessed at daac.ornl.gov.

R. A. Feagin, I. Forbrich, T.P. Huff, J.G. Barr, J. Ruiz-plancarte, J.D Fuentes, R.G. Najjar, R. Vargas, A. Vazquez-lule, L. Windham-Myers, K. Kroeger, E.J. Ward, G.W. Moore, M. Leclerc, K.W. Krauss, C.L. Stagg, M. Alber, S.H. Knox, K.V.R. Schafer, T.S., Bianchi, J.A. Hutchings, H.B. Nahrawi, A. Noormets, B. Mitra, A. Jaimes, A.L. Hinson, B. Bergamaschi, J. King, and G. Miao. (2020). Gross Primary Production Maps of Tidal Wetlands across Conterminous USA, 2000-2019 (Version 1). ORNL Distributed Active Archive Center. https://doi.org/10.3334/ORNLDAAC/1792. Accessed Nov. 2025.

The file named "SWR_90percentFinal.csv" contains the filtered SWR database used to calculate GPP. File named "temp_90percentFinal_rounded.csv" contains the filtered air temperature database used to calculate GPP. The file named "EVI(gapped_filled)_90percentFinal.csv" contains the final gap-filled EVI time series database used to calculate GPP. File named "QA_90%Final.csv" contains the filtered quality assurance (QA) values. Dates in the EVI database where QA = 3 were determined to be of poor quality, removed from the database, and replaced with "NA". Single NA gaps in the EVI database were gap-filled by taking the mean of the dates flanking the NA gap. File named "myGPP(gap_filled)_FINAL.csv" contains the calculated GPP estimates used throughout the study analysis. Code for implementing the data curation and processing workflow used to derive GPP can be found in the file: "CodingWorkflow_Lerner_et_al_2025_Ecology.R". This file also contains original functions for defining perturbation "events" from the GPP time series dataset, as well as code for producing figures found in the Main Text and Appendix S1.