Data and code from: Species compositions determine the ecosystem services of alternative forest transitions in the boreal-temperate ecotone
Data files
Jul 23, 2026 version files 48.02 KB
-
Ecotone_Data_ANA.r
19.07 KB
-
Ecotone_forest_multifunctionality_data_F.xlsx
25.03 KB
-
README.md
3.91 KB
Abstract
Boreal forests are integral for the maintenance of multiple ecosystem services to humanity, yet climate change, wildfire regimes, and land-use intensification are rapidly driving these conifer-dominated systems toward a boreal-temperate ecotone, characterized by a mixture of conifers and broadleaves, with uncertain implications for sustaining multiple ecosystem functions (i.e., multifunctionality). Here, we conducted a standardized 34-year field manipulative experiment in a boreal-temperate ecotone region of Northeast Asia, including four monocultures (conifers and broadleaves) and four two-species mixed forests (each combining coexisting broadleaf and conifer trees). We aimed to assess the ecological impacts of these forest transitions in the boreal-temperate ecotone. Our results showed that ecosystem multifunctionality depended strongly on tree species composition. Crucially, specific compositions generated distinct functional trade-offs rather than uniform changes in ecosystem multifunctionality. For example, transitioning to incompatible mixed-species stands maintained provisioning functions like wood production but resulted in significant reductions in critical regulating services, such as soil carbon sequestration, nutrient cycling, and water regulation. We further provide evidence that these functional declines were mechanistically driven by selection effects that superseded niche complementarity. Specific tree-compositions imposed constraints on soil resources and microbial properties, strongly influencing ecosystem services in the boreal-temperate ecotone. This research provides new insights into identifying species compositions that enhance ecosystem services in the boreal-temperate ecotone, thereby guiding climate change mitigation and forest management efforts in these sensitive regions.
Dataset DOI: 10.5061/dryad.z34tmpgwd
Description of the data and file structure
These data and code come from a standardized 34-year field manipulative experiment in the boreal-temperate ecotone of Northeast Asia, established in 1987, including four monocultures (conifers and broadleaves) and four two-species mixed forests (each combining coexisting broadleaf and conifer trees). The experiment was designed to assess the ecological impacts of these forest transitions in the boreal-temperate ecotone.
Files and variables
File: Ecotone_Data_ANA.r
Description: Data analysis of the R code
File: Ecotone_forest_multifunctionality_data_F.xlsx
Description: The dataset. Two sheets: Meta-data holds the data, one row per plot (24 plots x 42 variables); Notes holds the definition and unit of every column. Note: NA is used to represent missing data or "Not Applicable."
Variables
-
Identifiers
Observation— plot identifier (1-24)Forest— forest composition.JM= Juglans mandshurica,FM= Fraxinus mandshurica (broadleaf monocultures);PK= Pinus koraiensis,LG= Larix gmelinii (conifer monocultures);JM*PK,JM*LG,FM*PK,FM*LG(two-species broadleaf x conifer mixtures)Forest_type—B= broadleaf monoculture,C= conifer monoculture,M= two-species mixtureSample— replicate plot (1-3)Depth— soil layer sampled (cm); 0-10 for all soil variables -
Ecosystem services (min-max standardized, unitless, 0-1)
Soil carbon sequestration,Wood production,Nutrient cycling,Organic matter decomposition,Microbial habitat,Water regulationEMF— ecosystem multifunctionality; mean of the six services aboveWeighted EMF— weighted-average ecosystem multifunctionality -
Measured ecosystem functions
SOC— soil organic carbon (g kg-1)Plant tree area— stand basal area (m2 ha-1)TN— soil total nitrogen (g kg-1)TP— soil total phosphorus (g kg-1)AP— available phosphorus (mg kg-1)NO3— nitrate nitrogen (mg kg-1)NH4— ammonium nitrogen (mg kg-1)AcidP— acid phosphatase activity (umol g-1 h-1)BG— beta-1,4-glucosidase activity (umol g-1 h-1)BX— beta-1,4-xylosidase activity (umol g-1 h-1)NAG— N-acetyl-glucosaminidase activity (umol g-1 h-1)PPO— phenol oxidase activity (umol g-1 h-1)PER— peroxidase activity (umol g-1 h-1)MBC— microbial biomass carbon (mg kg-1)MBN— microbial biomass nitrogen (mg kg-1)MBP— microbial biomass phosphorus (mg kg-1)Max MOI— maximum water-holding capacity (%) -
Explanatory variables
MOI— soil moisture (%)BD— soil bulk density (g cm-3)pH— soil pH in water (unitless)DOC/SOC— dissolved organic carbon : soil organic carbon (ratio)C/N— soil carbon : nitrogen (ratio)N/P— soil nitrogen : phosphorus (ratio)MBC/N— microbial biomass carbon : nitrogen (ratio)MBN/MBP— microbial biomass nitrogen : phosphorus, computed as MBN / MBP (ratio)QCO2— microbial metabolic quotient (qCO2); basal respiration per unit microbial biomass carbon (ug CO2-C mg-1 MBC h-1)litter— annual litterfall, i.e. litter production (g m-2 yr-1)root— fine-root production measured with ingrowth cores (6.75 cm diameter, 30 cm long, filled with root-free sand; ten cores per plot, installed March 2021 and retrieved March 2022). Integrates 0-30 cm, unlike the soil variables (g m-2 yr-1)above litter— forest-floor litter standing mass (g m-2)
Code/software
My submission included the R code for this analysis.
