Unintended biological invasions: does risk vary by trading partner?

1. Unintended biological invasions: does risk vary by trading partner? C. Costello, M. Springborn, C. McAusland, & A. Solow JEEM 2007

2. The non-indigenous species problem Pecuniary damage • ~50,000 non-indigenous species introduced to the U.S. (Pimentel, 2005) • Estimation of yearly U.S. monetary losses: • $4.7 – 6.5 billion (Office of Technology Assessment, 1993) • $120 billion (Pimentel, 2005) Ecological damage • ~400 of the 958 species on the Endangered Species list are there primarily because of NIS (Wilcove et al. 1998) • SF Bay: no shallow water habitat remains uninvaded Asian Tiger Mosquito – introduced via tire shipments The Atlantic Green Crab, first seen in S.F. Bay in 1989, preys on juvenile Dungeness Crab.

3. Invasive Species - Rules of Thumb • Disturbed land is more susceptible to invasion (e.g. agricultural use versus primeval forest) • Trade in goods and services provides platform for unintentional introductions of non-native species (esp. agricultural imports, shipping and packing materials, ballast water, tourism) • Successful introductions are facilitated by bio-geographic similarities between host and source region

4. Invasive Species - Rules of Thumb continued. • Likelihood that an “arrival” will become established is increasing in the number of times the species is “exposed” to host region • “tens rule”: 10% of introduced species become casual , 10% of these become established (10% of these become a pest) • A Source’s potential pool of exotics is finite---sampling without replacement • Newly arrived exotics aren’t usually discovered for quite some time (chance, damages high, systematic species survey)

5. Unintentional introductions of NIS • Market failure stemming from foreign trade • Ballast water, Packing materials, Ship hulls, Hitchhiking with traded goods • Policy response depends on whether: • Most NIS already here (reactive policies) • Lots of new NIS likely to arrive (proactive policies)

6. Blunt Policy Response • Jenkins (1996) • “broad tools such as bans or restrictions on imports may be necessary to protect biodiversity”

7. Advisable?

8. Costello, Christopher and Carol McAusland. 2003. "Protectionism, Trade, and Measures of Damage from Exotic Species Introductions" American Journal of Agricultural Economics, 85(4) 2003: 964-975. • Build a dynamic general equilibrium model with stochastic introductions and damages • Arrivals increasing in volume of imports • Exogenous probability that a new arrival will become established • Damage from a newly established exotic is also a random variable

9. Show that import tariffs have two effects on damages for invasives • Tariffs reduce volume of imports • Shrinks the platform for new arrivals • Tariffs cause protected sector to expand • If country is a net importer of agricultural goods, then tariffs cause agriculture to expand • More agricultural activity • -> more disturbed land on which invasive species can get a foothold • -> more crops that can be damaged

10. Possibility: tariffs on agricultural goods can lead to higher estimates of agricultural damage

11. Shortcomings • Ignores attenuation • Treats all trade partners as identical • No data • Referee: “prove this has happened even once”

12. Example - Sugar • US support for sugar generates US price = 2 x ROW price • is similar to a 100% tariff except no beneficial tariff rent • Since 1934 harvested acreage for all crops in US fell by .1%/annum • over same period, land under sugarcane production grew at average annual rate of 1.6%

13. currently infests 20% of Texan sugarcane believed to have come in on imported goods detected on sugarcane, lemon grass, sorghum and broomcorn imports Texas damages estimated at $10 -$20 million (/yr?) while harvest valued only at $64 million Mexican rice borer Source: Texas A&M U http://insects.tamu.edu/images/insects/color/sorghum/sor069.jpg

14. Costello, Sprigborn, McAusland and Solow (2007) • Use data from San Francisco Bay to address role of attenuation and biogeographic similarity of partners • Asks: How risky is future trade?

15. Complications • Introduction rates may differ by region based on biogeographic differences • Region-specific information could be used to tailor policy • Attenuation: how might the introduction rate change over time and by partner? • Each region has a finite number of species to “contribute” to host region • Once an NIS is established, future introductions pose no risk • The more trade we’ve had with a particular partner, the lower the marginal invasion risk.

