| Title: | Reconstruct Historical Forest Stand Conditions |
|---|---|
| Description: | Reconstructs forest stand basal area and stem density at user-specified reference years using tree-level inventory data. The method estimates missing tree ages, back-calculates diameters using species-specific growth rates, and incorporates decay-class-based decomposition to infer mortality timing for dead trees. Results are returned in a tidy long format suitable for analysis and visualization. |
| Authors: | Kurt Riggin [aut, cre, cph] (ORCID: <https://orcid.org/0009-0004-4700-1251>) |
| Maintainer: | Kurt Riggin <[email protected]> |
| License: | MIT + file LICENSE |
| Version: | 0.1.0.9000 |
| Built: | 2026-05-14 09:26:36 UTC |
| Source: | https://github.com/kriggithub/standrecon |
Reconstruct stand conditions from a given reference year
standrecon( data, meas_year, ref_year, avg_inc_vec, plot_size, n_plots = 1, bark_eq_list = list(), species_col = "Species", age_col = "Age", dbh_col = "DBH", status_col = "Status", decay_col = "Decay", percentiles = c(0.25, 0.5, 0.75), min_dbh = 5 )standrecon( data, meas_year, ref_year, avg_inc_vec, plot_size, n_plots = 1, bark_eq_list = list(), species_col = "Species", age_col = "Age", dbh_col = "DBH", status_col = "Status", decay_col = "Decay", percentiles = c(0.25, 0.5, 0.75), min_dbh = 5 )
data |
Data with species, age, DBH, tree status, and tree decay columns. |
meas_year |
Year that data was measured. |
ref_year |
Reference years to reconstruct stand conditions. |
avg_inc_vec |
Vector of average increment growth per species with units mm/year. (PIEN = 0.5, ABBI = 0.3) |
plot_size |
Size of plots in m². |
n_plots |
Number of plots with plot_size. |
bark_eq_list |
List of bark correction equations. Default bark corrections for species codes used otherwise. list(PIPO = function(x) x*1.1029 + 0.7162) |
species_col |
Species column name in data as a string. |
age_col |
Species column name in data as a string. |
dbh_col |
DBH column name in data as a string. |
status_col |
Species status column name in data as a string. |
decay_col |
Species decay column name in data as a string. |
percentiles |
Percentiles to be calculated for sensitivity analysis. |
min_dbh |
Minimum DBH threshold for stand calculations in cm. |
A data frame with one row per species, reference year, and percentile. Columns include output type (reconstructed or measured), reference year, percentile, species code, basal area per hectare, and stem density per hectare.
data(standrecon_example_data) out <- standrecon( data = standrecon_example_data, meas_year = 2025, ref_year = c(1950, 1975), avg_inc_vec = c(PIEN = 0.5, ABBI = 0.3, PIPO = 0.4), plot_size = 1000) outdata(standrecon_example_data) out <- standrecon( data = standrecon_example_data, meas_year = 2025, ref_year = c(1950, 1975), avg_inc_vec = c(PIEN = 0.5, ABBI = 0.3, PIPO = 0.4), plot_size = 1000) out
A small synthetic data set for examples and tests.
standrecon_example_datastandrecon_example_data
A data frame with 200 rows and 5 variables:
Species code.
Tree age in years. Values that are known typically come from cored trees.
Tree status code.
Decay class for dead trees.
Diameter at breast height in cm.
Simulated data set.