Lab 6 Gypsy Moth Flight

Lab. 6. Phenology of Gypsy Moth Flight
The goal of this lab is to develop a model that can be used to predict the timing of the flight of gypsy moth males. Male moths numbers are monitored using pheromone traps. Pictures below show the trap and dead males removed from the trap.
Predicting the dates of male moth flight is important for selecting time when the traps should be placed and removed at various geographic locations
Save file MOTH CATCH Traps
Jul. Day No. 1 No. 2 No. 3 No. 4 No. 5 No. 6 No. 7 No. 8 No. 9 No. 10 No. 11 No. 12 No. 13 No. 14 No. 15 No. 16
6/10/96 162 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
6/13/96 165 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 0
6/17/96 169 0 0 0 0 0 0 0 0 0 0 0 0 0 0 175 0
6/21/96 173 0 0 0 0 0 1 3 0 0 0 0 0 0 0 240 2
6/24/96 176 0 0 0 0 0 0 2 0 0 0 0 0 0 0 233 42
6/28/96 180 0 0 0 0 26 1 9 0 0 0 0 0 0 0 208 75
7/1/96 183 0 0 0 0 10 0 16 3 1 0 0 1 6 0 102 40
7/5/96 187 9 1 3 2 10 15 5 4 13 8 7 8 13 2 71 29
7/8/96 190 19 0 0 3 6 8 16 36 2 3 3 0 8 3 50 20
7/12/96 194 51 0 2 0 2 1 12 0 1 0 12 0 1 2 81 7
7/15/96 197 20 1 0 0 1 6 4 0 4 8 65 1 11 23 73 4
7/19/96 201 14 6 17 1 0 0 0 0 14 15 67 16 27 67 1 0
7/22/96 204 17 21 10 15 0 0 6 3 78 99 91 37 33 13 0 0
7/26/96 208 16 19 16 47 0 0 0 5 60 64 62 24 33 54 0 0
7/29/96 211 2 4 4 48 0 0 0 26 236 265 142 42 52 14 0 0
8/2/96 215 0 3 1 74 0 0 0 32 135 59 29 18 23 9 0 0
8/5/96 218 2 2 0 49 0 0 0 124 335 129 75 20 10 15 0 0
8/9/96 222 1 1 0 48 0 0 0 141 216 70 17 10 7 3 0 0
8/12/96 225 0 1 0 16 0 0 0 91 124 36 11 3 3 1 0 0
8/16/96 229 0 0 0 6 0 0 0 138 73 14 1 0 2 0 0 0
8/19/96 232 0 0 0 3 0 0 0 241 87 14 1 0 0 0 0 0
8/23/96 236 0 0 0 6 0 0 0 65 19 5 1 0 0 0 0 0
8/26/96 239 0 0 0 1 0 0 0 10 0 1 0 0 0 0 0 0
8/29/96 242 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 0
* ELEV: 592 666 575 696 61 98 212 968 816 627 586 456 429 421 0 1
* LONG: -79.31291 -79.60217 -79.73571 -79.77274 -77.5196 -78.43374 -79.10205 -80.09523 -80.15904 -80.22949 -80.29023 -80.19657 -80.15506 -80.38629 -75.91171 -75.71411
* LAT: 37.73114 37.88688 38.05035 38.24153 37.46327 37.33706 37.38915 38.72598 38.77737 38.75695 38.81726 38.97192 39.02687 39.01885 36.27126 36.08459

Weather Stations
Name Long: Lat: Elev:
ELIZABETH_CITY -76.18 36.32 0
ASHLAND_1_SW -77.47 37.75 61
FARMVILLE_2_N -78.37 37.32 105
GATHRIGHT_DAM -79.93 37.95 574
LEXINGTON -79.42 37.77 325
LYNCHBURG_WSO_AIRPORT -79.18 37.32 267

