DOI: 10.5061/dryad.qnk98sfgj
title: Data for: Freeze tolerance influenced forest cover and hydrology during the Pennsylvanian

description: Global forest cover affects the Earth system by altering surface mass and
 energy exchange. Physiology determines plant environmental limits and
 influences geographical vegetation distribution. Ancient plant physiology,
 therefore, likely affected vegetation-climate feedbacks. We combine
 climate modeling and ecosystem-process modeling to simulate arboreal
 vegetation in the late Paleozoic ice age. Using GENESIS V3 GCM
 simulations, varying pCO2, pO2, and ice extent for the
 Pennsylvanian, and fossil-derived leaf C:N, maximum stomatal conductance,
 and specific conductivity for several major Carboniferous plant groups, we
 simulated global ecosystem processes at a 2-degree (longitude,
 latitude) resolution with Paleo-BGC. Based on leaf water
 constraints, Pangaea could have supported widespread arboreal plant growth
 and forest cover. However, these models do not account for the impacts of
 freezing on plants. According to our interpretation, freezing would have
 affected plants in 89% of unglaciated land during peak glacial periods,
 and 65% during the warmer interglacials. Comparing forest cover, minimum
 temperatures, and paleo-locations of Pennsylvanian-aged plant fossils from
 the Paleobiology Database supports restriction of global forest extent due
 to freezing. Many genera were limited to 25% of unglaciated land where
 temperatures remained above −4°C. Freeze-intolerance of Pennsylvanian
 arboreal vegetation had the potential to alter surface runoff, silicate
 weathering, CO2­ levels, and climate forcing. As a bounding case, we
 assume total plant mortality at −4°C and estimate that contracting forest
 cover increased net global surface runoff by up to 6.1%. Repeated freezing
 likely influenced freeze- and drought-tolerance evolution in lineages like
 the coniferophytes, which became increasingly dominant in the Permian and
 early Mesozoic.
