• Energy Research
• 2021-2025close
• CA close
• ZENODO close

These datasets are associated with the following paper: Metsaranta, J.M., Mamet, S.D., Maillett, J., Barr, A.G. (2021). Comparison of tree-ring and eddy covariance derived annual ecosystem production estimates for jack pine and trembling aspen forests in Saskatchewan, Canada. Agricultural and Forest Meteorology. There are two files: (1) CBMOutput.zip. This contains the hybrid biometric modelled ecosystem C stock and flux estimates. (2) StandReconstructionData.zip. This contains the field measurement data and the tree level biomass and wood volume data for the Stand Reconstruction plots used to develop the hybrid biometric modelled estimates. The data are formatted as .csv files, and an associated Microsoft Excel spreadsheet explains the data columns and provides information on the associated units of measure.

Most research on boreal populations of Woodland caribou (Rangifer tarandus caribou) has been conducted in areas of high anthropogenic disturbance. However, a large portion of the species’ range overlaps relatively pristine areas primarily disturbed by natural disturbances, such as wildfire. Climate-driven habitat change is a key concern for the conservation of boreal-dependent species, where management decisions have yet to consider knowledge from multiple ecological domains integrated into a cohesive and spatially explicit forecast of species-specific habitat and demography. We used a novel ecological forecasting framework to provide climate-sensitive projections of habitat and demography for five boreal caribou monitoring areas within the Northwest Territories (NWT), Canada, over 90 years. Importantly, we quantify uncertainty around forecasted mean values. Our results suggest habitat suitability may increase in central and southwest regions of the NWT’s Taiga Plains ecozone but decrease in southern and northwestern regions driven by conversion of coniferous to deciduous forests. We do not project boreal caribou population growth rates to change despite forecasted changes to habitat suitability. Our results emphasize the importance of efforts to protect and restore northern boreal caribou habitat despite climate uncertainty while highlighting expected spatial variations that are important considerations for local people who rely on them. An ability to reproduce previous work, and critical thought when incorporating sources of uncertainty, will be important to refine forecasts, derive management decisions, and improve conservation efficacy for northern species at risk. Please see the README document ("README.md") and the accompanying published article: Stewart, Micheletti et al. 2023. Climatepinformed forecasts reveal dramatic local habitat shifts and population uncertainty for nothern boreal caribou. Ecological Applications.

A lack of adequate high resolution climate proxy records for the Last Glacial Maximum (LGM) has prevented the extrapolation of climate–solar linkages on centennial time scales prior of the Holocene. Therefore, it is still unknown whether centennial climate variations of the last ten thousand years convey a universal climate change or merely represent a characteristic of the Holocene. Recently published high resolution climate proxy records for the LGM allowed us to extrapolate climate–solar linkages on centennial time scales ahead of the Holocene. Here we present the analysis of a high resolution pollen concentration record from Lake Kotokel in southern Siberia, Russia, during the LGM. The record reflects the dynamics of vegetation zones and temperature change with a resolution of ~ 40 years in the continental climate of north-eastern Asia. We demonstrate that our pollen concentration record, the oxygen isotope δ18O record from the Greenland ice core project NGRIP (NorthGRIP), the dust-fall contributions in Lake Qinghai, China, grain size in the Gulang and Jingyuan loess deposits, China, and the composite oxygen isotope δ18O record from the Alpine cave system 7H reveal cooler to warmer climate fluctuations between ~ 20.6 and 26 ka. Such fluctuations correspond to the ~ 1000-yr, 500-600-yr and 210-250-yr cycles possibly linked to the solar activity variations and recognized in high resolution Holocene proxies all over the world. We further show that climate fluctuations in the LGM and Holocene are spectrally similar suggesting that linkages between climate proxies and solar activity at the centennial time scale in the Holocene can be extended to the LGM. {"references": ["Vadim A. Kravchinsky, Rui Zhang, Ryan Borowiecki, Pavel E. Tarasov, Mirko van der Baan, Taslima Anwar, Avto Goguitchaichvili, Stefanie M\u00fcller, 2021. Centennial scale climate oscillations from southern Siberia in the Last Glacial Maximum. Quaternary Science Reviews, in press."]}