TEMPERATURE DATA
ELIZABETH_CITY ASHLAND_1_SW FARMVILLE_2_N GATHRIGHT_DAM LEXINGTON LYNCHBURG_WSO_AIRPORT
Jul.Day T-MIN T-MAX T-MIN T-MAX T-MIN T-MAX T-MIN T-MAX T-MIN T-MAX T-MIN T-MAX
1 5.6 10.6 2.2 3.9 2.2 11.7 -2.2 7.2 3.3 5.6 3.9 6.7
2 5.6 13.3 2.2 4.4 2.8 5 1.7 5.6 3.3 7.8 3.3 4.4
3 6.1 15 0.6 12.8 2.8 3.9 2.8 5.6 0 8.3 -2.2 10.6
4 -3.3 11.7 -7.2 2.8 -6.7 12.2 -8.9 5 -7.8 0 -6.1 2.2
5 -3.9 7.2 -3.9 3.9 -6.7 5.6 -8.9 -1.7 -3.9 5.6 -2.8 6.1
6 -3.9 3.3 -6.1 -2.8 -3.9 7.2 -6.1 3.9 -6.1 1.1 -7.8 -2.2
7 -3.9 7.2 -9.4 -6.1 -8.3 -2.2 -7.8 -0.6 -8.9 -5 -8.3 -5.6
8 -6.7 2.8 -8.9 -0.6 -11.1 -6.1 -10.6 -6.1 -7.8 0 -7.8 0.6
9 -7.2 6.1 -12.2 -1.1 -7.8 2.8 -12.2 -2.2 -12.8 4.4 -6.1 5
10 -2.8 6.1 -7.2 3.9 -5 3.3 -9.4 2.8 -2.8 4.4 -6.1 2.8
11 -7.2 3.9 -15 -2.2 -14.4 7.8 -14.4 -0.6 -
……………………..
http://home.comcast.net/~sharov/PopEcol/lab6/phenol.txt
………………………
348 -1.7 3.3 -1.7 3.3 -1.7 3.3 -1.7 3.3 -1.7 3.3 -1.7 3.3
349 0 11.1 0 11.1 0 11.1 0 11.1 0 11.1 0 11.1
350 8.3 16.1 8.3 16.1 8.3 16.1 8.3 16.1 8.3 16.1 8.3 16.1
351 -5.6 11.7 -5.6 11.7 -5.6 11.7 -5.6 11.7 -5.6 11.7 -5.6 11.7
352 1.1 11.1 1.1 11.1 1.1 11.1 1.1 11.1 1.1 11.1 1.1 11.1
353 -6.1 3.3 -6.1 3.3 -6.1 3.3 -6.1 3.3 -6.1 3.3 -6.1 3.3
354 0.4 1.4 0.4 1.4 0.4 1.4 0.4 1.4 0.4 1.4 0.4 1.4
355 -5 3.9 -5 3.9 -5 3.9 -5 3.9 -5 3.9 -5 3.9
356 -5 2.8 -5 2.8 -5 2.8 -5 2.8 -5 2.8 -5 2.8
357 -12.8 1.1 -12.8 1.1 -12.8 1.1 -12.8 1.1 -12.8 1.1 -12.8 1.1
358 -5 2.8 -5 2.8 -5 2.8 -5 2.8 -5 2.8 -5 2.8
359 -5.6 2.8 -5.6 2.8 -5.6 2.8 -5.6 2.8 -5.6 2.8 -5.6 2.8
360 -3.3 3.3 -3.3 3.3 -3.3 3.3 -3.3 3.3 -3.3 3.3 -3.3 3.3
361 -5.6 1.7 -5.6 1.7 -5.6 1.7 -5.6 1.7 -5.6 1.7 -5.6 1.7
362 -9.4 2.8 -9.4 2.8 -9.4 2.8 -9.4 2.8 -9.4 2.8 -9.4 2.8
363 -8.9 3.9 -8.9 3.9 -8.9 3.9 -8.9 3.9 -8.9 3.9 -8.9 3.9
364 -8.3 8.3 -8.3 8.3 -8.3 8.3 -8.3 8.3 -8.3 8.3 -8.3 8.3
365 -1.7 12.8 -1.7 12.8 -1.7 12.8 -1.7 12.8 -1.7 12.8 -1.7 12.8
and open it from the Microsoft Excel. This file contains the following information:
• Moth catches in pheromone traps located in Virginia, West Virginia, and North Carolina in 1996. Traps were checked twice a week and moths were removed from the trap each time.
• Geographic coordinated of traps: longitude and latitude in decimal degrees, and elevation at trap location.
• Minimum and maximum daily temperature at 6 weather stations in Virginia, West Virginia, and North Carolina in 1996.
• Geographic coordinated of weather stations: longitude and latitude in decimal degrees, and elevation.
You need to build the model of gypsy moth flight from literature data and validate it using the given data set. The lab include the following steps: build the model; simulate the dates of gypsy moth flight that are expected at locations of weather stations; using the regression of moth flight dates vs. elevation, latitude, and longitude predict moth flight at trap locations; compare actual and predicted moth flight dates.
1. Use data from Casagrande et al. (1987) on gypsy moth larval and pupal development (table below). Plot development rate vs. temperature and estimate the lower temperature threshold and degree-days separately for larval and pupal periods.
Table. Development time of gypsy moth males
Temperature Larval Pupal
C development development
======================================
13 93 -
18 55 27
20 44 18
25 30 16
28 28 12
30 27 12
======================================
2. Use the model of Johnson et al. (1983) for egg development: eggs start hatching after accumulating of 282 degree-days since January 1 (lower temperature threshold = 3 degrees C).
3. Build the model as follows. Put weather data from any station in the first 3 columns of a new worksheet: Julian day which is the day number since January 1, minimum temperature, and maximum temperature. The next column will be used for accumulating degree-days for the egg stage. Put zero in the upper cell (D2) of this column because initially no degree days have been accumulated. In the cell below (D3) put the following equation: =D2+0.5*(MAX(0,B2-3)+MAX(0,C2-3)). This equation implements the "rectangular method for accumulating degree days: it assumes that the minimum temperature (B2) has lasted for half day, and then the maximum temperature (C2) has lasted for another half day. Function MAX(0,B2-3) gives the effective temperature: if the temperature (B2) is greater than 3 (the lower temperature threshold) then the effective temperature = B2-3; if B2 < 3, then the effective temperature = 0 (because 0 is greater than B2-3).
4. In column E estimate degree-days for larvae. Degree-days for larvae start accumulating only when accumulated egg degree-days exceed 282 (model of Johnson et al. 1983). Thus, use the following equation for cell E3:=E2+(C2>282)*(effective temperature for larvae). The term (C2>282) equals 1 if C2 > 282, otherwise it equals 0. The equation for the effective temperature is the same as for the egg stage (write it yourself!). Use the lower temperature threshold for larval stage that you have estimated before (section 1).
5. In column F estimate degree-days for pupae. Degree-days for pupae start accumulating only when accumulated larval degree-days exceed the required amount of degree days which you have estimated before (section 1). You the equation similar to that in column E (write it yourself!). Use the lower temperature threshold for pupal stage that you have estimated before (section 1).
6. Adult males start emerging on the day when accumulated pupal degree-days exceed the amount required for pupal development, which you have estimated before (section 1). However, they continue emerging later for about 1 week. The model for egg phenology (Johnson et al. 1983) simulates the start of egg hatch, but many eggs hatch later than predicted by this model. Individuals that hatch later will have a delay in their development through all stages. This delay can be simulated by tracing the phenology of several cohorts that start their larval development at varying dates. However, we will use a simplified approach. We already found the date of the beginning of male moth flight. Now we will assume that emergence period is extended from this date till additional 150 degree-days are accumulated (for pupae). The cumulative proportion of emerged males, P, is estimated as P = (D - S)/150, where D = degree-days currently accumulated by pupae and S = degree required for the start of moth emergence. Also, P = 0 if D < S, and P = 1 if D > S+150. Combining all these conditions we get the equation that you need to put in cell G2: =MIN(150,MAX(0,F2-300))/150. Autofill it to the entire column.
7. The number of males caught in pheromone traps is assumed to be proportional to the density of males that are alive. We will assume that all males live for 10 days. Then the density of males that are alive is proportional to the difference between the cumulative proportion of emerged males, P(t), and the cumulative proportion of males emerged 10 days ago, P(t-10). Thus, estimate the value N(t) = P(t) - P(t-10) in column H.
8. Estimate the mean date (Julian day) of moth flight using the equation F1:
where N(t) was estimated in section 7. Use function SUMPRODUCT(array1, array2) to estimate the numerator. This function returns the sum of products of all elements of array1 and array2.
9. Estimate the mean moth flight date using weather data from all 6 weather stations. Estimate the regression of the date from elevation, latitude and longitude of these weather stations. Use this regression equation to predict the mean moth flight date at all trap locations.
10. Estimate the mean date of moth flight recorded by all traps (use the same equation for the mean; use the middle of the period between trap checks as t). Plot actual mean date of moth flight versus predicted. Use the diagonal line to check the accuracy of predictions. Does the model fit the data? Is there any bias in model predictions?
11. (Optional) How flight periods will change with global warming by 2 degrees?