--------------------------------------------- # Data from: Expanded potential growing region and yield increase for Agave americana with future climate Preferred citation (DataCite format): LeBauer, David Shaner; Davis, Sarah; Abatzoglou, John (2022). Data from: Expanded potential growing region and yield increase for Agave americana with future climate. University of Arizona Research Data Repository. Dataset. https://doi.org/10.25422/azu.data.16828279 Corresponding Author: Sarah Davis, Voinovich School of Leadership and Public Service, Ohio University, daviss6@ohio.edu License: CC0 (data) BSD 3-Clause (code) DOI: https://doi.org/10.25422/azu.data.16828279 --------------------------------------------- ## Summary In Davis et al (2021) we simulated the potential productivity of _Agave americana_ under current and future climates, with and without water management using either irrigation or rock mulch. This repository contains data generated for the analyses in Davis et al. (2021). The code and code documentation can be more easily browsed on GitHub at https://github.com/cct-datascience/agave-prediction, but are also archived here. The key dataset of interest is allbiomass.nc, which provides predicted annual biomass increments for Agave americana under present and future climates and under rainfed, irrigated, and rock mulch + irrigation scenarios. It contains the model predictions of annual biomass production with units of Mg/ha/y for six scenarios on the same latitude and longitude dimensions as the TERRACLIMATE dataset: - biomass: rainfed, historical climate - biomass_irrig: irrigated, historical climate - biomass_rockmulch: rainfed with rock mulch, historical climate - biomass4C: rainfed, +4C climate - biomass_irrig: irrigated, +4C climate - biomass_rockmulch: rainfed with rock mulch, +4C climate Other files likely of interest generated as intermediate steps in the analysis: absmin19812010.nc and absmin4C.nc These files contain the variable absmin, the absolute minimum temperature. These files also contain tmin, adjustC, and cadj used to compute absmin as described in the text and Calculating Absolute Minimum Annual Temperature section of the README. - tmin: Average minimum 2-m air temperature, from TERRACLIMATE - absmin: Absolute minimum 2-m temperature, calculated for this study to determine viable range of Agave americana coefs19812010.nc and coefs4C.nc These files provide calculated values of environmental productivity index model parameters alpha, beta, and gamma. These coefficients have values in the range [0,1], with 0 indicating that the environment does not allow growth and 1 indicating that the environmental factor does not limit growth. These coefficients are calculated monthly for each grid cell under 1981-2010 and +4C climate normals, and beta is calculated for the three water availability scenarios. These coefficients have values in the range [0,1], with 0 indicating that the environment does not allow growth and 1 indicating that the environmental factor does not limit growth. - alpha: light limitation coefficient - beta: water limitation coefficient - beta2: water limitation coefficient to simulate irrigation - beta3: water limitation coefficient to simulate rainfed with rock mulch - gamma: minimum temperature coefficient Files provided for reference, not intended for reuse: - README.pdf contains the README that explains the analysis steps, data sources, and outputs in more detail. - The file not_suitable.gpkg contains polygons representing the union of land that is either protected (UNEP-WCMC and IUCN 2021) or urbanized (Kelso and Patterson 2010). As UNEP-WCMC is updated monthly, researchers are advised to download the current version. - figures.zip contains all map figures generated using Panoply, as well as Panoply settings files with the extension .pcl. - agave-prediction-master_20220511.zip contains a snapshot of code in the GitHub repository. - The original files containing adjustC and cadj are archived in the file absmin_adjustments.zip because these are the inputs to the analysis. The same data layers are also provided alongside calculated absmin in the files absmin19802010.nc and absmin4C.nc. All NetCDF files (those ending in .nc) have additional metadata. The paper and software repository provide additional details. To cite the research: Davis, S.C.; Abatzoglou, J.T.; LeBauer, D.S. Expanded Potential Growing Region and Yield Increase for Agave americana with Future Climate. Agronomy 2021, 11, 2109. https://doi.org/10.3390/agronomy11112109 To cite these datasets (this archive): LeBauer, David Shaner; Davis, Sarah; Abatzoglou, John (2022): Data from: Expanded potential growing region and yield increase for Agave americana with future climate. University of Arizona Research Data Repository. Dataset. https://doi.org/10.25422/azu.data.16828279 To cite the code: LeBauer, D., Davis, S., & Abatzoglou, J. (2021). Data and Code from: Expanded potential growing region and yield increase for Agave americana with future climate (Version davis_etal_2021) [Computer software]. https://github.com/cct-datascience/agave-prediction --------------------------------------------- ## Files and Folders See Summary for more detailed descriptions of some files. - `allbiomass.nc`: Predicted annual biomass increments for _Agave americana_ under all climates and scenarios - `absmin19812010.nc`: Contains the variable _absmin_, the absolute minimum temperature, under 1981-2010 climate normals - `absmin4C.nc`: Contains the variable _absmin_, the absolute minimum temperature, under +4C climate normals - `coefs19812010.nc`: Calculated values of α, β, and Γ values for each grid cell under 1981-2010 climate normals for three water availability scenarios - `coefs4C.nc`: Calculated values of α, β, and Γ values for each grid cell under +4C climate normals - `README.pdf`: README that explains the analysis steps, data sources, and outputs in more detail - `not_suitable.gpkg`: Polygons representing the union of land that is either protected or urbanized - `figures.zip`: All map figures generated using Panoply and .pcl Panoply settings files - `agave-prediction-master.zip`: Snapshot of code in the GitHub repository - `01-download_climate.R` downloads the TERRA-Climate data to the `raw_data` directory. - `01.1-combine_masks.R` combines urban and protected area polygons and generates polygons stored in `/data/derived_data/not_suitable.gpkg` - `0.2-append_terraclimate.sh` combines TERRA-CLIMATE data into a single file for each scenario in the derived_data folder. `derived_data` named `TerraClimate4C.nc` and `TerraClimate19812010.nc` - `03-compute.biomass.sh` implements the model under different climate scenarios with rainfed, irrigation, and rock mulch scenarios (described below) - `04-update_nc_metadata.sh` corrects metadata in output netcdf files - `05-analysis.R` applies the urban and protected area masks, generates masked versions of the maps as netCDF files, and calculates the potential area and total biomass for each scenario. - `06-figures.Rmd` generates the tables using Panoply and also provides alternatives to generating maps in R, but these were not used. - `absmin_adjustments.zip`: Original files containing _adjustC_ and _cadj_ --------------------------------------------- ## Materials & Methods ## Software Used **NCO: NetCDF Operators (Zender 2014)** Zender, C. S. (2014), netCDF Operator (NCO) User Guide, Version 4.4.3, http://nco.sf.net/nco.pdf. **R** R Core Team (2018). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/. **fasterize** Noam Ross (2020). fasterize: Fast Polygon to Raster Conversion. R package version 1.0.3. https://CRAN.R-project.org/package=fasterize **rnaturalearth** Andy South (2017). rnaturalearth: World Map Data from Natural Earth. R package version 0.1.0. https://CRAN.R-project.org/package=rnaturalearth **sf** Pebesma, E., 2018. Simple Features for R: Standardized Support for Spatial Vector Data. The R Journal 10 (1), 439-446, https://doi.org/10.32614/RJ-2018-009 **exactextractr** Daniel Baston (2021). exactextractr: Fast Extraction from Raster Datasets using Polygons. R package version 0.6.0. https://CRAN.R-project.org/package=exactextractr **Panoply** PanoplyCL v 4.12.8 b0071 downloaded from https://www.giss.nasa.gov/tools/panoply/download/ can be used to recreate figures, see `.pcl` files in figures directory Full details on software and data used in the analysis are provided here: https://github.com/az-digitalag/agave-prediction/blob/davis_etal_2021/README.md --------------------------------------------- ## Contributor Roles The roles are defined by the CRediT taxonomy http://credit.niso.org/ - Sarah Davis, Ohio University: Conceptualization, Methodology, Formal Analysis, Resources, Data Curation, Writing - review & editing, Visualization, Project Administration, Funding Acquisition - John Abatzoglou, University of California at Merced: Methodology, Validation, Writing – review & editing - David LeBauer, University of Arizona: Methodology, Software, Validation, Formal Analysis, Resources, Data Curation, Writing – review & editing, Visualization --------------------------------------------- ## Additional Notes Links: - https://github.com/az-digitalag/agave-prediction