Estimating the Greening Effect on Florida Citrus Charles B. Moss 1 and Maria Bampasidou 1 1 University of Florida March 26, 2014
1 Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 2 Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation 3 4
Citrus Maladies Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 Agronomic/Weather Freezes January 1981, December 1983, January 1985, December 1989 Citrus Canker Most recently beginning in 2000 (declared uncontrollable in 2006) Huanglongbing or citrus greening Economics Increased trade Diet fads Decreased housing prices in Florida
Citrus Greening - The Disease Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 Bacterium Candidatus Liberibacter asiaticus Asiatic citrus psyllid (Diaphorina citri Kywayama) This insect was first found in Florida in 1998, and at that time was considered to be a pest of minor importance since the HLB pathogen was not known to be present. The 2005 discovery of HLB in Florida changed the status of this insect to a pest of great importance (Spann et al. 2014, p. 1). Effect of disease Phloem blockage typically leading to loss of root structure Misshapen and bitter fruit Tree decline and death
Canker - The Other Citrus Disease Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 Causes problems with the marketability of fresh fruit Florida experienced periodic bouts Citrus Health Response Program (CHRP) Handling recommendations 1900 foot kill zones Difficulties Florida hurricanes (2004 Charley, Frances, Ivan and Jeanne) Emergence of backyard citrus CHRP abandoned in 2006.
Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 Age Structure of the Citrus Grove 2011-2012 Year Orange Colored Set Early Midseason Valencia Grapefruit Pre 1968 661.3 348.5 1,073.7 47.4 1968-1977 282.7 186.7 620.5 360.1 1978-1987 4,254.5 427.0 4,531.1 827.8 1988-1990 3,808.1 390.4 6,697.7 717.7 1991-1993 3,283.4 558.5 5,873.9 784.6 1994-1996 1,239.6 175.6 2,430.9 108.3 1997-1999 1,658.5 341.7 3,975.5 108.3 2000-2002 2,225.4 368.3 2,919.6 150.3 2003-2005 2,130.5 230.0 2,188.0 255.2 2006-2008 1,834.0 214.8 2,128.9 194.8 Bearing 21,378.1 3,241.5 32,439.0 3,654.5 2009 669.8 76.0 668.9 114.0 2010 477.9 50.7 641.9 129.9 2011 460.5 54.1 499.4 40.0 Non-bearing 1,608.2 180.8 1,810.2 283.9 Total 22,986.3 3,422.3 32,250.0 3,838.4
Yield by Age in Boxes Citrus Maladies Citrus Greening - The Disease Canker - The Other Citrus Disease Age Structure of the Citrus Grove 2011-2012 Orange Colored Tree Age Non-Valencia Valencia Grapefruit 3-5 0.92 0.81 1.51 6-8 1.77 1.63 2.13 9-13 2.46 2.15 3.02 14-23 3.17 2.22 3.83 24 and over 4.37 3.27 5.21
Yield per Tree (Boxes) Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation We start with the basic hyperbolic tangent model of Zanzig, Moss, and Schmitz z (a it ) = β 0 2 [1 + tanh (β 1 + β 2 a it )] (1) 6 5 4 3 2 1 0 0 5 10 15 20 25 30 35 Tree Age (Years) Predicted Summary
Aggregate Yield Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation Yield if age components are known z T t = N i=3 D it β 0 2 [1 + tanh (β 1 + β 2 i)] (2) where D it is the number of trees in each age group.
Impact of Greening Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation Integrating the greening impact z T t = N i=3 θ it D it β 0 2 [1 + tanh (β 1 + β 2 i)] (3) θ it = g (i) (θ 0 + θ 1 t) where g (i) is a valid probability density function (generalized density function) N g (i) = 1 i=1 (4) g (i) 0 i
Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation Econometric Model Tree Age Structure First defining the total number of trees based on tree ages Taking the discrete change D T t = N D it (5) i=1 N Dt T = [D it D i,t 1 ] i=1 D T t = [D 1t D 1,t 1 ] + [D 2t D 2,t 1 ] + (6) D T t = D 1t + [D 2t D 1,t 1 ] + [D 3t D 2,t 1 ] +
Tree Age Structure Continued Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation Note that D 2t D 1,t 1 are the same treee cohorts We replace this with a random tree loss coefficient D T t = D 1t + κ 2t D 1,t 1 + κ 3t D 2,t 1 + (7) where κ it is a random fraction of one. Dividing through by the total tree numbers at time t yields ( d ln Dt T D 1t D T t ) D 1t D T t ( = d ln Dt T + ) N i=2 N i=2 ( ) Dit κ it ln D i 1,t 1 ( ) Dit κ it ln. D i 1,t 1 (8)
Bayesian Estimation Aggregate Yield Impact of Greening Econometric Model Bayesian Estimation α N ( α, Σ α ) (9) ˆα = L f (α l α, Σ α ) α l l=1 L f (α l α, Σ α ) l=1 θ (t) = θ 1 + θ 2 (t 2006) + θ 3 (t 2006) 2 (10)
Parameter No-Greening Greening β 0 3.877 3.874 β 1-0.153-0.153 β 2 0.659 0.671 θ 1-0.003 θ 2-0.010 θ 3 0.001
The preliminary estimates are showing little aggregate effect of greening on tree yields. Greening may affect tree age cohorts (maybe they are taking diseased trees out). We need to add some things to the model: Tree survival probabilities possibly a function of greening. Possibility that greening affects different tree ages differently.