Modeling Differential Equations in Biology

Modeling Differential Equations in Biology

Volumn , Issue , 2008, Pages 1-500

  Taubes, Clifford Henry ^a  

^a HARVARD UNIVERSITY   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

EID: 77951456296     PISSN: None     EISSN: None     Source Type: Book    
DOI: 10.1017/CBO9780511811364     Document Type: Book

References (621)

1. Nonaka, S. et al. Cell 95, 829-837 (1998).
2. Afzelius, B. A. Science 193, 317-319 (1976).
3. Supp, D. M., Witte, D. P., Potter, S. S. & Brueckner, M. Nature 389, 963-966 (1997).
4. Levin, M., Johnson, R. L., Stern, C. D., Kuehn, M. & Tabin, C. Cell 82, 803-814 (1995).
5. Levin, M. et al. Dev. Biol. 189, 57-67 (1997).
6. Danos, M. C. & Yost, H. J. Development 121, 1467-1474 (1995).
7. Levin, M, Roberts, D. J., Holmes, L. B. & Tabin, C. Nature 384, 321 (1996).
8. Sulik, K. et al. Dev. Dyn. 201, 260-278 (1994).
9. Bellomo, D., Lander, A., Harragan, I. & Brown, N. A. Dev. Dyn. 205, 471-485 (1996).
10. Brown, N. A. & Wolpert, L. Development 109, 1-9 (1990).
11. Pagan-Westphal, S. M. & Tabin, C. J. Cell 93, 25-35 (1998).
12. Harvey, R. P. Cell 941a, 273-276 (1998).
13. C. W. Clark, The Optimal Management of Renewable Resources (Wiley, New York, ed. 2, 1990), pp. 16-18
14. B. Dennis, Nat Resour. Model 3, 481 (1989)
15. C.W. Clark, J. Cons. Int. Explor. Mer. 36, 7 (1974)
16. J.A. Gulland, J. Cons. Int. Explor. Mer. 37, 199 (1977)
17. R.M. Peterman, J. Fish. Res. Board Can. 34, 1130 (1977).
18. J. R. Beddington, in Variability and Management of Large Marine Ecosystems, K. Sherman and L. M. Alexander, Eds. (Selected Symposium 99, American Association for the Advancement of Science, Washington, DC, 1986)
19. C.W. Clark, Mathematical Bioeconomics (Wiley-Interscience, New York, ed. 2, 1990)
20. U.S. Department of Commerce, Our Living Oceans (NOAA Technical Memo NMFSF/SPO-2, Washington, DC, 1992).
21. W. C. Allee, The Social Life of Animals (Norton, New York, 1938), p. 239
22. P. A. Larkin, R. F. Raleigh, N. J. Wilimovsky, J. Fish. Res. Board Can. 21, 477 (1964)
23. G.C. Varley, G. R. Gradwell, M. P. Hassell, Insect Population Ecology: An Analytical Approach (Blackwell, Oxford, 1973).
24. R. A. Myers, J. Bridson, N. J. Barrowman, Can. Tech. Rep. Fish. Aquat. Sci. 2024 (1995).
25. R. J. H. Beverton and S. J. Holt, On the Dynamics of Exploited Fish Populations (U.K. Ministry of Agriculture, Fisheries and Food Fishery Investment Series 2, no. 19, 1957)
26. G.G. Thompson, Can. Spec. Publ. Fish. Aquat. Sci. 120, 303 (1993).
27. J. Jakobssen, Rapp. P.-V. Reun. Cons. Int. Explor. Mer. 177, 23 (1980).
28. For the salmonids, there were two statistically significant estimates of the depensation parameter that were greater than 1 and two that were less than 1. The most convincing case of depensation is that of pink salmon in Sashin Creek, AK, at a population abundance of less than 100 females; it is perhaps at this level that depensation is expected to occur for salmonids.
29. G. I. Murphy, Proc. Calif. Acad. Sci. 34, 1 (1966)
30. V.C. Anthony and G. Waring, Rapp. P.-V. Reun. Cons. Int. Explor. Mer. 177, 72 (1980).
31. R. M. Peterman, Can. J. Fish. Aquat. Sci. 46, 2 (1990).
32. We repeated the analyses assuming gammainstead of log-normally distributed residuals; the results were almost identical. Robustness to estimation error in spawners and serial correlation in recruitment were investigated by introduction of these effects into the procedure used to estimate power. Lognormal errors in the estimation of spawners (cr = 0.2) and first-order autocorrelation of 0.4 in recruitment did not increase type 1 errors if 5 was held at 1. As expected, the power was reduced if depensation was present (<5 = 2); for the 26 high-power stocks, with errors in the estimation of spawners, the power was reduced by approximately 3% on average, whereas with autocorrelation, the power was reduced by approximately 1% on average. In addition, we tested the adequacy of the chi-square approximation to the distribution of the likelihood ratio statistic by calculating the type 1 error rate when 6 = 1 (approximately 3% on average for the 26 highpower stocks).
33. F. Neave, J. Fish. Res. Board Can. 9, 450 (1953)
34. J.G. Hunter, J. Fish. Res. Board Can. 16, 835 (1959)
35. O. Ulltang, Rapp. P.-V. Reun. Cons. Int. Explor. Mer. 177, 489 (1980).
36. R. A. Myers and N. G. Cadigan, Can. J. Fish. Aquat. Sci. 50, 8 (1993).
37. We thank the many fish population biologists who generously provided their data and the Canadian Department of Fisheries and Oceans Northern Cod Science Program for financial assistance.
38. J. Klicka and R. M. Zink, Science 277, 1666 (1997).
39. M. L. Rosenzweig, Species Diversity in Space and Time (Cambridge Univ. Press, Cambridge, 1995).
40. M. L. Rosenzweig and R. McCord, Paleobiology 17, 202 (1991).
41. P. A. Johnson, F. C. Hoppensteadt, J. J. Smith, G. L. Bush, Evol. Ecol. 10, 187 (1996).
42. M. L. Rosenzweig, in Biodiversity Dynamics; Turnover of Populations, Taxa and Communities, M. L. McKinney, Ed. (Columbia Univ. Press, New York, in press).
43. R. A. Martin and K. B. Fairbanks, Evol Ecol, in press.
44. J. K. McKee, J. Theor. Biol 172, 141 (1992).
45. E. S. Vrba, in Macroevolutionary Patterns in the Plio-Pleistocene of Africa, T. G. Brommage and F. Shrenk, Eds. (Oxford Univ. Press, Oxford, in press).
46. G. H. Denton, Macroevolutionary Patterns in the Plio-Pleistocene of Africa.
47. E. S. Vrba, in The Evolutionary History of Robust Australopithecines, F. E. Grine, Ed. (Al- dine Publishing, New York, 1988), pp. 405-426.
48. B. Wood and M. Collard, in [8].
49. J. J. Sepkoski, J. Geol. Soc. London 146, 7 (1989)
50. D. Jablonski, Science 253, 754 (1991)
51. G. J. Retallack, Science 267, 77 (1995).
52. Data for the Figure came from the following sources: Ages and durations from [W. Harland et al, A Geologic Time Scale, 1989 (Cambridge Univ. Press, Cambridge, 1990)] with the help of J. J. Sepkoski Jr. to adjust the beginning times of the Cambrian and Ordovician periods.
53. Rock areas and unadjusted diversities came from [D. M. Raup, Paleobiology 279 (1976a)
54. J. J. Sepkoski, Acta Palaeon- tol. Pol. 38, 175 (1993)].
55. Raup counted 71, 112 species entries in the zoological record from 1900 to 1970. He extrapolated 73, 139 more from these 71, 112. of these 144, 251 species, 7416 were insects, or uncertain, or lived in more than one period, or lived in the Pre-Cambrian; I excluded the latter. The remaining 136, 835 species represent a substantial fraction of all known fossil forms. The adjusted number of species in a period is the number (of the 136, 835) attributed to it divided by A0383, where A is the period's rock area divided by that of the Cenozoic.
56. P. W. Bretsky and S. W. Bretsky, Lethaia 9, 223 (1976)
57. W. D. Allmon, G. Rosenberg, R. W. Portell, K. S. Schindler, Science 260, 1626 (1993).
58. First, arrange the data in order of declining accumulation rates, so that the “first” period has the highest rate and the “last” has the lowest. This imposes the largest possible negative second derivative on the data. Second, reverse the order (so that the one with the lowest rate is first). This imposes the largest possible positive second derivative on the data. The results are, in the first case, an exponent of 0.866 in a power fit, and in the second, an exponent of 1.204. The Figure 1 shows the polynomial approximations of these curves.
59. Title respectfully purloined from my late colleague, G. G. Simpson, W. DiMichele, D. Jablonski, K. Flessa, M. McKinney, A. Miller, K. Niklas, S. Pimm, J. Sepkoski, and D. A. Thomson for stimulating ideas, advice, and feedback.
60. M. D. Hunger and P. W. Price, Ecology 73, 724 (1992).
61. N. G. Hairston, F. E. Smith, L. B. Slobodkin, Am. Nat. 94, 421 (1960).
62. S. L. Pimm, The Balance of Nature (Univ. of Chicago Press, Chicago, IL, 1992), chap. 14.
63. We considered the model in G. Caughley, Analysis of Vertebrate Populations (Wiley, New York, 1977), pp. 130-132.
64. A. T. Bergerud, TrendsEcol. Evol. 3, 68 (1988)
65. A. T. Bergerud, Can. J. Zool. 64, 1515 (1986).
66. R. O. Peterson, R. E. Page, K. M. Dodge, Science 224, 1350 (1984).
67. R. O. Peterson and R. E. Page, J. Mammal. 69, 89 (1988)
68. R. O. Peterson and R. E. Page, Acta Zool. Fenn. 174, 251 (1983).
69. K. L. Risenhoover, thesis, Michigan Technological University (1987), p. 49.
70. R. A. Janke, D. McKaig, R. Raymond, For. Sci. 24, 115 (1978)
71. B.E. McLaren, unpublished data (figures are for the west end of Isle Roy ale).
72. J. D. Snyder and R. A. Janke, Am. Midi. Nat. 95, 79 (1976).
73. Suppression of balsam fir growth on Isle Royale by moose foraging occurs because of continuous removal of upper branches in winter. Trees remain alive in the forest understory for several decades, because lateral branches are protected from browsing by snow cover and because balsam fir is tolerant of shade. Trees in this study (except those represented in Fig. 4A) ranged from 34 to 68 years of age but came from a consistent height interval of 1 to 3 m at the time of sampling. The age variation made no contribution to variation in growth rates.
74. T. A. Brandner, R. O. Peterson, K. L. Risenhoover, Ecology 71, 155 (1991).
75. G. Caughley and A. Gunn, Oikos 67, 47 (1993)
76. S.L. Pimm and A. Redfearn, Nature 334, 613 (1984)
77. M.H. Williamson, The Analysis of Biological Populations (Arnold, London, 1982), chap. 5.
78. A harmonic function fitted with an ordinary least squares approach, period 17.8 years, matched the west-end chronology with R2 = 0.38; a function with period 16.1 years matched the east-end chronology with R2 = 0.70. Following M. G. Bulmer Anim. Ecol. 43, 701 (1978)], functions incorporating both a harmonic term and a first-order autoregressive term were matched to the series, with improved fits of R2^= 0.68 and R2 = 0.81, respectively.
79. The cycle period inferred from empirical consideration of the three-trophic-level system is shorter than a period cited earlier, which was derived from a body-mass regression [6]. However, the confidence interval for the allometric Equation also included the shorter period of 16 to 18 years.
80. Ring-width suppression followed heightgrowth suppression in both samples by about 3 years.
81. G. Caughley, Oecologia 54, 309 (1976). The plant with the inferior growth rate lags behind herbivore density more than does the plant with the superior growth rate.
82. B. E. McLaren, unpublished data; R. M. Linn, thesis, Duke University (1957).
83. R. O. Peterson, unpublished data.
84. Our result is new primarily because studies of long-term processes in terrestrial food chains are rare, especially studies of those in which top carnivores still exist. A. R. E. Sinclair and colleagues [Am. Nat. 141, 173 (1993)] recently described the well-known snowshoe hare (Lepus amehcana) cycle in the Yukon Territory, Canada, as generating cyclic suppression of white spruce (Picea glauca). These authors, however, favored an explanation in meteorological driving forces. Other studies [for example, L. B. Keith, A. W. Todd, C. J. Brand, R. S. Adamcik, D. H. Rusch, Int. Congr. Game Biol. 13, 151 (1976)
85. L.B. Keith and L. A. Windberg, Midi. Monogr. 58, 70 (1978)] have considered plant-herbivore and herbivore-carnivore relations for the snowshoe hare, but in different study areas. The third trophic level in the fox-prey system studied by E. R. Lindstroem and colleagues [Ecology 75, 1042 (1994)] was a parasite.
86. J. M. Thurber and R. O. Peterson, J. Mammal. 74, 879 (1993).
87. D. J. McQueen, M. R. S. Johannes, J. R. Post, T. J. Stewart, D. R. S. Lean, Ecol. Monogr. 59, 289 (1989)
88. S.R. Carpenter et al., Ecology 68, 1863 (1987)
89. M. Lynch and J. Shapiro, Limnol Oceanogr. 26, 86 (1981)
90. J.L. Brooks and S. I. Dodson, Science 150, 28 (1965).
91. L. Oksanen, S. D. Fretwell, J. Arruda, P. Niemela, Am. Nat. 118, 240 (1981)
92. for other contrasting published ideas concerning predator-prey cycles on Isle Royale, see [6] and J. Pastor and R. J. Naiman, Am. Nat. 139, 690(1992).
93. M. L. Rosenzweig, Am. Nat. 102, 67 (1968). The index was calculated with the use of daily records from Thunder Bay, Ontario.
94. Understory fir density exceeds 250 ha"1 over about 70% of Isle Royale's land area and exceeds 10, 000 ha"1 in many shoreline forests. Values were obtained by averaging density estimates for stems < 2 m in 0.01-ha plots at the corners of each (square-mile) township section.
95. The ring-width index is the ratio of increments calculated as a summation of volume differences for a series of stacked conic sections representing the stems of sampled trees, divided by the cambial surface area at the beginning of each growing season for each year. An aggregate index, comprising ring-width measurements (±10~2 mm) averaged for four radii at 5-to 10-cm increments along stems, is presented. This intensity of sampling permits an accurate height-growth reconstruction of the trees and measurement of the wood volume increment throughout the trees' stems.
96. We acknowledge the National Park Service (NPS) for permission to sample trees; V. G. Smith, University of Toronto, for the design of the ring-width index; R. J. Miller, Ontario Ministry of Natural Resources, for the development of analytic software and loan of measurement equipment; Environment Canada for supplying weather records; and the following for financial support: National Geographic Society, NPS Earthwatch, the Boone and Crockett Club, and NSF grant DEB-9317401. W. C. Kerfoot and F. H. Wagner kindly reviewed this manuscript and offered many helpful suggestions.
97. R. Brookmeyer, Science 253, 37 (1991).
98. Centers for Disease Control, Morb. Mortal. Wkly. Rep. 43 826 (1994). This plateau was obscured after January 1993 by the expansion of case definition to include HIV-infected persons with severe CD4 depletion [Centers for Disease Control, 41 (no. RR-17) (1992)].
99. M. H. Gail, P. S. Rosenberg, J. J. Goedert, J. AIDS 3, 296 (1990).
100. N. G. Becker and I. C. Marschner, Biometrika 80, 165 (1993)
101. A. Verdecchia et ah, Int. J. Epidemiol. 23, 1027 (1994).
102. P. S. Rosenberg, Stat. Med. 13, 1975 (1994).
103. Centers for Disease Control, Morb. Mortal. Wkly. Rep. 36 (suppl. IS) (1987).
104. The analysis was based on all cases in the United States and its territories, including Puerto Rico, who were older than 13 years at diagnosis and who were reported to CDC as of March 1994. The AIDS incidence and mortality data for individuals who met the 1987 case definition [6] were adjusted for reporting delays with standard CDC adjustment weights. Incidence counts were inflated by 18%, an approximate figure, to account for unreported cases and for infected individuals who died before meeting the case definition. Adjusted totals of 336, 000 incident cases of AIDS and 227, 000 AIDS deaths as of 1 January 1993 were stratified by month of diagnosis, single year of age at diagnosis, gender, and race and ethnicity (non-Hispanic whites, non-Hispanic blacks, Hispanics, and others). Use of the 1987 definition maintains consistency between the AIDS data and available estimates of the incubation distribution.
105. M. E. Eyster et al, Ann. Intern. Med. 107, 1 (1987)
106. A.R. Moss et al, Br. Med. J. 296, 745 (1988)
107. A.B. Mariotto et al, Am. J. Epidemiol. 135, 428 (1992)
108. P.S. Rosenberg, J. J. Goedert, R. J. Biggar, AIDS 8, 803 (1994).
109. Technical details of the fitting procedures are described elsewhere [5].
110. This fast model may yield conservative estimates because the natural history hazard function may level off around 7 years after infection [J. C. M. Hendriks et al., AIDS 7, 231 (1993)]. The Weibull hazard, in contrast, continues to increase. To reflect that progression rates appear to be slower with infection at younger ages, the age-specific hazard of progression for the fast model was assumed to increase (decrease) by the factor 1.042 (0.960) per year increase (decrease) in the age at infection compared to a person aged 30 years at infection, as derived in [9].
111. The estimated linear slope effect was 0.01 per year at 2.6 years after infection and the corresponding age effect was 1.037 [9]. Both the fast and slow estimates of the incubation distribution were derived from a large cohort study of HIV-infected homosexual men [R. J. Biggar and the International Registry of Seroconverters, AIDS 4, 1059 (1990)].
112. Therapy was assumed to reduce the hazard of AIDS among treated individuals by the factor 0.50 for persons infected for about 5 years, and then to wear off [15]. The estimated cumulative proportion in treatment before AIDS was 40% for white males and 20% for other demographic groups. Therapy was assumed to have been introduced in April 1987 among white men and women and in April 1988 among other groups. These parameters are broadly consistent with clinical trial results [Concorde Coordinating Committee, Lancet 343, 871 (1994)] and with estimates of the extent of therapy use in different groups [N. M. H. Graham et al, J-AIDS 4, 267 (1991)
113. P.S. Rosenberg et al, J-AIDS 392
114. W. Lang et al, J-AIDS. 6, 191 (1993)].
115. To estimate age-specific prevalence at time T, note that individuals aged k years represent the surviving members of the cohort born at time T-k. Prevalence is estimated as the cumulative probability that a member of this cohort becomes infected as of T minus the cumulative probability that a member dies of AIDS as of r. These quantities can be estimated with actuarial formulas applied to census and AIDS mortality data and the results from the backcalculations.
116. Centers for Disease Control, Morb. Mortal. Wkly.Rep. 39, 110(1990).
117. P. S. Rosenberg, M. H. Gail, R. J. Carroll, Stat. Med. 11, 1633 (1992).
118. G. M. McQuillan et al, J. AIDS 7, 1195 (1994).
119. M. Gwinn et al, J. Am. Med. Assoc. 265, 1704 (1991).
120. J. M. Karon, personal communication. This estimate is based on an analysis of unpublished CDC data.
121. W. S. Cleveland, J. Am. Stat. Assoc. 74, 829 (1979).
122. I gratefully acknowledge M. Morgan for providing national AIDS surveillance data, T. Green for calculating the reporting delays, J. Smith of Information Management Systems for data management, and J. M. Karon and R. J. Biggar for many helpful discussions.
123. Wuebbles, D. J. J. geophys. Res. 88, 1433-1443 (1983).
124. Fisher, D. A. et al, Nature 344, 508-512 (1990).
125. Solomon, S. et al, J. geophys. Res. 97, 825-842 (1992).
126. Solomon, S. & Albritton, D. L. Nature 357, 33-37(1992).
127. Scientific Assessment of Ozone Depletion: 1991 Ch. 6 Global Ozone Research and Monitoring Project Rap. No. 25 (World Meteorological Organization, Geneva, 1992).
128. Scientific Assessment of Ozone Depletion: 1991 Ch. 10 Global Ozone Research and Monitoring Project Rep. No. 25 (World Meteorological Organization, Geneva, 1992).
129. Scientific Assessment of Ozone Depletion: 1991 Ch. 9 Global Ozone Research and Monitoring Project Rep. No. 25 (World Meteorological Organization, Geneva, 1992).
130. Methyl Bromide: Its Atmospheric Science, Technology and Economics Montreal Protocol Assessment Supplement (United Nations Environment Programme, Nairobi, 1992).
131. Prather, M. J. & Watson R. T. Nature 244, 729-734 (1990).
132. Scientific Assessment of Ozone Depletion: 1991 Ch. 8 Global Ozone Research and Monitoring Project Rep. No. 25 (World Meteorological Organization, Geneva, 1992).
133. Holton, J. R. J. atmos. Sci. 47, 392-395 (1990).
134. Ko, M.K.W. et al. J. geophys. Res. 96, 7547-7552(1991).
135. Chem. Engng. News p. 14(16 August 1993).
136. Prather, M. J. et al. J. geophys. Res. 95, 18583-18590(1990).
137. Science 260, 1884-1885 (1993).
138. Chem. Engng. News p. 18 (17 January 1994).
139. J. T. Wootton, M. E. Power, Proc. Natl. Acad. Sci. U.S.A. 90, 1384(1992).
140. M. L. Rosenzweig, Science 171, 385 (1971)
141. L. Oksanen, S. D. Fretwell, J. Arruda, P. Niemela, Am. Nat. 118, 240 (1981)
142. S.L. Pimm, Food Webs (Chapman and Hall, London, 1982)
143. D. Tilman, Am. Nat. 129, 769 (1987)
144. S.L. Pimm, R. L. Kitching, Oikos 50, 302 (1987)
145. P.A. Abrams, Am. Nat. 141, 351 (1993).
146. M. E. Power, Science 250, 811 (1990).
147. R. T. Paine, Am. Nat. 100, 65 (1966)
148. J.H. Connell, Science 199, 1302 (1978)
149. D.L. DeAngelis, Dynamics of Nutrient Cycling and Food Webs (Chapman and Hall, London, 1992)
150. J.T. Wootton, Anna. Rev. Ecol. Syst. 25, 443 (1994).
151. M. E. Power, M. S. Parker, J. T. Wootton, in Food Webs: Integration of Patterns and Dynamics, G. A. Polis and K. O. Winemiller, Eds. (Chapman and Hall, New York, 1995), pp. 286-297.
152. M. E. Power, Arch. Hydrobiol. 125, 385(1992).
153. Dicosmoecus exposed to seven predator species in 12.6-liter tubs (six replicates per predator species; 10 grazers and one predator per tub) experienced no mortality over 24 hours, whereas a representative grazer, the mayfly Centroptilum, experienced 50.7 ± 24.5% (SD) mortality due to predation pooled across all replicates of all predator species.
154. J. H. Connell, Am. Nat. 122, 661 (1983)
155. T. W. Schoener, Am. Nat., p. 240.
156. We visually counted all fish (Hesperoleucas symmetricus and Gasterosteus aculeatus) and odonates (primarily Aeshna californica and Archilestes californica) in the entire channel and counted all grazing insects [predominantly Centroptilum, Nixe, and Paroleptophlebia (Ephemeroptera); Mysticides, Lepidosoma, Tinodes, and Gumaga (Trichoptera); Ferrissia and Physella (Pulmonata); and Eubrianix (Psephenidae)], except for midge (Chironomidae) larvae, on the top and bottom of each tile. Midges were counted on only one tile because of their high abundance. We collected benthic algae (largely a mixture of Cladophora glomerata, Nostoc spp., and Epithemia spp.) from one randomly selected tile in each treatment by scraping the tile with a razor blade and measuring ashfree dry weight in the laboratory. All floating algal mats in a channel were collected with an aquarium net, spun in a salad spinner for a standard 50 turns to remove excess water, and weighed. Subsamples of floating algae were collected, dried, weighed, and analyzed through measurement of ash-free dry weight to calibrate wet mass with ashfree dry weight. We used a blocked, two-way MANOVA to test for community-wide differences among treatments, and then tested for the specific differences predicted a priori by the model, using one-tailed paired t tests.
157. L. B. Crowder and W. E. Cooper, Ecology 63, 1802 (1982)
158. M.E. Power, Oikos 58, 67 (1990).
159. J. T. Wootton, Am. Nat. 141, 71 (1993).
160. Methods are described in [6]. Data were taken from survey locations described in [6] and from surveys of reaches of the Mad River (regulated) and Van Duzen River (unregulated) in Six Rivers National Forest, CA, conducted throughout the summer of 1994. All variables except algal occurrence were logtransformed before analysis to stabilize the variance.
161. M. Waldichuk, Can. Bull. Fish. Aquat. Sci. 226, 295 (1993)
162. A.G. Maule, C. B. Schreck, C. S. Bradford, B. A. Barton, Trans. Am. Fish. Soc. 117, 245 (1988)
163. H.L. Raymond, Am. J. Fish. Manage. 8, 1 (1988).
164. R. W. Nelson, J. R. Dwyer, W. E. Greenberg, Environ. Manage. 11, 479 (1987).
165. F. K. Ligon, W. E. Dietrich, W. J. Trust. Bioscience 45, 183 (1995)
166. M.E. Power, A. Sun, G. Parker, W. E. Dietrich, J. T. Wootton, Bioscience., p. 159
167. C. S. Holling and G. K. Meffe, Conserv. Biol. 10, 328 (1996).
168. We thank B. Amerson, C. Bailey, J. Chase, M. Eskridge, D. Gordon, N. Guthrie, S. Kupferburg, S. Lane, M. Liu, J. Lyons, J. Marks, S. McGuire, K. Meier, E. Noonberg, C. Pfister, M. Pizer, W. Roberts, M. Salzer, A. Sun, C. Wang, and J. Wootton for field assistance and P. Steel for logistical support. Funded in part by NSF, the California State Water Resources Center, the Miller Institute for Basic Research, and the University of Chicago Block Fund.
169. C. B. Williams, Nature 152, 264 (1943)
170. F.W. Preston, Ecology 29, 254 (1948)
171. R.P. Mclntosh, The Background of Ecology (Cambridge Univ. Press, Cambridge, 1985).
172. R. H. MacArthur and E. O. Wilson, The Theory of Island Biogeography (Princeton Univ. Press, Princeton, NJ, 1967).
173. F. W. Preston, Ecology 41, 611 (1960)
174. E.F. Connor and E. D. McCoy, Am. Nat. 113, 791 (1979)
175. M.L. Rosenzweig, Species Diversity in Space and Time (Cambridge Univ. Press, Cambridge, 1995).
176. I. Hanski, Oikos 38, 210 (1982)
177. J. Kouki, A. Halkka, in Species Diversity in Ecological Communities: Historical and Geographical Perspectives, R. E. Ricklefs and D. Schluter, Eds. (Univ. of Chicago Press, Chicago, 1993), pp. 108-116
178. K. J. Gaston and J. H. Lawton, Oikos 58, 329 (1990)
179. J.H. Lawton, Trends Ecol. Evol. 8, 409 (1993)
180. K.J. Gaston, Rarity (Chapman and Hall, London, 1994)
181. T. M. Blackburn, J. H. Lawton, in preparation.
182. J. H. Brown, Am. Nat. 124, 255 (1984).
183. I. Hanski, Oikos 62, 88 (1991).
184. M. Gyllenberg and I. Hanski, Theor. Popul. Biol. 42, 35 (1992)
185. I. Hanski and M. Gyllenberg, Am. Nat. 142, 17 (1993).
186. S. Nee, R. D. Gregory, R. M. May, Oikos 62, 83 (1991)
187. L.A. Venier and L. Fahrig, Oikos. 76, 564 (1996).
188. To our knowledge, the only reference to a possible connection between the SA and DA curves is by I. Hanski, J. Kouki, and A. Halkka [in Species Diversity in Ecological Communities: Historical and Geographical Perspectives, R. E. Ricklefs and D. Schluter, Eds. (Univ. of Chicago Press, Chicago, 1993), p. 108.
189. We use the term "island" to refer both to true and habitat islands, and to what are called habitat patches, or fragments, in the metapopulation literature [7].
190. Without any loss of generality, we assume that mw equals unity.
191. M. E. Gilpin and J. M. Diamond, Proc. Natl. Acad. Sci. U.S.A. 73, 4130 (1976)
192. I. Hanski, Trends Ecol. Evol. 9, 131 (1994)
193. J.Anim. Ecol. 63, 151 (1994). In the latter two references, the more general relation m, = e/Kfj is used, but we assume here for simplicity that x = 1, which is a good approximation for many species. The rate parameter e has been here absorbed in the unit of island area, to give the per-year extinction probability 1-e~]/A for species with w = 1. A scaling constant d is given by d = A' I A, where A' is island area in (say) kilometers squared, d may be calculated from knowledge of the per-year extinction probability and A'.
194. T. W. Schoener, in 16th International Ornithological Congress, Canberra, Australia, 12 to 17 August 1974 (Australian Academy of Sciences, Canberra City, 1976), pp. 629-642.
195. Somewhat confusingly, the classical metapopulation scenario, where there is no external mainland, is referred to as the "mainland" regression in the species-area literature.
196. The equilibrium probability of species i occupying island j, p*-, which is called the incidence, Jij, is obtained from Equation 1 as [Equation Found] In the metapopulation model, the incidences can be calculated only numerically. Substituting Equation 2 into the expression C* = cWiXJijAj, which gives the equilibrium value of Ci(t) in the metapopulation model, we obtain[Equation Found] from which C* can be solved provided that cwf X A2-> 1, which is a necessary and sufficient condition for species i to persist in the network of islands. The incidences can then be calculated from Equation 2.
197. O. Arrhenius, J. Ecol. 9, 95 (1921)
198. H.A. Gleason, Ecology 3, 158 (1922)
199. F.W. Preston, Ecology. 41, 611 (1960).
200. M. H. Buys, J. S. Maritz, J. J. A. Van Der Walt, J. Veg. Sci. 5, 63 (1994)
201. M.R. Williams, Ecology 76, 2607(1995).
202. The nonlinear logistic model can be linearized with the logit-transformation, log[P/(l-P)] = a + blogiv, which we apply throughout this report.
203. Assuming that Q (species number in the pool) is large and that log if is uniformly distributed with zero mean, we obtain after some calculation the expected number of species on island j as[Equation Found] 20.
204. Denoting the per-year colonization and extinction probabilities by A and /u, we obtain[Equation Found] where E denotes the expected value. That the slope of the SA curve depends on the ratio of colonization to extinction probabilities has been suggested by R. E. Ricklefs and G. W. Cox [Am. Nat. 106, 195 (1972)]
205. M. P. Johnson and D. S. Simberloff [Biogeogr. 1, 149 (1974)]. [Equation Found] 23. We derived the slope values [19] also for backto-back exponential and lognormal distributions of w. The expressions are more complicated in these cases, but the results are very similar.
206. F. W. Preston, Ecology 43, 185 (1962)
207. R.M. May, in Ecology and Evolution of Communities, M. L. Cody and J. M. Diamond, Eds. (Belknap, Cambridge, MA, 1975), p. 81.
208. C. Wissel and B. Maier, J. Biogeogr. 19, 355 (1992).
209. B. D. Coleman, Math. Biosci. 54, 191 (1981)
210. H. Caswell and J. E. Cohen, in Species Diversity in Ecological Communities: Historical and Geographical Perspectives, R. E. Ricklefs and D. Schluter, Eds. (Univ. of Chicago Press, Chicago, 1993), p. 99.
211. B. D. Patterson, Conserv. Biol. 1, 323 (1987).
212. M. Nieminen, Oecologia 108, 643 (1996). The moth example includes 46 islands located about 15 km off the mainland in southwest Finland. Moths were trapped with sugar-bait traps in summer 1993. The parameter log w gives the abundance on mainland.
213. Z. Witkowski and P. Plonka, Bull. Pol. Acad. Sci. Biol. Sci. 32, 241 (1984)
214. W.B. Harms and P. Opdam, in Changing Landscapes: An Ecological Perspective, I. G. Zonneveld and R. T. T. Forman, Eds. (Springer-Verlag, New York, 1990), p.
215. The bird examples are from studies in Poland (Witkowski and Plonka) and The Netherlands (Harms and Opdam). "Islands" in this instance are habitat islands on mainland. The parameter logw in Figure 1H gives the estimated density of territories.
216. E. F. Connor and E. C. McCoy, Am. Nat. 113, 791 (1979).
217. We thank F. Adler, J. Brown, M. Camara, E. Connor, M. Kuussaari, J. Lawton, A. Moilanen, and T. Schoener for comments on the manuscript. W. W. Murdoch and National Center for Ecological Analysis and Synthesis (NCEAS) in Santa Barbara are thanked for supporting this research.
218. M. P. Hassell, R. M. May, Spatial heterogeneity and the dynamics of parasitoid-host systems, Ann. Zoo. Fennisi 25, (1998), pp. 55-61.
219. P. Kareiva, Population dynamics in spatially complex environments: theory and data, Phil. Trans. R. Soc. B 330 (1990), pp. 175-190.
220. R. M. May, in Large-scale Ecology and Conservssaon Biology (eds Edwards, P. I., May, R. M. and Webb. N. R.) (Blackwell Scientific, Oxford, 1994), pp. 1-17.
221. P. Kareiva, U. Wennergren, Connecting landscape pattern to ecosystem and population processes, Nature 373, (1995), pp. 299-302.
222. I. K. Waage, Aggregation in field parasitoid populations: foraging time allocation by a population of Diadegma (Hymenoptera, Ichneumonidae) Ecol Entomol. 8, (1963), pp. 447-453.
223. A. Kruess & T. Tscharntke, Habitat fragmentation, species loss, and biological control, Science 264, (1994), pp. 1581-1584.
224. P. A. Heads & J. H. Lawton, Studies on the natural enemy complex of the holly leaf-miner: the effects of scale on the detection of aggregative responses and the implications for biological control, Oikos 40, (1983), pp. 267-276.
225. L. D. Rothman, D. C. Darling, Spatial density dependence: effects of scale, host spatial pattern and parasitoid reproductive strategy, Oikos 62, (1991), pp. 221-230.
226. P. Kareiva, Habitat fragmentation and the stability of predator-prey interactions, Nature 326, (1987), pp. 388-390.
227. D. J. Rogers & B. G. Williams, in Large-scale Ecology and Conservation Biology (eds Edwards, P. J., May, R. M. and Webb. N. R.), (Blackwell Scientific, Oxtord, 1994), pp. 247-271.
228. P. C. Marino, D. A. Landis, Effects at landscape structure on parasitoid diversity and parasitism in agroecosystems, Ecol. Appl. 6, (1996), pp. 276-284.
229. W. L. Sippell, Outbreaks of the forest tent caterpillar, Maiacosoma disstria Hbn., periodic defoliator of broad-leaved trees in Ontario, Can. Entomol. 94, (1962), pp. 408-416.
230. J. A. Witter, H. M. Kulman, The parasite complex of the forest tent caterpillar in northern Minnesota, Env. Entomol, 8, (1979), pp. 723-731.
231. D. Parry, Larval and pupal parasitism of the forest-tent caterpillar Malacosoma disstria Hilbner (Lepidoptera, Lasiocampidae) in Alberta, Canada, Can. Entomol. 127, (1995), pp. 877-893.
232. H. C. J. Godfray, Parasitoids: Behavioral and Evolutionary Ecology (Princeton Univ. Press, Princeton, 1994).
233. J. Roland, Large-scale forest fragmentation increase the duration of forest tent caterpillar outbreak, Oecologia 93, (1993), pp. 25-301.
234. J. Roland, P. D. Taylor, in Population Dynamics: New Approaches and Synthesis (eds Cappuccino, N. and Price, P. W.) (Academic, San Diego, 1995, pp. 195-208).
235. R. M. Weseloh, Spatial distribution of gypsy moth (Lepidoptera: Lymantriidae) and some of its parisitoids within a forest environment, Entomophaga 17, (1972), pp. 339-351.
236. C. S. Holling, Cross-scale morphology, ge ometry and dynamics of ecosystems, Ecol. Monogr., (1992), pp. 447-502.
237. P. McCullagh, G. A. Nelder, Generalized Linear Models (Chapman and Hall, London, 1989).
238. L. Wilkinson, M. Hill, S. Miceli, G. Birkenbeuel, E. Vang, SYSTAT for Windows: V.5 (Evanston, Illinois, 1992).
239. H. E. Kubitschek, Nature 182, 234 (1958).
240. T. Allen, Particle Size Analysis 110-127 (Society for Anlytical Chemistry Publishers, London, 1967).
241. L. G. Bunville, Modern Methods of Particle Size Analysis (ed. H. G. Barth) 1-42 (Wiley, New York, 1984).
242. R. W. DeBlois et al., J. Colloid Interface Sci. 61, 323 (1977).
243. S. M. Bezrukov and I. Vodyanoy, Biophys. J. 64, 16 (1993).
244. M. S. P. Sansom, Eur. Biophys. J. 22, 105 (1993).
245. J. E. Hall et al, Biophys. J. 45, 233 (1984).
246. I. Vodyanoy et al, Biophys. J. 65, 2097 (1993). 9.
247. V. Krasilnikov et al, FEMS Microbiol. Immunol 105, 93 (1992).
248. S. M. Bezrukov and I. Vodyanoy, Membrane Electrochemistry (Advances in Chemistry Ser. No. 235, American Chemical Soc, Washington DC, in press).
249. G. Feher and M. Weissman, Proc. Nat. Acad. Sci. 70, 870 (1973).
250. S. J. Kuga, J. Chromatography 206, 449 (1981).
251. C. P. Bean, Membranes (ed. G. Eisenman) (Dekker, New York, 1972).
252. A. Couper and R. F. T. Stepto, Trans. Faraday Soc. 65, 2486 (1969).
253. P. Hess and R. W. Tsien, Nature 309, 453 (1984).
254. S. H. Heinemann and F. Sigworth, J. Biochim. Biophys. Acta 987, 8 (1989).
255. S. M. Simon and G. Blobel, Cell 65, 371 (1991).
256. M. Montal and P. Mueller, Proc. Nat. Acad. Sci. 69, 3561 (1972).
257. T. M. Balasubramanian et al, J. Am Chem. Soc. 103, 6127(1981).
258. L. G. M. Gordon and D. A. Haydon, Phil. Trans. R. Soc. B 270, 433 (1975).
259. Greenland Ice-Core Project (GRIP) members, Nature 364, 203 (1993).
260. W. Dansgaard et al, Nature., p. 218.
261. N. S. Gundestup, D. Dahl-Jensen, S. J. Johnsen, A. Rossi, Cold Reg. Set Technol. 21, 399 (1993).
262. G. D. Clow, R. W. Saltus, E. D. Waddington, J. Glaciol. 42, 576 (1996).
263. The deep borehole is located in a building, and the liquid surface in the borehole is found at a depth of 40 m. The temperatures measured in the top 40 m are very disturbed, so we used measurements from an air-filled shallow borehole (100 m) near the borehole.
264. N. S. Gundestrup and B. L. Hansen, J. Glaciol. 30, 282 (1984).
265. D. Dahl-Jensen and S. J. Johnsen, Nature 320, 250(1986).
266. D. Dahl-Jensen et al, I Glaciol 43, 300 (1997).
267. Between 50 and 20 ka, the ice thickness was 50 m less than at present, even though the ice sheet covered a larger area. The maximum ice thickness of 3230 m is found at 10 ka, after which the ice thickness gradually has decreased to the present 3028.6 m. The depression and uplift of the bedrock influences the elevation of the surface [S. J. Johnsen, D. Dahl-Jensen, W. Dansgaard, N. S. Gundestrup, Tellus 5 47, 624(1995)].
268. K. M. Cuffey and G. D. Clow, J. Geophys. Res. 102, 26383 (1997).
269. The past accumulation rates are determined by coupling them to the past (unknown) temperature through the relation A(T) = Ao exp[0.0467(T-To)-0.000227(7-7b)2], where A(T) is the accumulation rate at the surface temperature T, Ao is the present ice accumulation rate, which is 0.23 m/year at GRIP and 0.49 m/year at Dye 3, and To is the present surface temperatures at the sites:-31.7°C at GRIP and-20.1°C at Dye 3, respectively [9].
270. S. J. Johnsen, IAHS-AISH Publ. 118, 388 (1977).
271. S. J. Johnsen and W. Dansgaard, NATO ASI Ser. I Global Environ. Change 2, 13 (1992).
272. K. Mosegaard and A. Tarantola, J. Geophys. Res. 100, 12431 (1995).
273. K. Mosegaard, Inverse Problems 14, 405 (1998).
274. Our Monte Carlo scheme is a random walk in the high dimensional space of all possible models, m (temperature histories and geothermal heat flow densities). The temperature history has been divided in 125 intervals (interval length is 25 ky at 450 ka and 10 years at present). Including the geothermal heat flow density as an unknown the model space is 126-dimensional. In each step of the random walk, a perturbed model, wipert of the current model vector m1 is proposed. The next model becomes equal to mpert with an acceptance probability ^accept = min{l, exp(-[5(TOpert)-Sim1])}, where S(m) = Zj(gHm)-dJ obs)2, which is the misfit function measuring the difference between g{m), the calculated borehole temperatures, and dObS, the observed temperatures. If the perturbated model is rejected, the next model becomes equal to ml and a new perturbed model is proposed. To ensure an efficient sampling of all possible models, we developed ways of choosing the temperature histories and geothermal heat flow densities to be tested. The main scheme to perturb the models is to randomly select one of the 126 temperature/heat flow density parameters and change its value to a new value chosen uniformly at random within a given interval. A singular value decomposition (SVD) of the matrix G = {dgj/dnii}, evaluated in a near-optimal model, yields a set of eigenvectors in the model space whose orientations reveal efficient directions of perturbation for the random walk. The SVD method is included as a possible method of perturbing models especially in the start of the process as it speeds the Monte Carlo scheme significantly.
275. Of the 3.3 x 10(i models tested during the random walk 30% have been accepted by the Monte Carlo scheme [16]. Every 500 is chosen of those where the misfit function S [16] is less than the variance of the observations. The waiting time of 500 has been chosen to exceed the maximum correlation length of the output model parameters. This is a necessary condition for the 2000 models to be uncorrelated. To further ensure that the output models were uncorrelated, the random walk was frequently restarted at several random selected points in the model space.
276. The probabilistic formulation of the inverse problem leads to definition of a probability distribution in the model space, describing the likelihood of possible temperature histories and geothermal heat flow densities. The Monte Carlo scheme is constructed to sample according to this probability distribution. The histograms in Figure 2 describe the probability distribution of the geothermal heat flow density and temperatures at times before present. The maxima in the histograms thus describe the most likely values. The method does not constrain the distributions to have a single maximum, indeed there could be histograms with several maxima, reflecting that more than one value of the temperature at this time would give a good fit to the observed temperature in the borehole. The histograms however, are all seen to have a well-defined zone with most likely past temperatures. A soft curve is fitted to the histograms and the maximum value is taken as the most likely value. The standard deviations shown in Figure 3 are derived as deviations from the maximum value.
277. J. H. Sass, B. L. Nielsen, H. A. Wollenberg, R. J. Munroe, J. Geophys. Res. 77, 6435 (1972).
278. C. Clauser et al., J. Geophys. Res, 102, 18417 (1997)
279. L. Guillou-Frottier, J.-C. Marescal, J. Musset, J. Geophys. Res, 103, 7385 (1998)
280. H.N. Pollack, S. J. Hurter, J. R. Johnson, Rev. Geophys. 31, 267 (1993)
281. W.G. Powell, D. S. Chapman, N. Balling, A. E. Beck, in Handbook of Terrestrial Heat-Flow Density Determination (Kluwer Academic, New York, 1988), pp. 167-222.
282. In order to produce a past temperature record from the calculated past surface temperatures, the temperatures have been corrected to the present elevation of the GRIP site (and Dye 3 site respectively) using the surface elevation changes described in [9] and a lapse rate of 0.006 K/m.
283. K. M. Cuffey et al, Science 270, 455 (1995).
284. D. Dahl-Jensen, in Proceedings of the Interdisciplinary Inversion Workshop 2, Copenhagen, 19 May 1993, K. Mosegaard, Ed. (The Niels Bohr Institute for Astronomy, Physics and Geophysics, University of Copenhagen, Copenhagen, 1993), pp. 11-14.
285. C. U. Hammer et al, Report on the stratigraphic dating of the GRIP Ice Core. Special Report of the Geophysical Department (Niels Bohrs Institute for Astronomy, Physics and Geophysics, University of Copenhagen, Copenhagen, in press).
286. J. Firestone, J. Glaciol. 41, 39 (1995).
287. The amplitude of the warming at 1930 A.D. must be considered to be more uncertain. The information leading to this result are the measured temperatures in an open shallow borehole, where air movements could influence the measurements.
288. D. Fisher et al., NATO ASI Ser. I Global Environ. Change 41, 297 (1996)
289. J.W. C. White et al, J. Geophys. Res. 102, 26425 (1997)
290. J.W. Hurrel, Science 269, 676 (1995)
291. P. Frich et al, in DMI Scientific Report 96-1 (Danish Meteorological Institute, Copenhagen, 1996).
292. R. G. Barry and R. J. Charley, Atmosphere, Weather & Climate (Routledge, London, ed. 6, 1992)
293. J. Overpeck et al, Science 278, 1251 (1997)
294. H.H. Lamb, Climate History and the Modern World (Routledge, London, ed. 2, 1995)
295. N.W. T. Brink and A. Weidick, Quat. Res. 4, 429(1974).
296. G. D. Clow, Palaeogeogr. Palaeoclimatol. Palaeoecol. 98, 81 (1992).
297. To comply with this resolution the time steps have been chosen with increasing length back in time. The increasing length of the time steps can be considered as an efficient way of calculating the mean temperatures in the intervals so full available resolution is kept but the calculations are rationalized.
298. In [7], it is argued that parameter combinations of mean glacial temperature, mean glacial accumulation, and geothermal heat flow density can be found that fit the Dye 3 measurements due to the reduced resolution of the climate history reaching further back than 7 ka. A combination with a geothermal heat flow density of 38.7mW/m2 was chosen corresponding to a mean glacial temperature 12 K colder than the present temperatures. If a value of 51 mW/m2 is chosen as that found for our inversion, the mean glacial temperature is 19 K colder than the present, which is well in agreement with the results found for the GRIP reconstruction. Comparison of the Dye 3 temperature history presented in [7] and that presented here shows a general good agreement for the last 7 ky. The history presented in [7] is more intuitive and less detailed, and the history has not been corrected for elevation changes. The ice thickness was assumed constant in this reconstruction. L. K. Barlov, J. C. Rogers, M. C. Serreze, R. C. Barry, J. Geophys. Res. 102, 26333 (1997).
299. R. A. Keen, Occas. Pap. 34 (Institute of Arctic and Alpine Research, University of Colorado, Boulder, 1980).
300. S. Shubert, W. Higgins, C. K. Park, S. Moorthi, M. Svarez, An Atlas of ECMWF Analyses (1980-87). Part II: Second Moment Quantities, NASA Tech. Memorandum 100762 (1990)
301. F. Rex, World Surv. Climatol. 4, 1 (1969).
302. This is a contribution to the Greenland Ice Core Project (GRIP), a European Science Foundation program with eight nations and the European Economic Commission collaborating to drill through the central part of the Greenland Ice Sheet. G.D.C. thanks the USGS Climate History Program and NSF for support. L. Wolpert, Development (suppl.) 3-12 (1989).
303. J. B. A. Green & J. C. Smith, Trends Genet. 7, 245(1991).
304. J. B. A. Green & J. C. Smith, Nature 347, 391 (1990)
305. J. B. A. Green & J. C. Smith, Cell 71, 731 (1992).
306. J. C. Smith, B. M. J. Price, J. B. A. Green, D. Weigel & B. G. Hermann, Cell 67, 79 (1991).
307. J. B. A. Green, J. C. Smith & J. C. Gerhart, Development 120, 2271 (1994).
308. P. A. Wilson & D. A. Melton, Curr. Biol. 4, 676 (1994)
309. J. B. Gurdon, Development 105, 27 (1989).
310. J. B. Gurdon, K. Kao, K. Kato & N.D. Hopwood, Development (suppl.) 137-142 (1992).
311. S. Cascio & J.B. Gurdon, Development 100, 297(1987).
312. F. M. Rosa, Cell 57, 965 (1989).
313. K. W. Y. Cho, B. Blumberg, H. Steinbeisser & E. M. De Robertis, Cell 67, 1111 (1991).
314. P. Lemaire & J. B. Gurdon, Development 120, 1191 (1994).
315. L. Dale & J. M. W. Slack, Development 99, 527 (1987).
316. S. A. Moody, Devi Biol. 122, 300 (1987).
317. J. B. Gurdon, T.J. Mohun, S. Fairman & S. Brennan, Proc. Natn. Acad. Sci. U.S.A. 82, 139 (1985).
318. P. D. Nieuwkoop, Wilhelm Roux Arch. Dev. Biol. 162, 341 (1985).
319. O. Nakamura, H. Takasaki & M. Ishihara, Proc. Jap. Acad. 47, 313 (1971).
320. J. C. Smith, EMBO J. 12, 4463 (1993).
321. J. M. W. Slack, Mech. Dev. 41, 91 (1993).
322. A. Hemmati-Brivanlou & D. A. Melton, Nature 359, 609 (1992).
323. G. Thomsen et al, Cell 63, 485 (1990).
324. R. A. Cornell & D. Kimmelman, Development 120, 453(1994).
325. E. A. Jones & H.R. Woodland, Development 101, 557(1987).
326. M. Asashima et al, Proc. Natn. Acad. Sci. U.S.A. 88, 6511 (1991).
327. J. C. Smith, B. M. J. Price, K. Van Nimmen & D. Huylebroeck, Nature 345, 729 (1990).
328. G. H. Thomsen & D. A. Melton, Cell 74, 433 (1993).
329. L. Dale, G. Matthews & A. Colman, EMBOJ. 12, 4471 (1993).
330. H. Spemann & H. Mangold, Roux's Arch. Entw. Mech. 100, 599 (1924).
331. A. P. McMahon, Trends Genet. 8, 236 (1992).
332. L.W. Burrus, BioEssays 16, 155 (1994).
333. F. Gonzalez, L. Swales, A. Bejsovic, H. Skaer & A. Martinez-Aria, Mech. Dev. 35, 43 (1991).
334. J. P. Couso, S.A. Bishop & A. Martinez-Arias, Development 120, 621 (1994).
335. W. Driever & C. Nusslein-Volhard, Cell 54, 83 (1988).
336. G. Struhl, K. Struhl & P. M. Macdonald, Cell 57, 1259(1989).
337. D. St Johnston & C. Nusslein-Volhard, Cell 68, 201 (1992).
338. E. L. Ferguson & K. V. Anderson, Cell 71, 451 (1992).
339. J. Heemskerk & S. DiNardo, Cell 76, 449 (1994).
340. K. Basler & G. Struhl, Nature 368, 208 (1994).
341. C. Tickle, B. Alberts, L. Wolpert & J. Lee, Nature 296, 564 (1982).
342. Tickle, J. Lee & G. Eichele, Devi Biol. 109, 82 (1985).
343. J. Brockes, Neuron 2, 1285 (1989).
344. C. Neihrs, H. Steinbeisser & E. M. De Robertis, Science 263, 817(1994).
345. J-P. Vincent & P. A. Lawrence, Cell 77, 909 (1994).
346. M. Placzek, M. Tessier-Lavigne, T. M. Jessell & J. Dodd, Development 110, 19 (1990).
347. P. A. Lawrence, Cell 54, 1 (1988).
348. N. D. Hopwood, A. Pluck, J. B. Gurdon & S. M. Dilworth, Development 114, 31 (1992).
349. Krauss, N. et al. Photosystem I at 4-A resolution represents the first structural model of a joint photosynthetic reaction centre and core antenna system. Nature Struct. Biol. 3, 965-973 (1996).
350. Sidler, W. A. in The Molecular Biology of Cyanobacteria (ed. Bryant, D. A.) 139-216 (Kluwer Academic, Dordrecht, 1994).
351. McDermott G. et al. Crystal structure of an integral membrane light-harvesting complex from photosynthetic bacteria. Nature 374, 517-521 (1995).
352. Kuhlbrandt, W., Wang, D. N. & Filiyoshi, Y. Atomic model of plant light-harvesting complex by electron crystallography. Nature 367, 614-621 (1994).
353. Holzenburg, A., Bewley, M. C, Wilson F. H., Nicholson W. V. & Ford, K. C. Threedimensional structure of photosystem II, Nature 363, 470-473 (1993).
354. Boekema, E. J. et al. Supramolecular structure of the photosystem II complex from green plants and cyanobacteria. Proc. Natl. Acad. Sci. USA 92, 175-179 (1995).
355. Thomas, J. & Webb, W. W. in Non-invasive Techniques in Cell Biology 129-152 (Wiley-Liss, New York, 1990).
356. Walsby, A. E., van Rijn, J. & Cohen, Y. The biology of a new gas-vacuolate cyanobacterium, Dactylococcopsis salina sp. nov. in Solar Lake. Proc. R. Soc. Lond. B, 217, 417-447 (1983).
357. Zhang, F., Lee, G. M. & Jacobson, K. Protein lateral mobility as a reflection of membrane microstructure. BioEssays 15, 579-588 (1993).
358. Kubitscheck, U., Wedekind, F. & Peters, R. Lateral diffusion measurement at high spatial resolution by scanning microphotolysis in a confocal microscope. Biophys. J. 67, 948-956 (1994).
359. Drepper, F., Carlberg, I., Andersson, B. & Haehnel, W. Lateral diffusion of an integral membrane protein: Monte Carlo analysis of the migration of phosphorylated lightharvesting complex II in the thylakoid membrane. Biochemistry 32, 11915-11922 (1993).
360. Barber, J. & Andersson, B. Revealing the blueprint of photosynthesis. Nature 370, 31-34 (1994).
361. Holzwarth, A. R. Fluorescence lifetimes in photosynthetic systems. Photochem. Photobiol. 43, 707-725(1986).
362. Bald, D., Kruip, J. & Rogner, M. Supramolecular architecture of cyanobacterial thylakoid membranes: how is the phycobilisome connected with the photosystems? Photosynth. Res. 49, 103-118(1996).
363. Giddings, T. H., Wasmann, C. & Staehelin, L. A. Structure of the thylakoids and envelope membranes of the cyanelles of Cyanophora paradoxa. Plant Physiol. 71, 409-419 (1983).
364. Mustardy, L., Cunningham F. X. & Gantt, E. Photosynthetic membrane topography: quantitative in situ localisation of photosystems I and II. Proc. Natl. Acad. Sci. USA 89, 10021-10025(1992).
365. Adir, N., Shochat, S. & Ohad, I. Lightdependent Dl protein synthesis and translocation is regulated by reaction centre II. J. Biol. Chem. 265, 12563-12368 (1990).
366. Allen, J. F. Protein phosphorylation in regulation of photosynthesis. Biochem. Biophys. Acta 1098, 275-335 (1992).
367. Soitamo, A. J. et al. Overproduction of the Dl:2 protein makes Synechococcus cells more tolerant to photoinhibition of photosystem II. Plant Mol. Biol. 30, 467-478 (1996).
368. Mullineaux, C. W., Bittersmann, E., Allen, J. F. & Holzwarth, A. R. Picosecond timeresolved fluorescence emission spectra indicate decreased excitation energy transfer from the phycobilisome to photosystems II in the cyanobacterium Synexhococcus 6301. Biochim. Biophys. Acta 1015, 231-242 (1990).
369. Mallineaux, C. W. Excitation energy transfer from phycobilisomes to photosystem I in a cyanobacterial mutant lacking photosystem II. Biochim. Biophys. Acta 1184, 71-77 (1994).
370. van der Oord, C. J. R. et al. High-resolution confocal microscopy using synchrotron radiation. J. Microsc. 182, 217-224 (1996).
371. J. M. Davidenko, A. V. Pertsov, R. Salomonsz, W. Baxter and J. Jalife, Nature 355, pp. 349-351 (1992).
372. F. Siegert and C. J. Weijer, CellSci 93, pp. 325-335 (1989).
373. N. A. Gorelova and J. J. Bures, Neurobiol 14, pp. 353-363 (1983).
374. S. Jakubith, H. H. Rotermund, W. Engel, A. von Oertzen and G. Erth, Phys. Rev. Lett 65, pp. 3013-3016(1990).
375. A. T. Winfree, Science 175, pp. 634-636 (1972).
376. R. J. Field and M. Burger (eds) Oscillations and Traveling Waves in Chemical Systems (Wiley, New York, 1985).
377. W. Jahnke, W. E. Skaggs and A. T. Winfree, J. phys. Chem 93, pp. 740-749 (1989).
378. Th. Plesser, S. C. Muller and B. Hess, J. phys. Chem. 94, pp. 7501-7507 (1990).
379. G. S. Skinner and H. L. Swinney, Physica D40, pp. 1-16(1991).
380. L. Kuhnert, Naturwissenschaften 73, pp. 96-97(1986).
381. H.-J. Krug, L. Pohlmann and L Kuhnert, J. phys. Chem. 94, pp. 4862-4866 (1990).
382. K. I. Agladze, V. A. Davydov and A. S. Mikhailov, JETP Lett. 45, pp. 767-770 (1987).
383. L. Kuhnert, K. I. Agladze, and V. I. Krinsky, Nature 337, pp. 244-247 (1989).
384. M. Markus, Nagy Ungverai and B. Hess, Science 257, pp. 225-226 (1992). 19.
385. M. Braune and H. Engel, Chem. Phys. Lett, (in press). 20.
386. R. J. Field and R. M. Noyes, J. phys. Chem. 60, pp. 1877-1884 (1974). 21.
387. W. Jahnke and A. T. Winfree, Int. J. Bifurcation and Chaos 1, pp. 445-466 (1991).
388. T. Yamaguchi, L. Kuhnert, Zs. Nagy-Ungvarai
389. S. C. Miiller and B. Hess, J. phys. Chem. 95, pp. 5831-5837(1991).
390. O. Steinbock and S. C. Miiller, Physica A188, pp. 61-67(1992)
391. H. G. Schuster, Deterministic Chaos (VCH, Weinheim, 1989).
392. G. Hildebraudt, C Gutenbrunne and R. Moog (eds) Chronobiology and Chronomedicine (Peter Lang, Frankfurt, 1922).
393. V. A. Davydov, V. S. Zykov and A. S. Mikhailov, Soviet. Phys. Usp. 34, pp. 665-684 (1991).
394. L. B. Keith, Wildlife's Ten Year Cycle (Univ. of Wisconsin Press, Madison, WI, 1963)
395. J. R. Cary, O. J. Rongstad, M. C. Brittingham, Wildl. Monogr. 90, 1 (1984)
396. L.B. Keith, Proceedings of the World Lagomorph Conference, Guelph, Ontario, 12 to 16 August 1979
397. K. Myers and C. D. Maclnnes, Eds. (Univ. of Guelph, Ontario, Canada, 1981), pp. 395-440
398. M. R. Vaughan and L. B. Keith, J. Wildl. Manage. 45, 354 (1981)
399. J.P. Finerty, The Population Ecology of Cycles in Small Mammals (Yale Univ. Press, New Haven, CT, 1980).
400. C. J. Krebs, Oikos 52, 143 (1988)
401. D. Chitty, Can. J. Zool. 38, 99 (1960).
402. C. J. Krebs, B. S. Gilbert, S. Boutin, A. R. E. Sinclair, J. N. M. Smith, J. Anim. Ecol. 55, 963 (1986)
403. A.R. E. Sinclair, C. J. Krebs, J. N. M. Smith, S. Boutin, J. Anim. Ecol. 57, 787 (1988)
404. J.N. M. Smith, C. J. Krebs, A. R. E. Sinclair, R. Boonstra, J. Anim. Ecol., p. 269
405. C. J. Krebs, S. Boutin, B. S. Gilbert, Oecologia 70, 194 (1985).
406. E. A. Desy and G. O. Batzli, Ecology 70, 411 (1989)
407. G.O. Batzli, in Wildlife 2001: Populations, D. R. McCullough and R. H. Barrett, Eds. (Elsevier, London, 1992), pp. 831-850.
408. Blocks were spaced at least 1 km apart. Within each block, we surveyed checkerboard grids of 20-by-20 points with 30.5-m spacing and used these grids for snowshoe hare live trapping. Two experimental areas were provided with supplemental food (commercial rabbit chow, 16% protein) year round. In the summer of 1987, we built one electric fence around 1 km2 to exclude mammalian predators, and over the following year we covered 10 ha with monofilament to reduce avian predation. The monofilament was never effective in preventing avian predation inside the electric fences, and consequently we did not rely on it as a part of the treatment. In the summer of 1988, we built a second electric fence around 1 km2 to use for the combined predator reduction-food addition treatment. We modified the design of the electric fences in 1988 to make them more effective, and since then they have worked effectively to prevent mammalian predators from entering the area. The fences are permeable to snowshoe hares. We could not replicate either the predator reduction or the predator reduction-food addition treatment because of maintenance costs and the difficulty of maintaining electric fences in the Yukon winter with-45°C temperatures. The fences had to be checked every day during winter. From 1976 to 1985, we trapped hares in six areas within the main study region and found that their population trajectories were very similar [3]. We thus have no reason to suspect strong area effects on the unreplicated predator reduction plots.
409. We used three control areas but were not able to trap hares in all of them every year. We have more detailed data on hares from control area 1. The three control areas had quite different histories during the increase phase from 1986 to 1988. Control area 3 reached its greatest hare density in 1988 and remained at a plateau until 1990. Control area 2 reached its peak density in 1990, and control area 1 reached its peak in 1989. By the late peak in 1990 and during the decline phase, the control areas were much more similar to each other in hare densities.
410. The electric fence was 10-stranded, 2.2 m in height, and carried 8600 V. Snow tracking of mammalian predators meeting the fence illustrated its effectiveness. We excluded mammalian predators virtually continuously from January 1989 onward. Our attempts to use monofilament fishing line as a deterrent to birds of prey was largely ineffective because ice formation and snow accumulation on the lines in winter caused them to break or collapse to the ground. We used mononlament on 10 ha of the predator exclosure but did not attempt to use it on the combination treatment area. The predator exclosures thus were mammalian predator exclosures and were still subject to avian predation.
411. We fertilized two 1-km2 blocks of forest with commercial fertilizer. In May 1987, we used ammonium nitrate at 25 g/m2. In May 1988, we switched to NPK fertilizer and used 17.5 g of N/m2, 5 g of P/m2, and 2.5 g of K/m2. In 1989, we used half this amount, and in the years 1990 to 1994 we used the full amount as in 1988. The fertilizer was spread aerially and we did ground checks to make sure it was uniformly spread. We do not present the data here to show the plant growth responses, but all elements of the flora responded dramatically to the added nutrients (C. J. Krebs et al., unpublished data).
412. D. Tilman, in Long-Term Studies in Ecology, G. E. Likens, Ed. (Springer-Verlag, New York, 1989), pp. 136-157
413. J. F. Franklin, Long-Term Studies in Ecology., pp. 3-19.
414. K. H. Pollock, J. D. Nichols, C. Brownie, J. E. Hines, Wildl Monogr. 107, 1 (1990)
415. K.H. Pollock, J. Wildl Manage. 46, 752 (1982). Confidence limits were typically plus minus 15% of the population estimates. Snow tracking in winter confirmed density estimates.
416. L. B. Keith, Curr. Mammal. 2, 119 (1990).
417. S. Boutin and C. J. Krebs, J. Wildl. Manage. 50, 592 (1986)
418. K.H. Pollock, S. R. Winterstein, C. M. Bunck, P. D. Curtis, J. Wildl. Manage. 53, 7 (1989).
419. D. Hik, Wildl. Res. 22, 115 (1995)
420. S. Boutin et al., unpublished data.
421. M. O'Donoghue and C. J. Krebs, J. Anim. Ecol. 61, 631 (1992)
422. J.R. Cary and L. B. Keith, Can. J.Zool. 57, 375 (1979).
423. D. Hik, unpublished data
424. R. Boonstra and G. R. Singleton, Gen. Comp. Endocrinol. 91, 126 (1993).
425. R. Boonstra, Evol. Ecol. 8, 196 (1994)
426. D. Hik, G. R. Singleton, in preparation.
427. In snowshoe hares, dispersal occurs mainly in the juvenile stage; adult hares are mostly sedentary [S. Boutin, B. S. Gilbert, C. J. Krebs, A. R. E. Sinclair, J. N. M. Smith, Can. J. Zool. 63, 106 (1985)].
428. Artificial food addition causes immigration of adult hares [S. Boutin, Oecologia 62, 393 (1984)], and population changes in food addition areas are more affected by movements than are changes at other treatment or control sites.
429. L. Oksanen, in Perspectives on Plant Consumption, D. Tilman and J. Grace, Eds. (Academic Press, New York, 1990), pp. 445-474
430. S. D. Fretwell, Oikos 50, 291 (1987)
431. S.R. Carpenter, J. F. Kitchell, J. R. Hodgson, Bioscience 35, 634(1985).
432. L. B. Keith, Oikos 40, 385 (1983).
433. J. O. Wolff, Ecol. Monogr. 50, 111 (1980).
434. We thank the Natural Sciences and Engineering Research Council of Canada for supporting this research program through a Collaborative Special Project grant; the Arctic Institute of North America, University of Calgary, for the use of the Kluane Lake research station; and all the graduate students, technicians, and undergraduates who have assisted in this project over the past 9 years, in particular V. Nams, S. Schweiger, and M. O'Donoghue.
435. S. A. Nadler, M. S. Hafner, J. C. Hafner, D. J. Hafner, Evolution 44, 942 (1990).
436. R. D. Price and K. C. Emerson, J. Med. Entomol 8, 228(1971).
437. M. S. Hafner and S. A. Nadler, Syst. Zool 39, 192 (1990)
438. M. S. Hafner and S. A. Nadler, Nature 332, 258 (1988).
439. R. D. M. Page, Syst. Zool. 39, 205 (1990).
440. J. W. Demastes and M. S. Hafner, J. Mammal. 74, 521 (1993).
441. R. D. M. Page, Int. J. Parasitol. 23, 499 (1993).
442. Voucher specimens are deposited in the Museum of Natural Science, Louisiana State University (LSUMZ) or the New Mexico Museum of Natural History (NMMNH) and are as follows: Orthogeomys underwoodi (LSUMZ 29493), O. hispidus (LSUMZ 29231), O. cavator (LSUMZ 29253), O. cherriei (LSUMZ 29539), O. heterodus (LSUMZ 29501), Geomys breviceps (LSUMZ 33940), G. personatus (LSUMZ 31460), G. bursarius halli (LSUMZ 31463; designated "a" in Fig. 2), G. b. majusculus (LSUMZ 31448; designated "b" in Fig. 2), Cratogeomys castanops (LSUMZ 31455), C. merriami (LSUMZ 34343), Pappogeomys bulleri (LSUMZ 34338), Zygogeomys trichopus (LSUMZ 34340), Thomomys bottae (LSUMZ 29320 and 29569), and T. talpoides (NMMNH 1634 and 1637).
443. Pocket gophers were trapped, killed, and immediately brushed to recover lice. Gopher tissues and lice were stored at negative 70°Celsius. DNA extractions from gophers followed the phenol-chloroform technique [D. M. Hillis, A. Larson, S. K. Davis, E. A. Zimmer, in Molecular Systematics, D. M. Hillis and C. Moritz, Eds. (Sinauer, Sunderland, MA, 1990), pp. 318-370].
444. DNA extractions from lice followed a modification of the protocol described by H. Liu and A. T. Beckenbach [Mol. Phylogenet. Evol. 1, 41 (1992)]
445. and used two lice from each gopher. One /JI of the extraction solution was used for DNA amplification in a 50-/J1 reaction. A 379-base pair (bp) region of the mitochondrial COI gene was amplified by polymerase chain reaction (PCR) with two degenerate primers: L6625 (5'-CCGGATCCTTYTGRTTYTTYGGNCAYCC-3') and H7005 (5'-CCGGATCC ACNACRTARTANGTRTCRTG-3')-Primer names refer to the 3' position of each primer relative to the human mitochondrial genome [S. Anderson et ah, Nature 290, 457 (1981)]. Double-stranded amplifications were done by four cycles of 1 min of denaturation (95"Celsius), 1 min of annealing (45"Celsius), and 1 min of extension (72°Celsius), followed by 30 cycles at reduced denaturation temperatures (93°Celsius) and increased annealing temperatures (60° Celsius). Methods for production of single-stranded templates and for sequencing of single-stranded products are given elsewhere [18]. Because no two sequences were identical either within or between lice and gophers, cross-contamination is not suspected. Gopher and louse sequences are available through GenBank (accession numbers L32682 to L32696 for gophers and L32665 to L32681 for lice).
446. D. L. Swofford and G. J. Olsen, in Molecular Systematics, D. M. Hillis and C. Moritz, Eds. (Sinauer, Sunderland, MA, 1990), pp. 411-501
447. D. M. Hillis, M. W. Allard, M. M. Miyamoto, Methods Enzymol. 224, 456 (1993).
448. J. P. Huelsenbeck and D. M. Hillis, Syst. Biol. 42, 247(1993).
449. We used the PAUP program (Phylogenetic Analysis Using Parsimony, 3.1.1, D. L. Swof ford, Illinois Natural History Survey) for phylogenetic analyses of the sequence data. of the 134 variable sites in pocket gophers (178 in lice), 113 were phylogenetically informative (142 in lice). To compensate for the effects of multiple mutations at nucleotide positions, we used a parsimony analysis in which transitions at first and third codon positions were designated as more likely than transversions (with a 10:1 step matrix for gophers and a 17:1 step matrix for lice). These empirical estimates of transition bias represent maximal values, which were observed between taxa with the most similar sequences (8.7% sequence divergence between the most similar gophers; 4.7% between the most similar lice). Because the effects of multiple mutations should be least for closely related taxa, we consider these to be reasonable approximations of the actual transition bias. For each database (gophers and lice) we conducted 100 heuristic searches, varying the addition order of taxa (using the random addition option of PAUP) and using alternate branch-swapping algorithms (using the TBR and SPR options). All gopher and louse parsimony trees contained significant phylogenetic signal (g\ <-0.54, P < 0.01) [D. M. Hillis and J. P. Huelsenbeck, J. Hered. 83, 189 (1992)]. Gophers in the genus Thomomys and lice of the genus Thomomydoecus were designated as outgroups. The outgroup status of these taxa is supported by previous morphological and molecular analyses [16, 17] [R. A. Hellenthal and R. D. Price, J. Kans. Entomol. Soc. 57, 231 (1984)
450. S.A. Nadler and M. S. Hafner, Ann. Entomol. Soc. Am. 82, 109 (1989)].
451. We used PHYLIP [Phylogeny Inference Package, 3. 5c, J. Felsenstein, Department of Genetics, University of Washington, Seattle] for maximum-likelihood, Fitch-Margoliash [W. M. Fitch and E. Margoliash, Science 155, 279 (1967)], and neighbor-joining [N. Saitou and M. Nei, Mol. Biol. Evol. 4, 406 (1987)] analyses of the sequence data. Empirical base frequencies and maximal observed transition biases (10:1 for gophers and 17:1 for lice) were used in the maximum-likelihood analysis. The Fitch-Margoliash and neighborjoining analyses were based on distance matrices corrected with the Kimura twoparameter model [substituting maximal observed transition bias values for the default (2:1) value M. Kimura, J. Mol. Evol. 16, 111 (1980)]. We repeated each analysis at least 10 times, varying the input order of taxa using the jumble option of PHYLIP.
452. In the alternative parsimony tree for pocket gophers, the phylogenetic positions of P. bulleri and Z. trichopus were reversed and the two species of Cratogeomys were not depicted as sister taxa.
453. The neighbor-joining tree differed from the Fitch-Margoliash tree (Fig. IB) only in the placement of P. bulleri (immediately basal to the five-taxon clade containing O. hispidus) and Z. trichopus (basal to the eight-taxon clade containing C. castanops, C. merriami, and P. bulleri).
454. H. Kishino and M. Hasegawa, J. Mol. Evol. 29, 170 (1989), as used in PHYLIP [12].
455. R. S. Russell, Univ. Kans. Publ. Mus. Nat. Hist. 16, 473 (1968)
456. R.L. Honeycutt and S. L. Williams, J. Mammal. 63, 208 (1982)
457. M.S. Hafner, J. Mammal. 72, 1 (1991).
458. M. S. Hafner, Z. Zool. Syst. Evolutionsforsch. 20, 118(1982).
459. P. D. Sudman and M. S. Hafner, Mol. Phylogenet. Evol. 1, 17 (1992).
460. In the second parsimony tree, the louse G. actuosi was linked with the oklahomensisgeomydis clade, rather than the texanusewingi clade (as in Fig. 2A). The third tree showed G. thomomyus near the root of the tree between the outgroups T. minor and T. barbarae.
461. The neighbor-joining tree was identical to the Fitch-Margoliash tree (Fig. 2B), except that G. thomomyus and G. perotensis were linked as sister taxa (as in the parsimony and maximum-likelihood trees, Fig. 2A) and positioned basal to the G. chapini branch.
462. S. A. Nadler and M. S. Hafner, Int. J. Parasitol. 23, 191 (1993).
463. COMPONENT 2.0, R. D. M. Page, Natural History Museum, London
464. R. D. M. Page, Syst. Biol. 43, 58 (1994).
465. The test was based on the criterion of minimum number of independent losses necessary to reconcile the host and parasite trees [22]. Sixteen of the 20 pairwise comparisons were also found to be statistically significant (P <-0.05) according to the criterion of "number of leaves added" [22]. Number of leaves added is equivalent to onehalf the number of "items of error" [G. Nel son and N. Platnick, Systematics and Biogeography: Cladistics and Vicariance (Columbia Univ. Press, New York, 1981)]. The four comparisons that were not statistically significant (according to the criterion of number of leaves added) all involved the louse parsimony tree that depicted the outgroup Thomomydoecus as polyphyletic. The branching structure of this tree is challenged by evidence from previous morphological and allozyme studies [11] supporting monophyly of the louse genus Thomomydoecus.
466. The trees (Fig. 2) were reduced to the largest number of taxa that showed identical branching patterns in the hosts and parasites. All possible combinations of cospeciating taxa were compared. For example, there are four possible combinations of nine cospeciating taxa in Fig. 2A [that is, the nine-taxon gopher tree can include either O. underwoodi or O. cavator (but not both) and G. breviceps or G. personatus (but not both)]. There are 2 possible combinations of 10 cospeciating taxa in Fig. 2B (that is, the tree can include either G. breviceps or G. personatus, but not both). Because most of the uncertainty in the phylogenetic analyses involved branches near the base of the trees, only terminal and subterminal branches were compared between gophers and lice.
467. N. Mantel, Cancer Res. 27, 209 (1967). We calculated distance matrices using the maximum-likelihood model [12]. For the 9 cospeciating taxa in Fig. 2A, t = 4.073, P < 0.01. For the 10 cospeciating taxa in Fig. 2B, t = 4.289, P < 0.01. Degrees of freedom were adjusted to n-1, where n = equals number of taxa compared [3].
468. Trees for cospeciating taxa (24) were input as user trees into PHYLIP [12] in order to estimate maximum-likelihood branch lengths. Empirical estimates of transition bias, nucleotide frequency bias, and positional bias (first:second:third codon position: 5:1:40 for gophers and 5:1:21 for lice) were incorporated into the maximum-likelihood model [12]. We recognize that the parameters used in maximum-likelihood models will influence rate estimates [K. Fukami-Kobayashi and Y. Tateno, J. Mol. Evol. 32, 79 (1991)]. Our models use empirical estimates of these parameters (rather than arbitrary or constant values) to obtain a more realistic approximation of the actual mutational dynamics of the gene region being surveyed [12, 28].
469. We used Model II regression analysis (reduced major axis model), which does not permit conventional tests of significance [R. R. Sokal and F. J. Rohlf, Biometry (W. H. Freeman, San Francisco, 1969)]. However, the nonparametric Mantel's test documents a significant relation between genetic distances in gophers and lice [25], and the non-parametric Wilcoxon sign-rank test shows that louse branches are significantly longer than gopher branches. When the regression analyses were not constrained through the origin, none of the y-intercepts was significantly different from zero.
470. R. Kondo, S. Horai, Y. Satta, N. Takahata, J. Mol. Evol. 36, 517 (1993)
471. R.A. Capaldi, Annu. Rev. Biochem. 59, 569 (1990). Total sequence differences among the five species of Orthogeomys range from 4 to 11.9%. Sequence differences among the parasites of Orthogeomys range from 11.9 to 20.3%.
472. Differences in nucleotide composition were tested using the B-statistic [A. P. Martin, G. J. P. Naylor, S. P. Palumbi, Nature 357, 153 (1992)].
473. M. Kimura, The Neutral Theory of Molecular Evolution (Cambridge Univ. Press, London, 1983)
474. K.H. Wolfe, P. M. Sharp, W.-H. Li, Nature 337, 283 (1989).
475. The slope of this relationship is affected very little by use of different nucleotide frequency estimates in the maximum-likelihood model. The slopes are more sensitive, however, to different estimates of transition:transversion bias. Our model used empirical estimates of transition:transversion bias (4:1 in gophers; 10:1 in lice) because such estimates are more likely to reflect the actual mutational dynamics of the gene regions being surveyed [12, 28].
476. G. G. Brown and M. V. Simpson, Proc. Natl. Acad. Sci. U.S.A. 79, 3246 (1982); [28].
477. A.P. Martin and S. R. Palumbi, Proc. Natl. Acad. Sci. U.S.A. 90, 4087 (1993)
478. C.-I. Wu and W.-H. Li, Proc. Natl. Acad. Sci. U.S.A. 82, 1741 (1985).
479. R. W. Rust, Oecologia 15, 287 (1974)
480. A.G. Marshall, The Ecology of Ectoparasitic Insects (Academic Press, London, 1981).
481. We thank D. R. Wolstenholme, R. Okimoto, and C. T. Beakley for providing PCR primer sequences and protocols. Photography was provided by R. J. Bouchard. D. H. Clayton, S. V. Edwards, D. A. Good, D. J. Hafner, B. D. Marx, R. D. M. Page, S. Scheiner, and two anonymous reviewers provided helpful advice and comments on this manuscript. Supported by NSF grant BSR-8817329 and NSFLouisiana Education Quality Support Fund grant ADP-02.
482. Adrian, R. H. (1956). J. Physiol. 133, 631.
483. Aikin, O. C. & Thomas, R. C. (1975). J. Physiol. 252, 803.
484. Crick, C. M. (1941). Patterns and Problems in Development. Chicago: University of Chicago Press.
485. Crick, F. H. C. (1970). Nature, Lond. 225, 420.
486. Crick, F. H. C. & Lawrence, P. A. (1975). Science, JV. Y. 189, 340.
487. Dalcq, A. & Pasteels, J. (1937). Arch. Blol. 48, 669.
488. Edelstein, B. B. (1972). J. theor. Biol. 37, 221.
489. Furshipan, E. J. & Potter, D. D. (1968). Curr. Top. Devel. Biol. 3, 95.
490. Garcia-Bellido, A. (1975). In Cell Patterning, Ciba Foundation Symp. 29, 161. Amsterdam: Associated Scientific Publishers.
491. Gehring, W. (1972). in Biology of Imaginal Disks (H. Ursprung & R. Nothiger, eds) p.35. Berlin: Springer.
492. Gierer, A. & Meinhardt, H. (1972). Kybernetik 12, 30.
493. Goodwin, B. & Cohen, M. H. (1969). J. theor. Biol. 25, 49.
494. Hodgkin, A. L. (1964). The Conduction of the Nervous Impulse. Liverpool: University Press.
495. Kacser, H. (1960). Symp. Soc. exp. Biol. 14, 13.
496. Kauffman, S. (1975). In Cell Patterning, Ciba Foundation Symp. 29, 201. Amsterdam: Associated Scientific Publishers.
497. Kennedy, E. P. (1970). In The Lactose Operon (J. Beckwith & D. Zipser, eds), Cold Spring Harbor Lab.
498. Keynes, R. D. (1963). J. Physiol. 169, 690.
499. Killander, D. & Zetterbero, A. (1965). Expl. Cell Res. 38, 272.
500. Lawrence, P. A. (1975). In Cell Patterning, Ciba Foundation Symp. 29, 2. Amsterdam: Associated Scientific Publishers.
501. Lawrence, P. A., Crick, F. H. C. & Munro, M. (1972). J. cell Set 11, 815.
502. Loewenstein, W. R. (1973). Fedn. Proc. Fedn. Am. Socs. exp. Biol. 32, 60.
503. Maynard-Smith, J. (1960). Proc. R. Soc. B. 152, 397.
504. Mohr, H. & Oelze-Karow, H. (1976). In Light and Plant Development (H. Smith, ed) Kent: Butterworth.
505. Newsholme, E. A. & Start, C. (1973). Regulation in Metabolism. London: John Wiley & Sons Ltd.
506. Prigogine, I. & Nicolis, G. (1971). Q. Rev. Biophys. 4, 107.
507. Sander, K. (1975). In Cell Patterning, Ciba Foundation Symp. 29, 241. Amsterdam: Associated Scientific Publishers.
508. Slack, C. & Warner, A. E. (1975). J. Physiol 248, 97.
509. Spemann, H. (1938). In Embryonic Development and Induction, p.331. Yale: Yale University Press.
510. Sternberg, M. S. (1970). J. expl. Zool. 173, 395.
511. Summerbell, D., Lewis, J. H. & Wolpert, L. (1973). Nature, Lond. 244, 492.
512. Thorn, R. (1975). Structural Stability and Morphogenesis. Reading, Mass.: Benjamin.
513. Thornley, J. H. M. (1972). An. Bot. 36, 861.
514. Tickle, C, Summerbell, D. & Wolpert, L (1975). Nature, Lond. 254, 199.
515. Wolpert, L. (1971). Curr. Top. devel. Biol 6, 183.
516. Wolpert, L, Clarke, M. R. B. & Hornbruch, A. (1972). Nature, Lond. 239, 101.
517. Wolpert, L. & Lewis, J. H. (1975). Fedn. Proc. Fedn. Am. Socs. exp. Biol. 34, 1
518. G. H. Denton and W. Karlen, Quat. Res. 3, 155 (1973).
519. R. B. Alley et al, Geology 25, 483 (1997).
520. W. S. Broecker, Nature 372, 421 (1994).
521. T. Fronval and E. jansen, Paleoceanography 12, 443 (1997).
522. S. R. O'Brien et al., Science 270, 1962 (1995).
523. Core samples were collected with precision cutting tools to remove material in narrow slots at centimeter intervals. We avoided sampling the core's outer rind and obvious disturbances such as burrowing. All samples (dry) were weighed, then washed to separate the fine (<62 ^jm) from the coarse sediment fraction. All petrologic analyses were done in the 63-to 150-jjm size range, which consistently contains at least a few hundred lithic grains. A portion of that grain fraction was placed on a glass slide (with glycerin as a mounting medium) and 300 to 500 grains were counted per sample, using a line rather than point counting method. A key factor in our ability to measure certain tracers accurately is a specially prepared microscope. We placed a white reflector on the condenser lens and illuminated the sample slide from above with a halogen light source. By moving the substage reflector up and down, a position can be found that creates a strong impression of relief and brings into striking view details of surface textures and coatings on grains. The technique is especially useful for identifying hematite-stained grains, even when the stain is quite small (Figure 3). Counting of lithic grain concentrations was done in the >150 £/m fraction. To avoid errors introduced by splitting small grain populations, we did not split samples for this procedure. Counting of planktonic foraminiferal concentrations and abundances was done on splits of the >150 ji/m fraction. For the abundance measurements, 300 to 500 individuals were counted. Species in the >150 ji/m fraction were selected for planktonic and benthic isotopic measurements. All isotopic measurements were made at North Carolina State University by W.S. The samples were crushed and washed with reversed-osmosis water and airdried at 65°C, and then 30-to 100-fjg sample aliquots were isotopically analyzed in a Kiel Autocarbonate device attached to a FMAT 251 RMS. NBS-19, NBS-18, and the NCSU CM-1 marble standard were analyzed with every sample run for standard calibration and sample correction. During the period these samples were analyzed, for samples in the 10-to 100-/L/g size range the standard reproducibility was 0.07 per mil for 13C and 0.08 per mil for 18O. Planktonic isotopic analyses for VM 28-14 were made on 18 to 20 individuals per sample and were replicated at a number of depths. For VM 29-191, 10 repli cate analyses were done on two or three individuals from each centimeter depth with the objective of investigating small-sample variability over time. Because those results showed no coherent relation to our other proxies, we combined all measurements from each centimeter with a mass balance Equation using the sample size (micromolar) and the sample isotopic composition per mil, the results of which are shown in Figure 2. The final result is equivalent to isotopically analyzing 20 planktonic individuals at once.
524. G. C. Bond and R. Lotti, Science 267, 1005 (1995).
525. The counting error is the 2a standard deviation for the mean of 20 replicate counts in each of two representative samples.
526. That we find any measurable record of drift ice as far south as 54°N is not as extraordinary as it might seem. Between about 1880 and the mid-1960s there were more than 20 iceberg sightings in the eastern North Atlantic between 45°and 60°N, and icebergs have been reported from as far south as Bermuda and the Azores [40]. Even so, keeping in mind that the lithic grain concentrations are very low and that each sample integrates 50 to 100 years of time, the amounts of drifting ice reaching either site, even at peak concentrations, must have been small, and certainly much less than during the glaciation.
527. T. Johannessen, E. Jansen, A. Flatoy, A. C. Raveo, in Carbon Cycling in the Glacial Ocean: Constraints on the Ocean's Role in Global Change, R. Zahn et ah, Eds. (NATO ASI Series, Vol. I 17, Springer-Verlag, Berlin, 1994), pp. 61-85.
528. S. R. Epstein, R. Buchsbaum, H. A. Lowenstam, H. C. Urey, Geol. Soc. Am. Bull. 64, 1315 (1953).
529. M. Sarnthein et ah, Paleoceanography 10, 1063 (1995)
530. G.Wu and C. Hillaire-Marcel, Geochim. Cosmochim. Acta 58, 1303 (1994).
531. D. W. Oppo, M. Horowitz, S. J. Lehman, Paleoceanography 12, 51 (1997).
532. T. B. Kellogg, in Climatic Changes on a Yearly to Millennial Basis, N.-A. Morner and W. Karlen, -Eds. (Reidel, Dordrecht, Netherlands, 1984), pp. 123-133.
533. L. Labeyrie et al, Philos. Trans. R. Soc. London 348, 255 (1995).
534. B. E. Viekman and K. D. Baumer, International Ice Patrol Technical Report 95-03 (Department of Transportation, U.S. Coast Guard, International Ice Patrol, 1995).
535. Most of the cores (indicated by solid dots in Figure 1) have Holocene sections, documented on the basis of radiocarbon dating, isotopic measurements, or the presence of ash layer 1 (-10, 000 years) at depths of tens of centimeters. The existence of Holocene sections in the remainder, indicated by plus signs, is inferred from distinct color patterns and bulk carbonate measurements [41]. Core top ages likely vary from close to the present to perhaps a few thousand years; on average, then, the core top petrologic data reflect distributions of the tracers over time as well as over a large geographic area. The petrographic analyses of each core top lithic sample were done in the same way as for VM 29-191 and VM 28-14.
536. There probably were no tidewater glaciers in Iceland, at least during the middle and late parts of the Holocene, ruling out transport of debris to the sites by Icelandic icebergs; direct transport in eruption clouds is also ruled out by the large size of clasts (coarsest grains typically are 0.5 to 1 mm) and by evidence that only two of the six largest Holocene eruptions in Iceland are correlative with the petrologic events [42]. The only plausible mechanism, therefore, is eruption of Icelandic volcanic material onto nearby sea ice (and probably glacier ice as well) and subsequent transport of that ice in surface currents to the coring sites. Explosive eruptions from volcanoes such as Hekla, Katla, and Grimsvotn occur once every 20 to 30 years or so in Iceland [43], and the Icelandic lithic grains at both sites are dominantly (70 to 95%) clear rhyolitic glass shards, the characteristic product of explosive eruptions in silicic volcanoes. The importance of our interpretation of the Icelandic IRD is that its increases at VM 28-14 require increases in drift ice within reach of eruption clouds to the east. Hence, changes in surface hydrography favoring preservation and transportation of drift ice in the western North Atlantic must have taken place, at the very least, from Denmark Strait to the vicinity of Iceland or Jan Mayen.
537. S. Sathyendranath, A. Longhurst, C. M. Caverhill, T. Platt, Deep Sea Res. Pt. 1, 42, 1773 (1995).
538. The carbon isotopic composition of C. wuellerstorfi appears to record the 13C of total CO2 of sea water, and therefore is taken as a proxy of the nutrient content of that water [44]. Today, VM 29-191 lies within NADW, which is nutrient-depleted and hence <513C-enriched; below the depth of the core are modified waters of Southern Ocean origin, which are nutrient-enriched and <513Cdepleted [45], Variations in the <513C of that foraminifera are related to shifts in the rate of production of NADW and hence in the rate of convective overturning and thermohaline circulation in the North Atlantic [34, 44].
539. E. Boyle, Philos. Trans. R. Soc. London 348, 243 (1995).
540. M. S. McCartney, Prog. Oceanogr. 29, 283 (1992).
541. P. A. Mayewski et al, Science 263, 1747 (1994).
542. The calendar age model in the composite record for the interval between the Younger Dryas and 18, 000 radiocarbon years was constructed by linear interpolation between radiocarbon dates from VM 23-81 calibrated following Bard et al [46] (Table 1). For the remainder of the record, we transferred the GISP2 ages of Heinrich events 2 and 3 (H2 and H3) [47] into their equivalent levels in VM 23-81 (Table 1) [48] and interpolated linearly between the two ages. Our estimate of the age of H3 is corroborated by calibrated radiocarbon measurements in Barbados corals (Table 1) and by paired U-Th-14C measurements in aragonite from Lake Lisan in the Dead Sea Rift [49]. Detrital carbonate at the Younger Dryas level was deposited rapidly and is widespread, extending from the eastern North Atlantic to sources in Hudson Strait and Cumberland Strait [50]; it serves as an excellent chronostratigraphic marker. Calibrations of radiocarbon ages in VM 19-30 (Table 1) at 100 and 150 cm are from [46]. Remaining calibrations are from the age model for the Holocene-glacial composite record as described above.
543. The Heinrich events are tied to welldocumented lowerings of atmospheric and ocean surface temperatures and to large reductions in the rate of NADW formation (Figure 6) [34]. The two distinct ice-rafting peaks between interstadial 1 and the Younger Dryas, dated at-13, 200 and-14, 000 calendar years (-11, 300 and 12, 00014C years), are accompanied by increases in N. pachyderma (s.) measured in the same core [51], and they are coincident with prominent decreases in <518O in Greenland ice (Figure 6). Evidence of coolings at about the same times have been found in other deep sea cores from the North Atlantic [52] and the Nordic Seas [53] and in pollen and glacial records in Europe and North America [53]. From 23, 000 calendar years to the end of the record, the IRDpetrologic episodes between Heinrich events also occur at times of maxima in abundances of N. pachyderma (s.) [7] and at times of prominent cold phases in the ice core <518O records (Figure 6). We also find distinct IRD events punctuating prolonged stadials, such as between interstadials 2 and 3 and interstadials 4 and 5. There, the IRD-petrologic peaks correspond to prominent increases in dissolved Ca (Figure 6) and in the polar circulation index [23], suggesting a close association of the marine events with expansion of polar circulation above the ice cap. Within the long interval between interstadials 1 and 2, we have identified four other IRD-petrologic events in addition to HI. The N. pachyderma (s.) abundances are saturated through the entire interval, preventing reliable estimates of changes in sea surface temperatures. The youngest of the four events, dated at-15, 200 calendar years (-12, 90014C years), however, is correlative with glacial advances in the U.S. midwest [for example, Port Huron Stade [54]], with a glacial readvance on the Scotian shelf [55], and with a Heinrich-like cooling and ice-rafting event identified on the Scotian Slope [56]. The other three events, a, b, and c, previously identified by Bond and Lotti [7] and dated at about 19, 000, 21, 000, and 22, 000 calendar years (-16, 500, 17, 500, and ~18, 60014C years), are present in IRD records from a number of other sites in the subpolar North Atlantic and GIN Seas [57], and they also appear in magnetic susceptibility records from the Labrador Sea [58]. Two of these widespread events, a and c, are correlative with prominent coolings identified in pollen records from a long lacustrine sequence in northern Norway [59]. There is evidence that the two older events, b and c, may even be correlative with glacial advances in the Southern Hemisphere, implying a bihemi
544. E. J. Brook, T. Sowers, J. Orchardo, Science 273, 1087 (1996).
545. We define event pacing as the time separating adjacent peaks placed at the midpoints between the peaks (Figure 6). We define peaks as maxima with at least two points and amplitudes equal to or exceeding the 2a counting error (4% [8]); the measure of separation is made between maxima in peak percentage. At the levels of HI and H2, the sharp decreases in percentages of the two petrologic tracers reflect dilution by massive amounts of IRD discharged from the Labrador Sea [7], and we place the peaks for those two events at maxima in the lithic grain concentrations (Figure 6).
546. P. A. Mayewski et al, J. Geophys. Res., in press.
547. D. J. Thompson, Proc. IEEE 70, 1055 (1982).
548. P. Yiou and J. Jouzel, Geophys. Res. Lett 22, 2179 (1995)
549. L.D. Keigwin and G. A. Jones, J. Geophys. Res. 99, 12397 (1994)
550. E. Cortijo, P. Yiou, L. Labeyrie, M. Cremer, Paleoceanography 10, 911 (1995).
551. J. C. Rogers, J. Clim. 3, 1364 (1990).
552. J. M. Grove, The Little Ice Age (Methuen, London, 1988).
553. L. D. Keigwin, Science 274, 1504 (1996).
554. D. W. Oppo and S. J. Lehman, Paleoceanography 10, 901 (1995).
555. S. Rahmstorf, Nature 378, 145 (1995)
556. S. Manabe and R. J. Stouffer, Paleoceanography 12, 321 (1997)
557. T.F. Stocker, D. G. Wright, W. S. Broecker, Paleoceanography. 7, 529 (1992).
558. K. Sakai and W. R. Peltier, J. Clim. 10, 949 (1997)
559. E. Tziperman, Nature 386, 592 (1997)
560. A.J. Weaver and T. M. C. Hughes, Nature. 367, 447(1994).
561. R. T. Sutton and M. R. Allen, Nature 388, 563 (1997)
562. Y. Kushnir, J. Clim. 7, 141 (1994).
563. P. Pestiaux, I. Van der Mersch, A. Berger, J. C. Duplessy, Clim. Change 12, 9 (1988)
564. P. Yiou et al, J-Geophys. Res. 96, 20365 (1991).
565. M. Stuiver and T. F. Braziunas, The Holocene 3, 289 (1993)
566. E.A. Frieman, Solar Influences on Global Change (National Academy Press, Washington, DC, 1994).
567. Oceanographic Atlas of the North Atlantic (U.S. Naval Oceanographic Office, Washington, DC, 1968)
568. J.E. Murry, in Ice Seminar, V. E. Bohme, Ed. (Canadian Institute of Mining and Metallurgy, Washington, DC, 1969), pp. 3-18.
569. T. B. Kellogg, thesis, Columbia University (1973).
570. A. J. Dugmore, G. Larsen, A. J. Newton, The Holocene 5, 257 (1995)
571. S. Thorarinsson, in Tephra Studies, S. Self and R. S. J. Sparks, Eds. (Reidel, Dordrecht, Netherlands, 1981), pp. 109-134
572. R. A. J. Cas and J. V. Wright, Volcanic Successions (Allen and Unwin, Winchester, MA, 1987)
573. H.J. Gudmundsson, Quat. Sci. Rev. 16, 81 (1997).
574. T. Simkin et al, Smithsonian Institution Volcanoes of the World (Hutchinson Ross, Stroudsburg, PA, 1981).
575. E. A. Boyle and L. D. Keigwin, Nature 330, 35 (1987).
576. P. Kroopnick, Earth Planet. Sci. Lett. 49, 469 (1980).
577. E. Bard, M. Arnold, R. G. Fairbanks, B. Hamelin, Radiocarbon 35, 191 (1993)
578. M. Stuiver and J. Reimer, Radiocarbon., p. 215.
579. M. J. Bender et al, Nature 372, 663 (1994).
580. G. C. Bond et al, Nature 365, 143 (1993).
581. S. Goldstein, personal communication.
582. J. T. Andrews et al., Paleoceanography 10, 943 (1995).
583. S. J. Lehman and L. D. Keigwin, Nature 356, 757(1992).
584. H. P. Sejrup, H. Haflidason, D. K. Kristensen, J. Quat. Sci. 10, 385 (1995).
585. N. Koc, E. Jansen, M. Hald, L. Labeyrie, in Late Quaternary Palaeoceanography of the North Atlantic Margins, J. T. Andrews, W. E. N. Austin, H. Bergsten, A. E. Jennings, Eds. (Geological Society Special Publication 111, 1996), pp. 177-185.
586. G. H. Denton and T. J. Hughes, in The Last Great Ice Sheets (Wiley, New York, 1981).
587. R. R. Stea, R. Boyd, O. Costello, G. B. J. Fader, D. B. Scott, in [53], pp. 77-101.
588. L. D. Keigwin and G. A. Jones, Paleoceanography 10, 973 (1995).
589. R. Stein, S.-I. Nam, H. Grobe, H. Hubberten, in [53], pp. 135-151
590. T. Fronval and E. Jansen, Nature 383, 806 (1996)
591. H. Haflidason, H. P. Sejrup, D. K. Kristensen, S. Johnsen, Geology 23, 1059(1995).
592. J. S. Stoner, J. E. T. Channell, C. Hillaire-Marcel, Paleoceanography 11, 309 (1996).
593. T. Aim, Boreas 22, 171 (1993).
594. T. V. Lowell et al, Science 269, 1541 (1995).
595. I. M. Lagerklint, thesis, University of Maine (1995).
596. E. Bard et al, Earth Planet. Sci. Lett. 126, 275 (1994).
597. S. D. Woodruff, R. J. Slutz, R. L. Jenne, P. M. Steurer, Bull Am. Meteorol. Soc. 68, 521 (1987).
598. A. Mclntyre and B. Molfino, Science 274, 1867 (1996).
599. W. Dansgaard et al, Nature 364, 218 (1993).
600. N. J. Shackleton, J. Imbrie, M. A. Hall, Earth Planet. Sci. Lett. 65, 233 (1983).
601. D. G. Martinson et al, Quat. Res. 27, 1 (1987).
602. We thank G. Kukla, J. Lynch-Stieglitz, and A. Van Geen for comments on the manuscript. Supported in part by grants from NSF and the National Oceanic and Atmospheric Administration. Support for the core curation facilities of the Lamont-Doherty Earth Observatory Deep-Sea Sample Repository is provided by NSF through grant OCE94-02150 and the Office of Naval Research through grant N00014-96-1-0186.This is L-DEO contribution 5714.
603. May, R. M. J. theor. Biol. 51, 511-524 (1975).
604. May, R. M. Nature 261, 459-467 (1976).
605. May, R. M. & Oster, G. F. Am. Nat. 110, 573-599 (1976).
606. Strong, D. R. in Ecological Theory and Integrated Pest Management (ed. Kogan, M.) 37-58 (Wiley, New York, 1986).
607. Kareiva, P. in Perspectives in Ecological Theory (ed. Roughgarden, J., May, R. M. & Levin, S. A.) 68-88 (Princeton University Press, Princeton, NJ, 1989).
608. Bartlett, M. S. J. R. Statist. Soc. A 52, 321-347 (1990).
609. Logan, J. A. & Hain, F. (eds) Chaos and Insect Ecology (Virginia Polytechnic Institute and State University, Blacksburg, VA, 1991).
610. Costantino, R. F. & Desharnais, R. A. Population Dynamics and the Tribolium Model: Ge netics and Demography (Springer, New York, 1991).
611. Logan, J. A. & Allen, J. C. A. Rev. Ent. 37, 455-477(1992).
612. Hastings, A., Horn, C. L, Ellner, S., Turchin, P. & Godfray, H. C. J. A. Rev. Ecol. Syst. 24, 1-33 (1993).
613. Murdock, W. W. & McCauley, E. Nature 316, 628-630(1985).
614. Grenfell, B. J., Price, O. J., Albon, S. D. & Clutton-Brock, T. H. Nature 355, 823-826 (1992).
615. Hanski, I., Turchin, P., Korplmakl, E. & Henttonen, H. Nature 364, 232-235 (1993).
616. Tilman, D. & Wedin, D. Nature 353, 653-655 (1991).
617. Dennis, B., Desharnais, R. A., Cushing, J. M. & Costantino, R. F. Ecol. Monogr. (in press).
618. Tong, H. Non-linear Time Series: a Dynamical System Approach (Oxford Univ. Press, Oxford, 1990).
619. McCaffrey, D. J., Ellner, S., Gallant, A. R. & Nychka, D. W. J. Am. statist. Ass. 87, 682-695 (1992).
620. Turchin, P. & Taylor, A. Ecology 73, 289-305 (1992).
621. Turchin, P. Oikos 68, 167-172 (1993).