Climate-informed conservation priorities in British Columbia (Version 1.0) Territorial acknowledgement: We respectfully acknowledge that we live and work across diverse unceded territories and treaty lands and pay our respects to the First Nations, Inuit and Métis ancestors of these places. We honour our connections to these lands and waters and reaffirm our relationships with one another. Suggested citation: Stolar, J., D. Stralberg, I. Naujokaitis-Lewis, S.E. Nielsen, and G. Kehm. 2023. Spatial priorities for climate-change refugia and connectivity for British Columbia (Version 1.0). Place of publication: University of Alberta, Edmonton, Canada. doi: 10.5281/zenodo.8333303 Corresponding author: stolar@ualberta.ca Summary: The purpose of this project is to identify spatial locations of (a) vulnerabilities within British Columbia’s current network of protected areas and (b) priorities for conservation and management of natural landscapes within British Columbia under a range of future climate-change scenarios. This involved adaptation and implementation of existing continental- and provincial-scale frameworks for identifying areas that have potential to serve as refugia from climate change or corridors for species migration. Outcomes of this work include the provision of practical guidance for protected areas network design and vulnerabilities identification under climate change, with application to other regions and jurisdictions. Project results, in the form of multiple spatial prioritization scenarios, may be used to evaluate the resilience of the existing protected area network and other conservation designations to better understand the risks to British Columbia’s biodiversity in our changing climate. Description: These raster layers represent different scenarios of Zonation rankings of conservation priorities for climate resilience and connectivity between current and 2080s conditions for a provincial-scale analysis. Input conservation features included metrics of macrorefugia (forward and backward climate velocity (km/year), overlapping future and current habitat suitability for ~900 rare species in BC), microrefugia (presence of old growth ecosystems, drought refugia, glaciers/cool slopes/wetlands, and geodiversity), and connectivity. Please see details in the accompanying report. File nomenclature: .zip folder (Stolar_et_al_2023_CiCP_Zenodo_upload_Version_1.0.zip): Contains the files listed below. Macrorefugia (2080s_macrorefugia.tif): Scenarios for each taxonomic group (equal weightings for all species) (Core-area Zonation Function) Climate-type velocity + species scenarios from above (Core-area Zonation; equal weightings) Microrefugia (microrefugia.tif): Scenario with old growth forest habitat, landscape geodiversity, wetlands/cool slopes/glaciers, drought refugia (Core-area Zonation; equal weightings) Overall scenario (2080s_macro_micro_connectivity.tif): Inputs from above (with equal weightings) + connectivity metrics (each weighted at 0.1) (Additive Benefit Function Zonation) Conservation priorities (Conservation_priorities_2080s.tif): Overall scenario from above extracted to regions of low human footprint. Restoration priorities (Restoration_priorities_2080s.tif): Overall scenario from above extracted to regions of high human footprint. Accompanying report (Stolar_et_al_2023_CiCP_Zenodo_upload_Version_1.0.pdf): Documentation of rationale, methods and interpretation. READ_ME file (READ_ME_PLEASE.txt): Metadata. Legend interpretation: Ranked Zonation priorities increase from 0 (lowest) to 1 (highest). Raster information: Columns and Rows: 1597, 1368 Number of Bands: 1 Cell Size (X, Y): 1000, 1000 Format: TIFF Pixel Type: floating point Compression: LZW Spatial reference: XY Coordinate System: NAD_1983_Albers Linear Unit: Meter (1.000000) Angular Unit: Degree (0.0174532925199433) false_easting: 1000000 false_northing: 0 central_meridian: -126 standard_parallel_1: 50 standard_parallel_2: 58.5 latitude_of_origin: 45 Datum: D_North_American_1983 Extent: West -139.061502 East -110.430823 North 60.605550 South 47.680823 Disclaimer: The University of Alberta (UofA) is furnishing this deliverable "as is". UofA does not provide any warranty of the contents of the deliverable whatsoever, whether express, implied, or statutory, including, but not limited to, any warranty of merchantability or fitness for a particular purpose or any warranty that the contents of the deliverable will be error-free. Funding: We gratefully acknowledge the financial support of Environment and Climate Change Canada, the Province of British Columbia through the Ministry of Water, Land and Resource Stewardship) and the Ministry of Environment and Climate Change Strategy, the BC Parks Living Lab for Climate Change and Conservation, and the Wilburforce Foundation. We gratefully acknowledge the financial support of Environment and Climate Change Canada, the Province of British Columbia through the Ministry of Water, Land and Resource Stewardship) and the Ministry of Environment and Climate Change Strategy, the BC Parks Living Lab for Climate Change and Conservation, and the Wilburforce Foundation.