description: Forest cover in the Pennsylvanian
 was simulated for a global grid consisting of 2° × 2° cells using the
 ecosystem process model <i>Paleo</i>-BGC (1), a modified
 version of the ecosystem model BIOME-BGC (2, 3), driven by daily data
 derived from GENESIS V3 (4, 5) Forest cover in the Pennsylvanian was
 simulated for a global grid consisting of 2° × 2° cells using the
 ecosystem process model <i>Paleo</i>-BGC (1), which is a
 modified version of the ecosystem model BIOME-BGC (2, 3) driven by
 daily-resolution output derived from GENESIS V3 (4, 5)( Global
 Environmental and Ecological Simulation of Interactive Systems; Thompson
 and Pollard 1997, Alder et al. 2011).  The <i>Paleo</i>-BGC
 model (1) includes variable atmospheric
 <i>p</i>O<sub>2</sub>, mesophyll conductance
 important for plants with thickened leaves, and leaf hydraulic
 conductivity to account for simple and complex vascular pathways.
 Measurements of leaf carbon to nitrogen ratio (C:N; <i>kg
 </i><i>kg</i><sup><i><sup>-1</sup></i></sup>) were derived from preserved cuticular material. Maximum stomatal conductance (g<sub>smax</sub>; <i>mol </i><i>m</i><sup><i>-2</i></sup><i>s</i><sup><i>-1</i></sup>) was determined from stomatal density and size of leaf fossil impressions. Fossil measurements were collected for Pennsylvanian representative plant-fossil taxa (described in (6) including lycopsids, sphenophytes, pteridosperms, stem-group marattialean tree ferns (tree ferns), and early-diverging coniferophytes (cordaitaleans; (6, 7). Specific leaf area (SLA; <i>m</i><sub>­­</sub><sup><i>2</i></sup><i> kg </i><i>C</i><sup><i>-1</i></sup>) and related leaf attributes, important for converting carbon allocated to leaves to leaf area, were estimated using measured leaf C:N to SLA relationships from extant relatives (8, 9). Boundary layer conductance (g<sub>b</sub>; <i>μmol </i><i>s</i><sup><i>-1</i></sup><i>Pa</i><sup><i>-1</i></sup><i>m</i><sup><i>-2</i></sup>) was estimated from the mean leaf or leaflet width (10).  For leaf mesophyll conductance (g<sub>m</sub>; <i>μmol </i><i>s</i><sup><i>-1</i></sup><i>Pa</i><sup><i>-1</i></sup><i>m</i><sup><i>-2</i></sup>), a single, mean value of 0.273 mol H<sub>2</sub>O m<sup>-2</sup>s<sup>-1</sup> was used as a basis for measurements of leaf air space, leaf thickness, mean mesophyll cell width, and cell-wall thickness of leaf fossil cross-sections of Pennsylvanian taxa (1).  Finally, leaf hydraulic conductivity values (K<sub>leaf</sub>; <i>mmol </i><i>m</i><sup><i>-2</i></sup><i>s</i><sup><i>-1</i></sup><i>MPa</i><sup><i>-1</i></sup>) were calculated from an empirical relationship (11) using minimum and maximum leaf mesophyll pathlength derived from leaf fossil cross-sections (11).     Data was generated by GENESIS V3 as described in (12–14). Minimum temperatures were averaged over terrestrial locations that were not covered by glacial ice. Landmasses and glacial extents were specified after the description of Ziegler, Hulver, and Rowley (15) for ~290 Ma (earliest Artinskian). Vegetated land area was mapped following a two-part categorization process for each un-glaciated terrestrial grid cell.  Inputs taken for <i>Paleo</i>-BGC included daily maximum, minimum, and average temperature (T<sub>min</sub>, T<sub>max</sub>, and T<sub>ave</sub>), precipitation, and shortwave radiation for ten years. Vapor pressure deficit (VPD), a required input for <i>Paleo</i>-BGC, was calculated using T<sub>min </sub>and T<sub>ave</sub> in the modified Tetens Equation (16) for each grid cell. Finally, daily daylength was derived for each grid cell based on latitude and year day (10).  Two climate scenarios were evaluated, intended to provide a characterization of the glacial and interglacial intervals. One climate scenario with glacial ice and CO<sub>2</sub> of 182 ppm is referred to as the glacial, and one with reduced glacial ice and CO<sub>2</sub> of 546 ppm is referred to as the interglacial.  For both scenarios, atmospheric oxygen was set to 28%, which is considered a maximum pO2 for this period based on recent modeling (17).    1.         J. D. White, <i>et al.</i>, A Process-Based Ecosystem Model (Paleo-BGC) to Simulate the Dynamic Response of Late Carboniferous Plants to Elevated O2 And Aridification. <i>American Journal of Science</i> <b>In Press</b> (2020). 2.         M. A. White, P. E. Thornton, S. W. Running, R. R. Nemani, Parameterization and sensitivity analysis of the BIOME-BGC terrestrial ecosystem model: net primary production controls. <i>Earth interactions</i> <b>4</b>, 1–85 (2000). 3.         J. S. Golinkoff, “Estimation and modeling of forest attributes across large spatial scales using BiomeBGC, high-resolution imagery, LiDAR data, and inventory data.,”  University of Montana. (2013). 4.         S. L. Thompson, D. Pollard, Greenland and Antarctic Mass Balances for Present and Doubled Atmospheric CO <sub>2</sub> from the GENESIS Version-2 Global Climate Model. <i>Journal of Climate</i> <b>10</b>, 871–900 (1997). 5.         J. R. Alder, S. W. Hostetler, D. Pollard, A. Schmittner, Evaluation of a present-day climate simulation with a new coupled atmosphere-ocean model GENMOM. <i>Geoscientific Model Development</i> <b>4</b>, 69–83 (2011). 6.         I. P. Montañez, <i>et al.</i>, Climate, pCO2 and terrestrial carbon cycle linkages during late Palaeozoic glacial–interglacial cycles. <i>Nature Geoscience</i> <b>9</b>, 824–828 (2016). 7.         J. D. Richey, <i>et al.</i>, Influence of temporally varying weatherability on CO2-climate coupling and ecosystem change in the late Paleozoic. <i>Clim. Past</i> <b>16</b>, 1759–1775 (2020). 8.         M. T. van Wijk, M. Williams, G. R. Shaver, Tight coupling between leaf area index and foliage N content in arctic plant communities. <i>Oecologia</i> <b>142</b>, 421–427 (2005). 9.         Joseph. D. White, Neal. A. Scott, Specific leaf area and nitrogen distribution in New Zealand forests: Species independently respond to intercepted light. <i>Forest Ecology and Management</i> <b>226</b>, 319–329 (2006). 10.       G. s. Campbell, J. M. Norman, <i>An introduction to environmental biophysics</i> (Springer, 1998). 11.       T. J. Brodribb, T. S. Feild, G. J. Jordan, Leaf Maximum Photosynthetic Rate and Venation Are Linked by Hydraulics. <i>PLANT PHYSIOLOGY</i> <b>144</b>, 1890–1898 (2007). 12.       D. E. Horton, C. J. Poulsen, I. P. Montañez, W. A. DiMichele, Eccentricity-paced late Paleozoic climate change. <i>Palaeogeography, Palaeoclimatology, Palaeoecology</i> <b>331–332</b>, 150–161 (2012). 13.       D. E. Horton, C. J. Poulsen, D. Pollard, Influence of high-latitude vegetation feedbacks on late Palaeozoic glacial cycles. <i>Nature Geoscience</i> <b>3</b>, 572–577 (2010). 14.       D. E. Horton, C. J. Poulsen, Paradox of late Paleozoic glacioeustasy. <i>Geology</i> <b>37</b>, 715–718 (2009). 15.       I. P. Martini, Ed., <i>Late glacial and postglacial environmental changes : Quaternary, Carboniferous-Permian, and Proterozoic</i> (Oxford University Press, 1997). 16.       F. W. Murray, On the Computation of Saturation Vapor Pressure. <i>Journal of Applied Meteorology</i> <b>6</b>, 203–204 (1967). 17.       A. J. Krause, <i>et al.</i>, Stepwise oxygenation of the Paleozoic atmosphere. <i>Nature Communications</i> <b>9</b> (2018).
description: Scripts for porting GENESIS V3 data
 to <em>Paleo</em>-BGC input can be found at github.com/wjmatthaeus/bgc_utils.  <em>Paleo</em>-BGC and its associated parameterizations can be found at github.com/josephdwhite/paleo-bgc. GENESIS V3 files are named for their history variable, CO<sub>2</sub> and O<sub>2</sub> levels, simulation year, and resolution (i.e. PRECIP_546.28_41.2x2.nc). <em>Paleo</em>-BGC output is named for the glacial scenario (GLAC/IGLAC), CO<sub>2</sub> and O<sub>2</sub> levels, longitude and latitude, plant type, and fossil parameter values. Filenames were used to parse <em>Paleo</em>-BGC output into summary data structures in R for analysis. Please note: <em>Paleo</em>-BGC output is in a file structure that associates each model run with its parameterization an initialization files. This output is best processed using unix based command line utilities provided on github above.