16. Complications cont. 3. Discovery lag • discoveries reflect both an introduction process and a discovery process. • Delay may be influenced by observational effort, species population growth rate, level of damage, etc. • Solow and Costello (2004) show that even if introduction rate is constant and detection effort is constant, it will appear as if species are being introduced at an increasing rate if don’t account for discovery lag

17. Model features • We would like to estimate the region-specific invasion risk by trade partner • Must account for: • “Baseline invasion risk” (differs by region) • “Attenuation rate” as function of trade volume (differs by region) • “Discovery lag” (same for all regions) • We will require data on: • Trade volume over time (by region of origin) • Species discovery dates and region of origin

18. Literature treatment of dynamics

19. Literature treatment of dynamics

20. Model: introductions Let Nit measure (unobserved) introductions from region j in year t. Assume Njt has Poisson distribution with mean (rate per unit of imports) (1) βj is region specific “intrinsic infectiousness” of imports to importer sjt=import volume (measured in short tons) from region j in year t Sjt=cumulative import volume from region j through year t γj measures rate at which introductions attenuate (if γ<0) with cumulative import volume ω reflects time trend (e.g. if shipping speeds improve over time, may see more hitchhikers surviving passage) λt λtattenuates at rate γ Cum. shipping, St

21. Discoveries • Yjt defined as number of NIS from source j discovered in year t. • Define put as probability that an NIS introduced in year u is discovered in year t • Yit will have Poisson distribution with mean djt=∑tu=0γjuput (2)

22. Problem: we don’t have data on the discovery process Assume post-introduction waiting time to discovery is geometrically distributed: where π= probability that a species is observed (not necessarily for the first time) in any given year

23. Assumptions • Intrinsic infectiousness (β) and attenuation (γ) are region specific • Time trend (ω) and observation rate (π) are universal

24. Combining introductions & discoveries • Expected number of NIS discovered (for the first time) in year t: • This is “thinning” of a Poisson process • Therefore, number of discoveries in year t is a non-homogeneous Poisson random variable with mean (rate) dt. Arrived at time 1, not discovered until time t Arrived at time 2, not discovered until time t Arrived at time t Discovered immediately

25. Estimator • Our data (by trade region) are • NIS discoveries from that region over time • Shipping volume over time • Given discoveries distributed Poisson with rate dt, the likelihood of observed discovery record is • To estimate, choose β, γ, ωto maximize integrated likelihood • Then estimate “nuisance parameter” π by maximizing

26. Empirical application to SF Bay NIS discovery data (1853-1994): Cohen and Carlton (1995). “Nonindigenous aquatic species in a United States estuary: A case study of the biological invasions of the San Francisco bay and delta” Import data (1856-2000): • 1856-1945: Foreign Commerce and Navigation of the United States • 1946-1967: Foreign Trade through the San Francisco Customs District • 1968-1989: Annual Import Data Bank Files (tape) • 1990-2000: U.S. Imports of Merchandise (CD-ROM) U.S. Department of Commerce, Bureau of the Census records

27. NIS and trade into San Francisco

28. NIS and trade into San Francisco Restrict attention to regions with more than one NIS introduction between 1856 and 1994

29. Results and Predictions Maximum likelihood estimates for three jointly estimated regions (90% Confidence Intervals via Parametric Bootstrap) Using the likelihood ratio test, we reject the hypotheses that • βATM = βWPC • γATM = γWPC = γION Fail to reject that βION is equal to either βATM or βWPC Subsequently we restrict attenuation in WPC infection rate (γWPC=0).

30. Observations • ω=0.15 implies introductions increase 1.5% per year (other things equal) • Implied π=0.048 <-> discovery lag of 13 years. • Attenuation • Even though ATM has higher “intrinsic” infectiousness, attenuation is quick (about 8% per million short tons of trade) • WPC has lower β but reject attenuation • ION attenuation almost instantaneous (106% per million short tons)