Aim Despite an increasing number of studies highlighting the impacts of climate change on boreal species, the main factors that will drive changes in species assemblages remain ambiguous. We study how species community composition would change following anthropogenic and natural disturbances. We determine the main drivers of assemblage dissimilarity for bird and beetle communities. Location Côte-Nord, Québec, Canada. Methods We quantify two climate-induced pathways based on direct and indirect effects on species occurrence under different harvest management scenarios. The direct climate effects illustrate the impact of climate variables while the indirect effects are reflected through habitat-based climate change. We develop empirical models to predict the distribution of more than 100 species over the next century. We analyze the regional and the latitudinal species assemblage dissimilarity by decomposing it into 'balanced variation in species occupancy and occurrence' and 'occupancy and occurrence gradient'. Results Both pathways increased dissimilarity in species assemblage. At the regional scale, both effects have an impact on decreasing the number of winning species. Yet, responses are much larger in magnitude under mixed climate effects (a mixture of direct and indirect effects). Regional assemblage dissimilarity reached 0.77 and 0.69 under mixed effects versus 0.09 and 0.10 under indirect effects for beetles and birds, respectively, between RCP8.5 and baseline climate scenarios when considering harvesting. Latitudinally, assemblage dissimilarity increased following the climate conditions pattern. Main conclusions The two pathways are complementary and alter biodiversity, mainly caused by species turnover. Yet, responses are much larger in magnitude under mixed climate effects. Therefore, the inclusion of climatic variables considers aspects other than just those related to forest landscapes, such as life cycles of animal species. Moreover, we expect differences in occupancy between the two studied taxa. This could indicate the potential range of change in boreal species concerning novel environmental conditions.

