Open Access

Spatial Bioinvasion Externalities with Heterogeneous Landowner Preferences: A Two-Agent Bioeconomic Model

Shady S. Atallah

Article Figures & Data

  • Figure 1

    Short-Distance Dispersal (SDD) and Long-Distance Dispersal (LDD1 and LDD2)

    Note: Shaded cells represent invaded cells on the grid. Lightning bolts represent random long-distance dispersal events that are exogenous to the model (LDD1).

  • Figure 2

    Optimal Bioinvasion Control in Forestlands GA and GB with Market Ecosystem Service Preferences: left, Bioinvasion Starts on GA, Landowner A Moves First; right, Bioinvasion Starts on GB, Landowner B Moves First

  • Figure 3

    Optimal Bioinvasion Control with Nonmarket Ecosystem Service Preferences for Landowner A and Market Ecosystem Service Preferences for Landowner B: (a) Bioinvasion Starts on GA, A Moves First; (b) Bioinvasion Starts on GB, B Moves First

  • Table 1

    Model Parameters

    ParameterDescriptionValueUnit
    Panel A. Objective Function Parameters
    a0Intercept of the consumer surplus (CS) function−2.97a$/day
    b0Intercept of the user days (UD) function9.32aDays/year
    a1Linear parameter of the CS function0.24a$/day/healthy tree
    b1Linear parameter of the UD function0.24aDays/healthy tree
    a2Quadratic parameter of CS function−0.00017aUnitless
    b2Quadratic parameter of UD function−0.0002aUnitless
    pTimber prices for young and mature trees0.11, 0.13b$/board feet (BF)
    ym,n,tTimber yields for mature trees297cBF/tree
    cA, cBBioinvasion treatment cost4d$/invaded tree
    ρDiscount factor0.961538Year−1
    TForward-looking years (> 5 stages)5, 10,…, 25eYears
    Panel B. Spatial-Dynamic Externality Parameters
    αShort distance H to Eu transition rate0.69Year−1
    Ni,j,tNumber of contiguous Im or Ih neighbors0,…,8Cells
    γLong-distance power-law parameter3Unitless
    L1Average waiting time between invaded-undetectable and invaded-detectablef1Years
    L2Average waiting time between invaded-detectable to invaded-moderateg1Years
    L3Average waiting time between invaded-moderate to invaded-highh2Years
    τyoungYears between juvenile and young if tree and its immediate neighbors are healthy10Years
    τmatureYears between young and mature10Years
    I×JGrid GA dimensions49 × 16 = 784Rows × columns
    M×NGrid GB dimensions49 × 16 = 784(Trees per acre)
    • a Values are from Rosenberger et al. (2013), adjusting for tree density per acre in the case of a1 and b1.

    • b New Hampshire Department of Revenue Administration (2016).

    • c Hepp et al. (2015); Smalley et al. (2016).

    • d Ending the analysis at 25 years is consistent with landowner demographics and plans (most private forest landowners are above 65 years of age and plan to sell or transfer their land; U.S. Forest Service 2015).

    • e Lee (2017). This assumes repeated cutting in the spring, summer, and fall to avoid regrowth.

    • f Embedded Image.

    • g Embedded Image.

    • h Embedded Image.

  • Table 2

    Landowners’ Payoffs and Optimal Amount of Control: Nonmarket Ecosystem Service Preferences Case

    DMCPSocial Cost of Externality
    PayoffsControlPayoffControlDM with Respect to CP
    $/AcreTreatments$/2 AcresTreatments$/2 Acres%
    First mover841 (77)02,203 (82)783 (28)n/an/a
    Second mover1,591 (123)02,509 (64)682 (25)n/an/a
    Total2,432 (145)04,712 (104)1,465 (38)−2,280***−48
    • Note: Standard deviations are in parentheses. CP = central planner; DM = decentralized management.

    • *** p < 0.01.

  • Table 3

    Landowners’ Payoff and Optimal Amount of Control: Market Ecosystem Service Preferences Case

    DMCPSocial Cost of Externality
    PayoffsControlPayoffControlDM with Respect to CP
    $/AcreTreatments$/2 AcresTreatments$/2 Acres%
    First mover12,491 (2,885)1,203 (12)15,467 (3,346)944 (33)n/an/a
    Second mover15,318 (3,114)962 (13)15,847 (3,318)744 (28)n/an/a
    Total27,808 (4,245)2,165 (18)31,314 (4,712)1,688 (43)−3,506***−11
    • Note: Standard deviations are in parentheses. CP = central planner; DM = decentralized management.

    • *** p < 0.01.

  • Table 4

    Landowners’ Payoffs ($) and Optimal Amount of Control (Treatments) under Decentralized Management (DM) and Central Planner (CP) Cases for Landowner A and Landowner B

    DMCPSocial Cost of Externality
    PayoffsControlPayoffControlDM with Respect to CP
    $/AcreTreatments$/2 AcresTreatments$/2 Acres%
    Panel A. GA and GB Ecologically Isolated (No Long-Distance Dispersal)
    GA (NMES)2,999 (42)498 (7)n/a498 (7)n/an/a
    GB (MES)16,731 (3,384)497 (9)n/a497 (9)n/an/a
    GA + GB19,731 (3,384)995 (11)19,731 (3,384)995 (11)00
    Panel B. Bioinvasion Starts on GA, Landowner A (NMES) Moves First
    GA (NMES)831 (73)0 (0)2,187 (152)836 (32)n/an/a
    GB (MES)10,432 (2,453)1,543 (12)16,667 (3,367)700 (29)n/an/a
    GA + GB11,263 (2,453)1,543 (12)18,854 (3,370)1,536 (43)−7,591***−40
    Panel C. Bioinvasion Starts on GB, Landowner B (MES) Moves First
    GA (NMES)1,569 (123)0 (0)2,119 (72)812 (26)n/an/a
    GB (MES)10,490 (2,121)785 (0)15,127 (3,357)1,025 (33)n/an/a
    GA + GB12,060 (2,125)785 (0)17,246 (3,358)1,837 (42)−5,187***−30
    • Note: Landowner A = nonmarket ecosystem service preferences (NMES); landowner B = market ecosystem service preferences (MES). Standard deviations are in parentheses.

    • *** p < 0.01.

  • Table 5

    Net Aggregate Payoffs ($) for Different Levels of Uniform and Nonuniform Subsidies for Landowner A and Landowner B

    Cost (cA, cB)Subsidy (PA, PB)Control in GAControl in GBAggregate ControlAggregate PayoffsTotal Subsidy CostaNet Aggregate Payoffsa,bChange Relative to Baseline
    $/TreatmentTreatments$/2 Acres
    Panel A. Bioinvasion Starts on GA, Landowner A (NMES) Moves First
    (4, 4)(0, 0)01,5431,54311,263011,2630%
    (1, 1)(3, 3)1,3581,1802,53720,9924,65316,33945%
    (2, 2)(2, 2)1,0418651,90519,5472,53217,01551%
    (1, 4)(3, 0)1,0426031,64419,5462,01617,52956%
    Panel B. Bioinvasion Starts on GB, Landowner B (MES) Moves First
    (4, 4)(0, 0)078578512,060012,0600%
    (2, 2)(2, 2)1,3901,6493,03817,7513,76513,98616%
    (1, 1)(3, 3)1,2251,6732,89819,8025,56214,24018%
    (1, 4)(3, 0)1,4801,2072,68617,5472,92114,62621%
    • Note: Landowner A = nonmarket ecosystem service preferences (NMES); landowner B = market ecosystem service preferences (MES). Standard deviations are not shown.

    • a In present value terms (25 years, 4% annual discount rate).

    • b This is defined as aggregate payoffs minus total subsidy.