31. Calculating Marginal Invasion Risk (MIR) • MIRjt=βjexp(γjSjt+ωt) Calculated MIR as of 1994

32. Alternative hypotheses • Attenuation is a result of global, not regional, trade • Test: re-estimate model using where St=∑jSjt is global trade. Results: fitted betas, gammas and omegas nearly identical, while fitted n is small (0.0011) and statistically insignificant (p= 0.62) (6)

33. Model fit (ATM region) Fitted introductions Fitted discoveries Discovery data Cumulative Import Volume (x106 mt) Estimated Number of Undiscovered Species Marginal Invasion Risk = 0.11

34. “Undiscovered” species • Our model includes a lag between species introduction and discovery • At any point in time, we can estimate the number of undiscovered species • Equals # introduced species - # discovered

35. Forecasting trade • Forecasts of future imports into San Francisco Customs District are taken from Haveman and Hummels (2004). • Their forecasted values are drawn from GTAP, the Global Trade Analysis Project using the Walmsley et al. (2000) extension of the basic GTAP model

36. Predictions of future introductions and discoveries using forecasted trade volumes fitted intro. ATM cumulative trade to 2000 fitted disc. cum. disc. WPC ION

37. Estimated new introductions Predicted number of new NIS

38. So should we restrict trade?

39. Thought experiment • Suppose the US used trade restrictions to reduce by one the expected number of NIS in 2020 originating from each region. • By how much would the US have to reduce imports? • What are the costs of these trade restrictions? • How do they measure against the benefits of avoiding one NIS?

40. Aside • Do Costello, Springborn, Solow and I really want to imply that trade volumes should be curtailed? • No. But the thought experiment puts the damages from NIS into an economic context — lets the reader judge how big a problem trade-facilitated introductions are.

41. Standard formula for DWL from trade restriction Feenstra (2004, p.217) • DWLt≈ -½ Mt2PtMt/ε where • DWLt is deadweight loss from year t trade restrictions, • M is percentage reduction in imports, • ε is elasticity of import demand • ε = -1.23 (Hooper and Marquez 1995) • PM is value of imports. • Use actual M for 1995-2000 and forecasted M for 2001-2020 • Calculate Pi so PM2002,forecasted exactly equals PM2002,actual • PATM=$1769 • PWPC=$3399 ^ ^

42. In order to reduce by one number of expected NIS from ATM by 2020 need to reduce imports from ATM by 90%! • …from WPC…by 2% • Using 5% discount rate, total discounted DWL from restrictions on imports from • ATM = $9,520 million • WPC = $44 million

43. Benefits (B) from reducing trade Bt=D∑ts=1995[λus-λrs] where • D=annual damage from an average NIS • λus=mean introductions in year s when trade is unrestricted • λrs=…restricted

44. In order for costs and benefits of trade restrictions to balance… • would need annual damages from (prevented) NIS to be about • $1,063 million/year (ATM) • $8 million/year (WPC)

45. Damages from an “average” NIS • Pimental et al (2005) • 50,000 NIS present in US • Annual damage from NIS = $120 billion → $2.4 million = Crude estimate of average annual damage

46. What if we knew we were avoiding one of the worst NIS?

47. Teredo Navalis (Atlantic Shipworm) • $205 million/year (structural damage)

48. Zebra Mussel • $700 million/year (clogs intake valves, alters filtration)

49. Asian Clams • $1billion/year • “These small freshwater mussels can be drawn into power plants along with coolant water and clog tubes and pipes, resulting in economic costs (Fuller & Benson, 2003). • “The asian clam will also compete with native clams and mussels for habitat and food, and change benthic substrates (USGS, 2001).” • http://www.biology.duke.edu/bio217/2005/cjc6/Asianclam.html

50. Conclusions • Key policy variable is marginal invasion risk • Possible difference across trade partners • Theoretical reasons to expect attenuation • We find both (to some extent) • Riskiest partners likely to be new partners – not those who have delivered most species in past • Expect: 1.4 (ATM), 52.4 (WPC), and 0 (ION) by 2020 • Crudely restricting trade with either WTC or ATM not advisable