# Dataset: Subglacial freshwater driven speedup of East Antarctic outlet glacier retreat [https://doi.org/10.5061/dryad.1vhhmgr0b](https://doi.org/10.5061/dryad.1vhhmgr0b) Journal: Journal of Geophysical Research: Earth Surface Principle Investigator: * Tyler Pelle, Scripps Institution of Oceanography, University of California San Diego, [tpelle@ucsd.edu](mailto:tpelle@ucsd.edu) Co-Authors: * Dr. Jamin Greenbaum, Scripps Institution of Oceanography, University of California San Diego * Dr. Shivani Ehrenfeucht, Department of Geography and Environmental Management, University of Waterloo * Prof. Christine Dow, Department of Geography and Environmental Management, University of Waterloo * Dr. Felicity S. McCormack, Securing Antarctica's Environmental Future, School of Earth, Atmosphere, & Environment, Monash University Created on October 4, 2023 ## Description of the data and file structure ### File description: 1. runme.m - MATLAB script used to run coupled ISSM-GlaDS SSP5-8.5_{F,M} simulation - includes melt rate parameterization. 2. ssp585.mat – Yearly ice sheet model output from 2017-2100 for SSP5-8.5 simulation. 3. ssp585_F.mat – Yearly ice sheet model output from 2017-2100 for SSP5-8.5_{F} simulation. 4. ssp585_M.mat – Yearly ice sheet model output from 2017-2100 for SSP5-8.5_{M} simulation. 5. ssp585_FM.mat – Yearly ice sheet model output from 2017-2100 for SSP5-8.5_{F,M} simulation. 6. ssp126.mat – Yearly ice sheet model output from 2017-2100 for SSP1-2.6 simulation. 7. ssp126_F.mat – Yearly ice sheet model output from 2017-2100 for SSP1-2.6_{F} simulation. 8. ssp126_M.mat – Yearly ice sheet model output from 2017-2100 for SSP1-2.6_{M} simulation. 9. ssp126_FM.mat – Yearly ice sheet model output from 2017-2100 for SSP1-2.6_{F,M} simulation. 10. ssp585_Totten_T.mat - Bi-weekly ocean temperature (Ta) for Totten Glacier from January 1, 2017 to December 31, 2099 (high emission). 11. ssp585_Moscow_T.mat - Bi-weekly ocean temperature (Ta) for Moscow University Glacier from January 1, 2017 to December 31, 2099 (high emission). 12. ssp585_Vander_T.mat - Bi-weekly ocean temperature (Ta) for Vander Glacier from January 1, 2017 to December 31, 2099 (high emission). 13. ssp585_Totten_S.mat - Bi-weekly ocean salinity (Sa) for Totten Glacier from January 1, 2017 to December 31, 2099 (high emission). 14. ssp585_Moscow_S.mat - Bi-weekly ocean salinity (Sa) for Moscow University Glacier from January 1, 2017 to December 31, 2099 (high emission). 15. ssp585_Vander_S.mat - Bi-weekly ocean salinity (Sa) for Vander Glacier from January 1, 2017 to December 31, 2099 (high emission). 16. ssp126_Totten_T.mat - Bi-weekly ocean temperature (Ta) for Totten Glacier from January 1, 2017 to December 31, 2099 (low emission). 17. ssp126_Moscow_T.mat - Bi-weekly ocean temperature (Ta) for Moscow University Glacier from January 1, 2017 to December 31, 2099 (low emission). 18. ssp126_Vander_T.mat - Bi-weekly ocean temperature (Ta) for Vander Glacier from January 1, 2017 to December 31, 2099 (low emission). 19. ssp126_Totten_S.mat - Bi-weekly ocean salinity (Sa) for Totten Glacier from January 1, 2017 to December 31, 2099 (low emission). 20. ssp126_Moscow_S.mat - Bi-weekly ocean salinity (Sa) for Moscow University Glacier from January 1, 2017 to December 31, 2099 (low emission). 21. ssp126_Vander_S.mat - Bi-weekly ocean salinity (Sa) for Vander Glacier from January 1, 2017 to December 31, 2099 (low emission). 22. TotBasin.exp - Polygon that contains Totten Glacier over which Totten's ocean temperature is applied. 23. MuisBasin.exp - Polygon that contains Moscow University Glacier over which Totten's ocean temperature is applied. 24. VandBasin.exp - Polygon that contains Vanderford Glacier over which Totten's ocean temperature is applied. ### File specific information: **ASB_IceHydroModel.mat**: All data associated with the ice sheet and subglacial hydrology model initial state is held in ASB_IceHydroModel.mat, which contains a MATLAB ‘model’ object (for more information, see [https://issm.jpl.nasa.gov/documentation/modelclass/](https://issm.jpl.nasa.gov/documentation/modelclass/). In MATLAB, the model can be loaded and displayed by running load(‘ASB_IceHydroModel.mat’), which will load in the model variable ‘md’. Of particular interest will be the following data contained in md: md.mesh (mesh information), md.geometry (initial ice sheet geometry, ice shelf geometry, and bed topography), md.hydrology (initial hydrology model fields), md.initialization (model initialization fields) and md.mask (ice mask and grounded ice mask). Note that all fields are defined on the mesh nodes, and one can plot a given field in MATLAB using the ISSM tool ‘plotmodel’ (e.g., plotmodel(md,'data',md.geometry.bed) will plot the model bed topography). For more information on plotting, please see [https://issm.jpl.nasa.gov/documentation/plotmatlab/](https://issm.jpl.nasa.gov/documentation/plotmatlab/). **Model output files (e.g. ssp585_FM.mat)**: Yearly ice sheet model results between 2017-2100 for all model simulations described in the paper. Fields appended with '*' are included in results with changing subglacial hydrology (ssp126_F, ssp126_M, ssp126_FM, ssp585_F, ssp585_M, ssp585_FM). Fields appended with '**' are included in results where ice shelf melt is enhanced by subglacial discharge (ssp126_M, ssp126_FM, ssp585_M, ssp585_FM). These files contain a MATLAB variable that is the same as the file name, which is a model object of size 1x83 that contains the following yearly variables: * \* Vel (velocity norm, m/yr) * \* Thickness (ice sheet thickness, m) * \* Surface (ice sheet surface elevation, m) * \* Base (ice sheet base elevation, m) * \* BasalforcingsFloatingiceMeltingRate (ice shelf basal melting rate field, m/yr) * \* MaskOceanLevelset (ground ice mask, grounded ice if > 0, grounding line position if = 0, floating ice if < 0) * \* IceVolume (total ice volume in the model domain, t) * \* IceVolumeAboveFloatation (total ice volume in the model domain that is above hydrostatic equilibrium, t) * \* TotalFloatingBmb (Total floating basal mass balance, Gt) * \* \\*ChannelDischarge\\_Node (GlaDS-computed channel discharge interpolated onto model node, m3/s) * \* \\*ChannelDiameter\\_Node (GlaDS-computed channel diameter interpolated onto model node, m) * \* \\*ChannelArea (GlaDS-computed channel area defined on model edges, m2) * \* \\*ChannelDischarge (GlaDS\\_computed channel discharge defined on model edges, m3/s) * \* \\*EffectivePressure (GlaDS-computed ice sheet effective pressure, Pa) * \* \\*HydraulicPotential (GlaDS computed hydraulic potential, - * \* \\*HydrologySheetThickness (GlaDS-computed after sheet thickness, m) * \* \\*GroundedIceMeltingRate (Grounded ice melting rate defined on all grounded nodes, m/yr) * \* \\*\\*melt\\_nodis (ice shelf basal melting rate computed when discharge is set to zero, m/yr) * \* \\*\\*zgl (grounding line height field, m) * \* \\*\\*glfw (grounding line fresh water flux field, m2/s) * \* \\*\\*chan\\_wid (Domain average subglacial discharge channel width, m) * \* \\*\\*maxdist (5L' length scale used in melt computation, m) * \* \\*\\*maxis (maximum discharge at each subglacial outflow location, m2/s) * \**\\*\\_T.mat**: Bi-weekly ocean temperature extracted from an East Antarctic configuration of the MITgcm (Pelle et al., 2021), where '\\*' ssp126 (low emission) or ssp585 (high emission). Ocean temperature was averaged adjacent to each target ice front in both depth and in the contours shown in figure 1b. * \**\\*\\_S.mat**: Same as above, but for salinity in units on the Practical Salinity Scale (PSU). * \***.exp**: Exp files that contain coordinates that outline a polygon for the drainage basins of each major glacier in this study (Vanderford Glacier contains the drainage basins for Adams, Bond, and Underwood Glaciers as well). Recent studies have revealed the presence of a complex freshwater system underlying the Aurora Subglacial Basin (ASB), a region of East Antarctica that contains ~7 m of global sea level potential in ice mainly grounded below sea level. Yet, the impact that subglacial freshwater has on driving the evolution of the dynamic outlet glaciers that drain this basin has yet to be tested in a coupled ice sheet-subglacial hydrology numerical modeling framework. Here, we project the evolution of the primary outlet glaciers draining the ASB (Moscow University Ice Shelf, Totten, Vanderford, and Adams Glaciers) in response to an evolving subglacial hydrology system and to ocean forcing through 2100, following low and high CMIP6 emission scenarios. By 2100, ice-hydrology feedbacks enhance the ASB’s 2100 sea level contribution by ~30% (7.50 mm to 9.80 mm) in high emission scenarios and accelerate retreat of Totten Glacier’s main ice stream by 25 years. Ice-hydrology feedbacks are particularly influential in the retreat of the Vanderford and Adams Glaciers, driving an additional 10 km of retreat in fully-coupled simulations relative to uncoupled simulations. Hydrology-driven ice shelf melt enhancements are the primary cause of domain-wide mass loss in low emission scenarios, but are secondary to ice sheet frictional feedbacks under high emission scenarios. The results presented here demonstrate that ice-subglacial hydrology interactions can significantly accelerate retreat of dynamic Antarctic glaciers and that future Antarctic sea level assessments that do not take these interactions into account might be severely underestimating Antarctic Ice Sheet mass loss. In this data publication, we present the model output and results associated with the following manuscript recently submitted to the Journal of Geophysical Research: Earth Surface: “Subglacial discharge accelerates ocean driven retreat of Aurora Subglacial Basin outlet glaciers over the 21st century”. We include yearly ice sheet model output between 2017-2100 for eight numerical ice-subglacial hydrology model runs. We also include the ice sheet and subglacial hydrology model initial states. In addition, we include all ocean forcing time-series (temperature and salinity for the low emission and high emission climate forcing scenarios for three glacial regions), which are used as input into the melt parameterization. Lastly, we include a MATLAB script that contains the code used to couple the ice-subglacial hydrology models as well as a "readme" file with further information on all data in this publication. Ice sheet model results: Direct results taken from the Ice-sheet and Sea-level System Model (ISSM, Larour et al. 2012) with no processing applied, provided yearly as *.mat files. Ice sheet and subglacial hydrology model initial states: Initial state of the ice sheet model (ice geometry, mesh information, inversion results, etc.) and subglacial hydrology model (steady-state water column thickness, effective pressure, channelized discharge, etc.) containing Aurora Subglacial Basin outlet glaciers with no processing applied, provided as a *.mat file. The contents of the *.mat file is a MATLAB variable of class "model", which is compatible with ISSM. Model coupling script: Documented MATLAB script ready to run with the provided data sets. Ocean temperature and salinity timeseries: Bottom ocean temperature (°C) and salinity (PSU) timeseries (January 1st, 2017 through December 31, 2099) extracted from an East Antarctic configuration of the ocean component of the MITgcm (Pelle et al., 2021). Temperature and salinity are provided bi-weekly and averged both in depth and along the ice fronts of Moscow University, Totten, and Vanderford Glaciers (see white dashed contour in figure 1b of the main manuscript text). Data are provided as *.mat files. Polygons that provide locaion to apply ocean temperature and salinity: Polygons provided as a list of x/y coordinates (meters) are provided in three *.exp files that cover the drainage basins of Moscow University, Totten, and Vanderford Glaciers (the polygon for Vanderford also includes the drainage basins of Adams, Bond, and Underwood Glaciers). 

