## Overview

We established the Forest Resilience Threshold Experiment (FoRTE) with the goal of understanding how and why disturbance severity and type affect key C cycling processes and ecosystem C balance.

## Site Description

FoRTE is conducted at the University of Michigan Biological Station (UMBS) in northern lower Michigan, USA (45.56 N, -84.67 W), where mean annual air temperature is 5.5°C and mean annual precipitation is 817 mm (Gough and others 2013). The vegetation communities included in our analysis are transitioning 100-yr-old middle successional forests, with a dominant upper canopy comprised mainly of pioneer bigtooth and trembling aspen (Populus grandidentata and P. tremuloides, respectively) and paper birch (Betula papyrifera) established following region-wide harvesting and fire in the early 20th Century. This cohort of early successional tree species is rapidly declining (Gough and others 2010), and giving way to later‐successional red oak (Quercus rubra), eastern white pine (Pinus strobus), sugar maple (Acer saccharum), red maple (Acer rubrum), and American beech (Fagus grandifolia).

Despite sharing a relatively uniform disturbance history, tree species’ dominance and plant community composition along with structure and complexity vary substantially among the “landscape ecosystems” of UMBS (Table 1; Figure 2; Scheuermann and others 2018). The term “landscape ecosystems” is used to describe land units forming a unique complex of climate, soils, biota and landform (Pearsall et al. 1995).

## Factorial Design, Disturbance Severities and Types

During Spring 2019, we randomly assigned factorial combinations of disturbance severity at 4 levels (0, 45, 65, 85%) targeted leaf area reductions) and two treatment types (top-down - starting at the largest diameter canopy tree; bottom-up, starting at the smallest diameter canopy tree–with the caveat there is an 8 cm DBH cut-off for a tree to be considered a canopy tree) within four separate landscape ecosystems representing a range in site productivity characteristic of the forests of the Upper Great Lakes Region (Lapin and Barnes 1995), here represented by the replicate names (A, B, C, D; see map below). The code used to generate the random assignments is illustrated below:

#set random seed
set.seed(9925)

#define parameters
w <- c("top", "bottom")
x <- c("E", "W")
y <- c(1, 2, 3, 4)
z <- c(0, 45, 65, 85)

# example treatment using group A
# all other groups simulated in identical fashion
# then merged at the end of file and written to .csv
plot.disturbance.a <- data.frame(group = "A", plot = sample(y), disturbance = sample(z))
plot.disturbance.a <- plot.disturbance.a[order(plot.disturbance.a\$plot),]

# randomization of how which side E = east or W = west is chosen for
# top = top-down or bottom - bottom-up treatment
plot.treatment <- data.frame(plot.side = sample(w), treatment = sample(x))

Circular, 0.1 ha sampling subplots were established within each disturbance severity x type treatment in each of the four replicates, resulting in a total of 32 subplots surrounded by a 5-m wide measurement-free treatment buffer. Within each subplot and following methods detailed below, we measured ANPPw, Rs, and Asat, and we used terrestrial (i.e., ground-based) lidar and hemispherical imagery along with stem census data to derive measures of diversity, structure and complexity.

## Stem-Girdling

Disturbance was implemented via stem girdling, which, like wood boring insects (e.g., emerald ash borer, mountain pine beetle), causes gradual defoliation by eliminating the transport of photosynthate to roots, killing trees once carbohydrate reserves are exhausted over a period of two to three years (Dietze and others 2014; Gough and others 2013). Stem-girdling occurred in late May 2019 (DOY 140 to 143), prior to leaf-out.

## Mortality Assignment

We selected approximately 3600 trees to for girdling. 1207 of these trees were located in censused and mapped subplots 0.1 ha), with an additional 2400 trees located in the plot areas outside of subplot locations as well as in the 5 m buffer zones around all plots. To clarify, as mentioned above, each plot is 0.5 ha and contains 2 seperate 0.1 ha subplots. Only trees in subplots are censused. The estimate of 2400 trees girdled outside of the censuseed areas is based on extrapolation from plot-specific stem densities.

Mortality was assigned based on disturbance severity–0, 45, 65, 85%–where the disturbance levels are targetted reductions in leaf area. Estimated leaf area per stem was derived from species- and site- or region-specific allometries (Gough and others 2008) relating DBH to leaf area for all canopy trees (trees >8 cm DBH). to target gross LAI reductions within each plot of 0% (control), 45%, 65%, or 85%. Severity levels were assigned at random to four, 0.5 ha circular whole plots, which were split into 0.25 ha halves and randomly designated a “top-down” or “bottom-up” disturbance type treatment (Figure 2B). For the “top-down” disturbance type, we girdled the largest trees first, irrespective of species, starting with the highest leaf area individual and sequentially girdling lower leaf area trees until the assigned plot disturbance severity was reached. For the “bottom-up” treatment, we stem girdled trees (>8 cm DBH) with the lowest individual leaf area first, progressively girdling larger trees up to the targeted disturbance severity. A 10-cm band of phloem tissue was girdled ~1 m above the forest floor and below DBH using pry bars.

Individual mortality assignments at the stem level can be accessed via fd_mortality()