Logging is the main human disturbance impacting biodiversity in forest ecosystems. However, the impact of forest harvesting on biodiversity is modulated by abiotic conditions through complex relationships that remain poorly documented. Therefore, the interplay between forest management and climate change can no longer be ignored. Our aim was to study the expected long-term variations in the assemblage of bird and beetle communities following modifications in forest management under different climate change scenarios. We developed species distribution models to predict the occurrence of 88 species of birds and beetles in eastern Canadian boreal forests over the next century. We simulated three climate scenarios (baseline, RCP4.5 and RCP8.5) under which we varied the level of harvesting. We also analyzed the regional assemblage dissimilarity by decomposing it into balanced variations in species occupancy and occupancy gradient. We predict that forest harvesting will alter the diversity by increasing assemblage dissimilarity under all the studied climate scenarios, mainly due to species turnover. Species turnover intensity was greater for ground-dwelling beetles, probably because they have lower dispersal capacity than flying beetles or birds. A good dispersal capacity allows species to travel more easily between ecosystems across the landscape when they search for suitable habitats after a disturbance. Regionally, an overall increase in the probability of occupancy is projected for bird species, whereas a decrease is predicted for beetles, a variation that could reflect differences in ecological traits between taxa. Our results further predict a decrease in the number of species that increase their occupancy after harvest under the most severe climatic scenario for both taxa. We anticipate that under severe climate change, increasing forest disturbance will be detrimental to beetles associated with old forests but also with young forests after disturbances.

Premise: Seed recruitment niches along estuarine elevation gradients are seldom experimentally field-tested under tidal regimes of the Pacific Northwest of North America. Addressing this knowledge gap is important to better understand estuary restoration and plant community response to sea level rise. Methods: Germination was tested in marsh organ mesocosms across an elevation gradient (0.5–1.7 m above mean sea level). Seeds were sown on sterile peat moss, and the tops of pipes were secured with horticultural “frost cloth” to ensure no experimental seeds were washed out and no new seeds were introduced. The trials tested artificial and overwinter chilling regimes, as well as the presence and/or absence of a near-neighbor transplant. Results: Carex lyngbyei had significant elevation-driven germination after overwinter and artificial chilling. Schoenoplectus tabernaemontani had near-significant germination across elevation after overwinter chilling, and germination in the absence of competition was significantly greater than with a near-neighbor transplant. Discussion: Carex lyngbyei had the highest germination rate at higher elevations, which suggests restricted seed recruitment potential, and required clonal expansion to extend into lower marsh elevations. Identifying species-specific recruitment niches provides insight for restoration opportunities or invasive species monitoring, as well as for estuary migration under sea level rise. Datasets contain temperature data automatically recorded by iButton data loggers (iButtonLink, LLC, Whitewater, WI, USA. Ibuttonlink.com), or manually observed/recorded percent seed germination in experimental mesocosm (marsh organ). Please see ReadMe, then metadata files for each data file. R code is annotated with notes about analysis process, and often links to resources used to produce analyses.