
small (250x250 max)
medium (500x500 max)
large ( > 500x500)
Full Resolution


Adaptive Harvest Management 2009 HHunttinggSeasson U.S. Fish & Wildlife Service Adaptive Harvest Management 2009 Hunting Season PREFACE The process of setting waterfowl hunting regulations is conducted annually in the United States (Blohm 1989). This process involves a number of meetings where the status of waterfowl is reviewed by the agencies responsible for setting hunting regulations. In addition, the U.S. Fish and Wildlife Service (USFWS) publishes proposed regulations in the Federal Register to allow public comment. This document is part of a series of reports intended to support development of harvest regulations for the 2009 hunting season. Specifically, this report is intended to provide waterfowl managers and the public with information about the use of adaptive harvest management (AHM) for setting waterfowl hunting regulations in the United States. This report provides the most current data, analyses, and decisionmaking protocols. However, adaptive management is a dynamic process and some information presented in this report will differ from that in previous reports. ACKNOWLEDGMENTS A working group consisting of representatives from the USFWS, the U.S. Geological Survey (USGS), the Canadian Wildlife Service (CWS), and the four Flyway Councils (Appendix 1) was established in 1992 to review the scientific basis for managing waterfowl harvests. The working group, supported by technical experts from the waterfowl management and research communities, subsequently proposed a framework for adaptive harvest management, which was first implemented in 1995. The USFWS expresses its gratitude to the AHM Working Group and to the many other individuals, organizations, and agencies that have contributed to the development and implementation of AHM. This report was prepared by the USFWS Division of Migratory Bird Management. G. S. Boomer and T. A. Sanders were the principal authors. Individuals that provided essential information or otherwise assisted with report preparation were G. Zimmerman, M. Koneff, K. Richkus, E. Silverman, N. Zimpfer, J. Klimstra, K. Magruder, and P. Garrettson. Comments regarding this document should be sent to the Chief, Division of Migratory Bird Management  USFWS, 4401 North Fairfax Drive, MS MSP4107, Arlington, VA 22203. We are grateful for the continuing technical support from F. A. Johnson, M. C. Runge, and J. A. Royle (USGS), and acknowledge that information provided by USGS in this report has not received the Director's approval and, as such, is provisional and subject to revision. This information is released on the condition that neither the USGS nor the United States Government may be held liable for any damages resulting from its authorized or unauthorized use. 1 Citation: U.S. Fish and Wildlife Service. 2009. Adaptive Harvest Management: 2009 Hunting Season. U.S. Dept. Interior, Washington, D.C. 52pp. Online: http://www.fws.gov/migratorybirds/CurrentBirdIssues/Management/AHM/AHMintro.htm U.S. Fish & Wildlife Service Cover art: Joshua Spies’s painting of a longtailed duck (Clangula hyemalis) that was selected for the 2009 federal “duck stamp.”TABLE OF CONTENTS Executive Summary.............................................................................................................3 Background.........................................................................................................................4 Mallard Stocks and Flyway Management............................................................................5 Mallard Population Dynamics..............................................................................................6 HarvestManagement Objectives.......................................................................................13 Regulatory Alternatives.....................................................................................................14 Optimal Regulatory Strategies...........................................................................................18 Application of AHM Concepts to Other Stocks................................................................21 Emerging Issues within AHM............................................................................................26 Literature Cited..................................................................................................................27 Appendix 1: AHM Working Group..................................................................................30 Appendix 2: Midcontinent Mallard Models....................................................................33 Appendix 3: Eastern Mallard Models...............................................................................36 Appendix 4: Western Mallard Models..............................................................................39 Appendix 5: Modeling Mallard Harvest Rates.................................................................44 Appendix 6: Scaup Model.................................................................................................49 2 3 EXECUTIVE SUMMARY In 1995 the U.S. Fish and Wildlife Service (USFWS) implemented the Adaptive Harvest Management (AHM) program for setting duck hunting regulations in the United States. The AHM approach provides a framework for making objective decisions in the face of incomplete knowledge concerning waterfowl population dynamics and regulatory impacts. The AHM protocol is based on the population dynamics and status of three mallard (Anas platyrhynchos) stocks. Midcontinent mallards are defined as those breeding in the Waterfowl Breeding Population and Habitat Survey (WBPHS) strata 13–18, 20–50, and 75–77 plus mallards breeding in the states of Michigan, Minnesota, and Wisconsin (state surveys). The prescribed regulatory alternative for the Mississippi and Central Flyways depends exclusively on the status of these mallards. Eastern mallards are defined as those breeding in WBPHS strata 51–54 and 56 and breeding in the states of Virginia northward into New Hampshire (Atlantic Flyway Breeding Waterfowl Survey [AFBWS]). The regulatory choice for the Atlantic Flyway depends exclusively on the status of these mallards. Western mallards are defined as those birds breeding in WBPHS strata 1�����12 (hereafter Alaska) and those birds breeding in the states of California and Oregon (state surveys). The regulatory choice for the Pacific Flyway depends exclusively on the status of these mallards. Mallard population models are based on the best available information and account for uncertainty in population dynamics and the impact of harvest. Modelspecific weights reflect the relative confidence in alternative hypotheses and are updated annually using comparisons of predicted and observed population sizes. For midcontinent mallards, current model weights favor the weakly densitydependent reproductive hypothesis (88%) and suggest some preference for the additivemortality hypothesis (62%). For eastern mallards, virtually all of the weight is on models that have corrections for bias in estimates of survival or reproductive rates. Model weights do not discriminate between the strongly densitydependent (47%) and weakly densitydependent (53%) reproductive hypotheses. By consensus, hunting mortality is assumed to be additive in eastern mallards. Unlike midcontinent and eastern mallards, we consider a single functional form to predict western mallard population dynamics but consider a wide range of parameter values each weighted relative to the support from the data. For the 2009 hunting season, the USFWS is considering the same regulatory alternatives as last year. The nature of the restrictive, moderate, and liberal alternatives has remained essentially unchanged since 1997, except that extended framework dates have been offered in the moderate and liberal alternatives since 2002. Harvest rates associated with each of the regulatory alternatives have been updated based on bandreporting rate studies conducted since 1998. The expected harvest rates of adult males under liberal hunting seasons are 0.119 (SD = 0.020), 0.149 (SD = 0.043), and 0.115 (SD = 0.032) for midcontinent, eastern, and western mallards, respectively. Optimal regulatory strategies for the 2009 hunting season were calculated using: (1) harvestmanagement objectives specific to each mallard stock; (2) the 2009 regulatory alternatives; and (3) current population models. Based on this year’s survey results of 8.71 million midcontinent mallards, 3.57 million ponds in Prairie Canada, 0.908 million eastern mallards, and 0.884 million western mallards in Alaska (0.503 million) and California–Oregon (0.381 million), the optimal choice for all four flyways is the liberal regulatory alternative. AHM concepts and tools are also being applied to help improve harvest management for several other waterfowl stocks. In the last year, progress has been made in understanding the harvest potential of northern pintails (Anas acuta) and scaup (Aythya affinis, A. marila). While these biological assessments are ongoing, they are already informing decision makers and proving valuable in helping focus debate on the social aspects of harvest policy, including management objectives and the nature of regulatory alternatives.4 BACKGROUND The annual process of setting duckhunting regulations in the United States is based on a system of resource monitoring, data analyses, and rulemaking (Blohm 1989). Each year, monitoring activities such as aerial surveys and hunter questionnaires provide information on population size, habitat conditions, and harvest levels. Data collected from this monitoring program are analyzed each year, and proposals for duckhunting regulations are developed by the Flyway Councils, States, and USFWS. After extensive public review, the USFWS announces regulatory guidelines within which States can set their hunting seasons. In 1995, the USFWS adopted the concept of adaptive resource management (Walters 1986) for regulating duck harvests in the United States. This approach explicitly recognizes that the consequences of hunting regulations cannot be predicted with certainty and provides a framework for making objective decisions in the face of that uncertainty (Williams and Johnson 1995). Inherent in the adaptive approach is an awareness that management performance can be maximized only if regulatory effects can be predicted reliably. Thus, adaptive management relies on an iterative cycle of monitoring, assessment, and decisionmaking to clarify the relationships among hunting regulations, harvests, and waterfowl abundance. In regulating waterfowl harvests, managers face four fundamental sources of uncertainty (Nichols et al. 1995a, Johnson et al. 1996, Williams et al. 1996): (1) environmental variation  the temporal and spatial variation in weather conditions and other key features of waterfowl habitat; an example is the annual change in the number of ponds in the Prairie Pothole Region, where water conditions influence duck reproductive success; (2) partial controllability  the ability of managers to control harvest only within limits; the harvest resulting from a particular set of hunting regulations cannot be predicted with certainty because of variation in weather conditions, timing of migration, hunter effort, and other factors; (3) partial observability  the ability to estimate key population attributes (e.g., population size, reproductive rate, harvest) only within the precision afforded by extant monitoring programs; and (4) structural uncertainty  an incomplete understanding of biological processes; a familiar example is the longstanding debate about whether harvest is additive to other sources of mortality or whether populations compensate for hunting losses through reduced natural mortality. Structural uncertainty increases contentiousness in the decisionmaking process and decreases the extent to which managers can meet longterm conservation goals. AHM was developed as a systematic process for dealing objectively with these uncertainties. The key components of AHM include (Johnson et al. 1993, Williams and Johnson 1995): (1) a limited number of regulatory alternatives, which describe Flywayspecific season lengths, bag limits, and framework dates; (2) a set of population models describing various hypotheses about the effects of harvest and environmental factors on waterfowl abundance; (3) a measure of reliability (probability or "weight") for each population model; and (4) a mathematical description of the objective(s) of harvest management (i.e., an "objective function"), by which alternative regulatory strategies can be compared. These components are used in a stochastic optimization procedure to derive a regulatory strategy. A regulatory strategy specifies the optimal regulatory choice, with respect to the stated management objectives, for each possible combination of breeding population size, environmental conditions, and model weights (Johnson et al. 1997). The setting of annual hunting regulations then involves an iterative process: (1) each year, an optimal regulatory choice is identified based on resource and environmental conditions, and on current model weights; 5 (2) after the regulatory decision is made, modelspecific predictions for subsequent breeding population size are determined; (3) when monitoring data become available, model weights are increased to the extent that observations of population size agree with predictions, and decreased to the extent that they disagree; and (4) the new model weights are used to start another iteration of the process. By iteratively updating model weights and optimizing regulatory choices, the process should eventually identify which model is the best overall predictor of changes in population abundance. The process is optimal in the sense that it provides the regulatory choice each year necessary to maximize management performance. It is adaptive in the sense that the harvest strategy “evolves” to account for new knowledge generated by a comparison of predicted and observed population sizes. MALLARD STOCKS AND FLYWAY MANAGEMENT Since its inception AHM has focused on the population dynamics and harvest potential of mallards, especially those breeding in midcontinent North America. Mallards constitute a large portion of the total U.S. duck harvest, and traditionally have been a reliable indicator of the status of many other species. As management capabilities have grown, there has been increasing interest in the ecology and management of breeding mallards that occur outside the midcontinent region. Geographic differences in the reproduction, mortality, and migrations of mallard stocks suggest that there may be corresponding differences in optimal levels of sport harvest. The ability to regulate harvests of mallards originating from various breeding areas is complicated, however, by the fact that a large degree of mixing occurs during the hunting season. The challenge for managers, then, is to vary hunting regulations among Flyways in a manner that recognizes each Flyway’s unique breedingground derivation of mallards. Of course, no Flyway receives mallards exclusively from one breeding area; therefore, Flywayspecific harvest strategies ideally should account for multiple breeding stocks that are exposed to a common harvest. The optimization procedures used in AHM can account for breeding populations of mallards beyond the midcontinent region, and for the manner in which these ducks distribute themselves among the Flyways during the hunting season. An optimal approach would allow for Flywayspecific regulatory strategies, which represent an average of the optimal harvest strategies for each contributing breeding stock weighted by the relative size of each stock in the fall flight. This joint optimization of multiple mallard stocks requires: (1) models of population dynamics for all recognized stocks of mallards; (2) an objective function that accounts for harvestmanagement goals for all mallard stocks in the aggregate; and (3) decision rules allowing Flywayspecific regulatory choices. Currently, three stocks of mallards are officially recognized for the purposes of AHM (Fig. 1). We use a constrained approach to the optimization of these stocks’ harvest, in which the Atlantic Flyway regulatory strategy is based exclusively on the status of eastern mallards, the regulatory strategy for the Mississippi and Central Flyways is based exclusively on the status of midcontinent mallards, and the Pacific Flyway regulatory strategy is based exclusively on the status of western mallards. This approach has been determined to perform nearly as well as a jointoptimization because mixing of the three stocks during the hunting season is limited and because of the constraints imposed by management objectives and regulatory alternatives.Fig. 1. Survey areas currently assigned to the midcontinent, eastern, and western stocks of mallards for the purposes of AHM. MALLARD POPULATION DYNAMICS MidContinent Stock In 2008, midcontinent mallards were redefined as those breeding in WBPHS strata 13–18, 20–50, and 75–77, and in the Great Lakes region (Michigan, Minnesota, and Wisconsin; see Fig. 1). Estimates of the size of this population are available since 1992, and have varied from 6.4 to 11.2 million (Table 1, Fig. 2). Estimated breedingpopulation size in 2009 was 8.71 million (SE = 0.25 million), including 8.01 million (SE = 0.24 million) from the WBPHS and 0.696 million (SE = 0.056 million) from the Great Lakes region. Details describing the set of population models for midcontinent mallards are provided in Appendix 2. The set consists of four alternatives, formed by the combination of two survival hypotheses (additive vs. compensatory hunting mortality) and two reproductive hypotheses (strongly vs. weakly density dependent). Relative weights for the alternative models of midcontinent mallards changed little until all models underpredicted the change in population size from 1998 to 1999, perhaps indicating there is a significant factor affecting population dynamics that is absent from all four models (Fig. 3). Updated model weights suggest some preference for the additive 6 Table 1. Estimates (N) and associated standard errors (SE) of midcontinent mallards (in millions) observed in the WBPHS (strata 13–18, 20–50, and 75–77) and the Great Lakes region (Michigan, Minnesota, and Wisconsin). WBPHS area Great Lakes region Total Year N SE N SE N SE 1992 5.6304 0.2379 0.9946 0.1597 6.6249 0.2865 1993 5.4253 0.2068 0.9347 0.1457 6.3600 0.2529 1994 6.6292 0.2803 1.1505 0.1163 7.7797 0.3035 1995 7.7452 0.2793 1.1214 0.1965 8.8666 0.3415 1996 7.4193 0.2593 1.0251 0.1443 8.4444 0.2967 1997 9.3554 0.3041 1.0777 0.1445 10.4331 0.3367 1998 8.8041 0.2940 1.1224 0.1792 9.9266 0.3443 1999 10.0926 0.3374 1.0591 0.2122 11.1518 0.3986 2000 8.6999 0.2855 1.2350 0.1761 9.9348 0.3354 2001 7.1857 0.2204 0.8622 0.1086 8.0479 0.2457 2002 6.8364 0.2412 1.0820 0.1152 7.9184 0.2673 2003 7.1062 0.2589 0.8360 0.0734 7.9422 0.2691 2004 6.6142 0.2746 0.9333 0.0748 7.5474 0.2847 2005 6.0521 0.2754 0.7862 0.0650 6.8383 0.2830 2006 6.7607 0.2187 0.5881 0.0465 7.3488 0.2236 2007 7.7258 0.2805 0.7677 0.0584 8.4935 0.2865 2008 7.1914 0.2525 0.6750 0.0478 7.8664 0.2570 2009 8.0094 0.2442 0.6958 0.0564 8.7052 0.2506 02468101219911995199920032007YearPopulation size (millions)024681012WBPHS surveyGreat LakesTotal 7 Fig. 2. Population estimates of midcontinent mallards observed in the WBPHS (strata: 13–18, 20–50, and 75–77) and the Great Lakes region (Michigan, Minnesota, and Wisconsin). Error bars represent one standard error. 00.10.20.30.40.50.619951997199920012003200520072009YearModel Weight00.10.20.30.40.50.6ScRwSaRsSaRwScRs Fig. 3. Weights for models of midcontinent mallards (ScRs = compensatory mortality and strongly densitydependent reproduction, ScRw = compensatory mortality and weakly densitydependent reproduction, SaRs = additive mortality and strongly densitydependent reproduction, and SaRw = additive mortality and weakly densitydependent reproduction). Model weights were assumed to be equal in 1995. mortality models (62%) over those describing hunting mortality as compensatory (38%). For most of the time frame, model weights have strongly favored the weakly densitydependent reproductive models over the strongly densitydependent ones, with current model weights of 88% and 12%, respectively. The reader is cautioned, however, that models can sometimes make reliable predictions of population size for reasons having little to do with the biological hypotheses expressed therein (Johnson et al. 2002b). Eastern Stock Eastern mallards are defined as those breeding in southern Ontario and Quebec (WBPHS strata 51–54 and 56) and in the northeastern U.S. (AFBWS; Heusman and Sauer 2000; see Fig. 1). Estimates of population size have varied from 0.815 to 1.1 million since 1990, with the majority of the population accounted for in the northeastern U.S. (Table 2, Fig. 4). For 2009, the estimated breedingpopulation size of eastern mallards was 0.908 million (SE = 0.063 million), including 0.667 million (SE = 0.046 million) from the northeastern U.S. and 0.241 million (SE = 0.043 million) from the WBPHS. The reader is cautioned that these estimates differ from those reported in the USFWS annual waterfowl trend and status reports, which include composite estimates based on more fixedwing strata in eastern Canada and helicopter surveys conducted by the Canadian Wildlife Service (CWS). Details concerning the set of population models for eastern mallards are provided in Appendix 3. The set consists of six alternatives, formed by the combination of two reproductive hypotheses (strongly vs. weakly density dependent) and three hypotheses concerning bias in estimates of survival and reproductive rates (no bias vs. biased survival rates vs. biased reproductive rates). With respect to model weights, there is no single model that is clearly favored over the others at the current time. Collectively, current model weights provide little discrimination between the weakly densitydependent or strongly density dependent reproductive hypotheses, with current model weights of 53% and 47%, respectively (Fig. 5). In addition, there is overwhelming evidence of bias in extant estimates of survival or reproductive rates, assuming that survey estimates are unbiased. 8 9 Table 2. Estimates (N) and associated standard errors (SE) of eastern mallards (in millions) observed in the northeastern U.S. (AFBWS) and southern Ontario and Quebec (WBPHS strata 51–54 and 56). Northeastern U.S. Canadian survey strata Total Year N SE N SE N SE 1990 0.6651 0.0783 0.1907 0.0472 0.8558 0.0914 1991 0.7792 0.0883 0.1528 0.0337 0.9320 0.0945 1992 0.5622 0.0479 0.3203 0.0530 0.8825 0.0715 1993 0.6866 0.0499 0.2921 0.0482 0.9786 0.0694 1994 0.8563 0.0628 0.2195 0.0282 1.0758 0.0688 1995 0.8641 0.0704 0.1844 0.0400 1.0486 0.0810 1996 0.8486 0.0611 0.2831 0.0557 1.1317 0.0826 1997 0.7952 0.0496 0.2121 0.0396 1.0073 0.0634 1998 0.7752 0.0497 0.2638 0.0672 1.0390 0.0836 1999 0.8800 0.0602 0.2125 0.0369 1.0924 0.0706 2000 0.7626 0.0487 0.1323 0.0264 0.8948 0.0554 2001 0.8094 0.0516 0.2002 0.0356 1.0097 0.0627 2002 0.8335 0.0562 0.1915 0.0319 1.0250 0.0647 2003 0.7319 0.0470 0.3083 0.0554 1.0402 0.0726 2004 0.8066 0.0517 0.3015 0.0533 1.1081 0.0743 2005 0.7536 0.0536 0.2934 0.0531 1.0470 0.0755 2006 0.7214 0.0476 0.1740 0.0284 0.8954 0.0555 2007 0.6876 0.0467 0.2193 0.0336 0.9069 0.0576 2008 0.6191 0.0407 0.1960 0.0300 0.8151 0.0505 2009 0.6668 0.0457 0.2411 0.0434 0.9078 0.0630 00.20.40.60.811.21.419891994199920042009YearPopulation size (millions)00.20.40.60.811.21.4AFBWSWBPHS Total Fig. 4. Population estimates of eastern mallards observed in the northeastern states (AFBWS) and in southern Ontario and Quebec (WBPHS strata 51–54 and 56). Error bars represent one standard error. 00.050.10.150.20.250.30.3519951997199920012003200520072009YearModel WeightRw0Rs0RwSRsSRwRRsR Fig. 5. Weights for models of eastern mallards (Rw0 = weak densitydependent reproduction and no model bias, Rs0 = strong dependent reproduction and no model bias, RwS = weak densitydependent reproduction and biased survival rates, RsS = strong densitydependent reproduction and biased survival rates, RwR = weak densitydependent reproduction and biased reproductive rates, and RsR = strong densitydependent reproduction and biased reproductive rates). Model weights were assumed to be equal in 1996. 10 11 Western Stock Western mallards consist of 2 substocks and are defined as those birds breeding in Alaska (WBPHS strata 1–12) and those birds breeding in California and Oregon (state surveys; see Fig. 1). Estimates of the size of these subpopulations have varied from 0.283 to 0.843 million in Alaska since 1990 and 0.355 to 0.694 million in California and Oregon since 1992 (Table 3, Fig. 6). The total population size of western mallards has ranged from 0.748 to 1.407 million. Ideally, the western mallard stock assessment would account for mallards breeding in the states of the Pacific Flyway (including Alaska), British Columbia, and the Yukon Territory. However, we have had continuing concerns about our ability to determine changes in population size based on the collection of surveys conducted independently by Pacific Flyway States and the CWS in British Columbia. These surveys tend to vary in design and intensity, and in some cases lack measures of precision. We reviewed extant surveys to determine their adequacy for supporting a westernmallard AHM protocol and selected Alaska, California, and Oregon for modeling purposes. These three states likely harbor about 75% of the westernmallard breeding population. Nonetheless, this geographic delineation is considered temporary until surveys in other areas can be brought up to similar standards and an adequate record of population estimates is available for analysis. Details concerning the set of population models for western mallards are provided in Appendix 4. To predict changes in abundance we relied on a discrete logistic model, which combines reproduction and natural mortality into a single parameter, r, the intrinsic rate of growth. This model assumes densitydependent growth, which is regulated by the ratio of population size, N, to the carrying capacity of the environment, K (i.e., equilibrium population size in the absence of harvest). In the traditional formulation of the logistic model, harvest mortality is completely additive and any compensation for hunting losses occurs as a result of densitydependent responses beginning in the subsequent breeding season. To increase the model’s generality we included a scaling parameter for harvest that allows for the possibility of compensation prior to the breeding season. It is important to note, however, that this parameterization does not incorporate any hypothesized mechanism for harvest compensation and, therefore, must be interpreted cautiously. We modeled Alaska mallards independently of those in California and Oregon because of differing population trajectories (see Fig. 6) and substantial differences in the distribution of band recoveries. We used Bayesian estimation methods in combination with a statespace model that accounts explicitly for both process and observation error in breeding population size (Meyer and Millar 1999). Breeding population estimates of mallards in Alaska are available since 1955, but we had to limit the time series to 1990–2008 because of changes in survey methodology and insufficient bandrecovery data. The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed time series of Alaska population estimates. The estimated mean carrying capacity was 1.1 million, the intrinsic rate of growth was 0.31, while the scaling parameter estimate suggests that harvest mortality may be additive. Breeding population and harvestrate data were available for California–Oregon mallards for the period 1992–2008. The logistic model also provided a reasonable fit to these data, suggesting a mean carrying capacity of 0.7 million, an intrinsic rate of growth of 0.33, while the scaling parameter estimate suggests that harvest mortality may be only partially additive. 12 Table 3. Estimates (N) and associated standard errors (SE) of mallards (in millions) observed in Alaska (WBPHS strata 1–12) and California and Oregon (state surveys) combined. Alaska California–Oregon Total Year N SE N SE N SE 1990 0.3669 0.0370 1991 0.3853 0.0363 1992 0.3457 0.0387 0.4835 0.0605 0.8292 0.0718 1993 0.2830 0.0295 0.4654 0.0510 0.7484 0.0589 1994 0.3509 0.0371 0.4367 0.0426 0.7876 0.0565 1995 0.5242 0.0680 0.4541 0.0428 0.9783 0.0803 1996 0.5220 0.0436 0.6451 0.0802 1.1671 0.0912 1997 0.5842 0.0520 0.6390 0.1043 1.2232 0.1166 1998 0.8362 0.0673 0.4868 0.0489 1.3230 0.0832 1999 0.7131 0.0696 0.6937 0.1066 1.4068 0.1273 2000 0.7703 0.0522 0.4639 0.0532 1.2342 0.0745 2001 0.7183 0.0541 0.4044 0.0451 1.1227 0.0705 2002 0.6673 0.0507 0.3775 0.0327 1.0449 0.0603 2003 0.8435 0.0668 0.4340 0.0501 1.2775 0.0835 2004 0.8111 0.0639 0.3547 0.0352 1.1658 0.0729 2005 0.7031 0.0547 0.4014 0.0474 1.1045 0.0724 2006 0.5158 0.0469 0.4879 0.0576 1.0037 0.0743 2007 0.5815 0.0551 0.4900 0.0546 1.0715 0.0775 2008 0.5324 0.0468 0.3814 0.0478 0.9138 0.0669 2009 0.5030 0.0449 0.3815 0.0639 0.8844 0.0781 Year19901992199419961998200020022004200620082010Population size (millions)0.20.40.60.81.01.21.41.60.20.40.60.81.01.21.41.6AKCAORTotal Fig. 6. Population estimates of western mallards observed in Alaska (WBPHS strata 1–12) and California and Oregon (state surveys) combined. Error bars represent one standard error. Ideally, the development of AHM protocols for mallards would consider how different breeding stocks distribute themselves among the four Flyways so that Flywayspecific harvest strategies could account for the mixing of birds during the hunting season. At present, however, a joint optimization of western, midcontinent, and eastern stocks is not feasible due to computational hurdles. However, our preliminary analyses suggest that the lack of a joint optimization does not result in a significant decrease in performance. Therefore, the AHM protocol for western mallards is structured similarly to that used for eastern mallards, in which an optimal harvest strategy is based on the status of a single breeding stock and harvest regulations in a single flyway. Although the contribution of midcontinent mallards to the Pacific Flyway harvest is significant, we believe an independent harvest strategy for western mallards poses little risk to the midcontinent stock. Further analyses will be needed to confirm this conclusion, and to better understand the potential effect of midcontinent mallard status on sustainable hunting opportunities in the Pacific Flyway. HARVESTMANAGEMENT OBJECTIVES The basic harvestmanagement objective for midcontinent mallards is to maximize cumulative harvest over the long term, which inherently requires perpetuation of a viable population. Moreover, this objective is constrained to avoid regulations that could be expected to result in a subsequent population size below the goal of the North American Waterfowl Management Plan (NAWMP). According to this constraint, the value of harvest decreases proportionally as the difference between the goal and expected population size increases. This balance of harvest and population objectives results in a regulatory strategy that is more conservative than that for maximizing long 13 14 term harvest, but more liberal than a strategy to attain the NAWMP goal (regardless of effects on hunting opportunity). The current objective for midcontinent mallards uses a population goal of 8.5 million birds, which consists of 7.9 million mallards from the WBPHS (strata 13–18, 20–50, and 75–77) corresponding to the mallard population goal in the 1998 update of the NAWMP (less the portion of the mallard goal comprised of birds breeding in Alaska) and a goal of 0.6 million for the combined states of Michigan, Minnesota, and Wisconsin. For eastern and western mallards, there is no NAWMP goal or other established target for desired population size. Accordingly, the management objective for eastern and western mallards is simply to maximize longterm cumulative (i.e., sustainable) harvest. Additionally for western mallards, maximum longterm cumulative harvest is subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes. REGULATORY ALTERNATIVES Evolution of Alternatives When AHM was first implemented in 1995, three regulatory alternatives characterized as liberal, moderate, and restrictive were defined based on regulations used during 1979–84, 1985–87, and 1988–93, respectively. These regulatory alternatives also were considered for the 1996 hunting season. In 1997, the regulatory alternatives were modified to include: (1) the addition of a veryrestrictive alternative; (2) additional days and a higher duck bag limit in the moderate and liberal alternatives; and (3) an increase in the bag limit of hen mallards in the moderate and liberal alternatives. In 2002 the USFWS further modified the moderate and liberal alternatives to include extensions of approximately one week in both the opening and closing framework dates. In 2003 the veryrestrictive alternative was eliminated at the request of the Flyway Councils. Expected harvest rates under the veryrestrictive alternative did not differ significantly from those under the restrictive alternative, and the veryrestrictive alternative was expected to be prescribed for < 5% of all hunting seasons. Also in 2003, at the request of the Flyway Councils the USFWS agreed to exclude closed duckhunting seasons from the AHM protocol when the population size of midcontinent mallards was ≥ 5.5 million (WBPHS strata 1–18, 20–50, and 75–77 plus the Great Lakes region). Based on our original assessment, closed hunting seasons did not appear to be necessary from the perspective of sustainable harvesting when the midcontinent mallard population exceeded this level. The impact of maintaining open seasons above this level also appeared negligible for other midcontinent duck species, as based on population models developed by Johnson (2003). In 2008, because of the redefinition of the midcontinent mallard stock that excludes mallards breeding in Alaska, we rescaled the closedseason constraint. Initially, we attempted to adjust the original 5.5 million closure threshold by subtracting out the 1985 Alaska breeding population estimate, which was the year upon which the original closed season constraint was based. Our initial rescaling resulted in a new threshold equal to 5.25 million. Simulations based on optimal policies using this revised closed season constraint suggested that the Mississippi and Central Flyways would experience a 70% increase in the frequency of closed seasons. At this time, we agreed to consider alternative rescalings in order to minimize the effects on the midcontinent mallard strategy and account for the increase in mean breeding population sizes in Alaska over the past several decades. Based on this assessment, we recommended a revised closed season constraint of 4.75 million which resulted in a strategy performance equivalent to the performance expected prior to the redefinition of the midcontinent mallard stock. Because the performance of the revised strategy is essentially unchanged from the original strategy, we believe it will have no greater impact on other duck stocks in the Mississippi and Central Flyways. However, complete or partialseason closures for particular species or populations could still be deemed necessary in some situations regardless of the status of midcontinent mallards. Details of the regulatory alternatives for each Flyway are provided in Table 4. 15 Table 4. Regulatory alternatives for the 2009 duckhunting season. Flyway Regulation Atlantica Mississippi Centralb Pacific Shooting hours onehalf hour before sunrise to sunset Framework dates Restrictive Oct 1 – Jan 20 Saturday nearest Oct 1to the Sunday nearest Jan 20 Moderate and Liberal Saturday nearest September 24 to the last Sunday in January Season length (days) Restrictive 30 30 39 60 Moderate 45 45 60 86 Liberal 60 60 74 107 Bag limit (total / mallard / female mallard) Restrictive 3 / 3 / 1 3 / 2 / 1 3 / 3 / 1 4 / 3 / 1 Moderate 6 / 4 / 2 6 / 4 / 1 6 / 5 / 1 7 / 5 / 2 Liberal 6 / 4 / 2 6 / 4 / 2 6 / 5 / 2 7 / 7 / 2 a The states of Maine, Massachusetts, Connecticut, Pennsylvania, New Jersey, Maryland, Delaware, West Virginia, Virginia, and North Carolina are permitted to exclude Sundays, which are closed to hunting, from their total allotment of season days. b The High Plains Mallard Management Unit is allowed 12, 23, and 23 extra days in the restrictive, moderate, and liberal alternatives, respectively. c The Columbia Basin Mallard Management Unit is allowed seven extra days in the restrictive and moderate alternatives. RegulationSpecific Harvest Rates Harvest rates of mallards associated with each of the openseason regulatory alternatives were initially predicted using harvestrate estimates from 1979–84, which were adjusted to reflect current hunter numbers and contemporary specifications of season lengths and bag limits. In the case of closed seasons in the U.S., we assumed rates of harvest would be similar to those observed in Canada during 1988–93, which was a period of restrictive regulations both in Canada and the U.S. All harvestrate predictions were based only in part on bandrecovery data, and relied heavily on models of hunting effort and success derived from hunter surveys (Appendix C in USFWS 2002). As such, these predictions had large sampling variances and their accuracy was uncertain. In 2002, we began relying on Bayesian statistical methods for improving regulationspecific predictions of harvest rates, including predictions of the effects of frameworkdate extensions. Essentially, the idea is to use existing (prior) information to develop initial harvestrate predictions (as above), to make regulatory decisions based on those predictions, and then to observe realized harvest rates. Those observed harvest rates, in turn, are treated as 16 new sources of information for calculating updated (posterior) predictions. Bayesian methods are attractive because they provide a quantitative, formal, and an intuitive approach to adaptive management. For midcontinent mallards, we have empirical estimates of harvest rate from the recent period of liberal hunting regulations (1998–2008). The Bayesian methods thus allow us to combine these estimates with our prior predictions to provide updated estimates of harvest rates expected under the liberal regulatory alternative. Moreover, in the absence of experience (so far) with the restrictive and moderate regulatory alternatives, we reasoned that our initial predictions of harvest rates associated with those alternatives should be rescaled based on a comparison of predicted and observed harvest rates under the liberal regulatory alternative. In other words, if observed harvest rates under the liberal alternative were 10% less than predicted, then we might also expect that the mean harvest rate under the moderate alternative would be 10% less than predicted. The appropriate scaling factors currently are based exclusively on prior beliefs about differences in mean harvest rate among regulatory alternatives, but they will be updated once we have experience with something other than the liberal alternative. A detailed description of the analytical framework for modeling mallard harvest rates is provided in Appendix 5. Our models of regulationspecific harvest rates also allow for the marginal effect of frameworkdate extensions in the moderate and liberal alternatives. A previous analysis by the USFWS (2001) suggested that implementation of frameworkdate extensions might be expected to increase the harvest rate of midcontinent mallards by about 15%, or in absolute terms by about 0.02 (SD = 0.01). Based on the observed harvest rates during the 2002–2008 hunting seasons, the updated (posterior) estimate of the marginal change in harvest rate attributable to the frameworkdate extension is 0.008 (SD = 0.007). The estimated effect of the frameworkdate extension has been to increase harvest rate of midcontinent mallards by about 7% over what would otherwise be expected in the liberal alternative. However, the reader is strongly cautioned that reliable inference about the marginal effect of frameworkdate extensions ultimately depends on a rigorous experimental design (including controls and random application of treatments). Current predictions of harvest rates of adultmale midcontinent mallards associated with each of the regulatory alternatives are provided in Table 5. Predictions of harvest rates for the other age–sex cohorts are based on the historical ratios of cohortspecific harvest rates to adultmale rates (Runge et al. 2002). These ratios are considered fixed at their longterm averages and are 1.5407, 0.7191, and 1.1175 for young males, adult females, and young females, respectively. We make the simplifying assumption that the harvest rates of midcontinent mallards depend solely on the regulatory choice in the Mississippi and Central Flyways. Table 5. Predictions of harvest rates of adultmale midcontinent mallards expected with application of the 2009 regulatory alternatives in the Mississippi and Central Flyways. Regulatory alternative Mean SD Closed (U.S.) 0.0088 0.0019 Restrictive 0.0563 0.0129 Moderate 0.1029 0.0215 Liberal 0.1188 0.0198 17 The predicted harvest rates of eastern mallards are updated in the same fashion as that for midcontinent mallards based on reward banding conducted in eastern Canada and the northeastern U.S. (Appendix 5). Like midcontinent mallards, harvest rates of age and sex cohorts other than adult male mallards are based on constant rates of differential vulnerability as derived from bandrecovery data. For eastern mallards, these constants are 1.153, 1.331, and 1.509 for adult females, young males, and young females, respectively (Johnson et al. 2002a). Regulationspecific predictions of harvest rates of adultmale eastern mallards are provided in Table 6. In contrast to midcontinent mallards, frameworkdate extensions were expected to increase the harvest rate of eastern mallards by only about 5% (USFWS 2001), or in absolute terms by about 0.01 (SD = 0.01). Based on the observed harvest rates during the 2002–2008 hunting seasons, the updated (posterior) estimate of the marginal change in harvest rate attributable to the frameworkdate extension is 0.003 (SD = 0.010). The estimated effect of the frameworkdate extension has been to increase harvest rate of eastern mallards by about 2% over what would otherwise be expected in the liberal alternative. Table 6. Predictions of harvest rates of adultmale eastern mallards expected with application of the 2009 regulatory alternatives in the Atlantic Flyway. Regulatory alternative Mean SD Closed (U.S.) 0.0797 0.0233 Restrictive 0.1108 0.0393 Moderate 0.1412 0.0473 Liberal 0.1492 0.0433 Based on available estimates of harvest rates of mallards banded in California and Oregon during 1990–1995 and 2002–2007, there was no apparent relationship between harvest rate and regulatory changes in the Pacific Flyway. This is unusual given our ability to document such a relationship in other mallard stocks and in other species. We note, however, that the period 2002–2007 was comprised of both stable and liberal regulations and harvest rate estimates were based solely on reward bands. Regulations were relatively restrictive during most of the earlier period and harvest rates were estimated based on standard bands using reporting rates estimated from reward banding during 1987–1988. Additionally, 1993–1995 were transition years in which fulladdress and tollfree bands were being introduced and information to assess their reporting rates (and their effects on reporting rates of standard bands) is limited. Thus, the two periods in which we wish to compare harvest rates are characterized not only by changes in regulations, but also in estimation methods. Consequently, we lack a sound empirical basis for predicting harvest rates of western mallards associated with current regulatory alternatives in the Pacific Flyway. This year, we applied Bayesian statistical methods for improving regulationspecific predictions of harvest rates (see Appendix 5). The methodology is analogous to that currently in use for midcontinent and eastern mallards except that the marginal effect of framework date extensions in moderate and liberal alternatives is inestimable because there are no data prior to implementation of extensions. We specified prior regulationspecific harvest rates of 0.01, 0.06, 0.09, and 0.11 with associated standard deviations of 0.003, 0.02, 0.03, and 0.03 for the closed, restrictive, moderate, and liberal alternatives, respectively. The harvest rates for the liberal alternative were based on empirical estimates realized under the current liberal alternative during 2002–2007 and determined from adultmale mallards banded with reward bands in California and Oregon. Harvest rates for the moderate and restrictive alternatives were based on the proportional (0.85 and 0.51) difference in harvest rates expected for midcontinent mallards under the respective alternatives. And finally, harvest rate for the closed alternative was based on what we might realize with a closed season in the U.S. (including Alaska) and a very restrictive season in Canada, similar to that for midcontinent mallards. A relatively large standard deviation (CV=0.3) was chosen to reflect greater uncertainty about the means than that for midcontinent mallards (CV=0.2). Current predictions of harvest rates of adultmale western mallards associated with each regulatory alternative are provided in Table 7. 18 Table 7. Predictions of harvest rates of adultmale western mallards expected with application of the 2009 regulatory alternatives in the Pacific Flyway. Regulatory alternative Mean SD Closed (U.S.) 0.010 0.018 Restrictive 0.059 0.016 Moderate 0.098 0.027 Liberal 0.115 0.032 OPTIMAL REGULATORY STRATEGIES We calculated optimal regulatory strategies using stochastic dynamic programming (Lubow 1995, Johnson and Williams 1999). For the Mississippi and Central Flyways, we based this optimization on: (1) the 2009 regulatory alternatives, including the closedseason constraint; (2) current population models and associated weights for midcontinent mallards; and (3) the dual objectives of maximizing longterm cumulative harvest and achieving a population goal of 8.5 million midcontinent mallards. The resulting regulatory strategy (Table 8) is similar to that used last year. Note that prescriptions for closed seasons in this strategy represent resource conditions that are insufficient to support one of the current regulatory alternatives, given current harvestmanagement objectives and constraints. However, closed seasons under all of these conditions are not necessarily required for longterm resource protection, and simply reflect the NAWMP population goal and the nature of the current regulatory alternatives. Assuming that regulatory choices adhered to this strategy (and that current model weights accurately reflect population dynamics), breedingpopulation size would be expected to average 6.83 million (SD = 1.83 million). Based on an estimated population size of 8.71 million midcontinent mallards and 3.57 million ponds in Prairie Canada, the optimal choice for the Mississippi and Central Flyways in 2009 is the liberal regulatory alternative. We calculated an optimal regulatory strategy for the Atlantic Flyway based on: (1) the 2009 regulatory alternatives; (2) current population models and associated weights for eastern mallards; and (3) an objective to maximize longterm cumulative harvest. The resulting strategy suggests liberal regulations for all population sizes of record, and is characterized by a lack of intermediate regulations (Table 9). We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model weights accurately reflect population dynamics), breedingpopulation size would be expected to average 0.911 million (SD = 0.173 million). Based on an estimated breeding population size of 0.908 million mallards, the optimal choice for the Atlantic Flyway in 2009 is the liberal regulatory alternative. We calculated an optimal regulatory strategy for the Pacific Flyway based on: (1) the 2009 regulatory alternatives, (2) current (1990–2008) population models and parameter estimates, and (3) an objective to maximize longterm cumulative harvest subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes (Table 10). We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model parameters accurately reflect population dynamics), breedingpopulation size would be expected to average 1.0 million (SD = 0.22 million) in Alaska and 0.46 million (SD = 0.03 million) in California and Oregon. Based on an estimated breeding population size of 0.503 million mallards in Alaska and 0.381 million in California and Oregon, the optimal choice for the Pacific Flyway in 2009 is the liberal regulatory alternative (see Table 10).19 Table 8. Optimal regulatory strategya for the Mississippi and Central Flyways for the 2009 hunting season. This strategy is based on current regulatory alternatives (including the closedseason constraint), on current midcontinent mallard models and weights, and on the dual objectives of maximizing longterm cumulative harvest and achieving a population goal of 8.5 million mallards. The shaded cell indicates the regulatory prescription for 2009. Pondsc Bpopb 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 ≤ 4.5 C C C C C C C C C C 4.75–5.75 R R R R R R R R R R 6.00 R R R R R R R R M M 6.25 R R R R R R M M M L 6.5 R R R R M M M L L L 6.75 R R R M L L L L L L 7.0 R M M L L L L L L L 7.25 M M L L L L L L L L 7.5 M L L L L L L L L L ≥ 7.75 L L L L L L L L L L a C = closed season, R = restrictive, M = moderate, L = liberal. b Mallard breeding population size (in millions) in the WBPHS (strata 13–18, 20–50, 75–77) and Michigan, Minnesota, and Wisconsin. c Ponds (in millions) in Prairie Canada in May. Table 9. Optimal regulatory strategya for the Atlantic Flyway for the 2009 hunting season. This strategy is based on current regulatory alternatives, on current eastern mallard models and weights, and on an objective to maximize longterm cumulative harvest. The shaded cell indicates the regulatory prescription for 2009. Mallardsb Regulation ≤ 0.250 C 0.275 R ≥ 0.300 L a C = closed season, R = restrictive, M = moderate, and L = liberal. b Estimated number of mallards in eastern Canada (WBPHS strata 51–54, 56) and the northeastern U.S. (AFBWS), in millions. 20 Table 10. Optimal regulatory strategya for the Pacific Flyway during the 2009 hunting season. This strategy is based on the 2009 regulatory alternatives, current (1990–2008) population models and parameter estimates, and an objective to maximize longterm cumulative harvest subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes. The shaded cell indicates the regulatory prescription for 2009. Alaska BPOPb CAOR BPOPb 0 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 ≥ 0.55 0 C C M L L L L L L L L L 0.05 C C R R R M M M L L L L 0.10 C R R M M M L L L L L L 0.15 R R M M L L L L L L L L 0.20 M R L L L L L L L L L L 0.25 M R L L L L L L L L L L 0.30 M M L L L L L L L L L L 0.35 M M L L L L L L L L L L 0.40 M M L L L L L L L L L L ≥ 0.45 L M L L L L L L L L L L a C = closed season, R = restrictive, M = moderate, and L = liberal. b Estimated number of mallards in millions for Alaska (WBPHS strata 1–12) and in California and Oregon.Application of AHM Concepts to Other Stocks The USFWS is striving to apply the principles and tools of AHM to improve decisionmaking for several other stocks of waterfowl. Over the last year, some progress has been made to develop AHM frameworks for American black ducks (Anas rubripes) and the Atlantic Population of Canada geese (Branta canadensis), but these results are not yet finalized for inclusion in this year’s report. As these frameworks are developed further, we will continue to describe this work in subsequent reports. Below, we provide the 2009 updates for two decisionmaking frameworks that are currently informing harvest management. Northern Pintails The Flyway Councils have long identified the northern pintail as a highpriority species for inclusion in the AHM process. In 1997, the USFWS adopted a pintail harvest strategy to help align harvest opportunity with population status, while providing a foundation upon which to develop a formal AHM framework. Since 1997, the harvest strategy has undergone a number of technical improvements and policy revisions. However, the strategy continues to be a set of regulatory prescriptions born out of consensus, rather than an optimal strategy derived from agreedupon population models, management objectives, regulatory alternatives, and measures of uncertainty. In 2007, the USFWS and Flyway Councils took a major step towards a truly adaptive approach by incorporating alternative models of population dynamics. Two models are considered: one in which harvest is additive to natural mortality, and another in which harvest losses are compensated for by reductions in natural mortality. In the additive model, winter survival rate is a constant, whereas winter survival is densitydependent in the compensatory model. We here provide a summary of these recent modeling efforts. A detailed progress report is available online (http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/NOPI%20Harvest%20Strategy%202007.pdf). The predicted in year t + 1 () for the additive harvest mortality model is calculated as tcBPOP1+tcBPOP {}wttRsttscHRscBPOPcBPOP)1/(ˆ)ˆ1(1−−+=+γ where is the latitudeadjusted breeding population size in year t, and are the summer and winter survival rates, respectively, tcBPOPsswsRγ is a biascorrection constant for the age ratio, c is the crippling loss rate, is the predicted age ratio, and is the predicted continental harvest. The model uses the following constants: = 0.70, = 0.93, tRˆstHˆswsRγ = 0.80, and c = 0.20. The compensatory harvest mortality model serves as a hypothesis that stands in contrast to the additive harvest mortality model, positing a strong but realistic degree of compensation. The compensatory model assumes that the mechanism for compensation is densitydependent postharvest (winter) survival. The form is a logistic relationship between winter survival and postharvest population size, with the relationship anchored around the historic mean values for each variable. For the compensatory model then, predicted winter survival rate in year t () is calculated as ts[]1))((0101)(−−+−+−+=PPbattessss, 21 where (upper asymptote) is 1.0, (lower asymptote) is 0.7, b (slope term) is 1.0, is the postharvest population size in year t (expressed in millions),1s0stPP is the mean postharvest population size (4.295 million from 1974 through 2005), and a = logit010ssss⎛⎞−⎜⎟−⎝⎠ or ⎭⎬⎫⎩⎨⎧⎟⎟⎠⎞⎜⎜⎝⎛−−−−⎟⎟⎠⎞⎜⎜⎝⎛−−=0100101loglogssssssssa, where s is 0.93 (mean winter survival rate). At moderate population size and latitude, the compensatory model allows for greater harvest (Fig. 7) than does the additive model (note especially that the size of the restrictive region [seasonwithinaseason] is smaller and is invoked when the latitude is higher). Also, 2 and 3bird bag limits are called for under more circumstances. But, at high population sizes, the higher bag limits are called for less often, because the compensatory model predicts that growth of the population will be slower (densitydependence). The fit to historic data was used to compare the additive and compensatory harvest models. From the , , and observed harvest () for the period 1974–through year t, the subsequent year’s breeding population size (on the latitudeadjusted scale) was predicted with both the additive and compensatory model, and compared to the observed breeding population size (on the latitudeadjusted scale). The meansquared error of the predictions from the additive model () was calculated as: tcBPOPtmLATtHaddMSE Σ=−+−=ttaddttaddcBPOPcBPOPtMSE19752)(1)1975(1 and the meansquared error of the predictions from the compensatory model were calculated in a similar manner. The model weights for the additive and compensatory models were calculated from their relative meansquared errors. The model weight for the additive model () was calculated as: addW compaddaddaddMSEMSEMSEW111+=. The model weight for the compensatory model was found in a corresponding manner, or by subtracting the additive model weight from 1.0. As of 2008, the compensatory model did not fit the historic data as well as the additive model; the model weights were 0.597 for the additive model and 0.403 for the compensatory model, unchanged from 2007. The 2008 average model calls for a strategy that is intermediate between the additive and compensatory models (see Fig. 7). The USFWS remains committed to the development of an AHM framework to inform pintail harvest management based on a formal, derived strategy and clearly articulated management objectives. We are committed to continue to work with the Flyways to finalize an AHM protocol with the intent of implementing a revised harvest strategy for pintails in the near future. 22 23 1234567515253545556575859Average Latitude of the BPOPClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L21234567515253545556575859Average Latitude of the BPOPClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L21234567515253545556575859Average Latitude of the BPOPPintail BPOP (millions)ClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L2(A) Additive Model(B) Compensatory Model(C) Fig. 7. Statedependent harvest policy for northern pintails with (A) additive, (B) compensatory, and (C) 2008 weighted models. In each case the strategy assumes that the general duck hunting season is that prescribed under the liberal regulatory alternative. 2006 200824 Scaup In 2008, the USFWS implemented a decisionmaking framework for establishing scaup harvest regulations that was initially proposed in 2007 (Boomer and Johnson 2007). In addition, the USFWS committed to continue working with the Flyways to develop an alternative scaup population model for inclusion in the current decisionmaking framework. This model will capture the belief that the scaup population will decline to and stabilize at some lower equilibrium level in response to a declining carrying capacity and that harvest at current levels is completely compensatory. We plan to report on our efforts to develop an alternative model at the 2009 AHM Working Group meeting. In 2007, the USFWS also outlined methods to facilitate the specification of regulatory alternatives for scaup harvest management (Boomer et al. 2007). We proposed harvest thresholds to be considered under regulatory alternatives based on a simulation of an optimal policy that was derived under an objective to achieve 95% of the longterm cumulative harvest. We used this objective because it results in a strategy less sensitive to small change in population size as compared to a strategy derived under an objective to achieve 100% of longterm cumulative harvest. In addition, the 95% objective allows for some harvest opportunity at relatively low population sizes. We have continued to work with the Flyways to determine what regulations would achieve the allowable harvest thresholds set forth in Boomer et al. (2007). In 2008 during deliberations over scaup regulatory alternatives, the USFWS also agreed to consider a “hybrid season” option that would be available to all Flyways for the restrictive and moderate alternatives. In 2008, initial Restrictive, Moderate, and Liberal scaup regulatory alternatives were defined and implemented in all four Flyways. Subsequent concerns and dialogue led the USFWS to further clarify criteria associated with the establishment of hybrid seasons and to allow additional modifications of the alternatives for each Flyway in 2009. Final scaup regulatory alternatives were adopted for each Flyway in 2009. These alternatives will remain in place for a period of three years and then revisited as the latest harvest information is evaluated. The USFWS will continue to work with the Flyways to determine acceptable harvestmanagement objectives and evaluate regulatory alternatives to be used in the evolving decisionmaking framework for scaup harvest management. Presently, the scaup harvest strategy prescribes optimal harvest levels, not optimal regulatory alternatives. It is important to note that we currently have limited ability to predict expected scaup harvest under the newlyestablished, Flywayspecific scaup regulatory alternatives. The initial regulatory alternatives adopted for each Flyway were based on relatively crude predictions from harvest models developed in Boomer et al. (2007) or alternative harvest models proposed by the Flyways. As we gain experience with scaup regulatory alternatives, we will refine predicted harvest distributions associated with the Flywayspecific alternatives with the ultimate goal being to use regulatory alternatives, as opposed to harvest, as the control variable in deriving future scaup harvest policies. The lack of scaup demographic information over a sufficient timeframe and at a continental scale precludes the use of a traditional balance equation to represent scaup population and harvest dynamics. As a result, we used a discretetime, stochastic, logisticgrowth population model to represent changes in scaup abundance, while explicitly accounting for scaling issues associated with the monitoring data. Details describing the modeling and assessment framework that has been developed for scaup can be found in Appendix 6 and in Boomer and Johnson (2007). We updated the scaup assessment based on the current model formulation and data extending from 1974 through 2008. As in past analyses, the state space formulation and Bayesian analysis framework provided reasonable fits to the observed breeding population and total harvest estimates with realistic measures of variation. The posterior mean estimate of the intrinsic rate of increase (r) is 0.106 while the posterior mean estimate of the carrying 25 capacity (K) is 8.44 million birds. The posterior mean estimate of the scaling parameter (q) is 0.552, ranging between 0.479 and 0.634 with 95% probability. We calculated an optimal harvest policy for scaup based on: (1) a control variable of total harvest (U.S. and Canada combined), (2) current population model and updated parameter estimates, and (3) an objective to achieve 95% of the longterm cumulative harvest. We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model parameters accurately reflect population dynamics), breedingpopulation size would be expected to average 4.64 million (SD = 0.86 million). With an estimated breeding population size of 4.2 million scaup, the optimal harvest level for scaup is 0.3 million (Table 11). Based on the harvest thresholds specified in Boomer et al. (2007), this year’s optimal harvest corresponds to the moderate regulatory alternative. Table 11. Optimal scaup harvest levels (observed scale in millions) and corresponding breeding population sizes (in millions). This strategy is based on the current scaup population model, and an objective to maximize 95% of longterm cumulative harvest. The shaded cell indicates the optimal harvest level for 2009. BPOP Optimal Harvest 0.0 – 1.8 0 2.0 – 2.4 0.05 2.6 – 2.8 0.1 3.0 – 3.2 0.15 3.4 – 3.6 0.2 3.8 – 4.0 0.25 4.2 – 4.4 0.3 4.6 0.35 4.8 – 5.2 0.4 5.4 – 5.6 0.5 26 Emerging Issues in AHM Under the existing AHM protocol learning occurs passively as annual comparisons of model predictions to observations from monitoring programs are used to update model weights and relative beliefs about system responses to management (Johnson et al. 2002b) or as model parameters are updated based on an assessment of the most recent monitoring data (Boomer and Johnson 2007, Johnson et al. 2007). However, additional learning can also occur as decisionmaking frameworks are evaluated to determine if objectives are being achieved, if objectives are adequately specified or have changed, or if other aspects of the decision problem are adequately being addressed. Often the feedback resulting from this process results in a form of “double loop” learning (Lee 1993) that offers the opportunity to adapt decisionmaking frameworks in response to a shifting decision context, novel or emerging management alternatives, or to a need to revise models that may perform poorly or may need to account for new information. Adaptive management depends on this iterative process to ensure that decisionmaking protocols remain relevant in evolving natural and social systems. As a byproduct of the adaptive management process, it is natural to think about when or how decisionmaking protocols should be revised to incorporate new information or to accommodate changes in the overall management context. Recent outcomes from the 2008 Future of Waterfowl Management Workshop and evaluations of waterfowl harvest and habitat management programs (e.g., Anderson et al. 2007, Assessment Steering Committee 2007) suggest compelling reasons for a reevaluation of the objectives of waterfowl management as well as the technical and institutional frameworks through which harvest and habitat management decisions are made. We view such introspection as completely consistent with the principles underlying AHM, and, in fact, critical to its longterm utility as a decision framework for waterfowl harvest regulation. The AHM Working Group has begun to consider the need to better integrate harvest and habitat management objectives (sensu Runge et al. 2005, Anderson et al. 2007). In addition, we have begun to explore the implications of largescale system changes (e.g., related to climate change) in habitats and the environment on the decisionmaking frameworks used to inform harvest management. Additionally, we continue to explore questions related to the appropriate taxonomic, spatial, and temporal resolution of waterfowl harvest management. In light of all these issues, we remain committed to the continued reassessment of waterfowl harvest decision frameworks. Ultimately, we will need to prioritize the scope of work required to revisit the AHM protocol in the face of limited resources. We look forward to working with the Flyways as we continue to refine the adaptive harvest management framework. 27 LITERATURE CITED Anderson, D. R., and K. P. Burnham. 1976. Population ecology of the mallard. VI. The effect of exploitation on survival. U.S. Fish and Wildlife Service Resource Publication No. 128. 66pp. Anderson, M. G., D. Caswell, J. M. Eadie, J. T. Herbert, M. Huang, D. D. Humburg, F. A. Johnson, M. D. Koneff, S. E. Mott, T. D. Nudds, E. T. Reed, J. K. Ringleman, M. C. Runge, and B. C. Wilson. 2007. Report from the Joint Task Group for clarifying North American waterfowl management plan population objectives and their use in harvest management. Unpublished report. Assessment Steering Committee. 2007. North American waterfowl management plan: continental progress report. Unpublished final report. Blohm, R. J. 1989. Introduction to harvest  understanding surveys and season setting. Proceedings of the International Waterfowl Symposium 6:118–133. Blohm, R. J., R. E. Reynolds, J. P. Bladen, J. D. Nichols, J. E. Hines, K. P. Pollock, and R. T. Eberhardt. 1987. Mallard mortality rates on key breeding and wintering areas. Transactions of the North American Wildlife and Natural Resources Conference 52:246–263. Boomer, G. S., and F. A. Johnson. 2007. A proposed assessment and decisionmaking framework to inform scaup harvest management. Unpublished Report. U. S. Fish and Wildlife Service, Laurel, MD. 26pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/SCAUP2007Report.pdf Boomer, G. S., F. A. Johnson, M. D. Koneff, T. A. Sanders, and R. E. Trost. 2007. A process to determine scaup regulatory alternatives. Unpublished Scoping Document. U. S. Fish and Wildlife Service, Laurel, MD. 20pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/scaup_regs_scoping_draftVI.pdf Brooks, S. P., and A. Gelman. 1998 Alternative methods for monitoring convergence of iterative simulations. Journal of Computational and Graphical Statistics 7:434–455. Burnham, K. P., G. C. White, and D. R. Anderson. 1984. Estimating the effect of hunting on annual survival rates of adult mallards. Journal of Wildlife Management 48:350–361. Henny, C. J., and Burnham K. P. 1976. A reward band study of mallards to estimate reporting rates. Journal of Wildlife Management. 40:1–14. Heusman, H W, and J. R. Sauer. 2000. The northeastern states’ waterfowl breeding population survey. Wildlife Society Bulletin 28:355–364. Hodges, J. L., J. G. King, B. Conant, and H. A. Hanson. 1996. Aerial surveys of waterbirds in Alaska 1975–94: population trends and observer variability. Information and Technology Report 4, National Biological Service, U.S. Dept. of the Interior, Washington, D.C. 25pp. Johnson, F. A. 2003. Population dynamics of ducks other than mallards in midcontinent North America. Draft. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 15pp. 28 Johnson, F. A., G. S. Boomer, and T. A. Sanders. 2007. A proposed protocol for the adaptive harvest management of mallards breeding in Western North America. Unpublished Report. U. S. Fish and Wildlife Service, Laurel, MD. 33pp. Johnson, F. A., J. A. Dubovsky, M. C. Runge, and D. R. Eggeman. 2002a. A revised protocol for the adaptive harvest management of eastern mallards. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 13pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/emalahm2002.pdf. Johnson, F. A., W. L. Kendall, and J. A. Dubovsky. 2002b. Conditions and limitations on learning in the adaptive management of mallard harvests. Wildlife Society Bulletin 30:176–185. Johnson, F. A., C. T. Moore, W. L. Kendall, J. A. Dubovsky, D. F. Caithamer, J. R. Kelley, Jr., and B. K. Williams. 1997. Uncertainty and the management of mallard harvests. Journal of Wildlife Management 61:202–216. Johnson, F. A., and B. K. Williams. 1999. Protocol and practice in the adaptive management of waterfowl harvests. Conservation Ecology 3(1): 8. [online] URL: http://www.consecol.org/vol3/iss1/art8. Johnson, F. A., B. K. Williams, J. D. Nichols, J. E. Hines, W. L. Kendall, G. W. Smith, and D. F. Caithamer. 1993. Developing an adaptive management strategy for harvesting waterfowl in North America. Transactions of the North American Wildlife and Natural Resources Conference 58:565–583. Johnson, F. A., B. K. Williams, and P. R. Schmidt. 1996. Adaptive decisionmaking in waterfowl harvest and habitat management. Proceedings of the International Waterfowl Symposium 7:26–33. Lubow, B. C. 1995. SDP: Generalized software for solving stochastic dynamic optimization problems. Wildlife Society Bulletin 23:738–742. Meyer, R., and R. B. Millar. 1999. BUGS in Bayesian stock assessments. Canadian Journal of Fisheries and Aquatic Sciences 56:1078–1086. Millar, R. B., and R. Meyer. 2000. Nonlinear state space modeling of fisheries biomass dynamics by using MetropolisHastings within Gibbs sampling. Applied Statistics 49: 327–342. Nichols, J. D., F. A. Johnson, and B. K. Williams. 1995a. Managing North American waterfowl in the face of uncertainty. Annual Review of Ecology and Systematics 26:177–199. Nichols, J. D., R. E. Reynolds, R. J. Blohm, R. E. Trost, J. E. Hines, and J. P. Bladen. 1995b. Geographic variation in band reporting rates for mallards based on reward banding. Journal of Wildlife Management 59:697–708. Runge, M. C., F. A. Johnson, J. A. Dubovsky, W. L. Kendall, J. Lawrence, and J. Gammonley. 2002. A revised protocol for the adaptive harvest management of midcontinent mallards. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 28pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/MCMrevise2002.pdf. Runge, M. C., F. A. Johnson, M. G. Anderson, M. D. Koneff, E. T. Reed, and S. E. Mott. 2005. The need for coherence between waterfowl harvest and habitat management. Wildlife Society Bulletin 34:1231–1237. Schaefer, M. B. 1954. Some aspects of the dynamics of populations important to the management of commercial marine fisheries. Bulletin of the InterAmerican Tropical Tuna Commission 1:25–56. 29 Smith, G. W. 1995. A critical review of the aerial and ground surveys of breeding waterfowl in North America. Biological Science Report 5, National Biological Service, U.S. Dept. of the Interior, Washington, D.C. 252pp. Spiegelhalter, D. J., A. Thomas, N. Best, and D. Lunn. 2003. WinBUGS 1.4 User manual. MRC Biostatistics Unit, Institute of Public Health, Cambridge, UK. U.S. Fish and Wildlife Service. 2000. Adaptive harvest management: 2000 duck hunting season. U.S. Dept. Interior, Washington. D.C. 43pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2000/ahm2000.pdf U.S. Fish and Wildlife Service. 2001. Frameworkdate extensions for duck hunting in the United States: projected impacts & coping with uncertainty. U.S. Dept. Interior, Washington, D.C. 8pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2001/ahm2001.PDF U.S. Fish and Wildlife Service. 2002. Adaptive harvest management: 2002 duck hunting season. U.S. Dept. Interior, Washington. D.C. 34pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/2002AHMreport.pdf Walters, C. J. 1986. Adaptive management of renewable resources. MacMillan Publ. Co., New York, N.Y. 374pp. Williams, B. K., and F. A. Johnson. 1995. Adaptive management and the regulation of waterfowl harvests. Wildlife Society Bulletin 23:430–436. Williams, B. K., F. A. Johnson, and K. Wilkins. 1996. Uncertainty and the adaptive management of waterfowl harvests. Journal of Wildlife Management 60:223–232. 30 APPENDIX 1: AHM Working Group (Note: This list includes only permanent members of the AHM Working Group. Not listed here are numerous persons from federal and state agencies that assist the Working Group on an adhoc basis.) Coordinator: Scott Boomer U.S. Fish & Wildlife Service 11510 American Holly Drive Laurel, Maryland 207084017 phone: 3014975684 fax: 3014975871 email: scott_boomer@fws.gov USFWS Representatives: Brad Bortner (Region 1) U.S. Fish and Wildlife Service 911 NE 11th Ave. Portland, OR 972324181 phone: 5032316164 fax: 5032312364 email: brad_bortner@fws.gov Dave Case (contractor) D.J. Case & Associates 607 Lincolnway West Mishawaka, IN 46544 phone: 5742580100 fax: 5742580189 email: dave@djcase.com Jim Dubovsky (Region 6) U.S. Fish and Wildlife Service P.O. Box 25486DFC Denver, CO 802250486 phone: 3032364403 fax: 3032368680 email:james_dubovsky@fws.gov Jeff Haskins (Region 2) U.S. Fish and Wildlife Service P.O. Box 1306 Albuquerque, NM 87103 phone: 5052486827 (ext 30) fax: 5052487885 email: jeff_haskins@fws.gov Diane Pence (Region 5) Jim Kelley (Region 9) U.S. Fish and Wildlife Service 1 Federal Drive Fort Snelling, MN 551110458 phone: 6127135409 fax: 6127135393 email: james_r_kelley@fws.gov Sean Kelly (Region 3) U.S. Fish and Wildlife Service 1 Federal Drive Ft. Snelling, MN 551114056 phone: 6127135470 fax: 6127135393 email: sean_kelly@fws.gov Mark Koneff (Region 9) U.S. Fish & Wildlife Service 11510 American Holly Drive Laurel, Maryland 207084017 phone: 3014975648 fax: 3014975871 email: mark_koneff@fws.gov Paul Padding (Region 9) U.S. Fish and Wildlife Service 11510 American Holly Drive Laurel, MD 20708 phone: 3014975851 fax: 3014975885 email: paul_padding@fws.gov 31 U.S. Fish and Wildlife Service 300 Westgate Center Drive Hadley, MA 010359589 phone: 4132538577 fax: 4132538424 email: diane_pence@fws.gov Russ Oates (Region 7) U.S. Fish and Wildlife Service 1011 East Tudor Road Anchorage, AK 995036119 phone: 9077863446 fax: 9077863641 email: russ_oates@fws.gov Dave Sharp (Region 9) U.S. Fish and Wildlife Service P.O. Box 25486, DFC Denver, CO 802250486 phone: 3032752386 fax: 3032752384 email: dave_sharp@fws.gov Bob Trost (Region 9) U.S. Fish and Wildlife Service 911 NE 11th Ave. Portland, OR 972324181 phone: 5032316162 fax: 5032316228 email: robert_trost@fws.gov David Viker (Region 4) U.S. Fish and Wildlife Service 1875 Century Blvd., Suite 345 Atlanta, GA 30345 phone: 4046797188 fax: 4046797285 email: david_viker@fws.gov Canadian Wildlife Service Representatives: Dale Caswell Canadian Wildlife Service 123 Main St. Suite 150 Winnipeg, Manitoba, Canada R3C 4W2 phone: 2049835260 fax: 2049835248 email: dale.caswell@ec.gc.ca Eric Reed Canadian Wildlife Service 351 St. Joseph Boulevard Hull, QC K1A OH3, Canada phone: 8199530294 fax: 8199536283 email: eric.reed@ec.gc.ca Flyway Council Representatives: Min Huang (Atlantic Flyway) CT Dept. of Environmental Protection Franklin Wildlife Mgmt. Area 391 Route 32 North Franklin, CT 06254, USA Phone: 860/6426528 fax: 860/6427964 email: min.huang@po.state.ct.us Mike Johnson (Central Flyway) North Dakota Game and Fish Department 100 North Bismarck Expressway Bismarck, ND 585015095 phone: 7013286319 fax: 7013286352 email: mjohnson@state.nd.us Bryan Swift (Atlantic Flyway) Dept. Environmental Conservation Larry Reynolds (Mississippi Flyway) LA Dept. of Wildlife & Fisheries PO Box 98000 Baton Rouge, LA 708989000, USA Phone: 225/7650456 Fax: 225/7635456 email: lreynolds@wlf.state.la.us Jon Runge (Pacific Flyway) Colorado Division of Wildlife 317 W. Prospect Fort Collins, CO 80526 Phone: 9704724365 email: Jon.Runge@state.co.us Dan Yparraguirre (Pacific Flyway) California Dept. of Fish and Game 32 625 Broadway Albany, NY 122334754 phone: 5184028866 fax: 5184029027 or 4028925 email: blswift@gw.dec.state.ny.us Mark Vrtiska (Central Flyway) Nebraska Game and Parks Commission P.O. Box 30370 2200 North 33rd Street Lincoln, NE 685031417 phone: 4024715437 fax: 4024715528 email: mvrtiska@ngpc.state.ne.us 1812 Ninth Street Sacramento, CA 95814 phone: 9164453685 email: dyparraguirre@dfg.ca.gov Guy Zenner (Mississippi Flyway) Iowa Dept. of Natural Resources 1203 North Shore Drive Clear Lake, IA 50428 phone: 5153573517, ext. 23 fax: 5153575523 email: gzenner@netins.net USGS Technical Consultants: Fred Johnson Florida Integrated Science Center U. S. Geological Survey P.O. Box 110485 Gainesville, FL 32611 phone: 3523925075 fax: 3528460841 email: fjohnson@usgs.gov Mike Runge Patuxent Wildlife Research Center U. S. Geological Survey 12100 Beech Forest Rd. Laurel, MD 20708 phone: 3014975748 fax: 3014975545 email: mrunge@usgs.gov Andy Royle Patuxent Wildlife Research Center U. S. Geological Survey 12100 Beech Forest Rd. Laurel, MD 20708 phone: 3014975846 fax: 3014975545 email: aroyle@usgs.gov APPENDIX 2: Midcontinent Mallard Models In 1995, we developed population models to predict changes in midcontinent mallards based on the traditional survey area which includes individuals from Alaska (Johnson et al. 1997). In 1997, we added mallards from the Great Lakes region (Michigan, Minnesota, and Wisconsin) to the midcontinent mallard stock, assuming their population dynamics were equivalent. In 2002, we made extensive revisions to the set of alternative models describing the population dynamics of midcontinent mallards (Runge et al. 2002, USFWS 2002). In 2008, we redefined the population of midcontinent mallards to account for the removal of Alaskan birds (WBPHS strata 1–12) that are now considered to be in the western mallard stock and have subsequently rescaled the model set appropriately. Model Structure Collectively, the models express uncertainty (or disagreement) about whether harvest is an additive or compensatory form of mortality (Burnham et al. 1984), and whether the reproductive process is weakly or strongly densitydependent (i.e., the degree to which reproductive rates decline with increasing population size). All population models for midcontinent mallards share a common “balance equation” to predict changes in breedingpopulation size as a function of annual survival and reproductive rates: ()()()()NNmSmSRSStttAMtAFttJFtJMFsumMsum+=+−++11,,,,φφ where: N = breeding population size, m = proportion of males in the breeding population, SAM, SAF, SJF, and SJM = survival rates of adult males, adult females, young females, and young males, respectively, R = reproductive rate, defined as the fall age ratio of females, φφFsumMsum= the ratio of female (F) to male (M) summer survival, and t = year. We assumed that m and φφFsumMsumare fixed and known. We also assumed, based in part on information provided by Blohm et al. (1987), the ratio of female to male summer survival was equivalent to the ratio of annual survival rates in the absence of harvest. Based on this assumption, we estimated φφFsumMsum = 0.897. To estimate m we expressed the balance equation in matrix form: NNSRSSRSNNtAMtAFAMJMFsumMsumAFJFtAMtAF++⎡⎣⎢⎤⎦⎥=+⎡⎣⎢⎤⎦⎥⎡⎣⎢⎤⎦⎥110,,,,φφ and substituted the constant ratio of summer survival and means of estimated survival and reproductive rates. The right eigenvector of the transition matrix is the stable sex structure that the breeding population eventually would attain with these constant demographic rates. This eigenvector yielded an estimate of m = 0.5246. Using estimates of annual survival and reproductive rates, the balance equation for midcontinent mallards overpredicted observed population sizes by 11.0% on average. The source of the bias is unknown, so we modified the balance equation to eliminate the bias by adjusting both survival and reproductive rates: 33 ()()()()NNmSmSRSStSttAMtAFRttJFtJMFsumMsum+=+−++11γγ,,,, where γ denotes the biascorrection factors for survival (S) and reproduction (R). We used a least squares approach to estimate γS = 0.9407 and γR = 0.8647. Survival Process We considered two alternative hypotheses for the relationship between annual survival and harvest rates. For both models, we assumed that survival in the absence of harvest was the same for adults and young of the same sex. In the model where harvest mortality is additive to natural mortality: ()SsKtsexagesexAtsexage,,,,,=−01 and in the model where changes in natural mortality compensate for harvest losses (up to some threshold): SsifKKifKstsexagesexCtsexagesexCtsexagetsexagesexC,,,,,,,,,,=≤−−>⎧⎨⎪⎩⎪000111 where s0 = survival in the absence of harvest under the additive (A) or compensatory (C) model, and K = harvest rate adjusted for crippling loss (20%, Anderson and Burnham 1976). We averaged estimates of s0 across banding reference areas by weighting by breedingpopulation size. For the additive model, s0 = 0.7896 and 0.6886 for males and females, respectively. For the compensatory model, s0 = 0.6467 and 0.5965 for males and females, respectively. These estimates may seem counterintuitive because survival in the absence of harvest should be the same for both models. However, estimating a common (but still sexspecific) s0 for both models leads to alternative models that do not fit available bandrecovery data equally well. More importantly, it suggests that the greatest uncertainty about survival rates is when harvest rate is within the realm of experience. By allowing s0 to differ between additive and compensatory models, we acknowledge that the greatest uncertainty about survival rate is its value in the absence of harvest (i.e., where we have no experience). Reproductive Process Annual reproductive rates were estimated from age ratios in the harvest of females, corrected using a constant estimate of differential vulnerability. Predictor variables were the number of ponds in May in Prairie Canada (P, in millions) and the size of the breeding population (N, in millions). We estimated the bestfitting linear model, and then calculated the 80% confidence ellipsoid for all model parameters. We chose the two points on this ellipsoid with the largest and smallest values for the effect of breedingpopulation size, and generated a weakly densitydependent model: RPtt=+−071660108300373... and a strongly densitydependent model: RPtt= +−113900137601131... Predicted recruitment was then rescaled to reflect the current definition of midcontinent mallards which now excludes birds from Alaska but includes mallards observed in the Great Lakes region. 34 Pond Dynamics We modeled annual variation in Canadian pond numbers as a firstorder autoregressive process. The estimated model was: 35 t PPtt+=++12212703420..ε where ponds are in millions and εtis normally distributed with mean = 0 and variance = 1.2567. Variance of Prediction Errors Using the balance equation and submodels described above, predictions of breedingpopulation size in year t+1 depend only on specification of population size, pond numbers, and harvest rate in year t. For the period in which comparisons were possible, we compared these predictions with observed population sizes. We estimated the predictionerror variance by setting: ()()()()()[]()eNNeNNNnttobstprettobstpret=−=−Σlnln~,$lnlnthen assumingand estimating012σσ22 where Nobs and Npre are observed and predicted population sizes (in millions), respectively, and n = the number of years being compared. We were concerned about a variance estimate that was too small, either by chance or because the number of years in which comparisons were possible was small. Therefore, we calculated the upper 80% confidence limit for σ2 based on a Chisquared distribution for each combination of the alternative survival and reproductive submodels, and then averaged them. The final estimate of σ2 was 0.0280, equivalent to a coefficient of variation of about 18%. Model Implications The population model with additive hunting mortality and weakly densitydependent recruitment (SaRw) leads to the most conservative harvest strategy, whereas the model with compensatory hunting mortality and strongly densitydependent recruitment (ScRs) leads to the most liberal strategy. The other two models (SaRs and ScRw) lead to strategies that are intermediate between these extremes. Under the models with compensatory hunting mortality (ScRs and ScRw), the optimal strategy is to have a liberal regulation regardless of population size or number of ponds because at harvest rates achieved under the liberal alternative, harvest has no effect on population size. Under the strongly densitydependent model (ScRs), the density dependence regulates the population and keeps it within narrow bounds. Under the weakly density dependent model (ScRw), the densitydependence does not exert as strong a regulatory effect, and the population size fluctuates more. Model Weights Model weights are calculated as Bayesian probabilities, reflecting the relative ability of the individual alternative models to predict observed changes in population size. The Bayesian probability for each model is a function of the model’s previous (or prior) weight and the likelihood of the observed population size under that model. We used Bayes’ theorem to calculate model weights from a comparison of predicted and observed population sizes for the years 1996–2009, starting with equal model weights in 1995. APPENDIX 3: Eastern Mallard Models Model Structure We also revised the population models for eastern mallards in 2002 (Johnson et al. 2002a, USFWS 2002). The current set of six models: (1) relies solely on federal and state waterfowl surveys (rather than the Breeding Bird Survey) to estimate abundance; (2) allows for the possibility of a positive bias in estimates of survival or reproductive rates; (3) incorporates competing hypotheses of strongly and weakly densitydependent reproduction; and (4) assumes that hunting mortality is additive to other sources of mortality. As with midcontinent mallards, all population models for eastern mallards share a common balance equation to predict changes in breedingpopulation size as a function of annual survival and reproductive rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅+−⋅+⋅⋅+⋅⋅⋅11ψ where: N = breedingpopulation size, p = proportion of males in the breeding population, Sam, Saf, Sym, and Syf = survival rates of adult males, adult females, young males, and young females, respectively, Am = ratio of young males to adult males in the harvest, d = ratio of young male to adult male direct recovery rates, ψ = the ratio of male to female summer survival, and t = year. In this balance equation, we assume that p, d, and ψ are fixed and known. The parameter ψ is necessary to account for the difference in anniversary date between the breedingpopulation survey (May) and the survival and reproductive rate estimates (August). This model also assumes that the sex ratio of fledged young is 1:1; hence Am/d appears twice in the balance equation. We estimated d = 1.043 as the median ratio of young:adult male bandrecovery rates in those states from which wing receipts were obtained. We estimated ψ = 1.216 by regressing through the origin estimates of male survival against female survival in the absence of harvest, assuming that differences in natural mortality between males and females occur principally in summer. To estimate p, we used a population projection matrix of the form: ()()⎥⎦⎤⎢⎣⎡⋅⎥⎦⎤⎢⎣⎡⋅⋅⋅+=⎥⎦⎤⎢⎣⎡++ttafyfmymmamttFMSSdASdASFMψ011 where M and F are the relative number of males and females in the breeding populations, respectively. To parameterize the projection matrix we used average annual survival rate and age ratio estimates, and the estimates of d and ψ provided above. The right eigenvector of the projection matrix is the stable proportion of males and females the breeding population eventually would attain in the face of constant demographic rates. This eigenvector yielded an estimate of p = 0.544. We also attempted to determine whether estimates of survival and reproductive rates were unbiased. We relied on the balance equation provided above, except that we included additional parameters to correct for any bias that might exist. Because we were unsure of the source(s) of potential bias, we alternatively assumed that any bias resided solely in survival rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅⋅+−⋅+⋅⋅+⋅⋅⋅11Ωψ 36 (where Ω is the biascorrection factor for survival rates), or solely in reproductive rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅+−⋅+⋅⋅⋅+⋅⋅⋅⋅11αα (where α is the biascorrection factor for reproductive rates). We estimated Ω and α by determining the values of these parameters that minimized the sum of squared differences between observed and predicted population sizes. Based on this analysis, Ω = 0.836 and α = 0.701, suggesting a positive bias in survival or reproductive rates. However, because of the limited number of years available for comparing observed and predicted population sizes, we also retained the balance equation that assumes estimates of survival and reproductive rates are unbiased. Survival Process For purposes of AHM, annual survival rates must be predicted based on the specification of regulationspecific harvest rates (and perhaps on other uncontrolled factors). Annual survival for each age (i) and sex (j) class under a given regulatory alternative is: ()()Shvctijjtamij,,=⋅−⋅−⎛⎝⎜⎜⎞⎠⎟⎟θ11 where: S = annual survival, jθ= mean survival from natural causes, ham = harvest rate of adult males, v = harvest vulnerability relative to adult males, and c = rate of crippling (unretrieved harvest). This model assumes that annual variation in survival is due solely to variation in harvest rates, that relative harvest vulnerability of the different agesex classes is fixed and known, and that survival from natural causes is fixed at its sample mean. We estimated jθ= 0.7307 and 0.5950 for males and females, respectively. Reproductive process As with survival, annual reproductive rates must be predicted in advance of setting regulations. We relied on the apparent relationship between breedingpopulation size and reproductive rates: ()RabNtt=⋅⋅exp where Rt is the reproductive rate (i.e., dAmt), Nt is breedingpopulation size in millions, and a and b are model parameters. The leastsquares parameter estimates were a = 2.508 and b = 0.875. Because of both the importance and uncertainty of the relationship between population size and reproduction, we specified two alternative models in which the slope (b) was fixed at the leastsquares estimate ± one standard error, and in which the intercepts (a) were subsequently reestimated. This provided alternative hypotheses of strongly densitydependent (a = 4.154, b = 1.377) and weakly densitydependent reproduction (a = 1.518, b = 0.373). 37 Variance of Prediction Errors Using the balance equations and submodels provided above, predictions of breedingpopulation size in year t+1 depend only on the specification of a regulatory alternative and on an estimate of population size in year t. For the period in which comparisons were possible (1991–96), we were interested in how well these predictions corresponded with observed population sizes. In making these comparisons, we were primarily concerned with how well the biascorrected balance equations and reproductive and survival submodels performed. Therefore, we relied on estimates of harvest rates rather than regulations as model inputs. We estimated the predictionerror variance by setting: ()()()()()[]eNNeNNNttobstprettobstpret=−=−Σlnln~,$lnlnthen assumingand estimating02σσ22 where Nobs and Npre are observed and predicted population sizes (in millions), respectively, and n = 6. Variance estimates were similar regardless of whether we assumed that the bias was in reproductive rates or in survival, or whether we assumed that reproduction was strongly or weakly densitydependent. Thus, we averaged variance estimates to provide a final estimate of σ2 = 0.006, which is equivalent to a coefficient of variation (CV) of 8.0%. We were concerned, however, about the small number of years available for estimating this variance. Therefore, we estimated an 80% confidence interval for σ2 based on a Chisquared distribution and used the upper limit for σ2 = 0.018 (i.e., CV = 14.5%) to express the additional uncertainty about the magnitude of prediction errors attributable to potentially important environmental effects not expressed by the models. Model Implications Modelspecific regulatory strategies based on the hypothesis of weakly densitydependent reproduction are considerably more conservative than those based on the hypothesis of strongly densitydependent reproduction. The three models with weakly densitydependent reproduction suggest a carrying capacity (i.e., average population size in the absence of harvest) >2.0 million mallards, and prescribe extremely restrictive regulations for population size <1.0 million. The three models with strongly densitydependent reproduction suggest a carrying capacity of about 1.5 million mallards, and prescribe liberal regulations for population sizes >300 thousand. Optimal regulatory strategies are relatively insensitive to whether models include a bias correction or not. All modelspecific regulatory strategies are “knifeedged,” meaning that large differences in the optimal regulatory choice can be precipitated by only small changes in breedingpopulation size. This result is at least partially due to the small differences in predicted harvest rates among the current regulatory alternatives (see the section on Regulatory Alternatives later in this report). Model Weights We used Bayes’ theorem to calculate model weights from a comparison of predicted and observed population sizes for the years 1996–2009. We calculated weights for the alternative models based on an assumption of equal model weights in 1996 (the last year data was used to develop most model components) and on estimates of yearspecific harvest rates (Appendix 5). 38 APPENDIX 4: Western Mallard Models In contrast to midcontinent and eastern mallards, we did not model changes in population size of western mallards as an explicit function of survival and reproductive rate estimates (which in turn may be functions of harvest and environmental covariates). We believed this socalled “balanceequation approach” was not viable for western mallards because of insufficient banding in Alaska to estimate survival rates, and because of the difficulty in estimating stockspecific fall age ratios from a sample of wings derived from a mix of breeding stocks. We therefore relied on a discrete logistic model (Schaefer 1954), which combines reproduction and natural mortality into a single parameter r, the intrinsic rate of growth. The model assumes densitydependent growth, which is regulated by the ratio of population size, N, to the carrying capacity of the environment, K (i.e., equilibrium population size in the absence of harvest). In the traditional formulation, harvest mortality is additive to other sources of mortality, but compensation for hunting losses can occur through subsequent increases in production. However, we parameterized the model in a way that also allows for compensation of harvest mortality between the hunting and breeding seasons. It is important to note that compensation modeled in this way is purely phenomenological, in the sense that there is no explicit ecological mechanism for compensation (e.g., densitydependent mortality after the hunting season). The basic model for both the Alaska and California–Oregon stocks had the form: ()AMttttttthdwhereKNrNNN⋅=−⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛−+=+αα111 and where t = year, hAM = the harvest rate of adult males, and d = a scaling factor. The scaling factor is used to account for a combination of unobservable effects, including unretrieved harvest (i.e., crippling loss), differential harvest mortality of cohorts other than adult males, and for the possibility that some harvest mortality may not affect subsequent breedingpopulation size (i.e., the compensatory mortality hypothesis). Estimation Framework We used Bayesian estimation methods in combination with a statespace model that accounts explicitly for both process and observation error in breeding population size. This combination of methods is becoming widely used in natural resource modeling, in part because it facilitates the fitting of nonlinear models that may have nonnormal errors (Meyer and Millar 1999). The Bayesian approach also provides a natural and intuitive way to portray uncertainty, allows one to incorporate prior information about model parameters, and permits the updating of parameter estimates as further information becomes available. We first scaled N by K as recommended by Meyer and Millar (1999), and assumed that process errors et were lognormally distributed with mean 0 and variance σ2. Thus, the process model had the form: ()()[](){}()21111,0~11loglogσNewhereehdPrPPPKNPttAMtttttttt+⋅−−+==−−−− The observation model related the unknown population sizes (PtK) to the population sizes (Nt) estimated from the breedingpopulation surveys in Alaska and CaliforniaOregon. We assumed that the observation process yielded additive, normally distributed errors, which were represented by: 39 BPOPtttKPNε+=, ),0(~2BPOPBPOPtNwhereσε. Use of the observation model allowed us to account for the sampling error in population estimates, while permitting us to estimate the process error, which reflects the inability of the model to completely describe changes in population size. The process error reflects the combined effect of misspecification of an appropriate model form, as well as any unmodeled environmental drivers. We initially examined a number of possible environmental covariates, including the Palmer Drought Index in California and Oregon, spring temperature in Alaska, and the El Niño Southern Oscillation Index (http://www.cdc.noaa.gov/people/klaus.wolter/MEI/mei.html). While the estimated effects of these covariates on r or K were generally what one would expect, they were never of sufficient magnitude to have a meaningful effect on optimal harvest strategies. We therefore chose not to further pursue an investigation of environmental covariates, and posited that the process error was a sufficient surrogate for these unmodeled effects. Parameterization of the models also required measures of harvest rate. Beginning in 2002, harvest rates of adult males were estimated directly from the recovery of reward bands. Prior to 1993, we used direct recoveries of standard bands, corrected for bandreporting rates provided by Nichols et al. (1995b). We also used the bandreporting rates provided by Nichols et al. (1995b) for estimating harvest rates in 1994 and 1995, except that we inflated the reporting rates of fulladdress and tollfree bands based on an unpublished analysis by Clint Moore and Jim Nichols (Patuxent Wildlife Research Center). We were unwilling to estimate harvest rates for the years 1996–2001 because of suspected, but unknown, increases in the reporting rates of all bands. For simplicity, harvest rate estimates were treated as known values in our analysis, although future analyses might benefit from an appropriate observation model for these data. In a Bayesian analysis, one is interested in making probabilistic statements about the model parameters (θ), conditioned on the observed data. Thus, we are interested in evaluating P(θdata), which requires the specification of prior distributions for all model parameters and unobserved system states (θ) and the sampling distribution (likelihood) of the observed data P(data θ). Using Bayes theorem, we can represent the posterior probability distribution of model parameters, conditioned on the data, as: )()()(θθθdataPPdataP×∝. Accordingly, we specified prior distributions for model parameters r, K, d, and P0, which is the initial population size relative to carrying capacity. For both stocks, we specified the following prior distributions for r, d, and σ2: ()()()001.0,001.0~2,0~4427.1,0397.1~2gammaInverseUniformdnormalLogr−−−σ The prior distribution for r is centered at 0.35, which we believe to be a reasonable value for mallards based on lifehistory characteristics and estimates for other avian species. Yet the distribution also admits considerable uncertainty as to the value of r within what we believe to be realistic biological bounds. As for the harvestrate scalar, we would expect d ≥ 1 under the additive hypothesis and d < 1 under the compensatory hypothesis. As we had no data to specify an informative prior distribution, we specified a vague prior in which d could take on a wide range of values with equal probability. We used a traditional, uninformative prior distribution for σ2. Prior distributions for K and P0 were stockspecific and are described in the following sections. 40 We used the publicdomain software WinBUGS (http://www.mrcbsu.cam.ac.uk/bugs/) to derive samples from the joint posterior distribution of model parameters via MarkovChain Monte Carlo (MCMC) simulations. We obtained 510,000 samples from the joint posterior distribution, discarded the first 10,000, and then thinned the remainder by 50, resulting in a final sample of 10,000. Alaska mallards Data selection.Breeding population estimates of mallards in Alaska (and the Old Crow Flats in Yukon) are available since 1955 in WBPHS strata 1–12 (Smith 1995). However, a change in survey aircraft in 1977 instantaneously increased the detectability of waterfowl, and thus population estimates (Hodges et al. 1996). Moreover, there was a rapid increase in average annual temperature in Alaska at the same time, apparently tied to changes in the frequency and intensity of El Niño events (http://www.cdc.noaa.gov/people/klaus.wolter/MEI/mei.html). This confounding of changes in climate and survey methods led us to truncate the years 1955–1977 from the time series of population estimates. Modeling of the Alaska stock also depended on the availability of harvestrate estimates derived from bandrecovery data. Unfortunately, sufficient numbers of mallards were not banded in Alaska prior to 1990. A search for covariates that would have allowed us to make harvestrate predictions for years in which bandrecovery data were not available was not fruitful, and we were thus forced to further restrict the time series to 1990–2005. Even so, harvest rate estimates were not available for the years 1996–2001 because of unknown changes in bandreporting rates. Because available estimates of harvest rate showed no apparent variation over time, we simply used the mean and standard deviation of the available estimates and generated independent samples of predictions for the missing years based on a logit transformation and an assumption of normality: ()200119960830.0,3265.2~1ln−=−⎟⎟⎠⎞⎜⎜⎝⎛−tforNormalhhtt Prior distributions for K and P0.—We believed that sufficient information was available to use mildly informative priors for K and P0. In recent years the Alaska stock has contained approximately 0.8 million mallards. If harvest rates have been comparable to that necessary to achieve maximum sustained yield (MSY) under the logistic model (i.e., r/2), then we would expect K ≈ 1.6 million. On the other hand, if harvest rates have been less than those associated with MSY, then we would expect K < 1.6 million. Because we believed it was not likely that harvest rates were >r/2, we believed the likely range of K to be 0.8–1.6 million. We therefore specified a prior distribution that had a mean of 1.4 million, but had a sufficiently large variance to admit a wide range of possible values: ()41224.0,13035.0~LognormalK Extending this line of reasoning, we specified a prior distribution that assumed the estimated population size of approximately 0.4 million at the start of the timeseries (i.e., 1990) was 20–60% of K. Thus on a log scale: ()5108.0,6094.1~−−UniformPo Parameter estimates.—The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed timeseries of population estimates. The posterior means of K and r were similar to their priors, although their variances were considerably smaller (Table 1). However, the posterior distribution of d was essentially the same as its prior, reflecting the absence of information in the data necessary to reliably estimate this parameter. Table 1. Estimates of model parameters resulting from fitting a discrete logistic model with MCMC to a time 41 series of estimated population sizes and harvest rates of mallards breeding in Alaska, 1990–2008. Parameter Mean SD 95% credibility interval K 1.125 0.323 0.659–1.890 P0 0.342 0.097 0.207–0.560 d 1.066 0.541 0.091–1.946 r 0.310 0.131 0.095–0.594 σ2 0.024 0.013 0.008–0.057 California–Oregon mallards Data selection.—Breedingpopulation estimates of mallards in California are available starting in 1992, but not until 1994 in Oregon. Also, Oregon did not conduct a survey in 2001. To avoid truncating the time series, we used the admittedly weak relationship (P = 0.11) between California and Oregon population estimates to predict population sizes in Oregon in 1992, 1993, and 2001. The fitted linear model was: ()CAtORtNN0848.074486+= To derive realistic standard errors, we assumed that the predictions had the same mean coefficient of variation as the years when surveys were conducted (n = 14, CV = 0.082). The estimated sizes and variances of the California–Oregon stock were calculated by simply summing the statespecific estimates. We pooled banding and recovery data for California and Oregon and estimated harvest rates in the same manner as that for Alaska mallards. Although banded sample sizes were sufficient in all years, harvest rates could not be estimated for the years 1996–2001 because of unknown changes in bandreporting rates. As with Alaska, available estimates of harvest rate showed no apparent trend over time, and we simply used the mean and standard deviation of the available estimates and generated independent samples of predictions for the missing years based on a logit transformation and an assumption of normality: ()200119960335.0,9633.1~1ln−=−⎟⎟⎠⎞⎜⎜⎝⎛−tforNormalhhtt Prior distributions for K and P0.— Unlike the Alaska stock, the California–Oregon population has been relatively stable with a mean of 0.48 million mallards. We believed K should be in the range 0.48 – 0.96 million, assuming the logistic model and that harvest rates were ≤ r/2. We therefore specified a prior distribution on K that had a mean of 0.7 million, but with a variance sufficiently large to admit a wide range of possible values: ()41224.0,5628.0~−LognormalK The estimated size of the California–Oregon stock was 0.48 million at the start of the timeseries (i.e., 1992). We used a similar line of reasoning as that for Alaska for specifying a prior distribution P0, positing that initial population size was 40–100% of K. Thus on a log scale: ()0.0,9163.0~−UniformPo Parameter estimates.—The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed time series of population estimates. The posterior means of K and r were similar to their priors, although the variances were considerably smaller (Table 2). Interestingly, the posterior mean of d was <1, 42 43 suggestive of a compensatory response to harvest; however the standard deviation of the estimate was large, with the upper 95% credibility limit >1. Table 2. Estimates of model parameters resulting from fitting a discrete logistic model with MCMC to a timeseries of estimated population sizes and harvest rates of mallards breeding in California and Oregon, 1992–2008. Parameter Mean SD 95% credibility interval K 0.655 0.178 0.445–1.117 P0 0.741 0.161 0.429–0.986 d 0.599 0.410 0.038–1.588 r 0.335 0.220 0.067–0.883 σ2 0.013 0.012 0.001–0.044 APPENDIX 5: Modeling Mallard Harvest Rates Midcontinent We modeled harvest rates of midcontinent mallards within a Bayesian hierarchical framework. We developed a set of models to predict harvest rates under each regulatory alternative as a function of the harvest rates observed under the liberal alternative, using historical information. We modeled the probability of regulationspecific harvest rates (h) based on normal distributions with the following parameterizations: ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LfLLMfLMMRLRRCCCNhpNhpNhpNhpυδμυδμγυμγυμ++ For the restrictive and moderate alternatives we introduced the parameter γ to represent the relative difference between the harvest rate observed under the liberal alternative and the moderate or restrictive alternatives. Based on this parameterization, we are making use of the information that has been gained (under the liberal alternative) and are modeling harvest rates for the restrictive and moderate alternatives as a function of the mean harvest rate observed under the liberal alternative. For the harvestrate distributions assumed under the restrictive and moderate regulatory alternatives, we specified that γR and γM are equal to the prior estimates of the predicted mean harvest rates under the restrictive and moderate alternatives divided by the prior estimates of the predicted mean harvest rates observed under the liberal alternative. Thus, these parameters act to scale the mean of the restrictive and moderate distributions in relation to the mean harvest rate observed under the liberal regulatory alternative. We also considered the marginal effect of frameworkdate extensions under the moderate and liberal alternatives by including the parameter δf. To update the probability distributions of harvest rates realized under each regulatory alternative, we first needed to specify a prior probability distribution for each of the model parameters. These distributions represent prior beliefs regarding the relationship between each regulatory alternative and the expected harvest rates. We used a normal distribution to represent the mean and a scaled inversechisquare distribution to represent the variance of the normal distribution of the likelihood. For the mean (μ) of each harvestrate distribution associated with each regulatory alternative, we use the predicted mean harvest rates provided in USFWS (2000:13–14), assuming uniformity of regulatory prescriptions across flyways. We set prior values of each standard deviation (ν) equal to 20% of the mean (CV = 0.2) based on an analysis by Johnson et al. (1997). We then specified the following prior distributions and parameter values under each regulatory package: Closed (in U.S. only): )0018.0,6(~)()60018.0,0088.0(~)(2222χνμ−InvScaledpNpCC These closedseason parameter values are based on observed harvest rates in Canada during the 1988–93 seasons, which was a period of restrictive regulations in both Canada and the United States. For the restrictive and moderate alternatives, we specified that the standard error of the normal distribution of the scaling parameter is based on a coefficient of variation for the mean equal to 0.3. The scale parameter of the inversechisquare distribution was set equal to the standard deviation of the harvest rate mean under the restrictive and moderate regulation alternatives (i.e., CV = 0.2). 44 Restrictive: )0133.0,6(~)()615.0,51.0(~)(2222χνγ−InvScaledpNpRR Moderate: )0222.0,6(~)()626.0,85.0(~)(2222χνγ−InvScaledpNpRM Liberal: )0261.0,6(~)()60261.0,1305.0(~)(2222χνμ−InvScaledpNpRL The prior distribution for the marginal effect of the frameworkdate extension was specified as: )01.0,02.0(~)(2Npfδ The prior distributions were multiplied by the likelihood functions based on the last seven years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided in Table 1. Table 1. Parameter estimates for predicting midcontinent mallard harvest rates resulting from a hierarchical, Bayesian analysis of midcontinent mallard banding and recovery information from 1998 to 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0088 0.0021 h1998 0.1021 0.0069 Cν 0.0019 0.0005 h1999 0.0981 0.0071 Rγ 0.5107 0.0622 h2000 0.1249 0.0083 Rν 0.0129 0.0033 h2001 0.0923 0.0088 Mγ 0.8525 0.1066 h2002 0.1128 0.0062 Mν 0.0215 0.0055 h2003 0.1131 0.0069 μL 0.1103 0.0067 h2004 0.1200 0.0109 Lν 0.0198 0.0035 h2005 0.1157 0.0085 fδ 0.0085 0.0075 h2006 0.1094 0.0078 h2007 0.1034 0.0066 h2008 0.1112 0.0065 45 Eastern We modeled harvest rates of eastern mallards using the same parameterizations as those for midcontinent mallards: ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LfLLMfLMMRLRRCCCNhpNhpNhpNhpυδμυδμγυμγυμ++ We set prior values of each standard deviation (ν) equal to 30% of the mean (CV = 0.3) to account for additional variation due to changes in regulations in the other Flyways and their unpredictable effects on the harvest rates of eastern mallards. We then specified the following prior distribution and parameter values for the liberal regulatory alternative: Closed (in U.S. only): )024.0,6(~)()6024.0,08.0(~)(2222χνμ−InvScaledpNpCC Restrictive: )0404.0,6(~)()6228.0,76.0(~)(2222χνγ−InvScaledpNpRR Moderate: )0488.0,6(~)()628.0,92.0(~)(2222χνγ−InvScaledpNpRM Liberal: )0531.0,6(~)()60531.0,1771.0(~)(2222χνμ−InvScaledpNpRL A previous analysis suggested that the effect of the frameworkdate extension on eastern mallards would be of lower magnitude and more variable than on midcontinent mallards (USFWS 2000). Therefore, we specified the following prior distribution for the marginal effect of the frameworkdate extension for eastern mallards as: )01.0,01.0(~)(2Npfδ 46 The prior distributions were multiplied by the likelihood functions based on the last four years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided in Table 2. Table 2. Parameter estimates for predicting eastern mallard harvest rates resulting from a hierarchical, Bayesian analysis of eastern mallard banding and recovery information from 2002 to 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0797 0.0257 h2002 0.1617 0.0125 Cν 0.0233 0.0061 h2003 0.1456 0.0105 Rγ 0.7619 0.0919 h2004 0.1360 0.0114 Rν 0.0393 0.0100 h2005 0.1308 0.0119 Mγ 0.9217 0.1145 h2006 0.1034 0.0132 Mν 0.0473 0.0118 h2007 0.1209 0.0134 μL 0.1459 0.0153 h2008 0.1206 0.0120 Lν 0.0433 0.0087 fδ 0.0034 0.0096 Western We modeled harvest rates of western mallards using a similar parameterization as that used for midcontinent and eastern mallards. However, we did not explicitly model the effect of the framework date extension because we did not use data observed prior to when framework date extensions were available. In the western mallard parameterization, the effect of the framework date extensions are implicit in the expected mean harvest rate expected under the liberal regulatory option. ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LLLMLMMRLRRCCCNhpNhpNhpNhpυμυμγυμγυμ We set prior values of each standard deviation (ν) equal to 30% of the mean (CV = 0.3) to account for additional variation due to changes in regulations in the other Flyways and their unpredictable effects on the harvest rates of western mallards. We then specified the following prior distribution and parameter values for the liberal regulatory alternative: Closed (in U.S. only): )00264.0,6(~)()600264.0,01.0(~)(2222χνμ−InvScaledpNpCC Restrictive: 47 )01683.0,6(~)()6153.0,51.0(~)(2222χνγ−InvScaledpNpRR Moderate: )02805.0,6(~)()6255.0,85.0(~)(2222χνγ−InvScaledpNpRM Liberal: )033.0,6(~)()6033.0,11.0(~)(2222χνμ−InvScaledpNpRL The prior distributions were multiplied by the likelihood functions based on the last four years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided Table 3. Table 3. Parameter estimates for predicting western mallard harvest rates resulting from a hierarchical, Bayesian analysis of western mallard banding and recovery information from 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0102 0.0187 h2008 0.1435 0.0243 Cν 0.0181 0.0045 Rγ 0.5103 0.0636 Rν 0.0164 0.0042 Mγ 0.8524 0.1025 Mν 0.0272 0.0069 μL 0.1155 0.0124 Lν 0.0318 0.0076 48 APPENDIX 6: Scaup Model We use a statespace formulation of scaup population and harvest dynamics within a Bayesian estimation framework (Meyer and Millar 1999, Millar and Meyer 2000). This analytical framework allows us to represent uncertainty associated with the monitoring programs (observation error) and the ability of our model formulation to predict actual changes in the system (process error). 8.1 Process Model Given a logistic growth population model that includes harvest (Schaefer 1954), scaup population and harvest dynamics are calculated as a function of the intrinsic rate of increase (r), carrying capacity (K), and harvest (Ht). Following Meyer and Millar (1999), we scaled population sizes by K (i.e., Pt = Nt/K) and assumed that process errors (εt) are lognormally distributed with a mean of 0 and variance. The state dynamics can be expressed as 2Processσ 197401974εePP= ()008,1975,...,2t,/)1(1111t=−−+=−−−−teKHPrPPPttttε where P0 is the initial ratio of population size to carrying capacity. To predict total scaup harvest levels, we modeled scaup harvest rates (ht) as a function of the pooled direct recovery rate (ft) observed each year with ./tttfhλ= We specified reporting rate (λt) distributions based on estimates for mallards (Anas platyrhynchos) from large scale historical and existing reward banding studies (Henny and Burham 1976, Nichols et. al. 1995b, P. Garrettson unpublished data). We accounted for increases in reporting rate believed to be associated with changes in band type (e.g., from AVISE and new address bands to 1800 toll free bands) by specifying year specific reporting rates according to 9961974,...,1 )04.0,38.0(~=tNormaltλ . 0081997,...,2 )04.0,70.0(~=tNormaltλ We then predicted total scaup harvest (Ht) with ()[].0081974,...,2t,1t=−+=KPrPPhHttt 8.2 Observation Model We compared our predictions of population and harvest numbers from our process model to the observations collected by the Waterfowl and Breeding Habitat Survey (WBPHS) and the Harvest Survey programs with the following relationships, assuming that the population and harvest observation errors were additive and normally distributed. May breeding population estimates were related to model predictions by ,2008,...,1974),,0(~where,2,==−tNKPNBPOPtBPOPtBPOPttObservedtσεε where is specified for each year with the variance estimates resulting from the WBPHS. 2,BPOPtσ 49 We adjusted our harvest predictions to the observed harvest data estimates with a scaling parameter (q) according to ()[]()).,0(~ where008,1974,...,2t,12,HarvesttHtHttttObservedNε/qKPrPPhHtσε==−+− We assumed that appropriate measures of the harvest observation error could be approximated by assuming a coefficient of variation for each annual harvest estimate equal to 0.15 (Paul Padding pers. comm.). The final component of the likelihood included the year specific direct recovery rates that were represented by the rate parameter (ft) of a Binomial distribution indexed by the total number of birds banded preseason and estimated with. 2,Harvesttσ ),(~,/ttttttfMBinomialmMmf= where mt is the total number of scaup banded preseason in year t and recovered during the hunting season in year t and Mt is the total number of scaup banded preseason in year t. 8.3 Bayesian Analysis Following Meyer and Millar (1999), we developed a fully conditional joint probability model, by first proposing prior distributions for all model parameters and unobserved system states and secondly by developing a fully conditional likelihood for each sampling distribution. Prior Distributions For this analysis, a joint prior distribution is required because the unknown system states P are assumed to be conditionally independent (Meyer and Millar 1999). This leads to the following joint prior distribution for the model parameters and unobserved system states ).,,,,,(),()()()()()()()(),,,,,,,(2Pr21112Pr0102Pr,...,102ProcessntttttocessocessttTocessttfKrPPpPPpPpppfpqpKprpPPfqKrPσλσσλσλΠ=−−−×= In general, we chose noninformative priors to represent the uncertainty we have in specifying the value of the parameters used in our assessment. However, we were required to use existing information to specify informative priors for the initial ratio of population size to carrying capacity (P0) as well as the reporting rate values (λt) specified above that were used to adjust the direct recovery rate estimates to harvest rates. We specified that the value of P0, ranged from the population size at maximum sustained yield (P0 = NMSY/K = (K/2)/K = 0.5) to the carrying capacity (P0= N/K= 1), using a uniform distribution on the log scale to represent this range of values. We assumed that the exploitation experienced at this population state was somewhere on the righthand shoulder of a sustained yield curve (i.e., between MSY and K). Given that we have very little evidence to suggest that historical scaup harvest levels were limiting scaup population growth, this seems like a reasonable prior distribution. We used noninformative prior distributions to represent the variance and scaling terms, while the priors for the 50 population parameters r and K were chosen to be vague but within biological bounds. These distributions were specified according to P0 ~ Uniform(ln(0.5),0), K ~ Lognormal (2.17, 0.667), r ~ Uniform (0.00001, 2), ft ~ Beta(0.5,0.5), q~ Uniform(0.0, 2), ~ Inverse Gamma (0.001, 0.001). 2Processσ Likelihood We related the observed population, total harvest estimates, and observed direct recoveries to the model parameters and unobserved system states with the following likelihood function: ).,(),,,,,,(),,(),,,,,,,,,,(12121,...,122Pr,...,1,...,1,...,1,...,1ΠΠΠ===××=TtttHarvestTttttBPOPTtttTHarvestocessttTTTTfMmpqfKrPHpKPNpPqfKrMmHNPttσλσσσλ Posterior Evaluation Using Bayes theorem we then specified a posterior distribution for the fully conditional joint probability distribution of the parameters given the observed information according to ).,,(),,,,,,(),,(),,,,,(),()()()()()()()(),,,,,,,,,,(121212Pr21112Pr0102Pr,...,1,...,1,...,1,...,1,...,102PrtTtttHarvestTttttBPOPTtttocessntttttocessocessttTTTTTocess
Click tabs to swap between content that is broken into logical sections.
Rating  
Title  Adaptive harvest management 2009 duck hunting season 
Description  adaptharvestmallards2009.pdf 
FWS Resource Links  http://library.fws.gov 
Subject 
Document Birds 
Publisher  U.S. Fish and Wildlife Service 
Date of Original  2009 
Type  Text 
Format  
Source  NCTC Conservation Library 
Rights  Public domain 
File Size  681325 Bytes 
Original Format  Document 
Length  53 
Full Resolution File Size  681325 Bytes 
Transcript  Adaptive Harvest Management 2009 HHunttinggSeasson U.S. Fish & Wildlife Service Adaptive Harvest Management 2009 Hunting Season PREFACE The process of setting waterfowl hunting regulations is conducted annually in the United States (Blohm 1989). This process involves a number of meetings where the status of waterfowl is reviewed by the agencies responsible for setting hunting regulations. In addition, the U.S. Fish and Wildlife Service (USFWS) publishes proposed regulations in the Federal Register to allow public comment. This document is part of a series of reports intended to support development of harvest regulations for the 2009 hunting season. Specifically, this report is intended to provide waterfowl managers and the public with information about the use of adaptive harvest management (AHM) for setting waterfowl hunting regulations in the United States. This report provides the most current data, analyses, and decisionmaking protocols. However, adaptive management is a dynamic process and some information presented in this report will differ from that in previous reports. ACKNOWLEDGMENTS A working group consisting of representatives from the USFWS, the U.S. Geological Survey (USGS), the Canadian Wildlife Service (CWS), and the four Flyway Councils (Appendix 1) was established in 1992 to review the scientific basis for managing waterfowl harvests. The working group, supported by technical experts from the waterfowl management and research communities, subsequently proposed a framework for adaptive harvest management, which was first implemented in 1995. The USFWS expresses its gratitude to the AHM Working Group and to the many other individuals, organizations, and agencies that have contributed to the development and implementation of AHM. This report was prepared by the USFWS Division of Migratory Bird Management. G. S. Boomer and T. A. Sanders were the principal authors. Individuals that provided essential information or otherwise assisted with report preparation were G. Zimmerman, M. Koneff, K. Richkus, E. Silverman, N. Zimpfer, J. Klimstra, K. Magruder, and P. Garrettson. Comments regarding this document should be sent to the Chief, Division of Migratory Bird Management  USFWS, 4401 North Fairfax Drive, MS MSP4107, Arlington, VA 22203. We are grateful for the continuing technical support from F. A. Johnson, M. C. Runge, and J. A. Royle (USGS), and acknowledge that information provided by USGS in this report has not received the Director's approval and, as such, is provisional and subject to revision. This information is released on the condition that neither the USGS nor the United States Government may be held liable for any damages resulting from its authorized or unauthorized use. 1 Citation: U.S. Fish and Wildlife Service. 2009. Adaptive Harvest Management: 2009 Hunting Season. U.S. Dept. Interior, Washington, D.C. 52pp. Online: http://www.fws.gov/migratorybirds/CurrentBirdIssues/Management/AHM/AHMintro.htm U.S. Fish & Wildlife Service Cover art: Joshua Spies’s painting of a longtailed duck (Clangula hyemalis) that was selected for the 2009 federal “duck stamp.”TABLE OF CONTENTS Executive Summary.............................................................................................................3 Background.........................................................................................................................4 Mallard Stocks and Flyway Management............................................................................5 Mallard Population Dynamics..............................................................................................6 HarvestManagement Objectives.......................................................................................13 Regulatory Alternatives.....................................................................................................14 Optimal Regulatory Strategies...........................................................................................18 Application of AHM Concepts to Other Stocks................................................................21 Emerging Issues within AHM............................................................................................26 Literature Cited..................................................................................................................27 Appendix 1: AHM Working Group..................................................................................30 Appendix 2: Midcontinent Mallard Models....................................................................33 Appendix 3: Eastern Mallard Models...............................................................................36 Appendix 4: Western Mallard Models..............................................................................39 Appendix 5: Modeling Mallard Harvest Rates.................................................................44 Appendix 6: Scaup Model.................................................................................................49 2 3 EXECUTIVE SUMMARY In 1995 the U.S. Fish and Wildlife Service (USFWS) implemented the Adaptive Harvest Management (AHM) program for setting duck hunting regulations in the United States. The AHM approach provides a framework for making objective decisions in the face of incomplete knowledge concerning waterfowl population dynamics and regulatory impacts. The AHM protocol is based on the population dynamics and status of three mallard (Anas platyrhynchos) stocks. Midcontinent mallards are defined as those breeding in the Waterfowl Breeding Population and Habitat Survey (WBPHS) strata 13–18, 20–50, and 75–77 plus mallards breeding in the states of Michigan, Minnesota, and Wisconsin (state surveys). The prescribed regulatory alternative for the Mississippi and Central Flyways depends exclusively on the status of these mallards. Eastern mallards are defined as those breeding in WBPHS strata 51–54 and 56 and breeding in the states of Virginia northward into New Hampshire (Atlantic Flyway Breeding Waterfowl Survey [AFBWS]). The regulatory choice for the Atlantic Flyway depends exclusively on the status of these mallards. Western mallards are defined as those birds breeding in WBPHS strata 1�����12 (hereafter Alaska) and those birds breeding in the states of California and Oregon (state surveys). The regulatory choice for the Pacific Flyway depends exclusively on the status of these mallards. Mallard population models are based on the best available information and account for uncertainty in population dynamics and the impact of harvest. Modelspecific weights reflect the relative confidence in alternative hypotheses and are updated annually using comparisons of predicted and observed population sizes. For midcontinent mallards, current model weights favor the weakly densitydependent reproductive hypothesis (88%) and suggest some preference for the additivemortality hypothesis (62%). For eastern mallards, virtually all of the weight is on models that have corrections for bias in estimates of survival or reproductive rates. Model weights do not discriminate between the strongly densitydependent (47%) and weakly densitydependent (53%) reproductive hypotheses. By consensus, hunting mortality is assumed to be additive in eastern mallards. Unlike midcontinent and eastern mallards, we consider a single functional form to predict western mallard population dynamics but consider a wide range of parameter values each weighted relative to the support from the data. For the 2009 hunting season, the USFWS is considering the same regulatory alternatives as last year. The nature of the restrictive, moderate, and liberal alternatives has remained essentially unchanged since 1997, except that extended framework dates have been offered in the moderate and liberal alternatives since 2002. Harvest rates associated with each of the regulatory alternatives have been updated based on bandreporting rate studies conducted since 1998. The expected harvest rates of adult males under liberal hunting seasons are 0.119 (SD = 0.020), 0.149 (SD = 0.043), and 0.115 (SD = 0.032) for midcontinent, eastern, and western mallards, respectively. Optimal regulatory strategies for the 2009 hunting season were calculated using: (1) harvestmanagement objectives specific to each mallard stock; (2) the 2009 regulatory alternatives; and (3) current population models. Based on this year’s survey results of 8.71 million midcontinent mallards, 3.57 million ponds in Prairie Canada, 0.908 million eastern mallards, and 0.884 million western mallards in Alaska (0.503 million) and California–Oregon (0.381 million), the optimal choice for all four flyways is the liberal regulatory alternative. AHM concepts and tools are also being applied to help improve harvest management for several other waterfowl stocks. In the last year, progress has been made in understanding the harvest potential of northern pintails (Anas acuta) and scaup (Aythya affinis, A. marila). While these biological assessments are ongoing, they are already informing decision makers and proving valuable in helping focus debate on the social aspects of harvest policy, including management objectives and the nature of regulatory alternatives.4 BACKGROUND The annual process of setting duckhunting regulations in the United States is based on a system of resource monitoring, data analyses, and rulemaking (Blohm 1989). Each year, monitoring activities such as aerial surveys and hunter questionnaires provide information on population size, habitat conditions, and harvest levels. Data collected from this monitoring program are analyzed each year, and proposals for duckhunting regulations are developed by the Flyway Councils, States, and USFWS. After extensive public review, the USFWS announces regulatory guidelines within which States can set their hunting seasons. In 1995, the USFWS adopted the concept of adaptive resource management (Walters 1986) for regulating duck harvests in the United States. This approach explicitly recognizes that the consequences of hunting regulations cannot be predicted with certainty and provides a framework for making objective decisions in the face of that uncertainty (Williams and Johnson 1995). Inherent in the adaptive approach is an awareness that management performance can be maximized only if regulatory effects can be predicted reliably. Thus, adaptive management relies on an iterative cycle of monitoring, assessment, and decisionmaking to clarify the relationships among hunting regulations, harvests, and waterfowl abundance. In regulating waterfowl harvests, managers face four fundamental sources of uncertainty (Nichols et al. 1995a, Johnson et al. 1996, Williams et al. 1996): (1) environmental variation  the temporal and spatial variation in weather conditions and other key features of waterfowl habitat; an example is the annual change in the number of ponds in the Prairie Pothole Region, where water conditions influence duck reproductive success; (2) partial controllability  the ability of managers to control harvest only within limits; the harvest resulting from a particular set of hunting regulations cannot be predicted with certainty because of variation in weather conditions, timing of migration, hunter effort, and other factors; (3) partial observability  the ability to estimate key population attributes (e.g., population size, reproductive rate, harvest) only within the precision afforded by extant monitoring programs; and (4) structural uncertainty  an incomplete understanding of biological processes; a familiar example is the longstanding debate about whether harvest is additive to other sources of mortality or whether populations compensate for hunting losses through reduced natural mortality. Structural uncertainty increases contentiousness in the decisionmaking process and decreases the extent to which managers can meet longterm conservation goals. AHM was developed as a systematic process for dealing objectively with these uncertainties. The key components of AHM include (Johnson et al. 1993, Williams and Johnson 1995): (1) a limited number of regulatory alternatives, which describe Flywayspecific season lengths, bag limits, and framework dates; (2) a set of population models describing various hypotheses about the effects of harvest and environmental factors on waterfowl abundance; (3) a measure of reliability (probability or "weight") for each population model; and (4) a mathematical description of the objective(s) of harvest management (i.e., an "objective function"), by which alternative regulatory strategies can be compared. These components are used in a stochastic optimization procedure to derive a regulatory strategy. A regulatory strategy specifies the optimal regulatory choice, with respect to the stated management objectives, for each possible combination of breeding population size, environmental conditions, and model weights (Johnson et al. 1997). The setting of annual hunting regulations then involves an iterative process: (1) each year, an optimal regulatory choice is identified based on resource and environmental conditions, and on current model weights; 5 (2) after the regulatory decision is made, modelspecific predictions for subsequent breeding population size are determined; (3) when monitoring data become available, model weights are increased to the extent that observations of population size agree with predictions, and decreased to the extent that they disagree; and (4) the new model weights are used to start another iteration of the process. By iteratively updating model weights and optimizing regulatory choices, the process should eventually identify which model is the best overall predictor of changes in population abundance. The process is optimal in the sense that it provides the regulatory choice each year necessary to maximize management performance. It is adaptive in the sense that the harvest strategy “evolves” to account for new knowledge generated by a comparison of predicted and observed population sizes. MALLARD STOCKS AND FLYWAY MANAGEMENT Since its inception AHM has focused on the population dynamics and harvest potential of mallards, especially those breeding in midcontinent North America. Mallards constitute a large portion of the total U.S. duck harvest, and traditionally have been a reliable indicator of the status of many other species. As management capabilities have grown, there has been increasing interest in the ecology and management of breeding mallards that occur outside the midcontinent region. Geographic differences in the reproduction, mortality, and migrations of mallard stocks suggest that there may be corresponding differences in optimal levels of sport harvest. The ability to regulate harvests of mallards originating from various breeding areas is complicated, however, by the fact that a large degree of mixing occurs during the hunting season. The challenge for managers, then, is to vary hunting regulations among Flyways in a manner that recognizes each Flyway’s unique breedingground derivation of mallards. Of course, no Flyway receives mallards exclusively from one breeding area; therefore, Flywayspecific harvest strategies ideally should account for multiple breeding stocks that are exposed to a common harvest. The optimization procedures used in AHM can account for breeding populations of mallards beyond the midcontinent region, and for the manner in which these ducks distribute themselves among the Flyways during the hunting season. An optimal approach would allow for Flywayspecific regulatory strategies, which represent an average of the optimal harvest strategies for each contributing breeding stock weighted by the relative size of each stock in the fall flight. This joint optimization of multiple mallard stocks requires: (1) models of population dynamics for all recognized stocks of mallards; (2) an objective function that accounts for harvestmanagement goals for all mallard stocks in the aggregate; and (3) decision rules allowing Flywayspecific regulatory choices. Currently, three stocks of mallards are officially recognized for the purposes of AHM (Fig. 1). We use a constrained approach to the optimization of these stocks’ harvest, in which the Atlantic Flyway regulatory strategy is based exclusively on the status of eastern mallards, the regulatory strategy for the Mississippi and Central Flyways is based exclusively on the status of midcontinent mallards, and the Pacific Flyway regulatory strategy is based exclusively on the status of western mallards. This approach has been determined to perform nearly as well as a jointoptimization because mixing of the three stocks during the hunting season is limited and because of the constraints imposed by management objectives and regulatory alternatives.Fig. 1. Survey areas currently assigned to the midcontinent, eastern, and western stocks of mallards for the purposes of AHM. MALLARD POPULATION DYNAMICS MidContinent Stock In 2008, midcontinent mallards were redefined as those breeding in WBPHS strata 13–18, 20–50, and 75–77, and in the Great Lakes region (Michigan, Minnesota, and Wisconsin; see Fig. 1). Estimates of the size of this population are available since 1992, and have varied from 6.4 to 11.2 million (Table 1, Fig. 2). Estimated breedingpopulation size in 2009 was 8.71 million (SE = 0.25 million), including 8.01 million (SE = 0.24 million) from the WBPHS and 0.696 million (SE = 0.056 million) from the Great Lakes region. Details describing the set of population models for midcontinent mallards are provided in Appendix 2. The set consists of four alternatives, formed by the combination of two survival hypotheses (additive vs. compensatory hunting mortality) and two reproductive hypotheses (strongly vs. weakly density dependent). Relative weights for the alternative models of midcontinent mallards changed little until all models underpredicted the change in population size from 1998 to 1999, perhaps indicating there is a significant factor affecting population dynamics that is absent from all four models (Fig. 3). Updated model weights suggest some preference for the additive 6 Table 1. Estimates (N) and associated standard errors (SE) of midcontinent mallards (in millions) observed in the WBPHS (strata 13–18, 20–50, and 75–77) and the Great Lakes region (Michigan, Minnesota, and Wisconsin). WBPHS area Great Lakes region Total Year N SE N SE N SE 1992 5.6304 0.2379 0.9946 0.1597 6.6249 0.2865 1993 5.4253 0.2068 0.9347 0.1457 6.3600 0.2529 1994 6.6292 0.2803 1.1505 0.1163 7.7797 0.3035 1995 7.7452 0.2793 1.1214 0.1965 8.8666 0.3415 1996 7.4193 0.2593 1.0251 0.1443 8.4444 0.2967 1997 9.3554 0.3041 1.0777 0.1445 10.4331 0.3367 1998 8.8041 0.2940 1.1224 0.1792 9.9266 0.3443 1999 10.0926 0.3374 1.0591 0.2122 11.1518 0.3986 2000 8.6999 0.2855 1.2350 0.1761 9.9348 0.3354 2001 7.1857 0.2204 0.8622 0.1086 8.0479 0.2457 2002 6.8364 0.2412 1.0820 0.1152 7.9184 0.2673 2003 7.1062 0.2589 0.8360 0.0734 7.9422 0.2691 2004 6.6142 0.2746 0.9333 0.0748 7.5474 0.2847 2005 6.0521 0.2754 0.7862 0.0650 6.8383 0.2830 2006 6.7607 0.2187 0.5881 0.0465 7.3488 0.2236 2007 7.7258 0.2805 0.7677 0.0584 8.4935 0.2865 2008 7.1914 0.2525 0.6750 0.0478 7.8664 0.2570 2009 8.0094 0.2442 0.6958 0.0564 8.7052 0.2506 02468101219911995199920032007YearPopulation size (millions)024681012WBPHS surveyGreat LakesTotal 7 Fig. 2. Population estimates of midcontinent mallards observed in the WBPHS (strata: 13–18, 20–50, and 75–77) and the Great Lakes region (Michigan, Minnesota, and Wisconsin). Error bars represent one standard error. 00.10.20.30.40.50.619951997199920012003200520072009YearModel Weight00.10.20.30.40.50.6ScRwSaRsSaRwScRs Fig. 3. Weights for models of midcontinent mallards (ScRs = compensatory mortality and strongly densitydependent reproduction, ScRw = compensatory mortality and weakly densitydependent reproduction, SaRs = additive mortality and strongly densitydependent reproduction, and SaRw = additive mortality and weakly densitydependent reproduction). Model weights were assumed to be equal in 1995. mortality models (62%) over those describing hunting mortality as compensatory (38%). For most of the time frame, model weights have strongly favored the weakly densitydependent reproductive models over the strongly densitydependent ones, with current model weights of 88% and 12%, respectively. The reader is cautioned, however, that models can sometimes make reliable predictions of population size for reasons having little to do with the biological hypotheses expressed therein (Johnson et al. 2002b). Eastern Stock Eastern mallards are defined as those breeding in southern Ontario and Quebec (WBPHS strata 51–54 and 56) and in the northeastern U.S. (AFBWS; Heusman and Sauer 2000; see Fig. 1). Estimates of population size have varied from 0.815 to 1.1 million since 1990, with the majority of the population accounted for in the northeastern U.S. (Table 2, Fig. 4). For 2009, the estimated breedingpopulation size of eastern mallards was 0.908 million (SE = 0.063 million), including 0.667 million (SE = 0.046 million) from the northeastern U.S. and 0.241 million (SE = 0.043 million) from the WBPHS. The reader is cautioned that these estimates differ from those reported in the USFWS annual waterfowl trend and status reports, which include composite estimates based on more fixedwing strata in eastern Canada and helicopter surveys conducted by the Canadian Wildlife Service (CWS). Details concerning the set of population models for eastern mallards are provided in Appendix 3. The set consists of six alternatives, formed by the combination of two reproductive hypotheses (strongly vs. weakly density dependent) and three hypotheses concerning bias in estimates of survival and reproductive rates (no bias vs. biased survival rates vs. biased reproductive rates). With respect to model weights, there is no single model that is clearly favored over the others at the current time. Collectively, current model weights provide little discrimination between the weakly densitydependent or strongly density dependent reproductive hypotheses, with current model weights of 53% and 47%, respectively (Fig. 5). In addition, there is overwhelming evidence of bias in extant estimates of survival or reproductive rates, assuming that survey estimates are unbiased. 8 9 Table 2. Estimates (N) and associated standard errors (SE) of eastern mallards (in millions) observed in the northeastern U.S. (AFBWS) and southern Ontario and Quebec (WBPHS strata 51–54 and 56). Northeastern U.S. Canadian survey strata Total Year N SE N SE N SE 1990 0.6651 0.0783 0.1907 0.0472 0.8558 0.0914 1991 0.7792 0.0883 0.1528 0.0337 0.9320 0.0945 1992 0.5622 0.0479 0.3203 0.0530 0.8825 0.0715 1993 0.6866 0.0499 0.2921 0.0482 0.9786 0.0694 1994 0.8563 0.0628 0.2195 0.0282 1.0758 0.0688 1995 0.8641 0.0704 0.1844 0.0400 1.0486 0.0810 1996 0.8486 0.0611 0.2831 0.0557 1.1317 0.0826 1997 0.7952 0.0496 0.2121 0.0396 1.0073 0.0634 1998 0.7752 0.0497 0.2638 0.0672 1.0390 0.0836 1999 0.8800 0.0602 0.2125 0.0369 1.0924 0.0706 2000 0.7626 0.0487 0.1323 0.0264 0.8948 0.0554 2001 0.8094 0.0516 0.2002 0.0356 1.0097 0.0627 2002 0.8335 0.0562 0.1915 0.0319 1.0250 0.0647 2003 0.7319 0.0470 0.3083 0.0554 1.0402 0.0726 2004 0.8066 0.0517 0.3015 0.0533 1.1081 0.0743 2005 0.7536 0.0536 0.2934 0.0531 1.0470 0.0755 2006 0.7214 0.0476 0.1740 0.0284 0.8954 0.0555 2007 0.6876 0.0467 0.2193 0.0336 0.9069 0.0576 2008 0.6191 0.0407 0.1960 0.0300 0.8151 0.0505 2009 0.6668 0.0457 0.2411 0.0434 0.9078 0.0630 00.20.40.60.811.21.419891994199920042009YearPopulation size (millions)00.20.40.60.811.21.4AFBWSWBPHS Total Fig. 4. Population estimates of eastern mallards observed in the northeastern states (AFBWS) and in southern Ontario and Quebec (WBPHS strata 51–54 and 56). Error bars represent one standard error. 00.050.10.150.20.250.30.3519951997199920012003200520072009YearModel WeightRw0Rs0RwSRsSRwRRsR Fig. 5. Weights for models of eastern mallards (Rw0 = weak densitydependent reproduction and no model bias, Rs0 = strong dependent reproduction and no model bias, RwS = weak densitydependent reproduction and biased survival rates, RsS = strong densitydependent reproduction and biased survival rates, RwR = weak densitydependent reproduction and biased reproductive rates, and RsR = strong densitydependent reproduction and biased reproductive rates). Model weights were assumed to be equal in 1996. 10 11 Western Stock Western mallards consist of 2 substocks and are defined as those birds breeding in Alaska (WBPHS strata 1–12) and those birds breeding in California and Oregon (state surveys; see Fig. 1). Estimates of the size of these subpopulations have varied from 0.283 to 0.843 million in Alaska since 1990 and 0.355 to 0.694 million in California and Oregon since 1992 (Table 3, Fig. 6). The total population size of western mallards has ranged from 0.748 to 1.407 million. Ideally, the western mallard stock assessment would account for mallards breeding in the states of the Pacific Flyway (including Alaska), British Columbia, and the Yukon Territory. However, we have had continuing concerns about our ability to determine changes in population size based on the collection of surveys conducted independently by Pacific Flyway States and the CWS in British Columbia. These surveys tend to vary in design and intensity, and in some cases lack measures of precision. We reviewed extant surveys to determine their adequacy for supporting a westernmallard AHM protocol and selected Alaska, California, and Oregon for modeling purposes. These three states likely harbor about 75% of the westernmallard breeding population. Nonetheless, this geographic delineation is considered temporary until surveys in other areas can be brought up to similar standards and an adequate record of population estimates is available for analysis. Details concerning the set of population models for western mallards are provided in Appendix 4. To predict changes in abundance we relied on a discrete logistic model, which combines reproduction and natural mortality into a single parameter, r, the intrinsic rate of growth. This model assumes densitydependent growth, which is regulated by the ratio of population size, N, to the carrying capacity of the environment, K (i.e., equilibrium population size in the absence of harvest). In the traditional formulation of the logistic model, harvest mortality is completely additive and any compensation for hunting losses occurs as a result of densitydependent responses beginning in the subsequent breeding season. To increase the model’s generality we included a scaling parameter for harvest that allows for the possibility of compensation prior to the breeding season. It is important to note, however, that this parameterization does not incorporate any hypothesized mechanism for harvest compensation and, therefore, must be interpreted cautiously. We modeled Alaska mallards independently of those in California and Oregon because of differing population trajectories (see Fig. 6) and substantial differences in the distribution of band recoveries. We used Bayesian estimation methods in combination with a statespace model that accounts explicitly for both process and observation error in breeding population size (Meyer and Millar 1999). Breeding population estimates of mallards in Alaska are available since 1955, but we had to limit the time series to 1990–2008 because of changes in survey methodology and insufficient bandrecovery data. The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed time series of Alaska population estimates. The estimated mean carrying capacity was 1.1 million, the intrinsic rate of growth was 0.31, while the scaling parameter estimate suggests that harvest mortality may be additive. Breeding population and harvestrate data were available for California–Oregon mallards for the period 1992–2008. The logistic model also provided a reasonable fit to these data, suggesting a mean carrying capacity of 0.7 million, an intrinsic rate of growth of 0.33, while the scaling parameter estimate suggests that harvest mortality may be only partially additive. 12 Table 3. Estimates (N) and associated standard errors (SE) of mallards (in millions) observed in Alaska (WBPHS strata 1–12) and California and Oregon (state surveys) combined. Alaska California–Oregon Total Year N SE N SE N SE 1990 0.3669 0.0370 1991 0.3853 0.0363 1992 0.3457 0.0387 0.4835 0.0605 0.8292 0.0718 1993 0.2830 0.0295 0.4654 0.0510 0.7484 0.0589 1994 0.3509 0.0371 0.4367 0.0426 0.7876 0.0565 1995 0.5242 0.0680 0.4541 0.0428 0.9783 0.0803 1996 0.5220 0.0436 0.6451 0.0802 1.1671 0.0912 1997 0.5842 0.0520 0.6390 0.1043 1.2232 0.1166 1998 0.8362 0.0673 0.4868 0.0489 1.3230 0.0832 1999 0.7131 0.0696 0.6937 0.1066 1.4068 0.1273 2000 0.7703 0.0522 0.4639 0.0532 1.2342 0.0745 2001 0.7183 0.0541 0.4044 0.0451 1.1227 0.0705 2002 0.6673 0.0507 0.3775 0.0327 1.0449 0.0603 2003 0.8435 0.0668 0.4340 0.0501 1.2775 0.0835 2004 0.8111 0.0639 0.3547 0.0352 1.1658 0.0729 2005 0.7031 0.0547 0.4014 0.0474 1.1045 0.0724 2006 0.5158 0.0469 0.4879 0.0576 1.0037 0.0743 2007 0.5815 0.0551 0.4900 0.0546 1.0715 0.0775 2008 0.5324 0.0468 0.3814 0.0478 0.9138 0.0669 2009 0.5030 0.0449 0.3815 0.0639 0.8844 0.0781 Year19901992199419961998200020022004200620082010Population size (millions)0.20.40.60.81.01.21.41.60.20.40.60.81.01.21.41.6AKCAORTotal Fig. 6. Population estimates of western mallards observed in Alaska (WBPHS strata 1–12) and California and Oregon (state surveys) combined. Error bars represent one standard error. Ideally, the development of AHM protocols for mallards would consider how different breeding stocks distribute themselves among the four Flyways so that Flywayspecific harvest strategies could account for the mixing of birds during the hunting season. At present, however, a joint optimization of western, midcontinent, and eastern stocks is not feasible due to computational hurdles. However, our preliminary analyses suggest that the lack of a joint optimization does not result in a significant decrease in performance. Therefore, the AHM protocol for western mallards is structured similarly to that used for eastern mallards, in which an optimal harvest strategy is based on the status of a single breeding stock and harvest regulations in a single flyway. Although the contribution of midcontinent mallards to the Pacific Flyway harvest is significant, we believe an independent harvest strategy for western mallards poses little risk to the midcontinent stock. Further analyses will be needed to confirm this conclusion, and to better understand the potential effect of midcontinent mallard status on sustainable hunting opportunities in the Pacific Flyway. HARVESTMANAGEMENT OBJECTIVES The basic harvestmanagement objective for midcontinent mallards is to maximize cumulative harvest over the long term, which inherently requires perpetuation of a viable population. Moreover, this objective is constrained to avoid regulations that could be expected to result in a subsequent population size below the goal of the North American Waterfowl Management Plan (NAWMP). According to this constraint, the value of harvest decreases proportionally as the difference between the goal and expected population size increases. This balance of harvest and population objectives results in a regulatory strategy that is more conservative than that for maximizing long 13 14 term harvest, but more liberal than a strategy to attain the NAWMP goal (regardless of effects on hunting opportunity). The current objective for midcontinent mallards uses a population goal of 8.5 million birds, which consists of 7.9 million mallards from the WBPHS (strata 13–18, 20–50, and 75–77) corresponding to the mallard population goal in the 1998 update of the NAWMP (less the portion of the mallard goal comprised of birds breeding in Alaska) and a goal of 0.6 million for the combined states of Michigan, Minnesota, and Wisconsin. For eastern and western mallards, there is no NAWMP goal or other established target for desired population size. Accordingly, the management objective for eastern and western mallards is simply to maximize longterm cumulative (i.e., sustainable) harvest. Additionally for western mallards, maximum longterm cumulative harvest is subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes. REGULATORY ALTERNATIVES Evolution of Alternatives When AHM was first implemented in 1995, three regulatory alternatives characterized as liberal, moderate, and restrictive were defined based on regulations used during 1979–84, 1985–87, and 1988–93, respectively. These regulatory alternatives also were considered for the 1996 hunting season. In 1997, the regulatory alternatives were modified to include: (1) the addition of a veryrestrictive alternative; (2) additional days and a higher duck bag limit in the moderate and liberal alternatives; and (3) an increase in the bag limit of hen mallards in the moderate and liberal alternatives. In 2002 the USFWS further modified the moderate and liberal alternatives to include extensions of approximately one week in both the opening and closing framework dates. In 2003 the veryrestrictive alternative was eliminated at the request of the Flyway Councils. Expected harvest rates under the veryrestrictive alternative did not differ significantly from those under the restrictive alternative, and the veryrestrictive alternative was expected to be prescribed for < 5% of all hunting seasons. Also in 2003, at the request of the Flyway Councils the USFWS agreed to exclude closed duckhunting seasons from the AHM protocol when the population size of midcontinent mallards was ≥ 5.5 million (WBPHS strata 1–18, 20–50, and 75–77 plus the Great Lakes region). Based on our original assessment, closed hunting seasons did not appear to be necessary from the perspective of sustainable harvesting when the midcontinent mallard population exceeded this level. The impact of maintaining open seasons above this level also appeared negligible for other midcontinent duck species, as based on population models developed by Johnson (2003). In 2008, because of the redefinition of the midcontinent mallard stock that excludes mallards breeding in Alaska, we rescaled the closedseason constraint. Initially, we attempted to adjust the original 5.5 million closure threshold by subtracting out the 1985 Alaska breeding population estimate, which was the year upon which the original closed season constraint was based. Our initial rescaling resulted in a new threshold equal to 5.25 million. Simulations based on optimal policies using this revised closed season constraint suggested that the Mississippi and Central Flyways would experience a 70% increase in the frequency of closed seasons. At this time, we agreed to consider alternative rescalings in order to minimize the effects on the midcontinent mallard strategy and account for the increase in mean breeding population sizes in Alaska over the past several decades. Based on this assessment, we recommended a revised closed season constraint of 4.75 million which resulted in a strategy performance equivalent to the performance expected prior to the redefinition of the midcontinent mallard stock. Because the performance of the revised strategy is essentially unchanged from the original strategy, we believe it will have no greater impact on other duck stocks in the Mississippi and Central Flyways. However, complete or partialseason closures for particular species or populations could still be deemed necessary in some situations regardless of the status of midcontinent mallards. Details of the regulatory alternatives for each Flyway are provided in Table 4. 15 Table 4. Regulatory alternatives for the 2009 duckhunting season. Flyway Regulation Atlantica Mississippi Centralb Pacific Shooting hours onehalf hour before sunrise to sunset Framework dates Restrictive Oct 1 – Jan 20 Saturday nearest Oct 1to the Sunday nearest Jan 20 Moderate and Liberal Saturday nearest September 24 to the last Sunday in January Season length (days) Restrictive 30 30 39 60 Moderate 45 45 60 86 Liberal 60 60 74 107 Bag limit (total / mallard / female mallard) Restrictive 3 / 3 / 1 3 / 2 / 1 3 / 3 / 1 4 / 3 / 1 Moderate 6 / 4 / 2 6 / 4 / 1 6 / 5 / 1 7 / 5 / 2 Liberal 6 / 4 / 2 6 / 4 / 2 6 / 5 / 2 7 / 7 / 2 a The states of Maine, Massachusetts, Connecticut, Pennsylvania, New Jersey, Maryland, Delaware, West Virginia, Virginia, and North Carolina are permitted to exclude Sundays, which are closed to hunting, from their total allotment of season days. b The High Plains Mallard Management Unit is allowed 12, 23, and 23 extra days in the restrictive, moderate, and liberal alternatives, respectively. c The Columbia Basin Mallard Management Unit is allowed seven extra days in the restrictive and moderate alternatives. RegulationSpecific Harvest Rates Harvest rates of mallards associated with each of the openseason regulatory alternatives were initially predicted using harvestrate estimates from 1979–84, which were adjusted to reflect current hunter numbers and contemporary specifications of season lengths and bag limits. In the case of closed seasons in the U.S., we assumed rates of harvest would be similar to those observed in Canada during 1988–93, which was a period of restrictive regulations both in Canada and the U.S. All harvestrate predictions were based only in part on bandrecovery data, and relied heavily on models of hunting effort and success derived from hunter surveys (Appendix C in USFWS 2002). As such, these predictions had large sampling variances and their accuracy was uncertain. In 2002, we began relying on Bayesian statistical methods for improving regulationspecific predictions of harvest rates, including predictions of the effects of frameworkdate extensions. Essentially, the idea is to use existing (prior) information to develop initial harvestrate predictions (as above), to make regulatory decisions based on those predictions, and then to observe realized harvest rates. Those observed harvest rates, in turn, are treated as 16 new sources of information for calculating updated (posterior) predictions. Bayesian methods are attractive because they provide a quantitative, formal, and an intuitive approach to adaptive management. For midcontinent mallards, we have empirical estimates of harvest rate from the recent period of liberal hunting regulations (1998–2008). The Bayesian methods thus allow us to combine these estimates with our prior predictions to provide updated estimates of harvest rates expected under the liberal regulatory alternative. Moreover, in the absence of experience (so far) with the restrictive and moderate regulatory alternatives, we reasoned that our initial predictions of harvest rates associated with those alternatives should be rescaled based on a comparison of predicted and observed harvest rates under the liberal regulatory alternative. In other words, if observed harvest rates under the liberal alternative were 10% less than predicted, then we might also expect that the mean harvest rate under the moderate alternative would be 10% less than predicted. The appropriate scaling factors currently are based exclusively on prior beliefs about differences in mean harvest rate among regulatory alternatives, but they will be updated once we have experience with something other than the liberal alternative. A detailed description of the analytical framework for modeling mallard harvest rates is provided in Appendix 5. Our models of regulationspecific harvest rates also allow for the marginal effect of frameworkdate extensions in the moderate and liberal alternatives. A previous analysis by the USFWS (2001) suggested that implementation of frameworkdate extensions might be expected to increase the harvest rate of midcontinent mallards by about 15%, or in absolute terms by about 0.02 (SD = 0.01). Based on the observed harvest rates during the 2002–2008 hunting seasons, the updated (posterior) estimate of the marginal change in harvest rate attributable to the frameworkdate extension is 0.008 (SD = 0.007). The estimated effect of the frameworkdate extension has been to increase harvest rate of midcontinent mallards by about 7% over what would otherwise be expected in the liberal alternative. However, the reader is strongly cautioned that reliable inference about the marginal effect of frameworkdate extensions ultimately depends on a rigorous experimental design (including controls and random application of treatments). Current predictions of harvest rates of adultmale midcontinent mallards associated with each of the regulatory alternatives are provided in Table 5. Predictions of harvest rates for the other age–sex cohorts are based on the historical ratios of cohortspecific harvest rates to adultmale rates (Runge et al. 2002). These ratios are considered fixed at their longterm averages and are 1.5407, 0.7191, and 1.1175 for young males, adult females, and young females, respectively. We make the simplifying assumption that the harvest rates of midcontinent mallards depend solely on the regulatory choice in the Mississippi and Central Flyways. Table 5. Predictions of harvest rates of adultmale midcontinent mallards expected with application of the 2009 regulatory alternatives in the Mississippi and Central Flyways. Regulatory alternative Mean SD Closed (U.S.) 0.0088 0.0019 Restrictive 0.0563 0.0129 Moderate 0.1029 0.0215 Liberal 0.1188 0.0198 17 The predicted harvest rates of eastern mallards are updated in the same fashion as that for midcontinent mallards based on reward banding conducted in eastern Canada and the northeastern U.S. (Appendix 5). Like midcontinent mallards, harvest rates of age and sex cohorts other than adult male mallards are based on constant rates of differential vulnerability as derived from bandrecovery data. For eastern mallards, these constants are 1.153, 1.331, and 1.509 for adult females, young males, and young females, respectively (Johnson et al. 2002a). Regulationspecific predictions of harvest rates of adultmale eastern mallards are provided in Table 6. In contrast to midcontinent mallards, frameworkdate extensions were expected to increase the harvest rate of eastern mallards by only about 5% (USFWS 2001), or in absolute terms by about 0.01 (SD = 0.01). Based on the observed harvest rates during the 2002–2008 hunting seasons, the updated (posterior) estimate of the marginal change in harvest rate attributable to the frameworkdate extension is 0.003 (SD = 0.010). The estimated effect of the frameworkdate extension has been to increase harvest rate of eastern mallards by about 2% over what would otherwise be expected in the liberal alternative. Table 6. Predictions of harvest rates of adultmale eastern mallards expected with application of the 2009 regulatory alternatives in the Atlantic Flyway. Regulatory alternative Mean SD Closed (U.S.) 0.0797 0.0233 Restrictive 0.1108 0.0393 Moderate 0.1412 0.0473 Liberal 0.1492 0.0433 Based on available estimates of harvest rates of mallards banded in California and Oregon during 1990–1995 and 2002–2007, there was no apparent relationship between harvest rate and regulatory changes in the Pacific Flyway. This is unusual given our ability to document such a relationship in other mallard stocks and in other species. We note, however, that the period 2002–2007 was comprised of both stable and liberal regulations and harvest rate estimates were based solely on reward bands. Regulations were relatively restrictive during most of the earlier period and harvest rates were estimated based on standard bands using reporting rates estimated from reward banding during 1987–1988. Additionally, 1993–1995 were transition years in which fulladdress and tollfree bands were being introduced and information to assess their reporting rates (and their effects on reporting rates of standard bands) is limited. Thus, the two periods in which we wish to compare harvest rates are characterized not only by changes in regulations, but also in estimation methods. Consequently, we lack a sound empirical basis for predicting harvest rates of western mallards associated with current regulatory alternatives in the Pacific Flyway. This year, we applied Bayesian statistical methods for improving regulationspecific predictions of harvest rates (see Appendix 5). The methodology is analogous to that currently in use for midcontinent and eastern mallards except that the marginal effect of framework date extensions in moderate and liberal alternatives is inestimable because there are no data prior to implementation of extensions. We specified prior regulationspecific harvest rates of 0.01, 0.06, 0.09, and 0.11 with associated standard deviations of 0.003, 0.02, 0.03, and 0.03 for the closed, restrictive, moderate, and liberal alternatives, respectively. The harvest rates for the liberal alternative were based on empirical estimates realized under the current liberal alternative during 2002–2007 and determined from adultmale mallards banded with reward bands in California and Oregon. Harvest rates for the moderate and restrictive alternatives were based on the proportional (0.85 and 0.51) difference in harvest rates expected for midcontinent mallards under the respective alternatives. And finally, harvest rate for the closed alternative was based on what we might realize with a closed season in the U.S. (including Alaska) and a very restrictive season in Canada, similar to that for midcontinent mallards. A relatively large standard deviation (CV=0.3) was chosen to reflect greater uncertainty about the means than that for midcontinent mallards (CV=0.2). Current predictions of harvest rates of adultmale western mallards associated with each regulatory alternative are provided in Table 7. 18 Table 7. Predictions of harvest rates of adultmale western mallards expected with application of the 2009 regulatory alternatives in the Pacific Flyway. Regulatory alternative Mean SD Closed (U.S.) 0.010 0.018 Restrictive 0.059 0.016 Moderate 0.098 0.027 Liberal 0.115 0.032 OPTIMAL REGULATORY STRATEGIES We calculated optimal regulatory strategies using stochastic dynamic programming (Lubow 1995, Johnson and Williams 1999). For the Mississippi and Central Flyways, we based this optimization on: (1) the 2009 regulatory alternatives, including the closedseason constraint; (2) current population models and associated weights for midcontinent mallards; and (3) the dual objectives of maximizing longterm cumulative harvest and achieving a population goal of 8.5 million midcontinent mallards. The resulting regulatory strategy (Table 8) is similar to that used last year. Note that prescriptions for closed seasons in this strategy represent resource conditions that are insufficient to support one of the current regulatory alternatives, given current harvestmanagement objectives and constraints. However, closed seasons under all of these conditions are not necessarily required for longterm resource protection, and simply reflect the NAWMP population goal and the nature of the current regulatory alternatives. Assuming that regulatory choices adhered to this strategy (and that current model weights accurately reflect population dynamics), breedingpopulation size would be expected to average 6.83 million (SD = 1.83 million). Based on an estimated population size of 8.71 million midcontinent mallards and 3.57 million ponds in Prairie Canada, the optimal choice for the Mississippi and Central Flyways in 2009 is the liberal regulatory alternative. We calculated an optimal regulatory strategy for the Atlantic Flyway based on: (1) the 2009 regulatory alternatives; (2) current population models and associated weights for eastern mallards; and (3) an objective to maximize longterm cumulative harvest. The resulting strategy suggests liberal regulations for all population sizes of record, and is characterized by a lack of intermediate regulations (Table 9). We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model weights accurately reflect population dynamics), breedingpopulation size would be expected to average 0.911 million (SD = 0.173 million). Based on an estimated breeding population size of 0.908 million mallards, the optimal choice for the Atlantic Flyway in 2009 is the liberal regulatory alternative. We calculated an optimal regulatory strategy for the Pacific Flyway based on: (1) the 2009 regulatory alternatives, (2) current (1990–2008) population models and parameter estimates, and (3) an objective to maximize longterm cumulative harvest subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes (Table 10). We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model parameters accurately reflect population dynamics), breedingpopulation size would be expected to average 1.0 million (SD = 0.22 million) in Alaska and 0.46 million (SD = 0.03 million) in California and Oregon. Based on an estimated breeding population size of 0.503 million mallards in Alaska and 0.381 million in California and Oregon, the optimal choice for the Pacific Flyway in 2009 is the liberal regulatory alternative (see Table 10).19 Table 8. Optimal regulatory strategya for the Mississippi and Central Flyways for the 2009 hunting season. This strategy is based on current regulatory alternatives (including the closedseason constraint), on current midcontinent mallard models and weights, and on the dual objectives of maximizing longterm cumulative harvest and achieving a population goal of 8.5 million mallards. The shaded cell indicates the regulatory prescription for 2009. Pondsc Bpopb 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 ≤ 4.5 C C C C C C C C C C 4.75–5.75 R R R R R R R R R R 6.00 R R R R R R R R M M 6.25 R R R R R R M M M L 6.5 R R R R M M M L L L 6.75 R R R M L L L L L L 7.0 R M M L L L L L L L 7.25 M M L L L L L L L L 7.5 M L L L L L L L L L ≥ 7.75 L L L L L L L L L L a C = closed season, R = restrictive, M = moderate, L = liberal. b Mallard breeding population size (in millions) in the WBPHS (strata 13–18, 20–50, 75–77) and Michigan, Minnesota, and Wisconsin. c Ponds (in millions) in Prairie Canada in May. Table 9. Optimal regulatory strategya for the Atlantic Flyway for the 2009 hunting season. This strategy is based on current regulatory alternatives, on current eastern mallard models and weights, and on an objective to maximize longterm cumulative harvest. The shaded cell indicates the regulatory prescription for 2009. Mallardsb Regulation ≤ 0.250 C 0.275 R ≥ 0.300 L a C = closed season, R = restrictive, M = moderate, and L = liberal. b Estimated number of mallards in eastern Canada (WBPHS strata 51–54, 56) and the northeastern U.S. (AFBWS), in millions. 20 Table 10. Optimal regulatory strategya for the Pacific Flyway during the 2009 hunting season. This strategy is based on the 2009 regulatory alternatives, current (1990–2008) population models and parameter estimates, and an objective to maximize longterm cumulative harvest subject to a constraint intended to prevent extreme changes in regulations associated with relatively small changes in population sizes. The shaded cell indicates the regulatory prescription for 2009. Alaska BPOPb CAOR BPOPb 0 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 ≥ 0.55 0 C C M L L L L L L L L L 0.05 C C R R R M M M L L L L 0.10 C R R M M M L L L L L L 0.15 R R M M L L L L L L L L 0.20 M R L L L L L L L L L L 0.25 M R L L L L L L L L L L 0.30 M M L L L L L L L L L L 0.35 M M L L L L L L L L L L 0.40 M M L L L L L L L L L L ≥ 0.45 L M L L L L L L L L L L a C = closed season, R = restrictive, M = moderate, and L = liberal. b Estimated number of mallards in millions for Alaska (WBPHS strata 1–12) and in California and Oregon.Application of AHM Concepts to Other Stocks The USFWS is striving to apply the principles and tools of AHM to improve decisionmaking for several other stocks of waterfowl. Over the last year, some progress has been made to develop AHM frameworks for American black ducks (Anas rubripes) and the Atlantic Population of Canada geese (Branta canadensis), but these results are not yet finalized for inclusion in this year’s report. As these frameworks are developed further, we will continue to describe this work in subsequent reports. Below, we provide the 2009 updates for two decisionmaking frameworks that are currently informing harvest management. Northern Pintails The Flyway Councils have long identified the northern pintail as a highpriority species for inclusion in the AHM process. In 1997, the USFWS adopted a pintail harvest strategy to help align harvest opportunity with population status, while providing a foundation upon which to develop a formal AHM framework. Since 1997, the harvest strategy has undergone a number of technical improvements and policy revisions. However, the strategy continues to be a set of regulatory prescriptions born out of consensus, rather than an optimal strategy derived from agreedupon population models, management objectives, regulatory alternatives, and measures of uncertainty. In 2007, the USFWS and Flyway Councils took a major step towards a truly adaptive approach by incorporating alternative models of population dynamics. Two models are considered: one in which harvest is additive to natural mortality, and another in which harvest losses are compensated for by reductions in natural mortality. In the additive model, winter survival rate is a constant, whereas winter survival is densitydependent in the compensatory model. We here provide a summary of these recent modeling efforts. A detailed progress report is available online (http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/NOPI%20Harvest%20Strategy%202007.pdf). The predicted in year t + 1 () for the additive harvest mortality model is calculated as tcBPOP1+tcBPOP {}wttRsttscHRscBPOPcBPOP)1/(ˆ)ˆ1(1−−+=+γ where is the latitudeadjusted breeding population size in year t, and are the summer and winter survival rates, respectively, tcBPOPsswsRγ is a biascorrection constant for the age ratio, c is the crippling loss rate, is the predicted age ratio, and is the predicted continental harvest. The model uses the following constants: = 0.70, = 0.93, tRˆstHˆswsRγ = 0.80, and c = 0.20. The compensatory harvest mortality model serves as a hypothesis that stands in contrast to the additive harvest mortality model, positing a strong but realistic degree of compensation. The compensatory model assumes that the mechanism for compensation is densitydependent postharvest (winter) survival. The form is a logistic relationship between winter survival and postharvest population size, with the relationship anchored around the historic mean values for each variable. For the compensatory model then, predicted winter survival rate in year t () is calculated as ts[]1))((0101)(−−+−+−+=PPbattessss, 21 where (upper asymptote) is 1.0, (lower asymptote) is 0.7, b (slope term) is 1.0, is the postharvest population size in year t (expressed in millions),1s0stPP is the mean postharvest population size (4.295 million from 1974 through 2005), and a = logit010ssss⎛⎞−⎜⎟−⎝⎠ or ⎭⎬⎫⎩⎨⎧⎟⎟⎠⎞⎜⎜⎝⎛−−−−⎟⎟⎠⎞⎜⎜⎝⎛−−=0100101loglogssssssssa, where s is 0.93 (mean winter survival rate). At moderate population size and latitude, the compensatory model allows for greater harvest (Fig. 7) than does the additive model (note especially that the size of the restrictive region [seasonwithinaseason] is smaller and is invoked when the latitude is higher). Also, 2 and 3bird bag limits are called for under more circumstances. But, at high population sizes, the higher bag limits are called for less often, because the compensatory model predicts that growth of the population will be slower (densitydependence). The fit to historic data was used to compare the additive and compensatory harvest models. From the , , and observed harvest () for the period 1974–through year t, the subsequent year’s breeding population size (on the latitudeadjusted scale) was predicted with both the additive and compensatory model, and compared to the observed breeding population size (on the latitudeadjusted scale). The meansquared error of the predictions from the additive model () was calculated as: tcBPOPtmLATtHaddMSE Σ=−+−=ttaddttaddcBPOPcBPOPtMSE19752)(1)1975(1 and the meansquared error of the predictions from the compensatory model were calculated in a similar manner. The model weights for the additive and compensatory models were calculated from their relative meansquared errors. The model weight for the additive model () was calculated as: addW compaddaddaddMSEMSEMSEW111+=. The model weight for the compensatory model was found in a corresponding manner, or by subtracting the additive model weight from 1.0. As of 2008, the compensatory model did not fit the historic data as well as the additive model; the model weights were 0.597 for the additive model and 0.403 for the compensatory model, unchanged from 2007. The 2008 average model calls for a strategy that is intermediate between the additive and compensatory models (see Fig. 7). The USFWS remains committed to the development of an AHM framework to inform pintail harvest management based on a formal, derived strategy and clearly articulated management objectives. We are committed to continue to work with the Flyways to finalize an AHM protocol with the intent of implementing a revised harvest strategy for pintails in the near future. 22 23 1234567515253545556575859Average Latitude of the BPOPClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L21234567515253545556575859Average Latitude of the BPOPClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L21234567515253545556575859Average Latitude of the BPOPPintail BPOP (millions)ClosedRestrictiveLiberal (1bird)Liberal (3birds)L212345671234567515253545556575859515253545556575859Average L2(A) Additive Model(B) Compensatory Model(C) Fig. 7. Statedependent harvest policy for northern pintails with (A) additive, (B) compensatory, and (C) 2008 weighted models. In each case the strategy assumes that the general duck hunting season is that prescribed under the liberal regulatory alternative. 2006 200824 Scaup In 2008, the USFWS implemented a decisionmaking framework for establishing scaup harvest regulations that was initially proposed in 2007 (Boomer and Johnson 2007). In addition, the USFWS committed to continue working with the Flyways to develop an alternative scaup population model for inclusion in the current decisionmaking framework. This model will capture the belief that the scaup population will decline to and stabilize at some lower equilibrium level in response to a declining carrying capacity and that harvest at current levels is completely compensatory. We plan to report on our efforts to develop an alternative model at the 2009 AHM Working Group meeting. In 2007, the USFWS also outlined methods to facilitate the specification of regulatory alternatives for scaup harvest management (Boomer et al. 2007). We proposed harvest thresholds to be considered under regulatory alternatives based on a simulation of an optimal policy that was derived under an objective to achieve 95% of the longterm cumulative harvest. We used this objective because it results in a strategy less sensitive to small change in population size as compared to a strategy derived under an objective to achieve 100% of longterm cumulative harvest. In addition, the 95% objective allows for some harvest opportunity at relatively low population sizes. We have continued to work with the Flyways to determine what regulations would achieve the allowable harvest thresholds set forth in Boomer et al. (2007). In 2008 during deliberations over scaup regulatory alternatives, the USFWS also agreed to consider a “hybrid season” option that would be available to all Flyways for the restrictive and moderate alternatives. In 2008, initial Restrictive, Moderate, and Liberal scaup regulatory alternatives were defined and implemented in all four Flyways. Subsequent concerns and dialogue led the USFWS to further clarify criteria associated with the establishment of hybrid seasons and to allow additional modifications of the alternatives for each Flyway in 2009. Final scaup regulatory alternatives were adopted for each Flyway in 2009. These alternatives will remain in place for a period of three years and then revisited as the latest harvest information is evaluated. The USFWS will continue to work with the Flyways to determine acceptable harvestmanagement objectives and evaluate regulatory alternatives to be used in the evolving decisionmaking framework for scaup harvest management. Presently, the scaup harvest strategy prescribes optimal harvest levels, not optimal regulatory alternatives. It is important to note that we currently have limited ability to predict expected scaup harvest under the newlyestablished, Flywayspecific scaup regulatory alternatives. The initial regulatory alternatives adopted for each Flyway were based on relatively crude predictions from harvest models developed in Boomer et al. (2007) or alternative harvest models proposed by the Flyways. As we gain experience with scaup regulatory alternatives, we will refine predicted harvest distributions associated with the Flywayspecific alternatives with the ultimate goal being to use regulatory alternatives, as opposed to harvest, as the control variable in deriving future scaup harvest policies. The lack of scaup demographic information over a sufficient timeframe and at a continental scale precludes the use of a traditional balance equation to represent scaup population and harvest dynamics. As a result, we used a discretetime, stochastic, logisticgrowth population model to represent changes in scaup abundance, while explicitly accounting for scaling issues associated with the monitoring data. Details describing the modeling and assessment framework that has been developed for scaup can be found in Appendix 6 and in Boomer and Johnson (2007). We updated the scaup assessment based on the current model formulation and data extending from 1974 through 2008. As in past analyses, the state space formulation and Bayesian analysis framework provided reasonable fits to the observed breeding population and total harvest estimates with realistic measures of variation. The posterior mean estimate of the intrinsic rate of increase (r) is 0.106 while the posterior mean estimate of the carrying 25 capacity (K) is 8.44 million birds. The posterior mean estimate of the scaling parameter (q) is 0.552, ranging between 0.479 and 0.634 with 95% probability. We calculated an optimal harvest policy for scaup based on: (1) a control variable of total harvest (U.S. and Canada combined), (2) current population model and updated parameter estimates, and (3) an objective to achieve 95% of the longterm cumulative harvest. We simulated the use of this regulatory strategy to determine expected performance characteristics. Assuming that harvest management adhered to this strategy (and that current model parameters accurately reflect population dynamics), breedingpopulation size would be expected to average 4.64 million (SD = 0.86 million). With an estimated breeding population size of 4.2 million scaup, the optimal harvest level for scaup is 0.3 million (Table 11). Based on the harvest thresholds specified in Boomer et al. (2007), this year’s optimal harvest corresponds to the moderate regulatory alternative. Table 11. Optimal scaup harvest levels (observed scale in millions) and corresponding breeding population sizes (in millions). This strategy is based on the current scaup population model, and an objective to maximize 95% of longterm cumulative harvest. The shaded cell indicates the optimal harvest level for 2009. BPOP Optimal Harvest 0.0 – 1.8 0 2.0 – 2.4 0.05 2.6 – 2.8 0.1 3.0 – 3.2 0.15 3.4 – 3.6 0.2 3.8 – 4.0 0.25 4.2 – 4.4 0.3 4.6 0.35 4.8 – 5.2 0.4 5.4 – 5.6 0.5 26 Emerging Issues in AHM Under the existing AHM protocol learning occurs passively as annual comparisons of model predictions to observations from monitoring programs are used to update model weights and relative beliefs about system responses to management (Johnson et al. 2002b) or as model parameters are updated based on an assessment of the most recent monitoring data (Boomer and Johnson 2007, Johnson et al. 2007). However, additional learning can also occur as decisionmaking frameworks are evaluated to determine if objectives are being achieved, if objectives are adequately specified or have changed, or if other aspects of the decision problem are adequately being addressed. Often the feedback resulting from this process results in a form of “double loop” learning (Lee 1993) that offers the opportunity to adapt decisionmaking frameworks in response to a shifting decision context, novel or emerging management alternatives, or to a need to revise models that may perform poorly or may need to account for new information. Adaptive management depends on this iterative process to ensure that decisionmaking protocols remain relevant in evolving natural and social systems. As a byproduct of the adaptive management process, it is natural to think about when or how decisionmaking protocols should be revised to incorporate new information or to accommodate changes in the overall management context. Recent outcomes from the 2008 Future of Waterfowl Management Workshop and evaluations of waterfowl harvest and habitat management programs (e.g., Anderson et al. 2007, Assessment Steering Committee 2007) suggest compelling reasons for a reevaluation of the objectives of waterfowl management as well as the technical and institutional frameworks through which harvest and habitat management decisions are made. We view such introspection as completely consistent with the principles underlying AHM, and, in fact, critical to its longterm utility as a decision framework for waterfowl harvest regulation. The AHM Working Group has begun to consider the need to better integrate harvest and habitat management objectives (sensu Runge et al. 2005, Anderson et al. 2007). In addition, we have begun to explore the implications of largescale system changes (e.g., related to climate change) in habitats and the environment on the decisionmaking frameworks used to inform harvest management. Additionally, we continue to explore questions related to the appropriate taxonomic, spatial, and temporal resolution of waterfowl harvest management. In light of all these issues, we remain committed to the continued reassessment of waterfowl harvest decision frameworks. Ultimately, we will need to prioritize the scope of work required to revisit the AHM protocol in the face of limited resources. We look forward to working with the Flyways as we continue to refine the adaptive harvest management framework. 27 LITERATURE CITED Anderson, D. R., and K. P. Burnham. 1976. Population ecology of the mallard. VI. The effect of exploitation on survival. U.S. Fish and Wildlife Service Resource Publication No. 128. 66pp. Anderson, M. G., D. Caswell, J. M. Eadie, J. T. Herbert, M. Huang, D. D. Humburg, F. A. Johnson, M. D. Koneff, S. E. Mott, T. D. Nudds, E. T. Reed, J. K. Ringleman, M. C. Runge, and B. C. Wilson. 2007. Report from the Joint Task Group for clarifying North American waterfowl management plan population objectives and their use in harvest management. Unpublished report. Assessment Steering Committee. 2007. North American waterfowl management plan: continental progress report. Unpublished final report. Blohm, R. J. 1989. Introduction to harvest  understanding surveys and season setting. Proceedings of the International Waterfowl Symposium 6:118–133. Blohm, R. J., R. E. Reynolds, J. P. Bladen, J. D. Nichols, J. E. Hines, K. P. Pollock, and R. T. Eberhardt. 1987. Mallard mortality rates on key breeding and wintering areas. Transactions of the North American Wildlife and Natural Resources Conference 52:246–263. Boomer, G. S., and F. A. Johnson. 2007. A proposed assessment and decisionmaking framework to inform scaup harvest management. Unpublished Report. U. S. Fish and Wildlife Service, Laurel, MD. 26pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/SCAUP2007Report.pdf Boomer, G. S., F. A. Johnson, M. D. Koneff, T. A. Sanders, and R. E. Trost. 2007. A process to determine scaup regulatory alternatives. Unpublished Scoping Document. U. S. Fish and Wildlife Service, Laurel, MD. 20pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/SpecialTopics/BySpecies/scaup_regs_scoping_draftVI.pdf Brooks, S. P., and A. Gelman. 1998 Alternative methods for monitoring convergence of iterative simulations. Journal of Computational and Graphical Statistics 7:434–455. Burnham, K. P., G. C. White, and D. R. Anderson. 1984. Estimating the effect of hunting on annual survival rates of adult mallards. Journal of Wildlife Management 48:350–361. Henny, C. J., and Burnham K. P. 1976. A reward band study of mallards to estimate reporting rates. Journal of Wildlife Management. 40:1–14. Heusman, H W, and J. R. Sauer. 2000. The northeastern states’ waterfowl breeding population survey. Wildlife Society Bulletin 28:355–364. Hodges, J. L., J. G. King, B. Conant, and H. A. Hanson. 1996. Aerial surveys of waterbirds in Alaska 1975–94: population trends and observer variability. Information and Technology Report 4, National Biological Service, U.S. Dept. of the Interior, Washington, D.C. 25pp. Johnson, F. A. 2003. Population dynamics of ducks other than mallards in midcontinent North America. Draft. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 15pp. 28 Johnson, F. A., G. S. Boomer, and T. A. Sanders. 2007. A proposed protocol for the adaptive harvest management of mallards breeding in Western North America. Unpublished Report. U. S. Fish and Wildlife Service, Laurel, MD. 33pp. Johnson, F. A., J. A. Dubovsky, M. C. Runge, and D. R. Eggeman. 2002a. A revised protocol for the adaptive harvest management of eastern mallards. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 13pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/emalahm2002.pdf. Johnson, F. A., W. L. Kendall, and J. A. Dubovsky. 2002b. Conditions and limitations on learning in the adaptive management of mallard harvests. Wildlife Society Bulletin 30:176–185. Johnson, F. A., C. T. Moore, W. L. Kendall, J. A. Dubovsky, D. F. Caithamer, J. R. Kelley, Jr., and B. K. Williams. 1997. Uncertainty and the management of mallard harvests. Journal of Wildlife Management 61:202–216. Johnson, F. A., and B. K. Williams. 1999. Protocol and practice in the adaptive management of waterfowl harvests. Conservation Ecology 3(1): 8. [online] URL: http://www.consecol.org/vol3/iss1/art8. Johnson, F. A., B. K. Williams, J. D. Nichols, J. E. Hines, W. L. Kendall, G. W. Smith, and D. F. Caithamer. 1993. Developing an adaptive management strategy for harvesting waterfowl in North America. Transactions of the North American Wildlife and Natural Resources Conference 58:565–583. Johnson, F. A., B. K. Williams, and P. R. Schmidt. 1996. Adaptive decisionmaking in waterfowl harvest and habitat management. Proceedings of the International Waterfowl Symposium 7:26–33. Lubow, B. C. 1995. SDP: Generalized software for solving stochastic dynamic optimization problems. Wildlife Society Bulletin 23:738–742. Meyer, R., and R. B. Millar. 1999. BUGS in Bayesian stock assessments. Canadian Journal of Fisheries and Aquatic Sciences 56:1078–1086. Millar, R. B., and R. Meyer. 2000. Nonlinear state space modeling of fisheries biomass dynamics by using MetropolisHastings within Gibbs sampling. Applied Statistics 49: 327–342. Nichols, J. D., F. A. Johnson, and B. K. Williams. 1995a. Managing North American waterfowl in the face of uncertainty. Annual Review of Ecology and Systematics 26:177–199. Nichols, J. D., R. E. Reynolds, R. J. Blohm, R. E. Trost, J. E. Hines, and J. P. Bladen. 1995b. Geographic variation in band reporting rates for mallards based on reward banding. Journal of Wildlife Management 59:697–708. Runge, M. C., F. A. Johnson, J. A. Dubovsky, W. L. Kendall, J. Lawrence, and J. Gammonley. 2002. A revised protocol for the adaptive harvest management of midcontinent mallards. Fish and Wildlife Service, U.S. Dept. Interior, Washington, D.C. 28pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/MCMrevise2002.pdf. Runge, M. C., F. A. Johnson, M. G. Anderson, M. D. Koneff, E. T. Reed, and S. E. Mott. 2005. The need for coherence between waterfowl harvest and habitat management. Wildlife Society Bulletin 34:1231–1237. Schaefer, M. B. 1954. Some aspects of the dynamics of populations important to the management of commercial marine fisheries. Bulletin of the InterAmerican Tropical Tuna Commission 1:25–56. 29 Smith, G. W. 1995. A critical review of the aerial and ground surveys of breeding waterfowl in North America. Biological Science Report 5, National Biological Service, U.S. Dept. of the Interior, Washington, D.C. 252pp. Spiegelhalter, D. J., A. Thomas, N. Best, and D. Lunn. 2003. WinBUGS 1.4 User manual. MRC Biostatistics Unit, Institute of Public Health, Cambridge, UK. U.S. Fish and Wildlife Service. 2000. Adaptive harvest management: 2000 duck hunting season. U.S. Dept. Interior, Washington. D.C. 43pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2000/ahm2000.pdf U.S. Fish and Wildlife Service. 2001. Frameworkdate extensions for duck hunting in the United States: projected impacts & coping with uncertainty. U.S. Dept. Interior, Washington, D.C. 8pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2001/ahm2001.PDF U.S. Fish and Wildlife Service. 2002. Adaptive harvest management: 2002 duck hunting season. U.S. Dept. Interior, Washington. D.C. 34pp. [online] URL: http://www.fws.gov/migratorybirds/NewReportsPublications/AHM/Year2002/2002AHMreport.pdf Walters, C. J. 1986. Adaptive management of renewable resources. MacMillan Publ. Co., New York, N.Y. 374pp. Williams, B. K., and F. A. Johnson. 1995. Adaptive management and the regulation of waterfowl harvests. Wildlife Society Bulletin 23:430–436. Williams, B. K., F. A. Johnson, and K. Wilkins. 1996. Uncertainty and the adaptive management of waterfowl harvests. Journal of Wildlife Management 60:223–232. 30 APPENDIX 1: AHM Working Group (Note: This list includes only permanent members of the AHM Working Group. Not listed here are numerous persons from federal and state agencies that assist the Working Group on an adhoc basis.) Coordinator: Scott Boomer U.S. Fish & Wildlife Service 11510 American Holly Drive Laurel, Maryland 207084017 phone: 3014975684 fax: 3014975871 email: scott_boomer@fws.gov USFWS Representatives: Brad Bortner (Region 1) U.S. Fish and Wildlife Service 911 NE 11th Ave. Portland, OR 972324181 phone: 5032316164 fax: 5032312364 email: brad_bortner@fws.gov Dave Case (contractor) D.J. Case & Associates 607 Lincolnway West Mishawaka, IN 46544 phone: 5742580100 fax: 5742580189 email: dave@djcase.com Jim Dubovsky (Region 6) U.S. Fish and Wildlife Service P.O. Box 25486DFC Denver, CO 802250486 phone: 3032364403 fax: 3032368680 email:james_dubovsky@fws.gov Jeff Haskins (Region 2) U.S. Fish and Wildlife Service P.O. Box 1306 Albuquerque, NM 87103 phone: 5052486827 (ext 30) fax: 5052487885 email: jeff_haskins@fws.gov Diane Pence (Region 5) Jim Kelley (Region 9) U.S. Fish and Wildlife Service 1 Federal Drive Fort Snelling, MN 551110458 phone: 6127135409 fax: 6127135393 email: james_r_kelley@fws.gov Sean Kelly (Region 3) U.S. Fish and Wildlife Service 1 Federal Drive Ft. Snelling, MN 551114056 phone: 6127135470 fax: 6127135393 email: sean_kelly@fws.gov Mark Koneff (Region 9) U.S. Fish & Wildlife Service 11510 American Holly Drive Laurel, Maryland 207084017 phone: 3014975648 fax: 3014975871 email: mark_koneff@fws.gov Paul Padding (Region 9) U.S. Fish and Wildlife Service 11510 American Holly Drive Laurel, MD 20708 phone: 3014975851 fax: 3014975885 email: paul_padding@fws.gov 31 U.S. Fish and Wildlife Service 300 Westgate Center Drive Hadley, MA 010359589 phone: 4132538577 fax: 4132538424 email: diane_pence@fws.gov Russ Oates (Region 7) U.S. Fish and Wildlife Service 1011 East Tudor Road Anchorage, AK 995036119 phone: 9077863446 fax: 9077863641 email: russ_oates@fws.gov Dave Sharp (Region 9) U.S. Fish and Wildlife Service P.O. Box 25486, DFC Denver, CO 802250486 phone: 3032752386 fax: 3032752384 email: dave_sharp@fws.gov Bob Trost (Region 9) U.S. Fish and Wildlife Service 911 NE 11th Ave. Portland, OR 972324181 phone: 5032316162 fax: 5032316228 email: robert_trost@fws.gov David Viker (Region 4) U.S. Fish and Wildlife Service 1875 Century Blvd., Suite 345 Atlanta, GA 30345 phone: 4046797188 fax: 4046797285 email: david_viker@fws.gov Canadian Wildlife Service Representatives: Dale Caswell Canadian Wildlife Service 123 Main St. Suite 150 Winnipeg, Manitoba, Canada R3C 4W2 phone: 2049835260 fax: 2049835248 email: dale.caswell@ec.gc.ca Eric Reed Canadian Wildlife Service 351 St. Joseph Boulevard Hull, QC K1A OH3, Canada phone: 8199530294 fax: 8199536283 email: eric.reed@ec.gc.ca Flyway Council Representatives: Min Huang (Atlantic Flyway) CT Dept. of Environmental Protection Franklin Wildlife Mgmt. Area 391 Route 32 North Franklin, CT 06254, USA Phone: 860/6426528 fax: 860/6427964 email: min.huang@po.state.ct.us Mike Johnson (Central Flyway) North Dakota Game and Fish Department 100 North Bismarck Expressway Bismarck, ND 585015095 phone: 7013286319 fax: 7013286352 email: mjohnson@state.nd.us Bryan Swift (Atlantic Flyway) Dept. Environmental Conservation Larry Reynolds (Mississippi Flyway) LA Dept. of Wildlife & Fisheries PO Box 98000 Baton Rouge, LA 708989000, USA Phone: 225/7650456 Fax: 225/7635456 email: lreynolds@wlf.state.la.us Jon Runge (Pacific Flyway) Colorado Division of Wildlife 317 W. Prospect Fort Collins, CO 80526 Phone: 9704724365 email: Jon.Runge@state.co.us Dan Yparraguirre (Pacific Flyway) California Dept. of Fish and Game 32 625 Broadway Albany, NY 122334754 phone: 5184028866 fax: 5184029027 or 4028925 email: blswift@gw.dec.state.ny.us Mark Vrtiska (Central Flyway) Nebraska Game and Parks Commission P.O. Box 30370 2200 North 33rd Street Lincoln, NE 685031417 phone: 4024715437 fax: 4024715528 email: mvrtiska@ngpc.state.ne.us 1812 Ninth Street Sacramento, CA 95814 phone: 9164453685 email: dyparraguirre@dfg.ca.gov Guy Zenner (Mississippi Flyway) Iowa Dept. of Natural Resources 1203 North Shore Drive Clear Lake, IA 50428 phone: 5153573517, ext. 23 fax: 5153575523 email: gzenner@netins.net USGS Technical Consultants: Fred Johnson Florida Integrated Science Center U. S. Geological Survey P.O. Box 110485 Gainesville, FL 32611 phone: 3523925075 fax: 3528460841 email: fjohnson@usgs.gov Mike Runge Patuxent Wildlife Research Center U. S. Geological Survey 12100 Beech Forest Rd. Laurel, MD 20708 phone: 3014975748 fax: 3014975545 email: mrunge@usgs.gov Andy Royle Patuxent Wildlife Research Center U. S. Geological Survey 12100 Beech Forest Rd. Laurel, MD 20708 phone: 3014975846 fax: 3014975545 email: aroyle@usgs.gov APPENDIX 2: Midcontinent Mallard Models In 1995, we developed population models to predict changes in midcontinent mallards based on the traditional survey area which includes individuals from Alaska (Johnson et al. 1997). In 1997, we added mallards from the Great Lakes region (Michigan, Minnesota, and Wisconsin) to the midcontinent mallard stock, assuming their population dynamics were equivalent. In 2002, we made extensive revisions to the set of alternative models describing the population dynamics of midcontinent mallards (Runge et al. 2002, USFWS 2002). In 2008, we redefined the population of midcontinent mallards to account for the removal of Alaskan birds (WBPHS strata 1–12) that are now considered to be in the western mallard stock and have subsequently rescaled the model set appropriately. Model Structure Collectively, the models express uncertainty (or disagreement) about whether harvest is an additive or compensatory form of mortality (Burnham et al. 1984), and whether the reproductive process is weakly or strongly densitydependent (i.e., the degree to which reproductive rates decline with increasing population size). All population models for midcontinent mallards share a common “balance equation” to predict changes in breedingpopulation size as a function of annual survival and reproductive rates: ()()()()NNmSmSRSStttAMtAFttJFtJMFsumMsum+=+−++11,,,,φφ where: N = breeding population size, m = proportion of males in the breeding population, SAM, SAF, SJF, and SJM = survival rates of adult males, adult females, young females, and young males, respectively, R = reproductive rate, defined as the fall age ratio of females, φφFsumMsum= the ratio of female (F) to male (M) summer survival, and t = year. We assumed that m and φφFsumMsumare fixed and known. We also assumed, based in part on information provided by Blohm et al. (1987), the ratio of female to male summer survival was equivalent to the ratio of annual survival rates in the absence of harvest. Based on this assumption, we estimated φφFsumMsum = 0.897. To estimate m we expressed the balance equation in matrix form: NNSRSSRSNNtAMtAFAMJMFsumMsumAFJFtAMtAF++⎡⎣⎢⎤⎦⎥=+⎡⎣⎢⎤⎦⎥⎡⎣⎢⎤⎦⎥110,,,,φφ and substituted the constant ratio of summer survival and means of estimated survival and reproductive rates. The right eigenvector of the transition matrix is the stable sex structure that the breeding population eventually would attain with these constant demographic rates. This eigenvector yielded an estimate of m = 0.5246. Using estimates of annual survival and reproductive rates, the balance equation for midcontinent mallards overpredicted observed population sizes by 11.0% on average. The source of the bias is unknown, so we modified the balance equation to eliminate the bias by adjusting both survival and reproductive rates: 33 ()()()()NNmSmSRSStSttAMtAFRttJFtJMFsumMsum+=+−++11γγ,,,, where γ denotes the biascorrection factors for survival (S) and reproduction (R). We used a least squares approach to estimate γS = 0.9407 and γR = 0.8647. Survival Process We considered two alternative hypotheses for the relationship between annual survival and harvest rates. For both models, we assumed that survival in the absence of harvest was the same for adults and young of the same sex. In the model where harvest mortality is additive to natural mortality: ()SsKtsexagesexAtsexage,,,,,=−01 and in the model where changes in natural mortality compensate for harvest losses (up to some threshold): SsifKKifKstsexagesexCtsexagesexCtsexagetsexagesexC,,,,,,,,,,=≤−−>⎧⎨⎪⎩⎪000111 where s0 = survival in the absence of harvest under the additive (A) or compensatory (C) model, and K = harvest rate adjusted for crippling loss (20%, Anderson and Burnham 1976). We averaged estimates of s0 across banding reference areas by weighting by breedingpopulation size. For the additive model, s0 = 0.7896 and 0.6886 for males and females, respectively. For the compensatory model, s0 = 0.6467 and 0.5965 for males and females, respectively. These estimates may seem counterintuitive because survival in the absence of harvest should be the same for both models. However, estimating a common (but still sexspecific) s0 for both models leads to alternative models that do not fit available bandrecovery data equally well. More importantly, it suggests that the greatest uncertainty about survival rates is when harvest rate is within the realm of experience. By allowing s0 to differ between additive and compensatory models, we acknowledge that the greatest uncertainty about survival rate is its value in the absence of harvest (i.e., where we have no experience). Reproductive Process Annual reproductive rates were estimated from age ratios in the harvest of females, corrected using a constant estimate of differential vulnerability. Predictor variables were the number of ponds in May in Prairie Canada (P, in millions) and the size of the breeding population (N, in millions). We estimated the bestfitting linear model, and then calculated the 80% confidence ellipsoid for all model parameters. We chose the two points on this ellipsoid with the largest and smallest values for the effect of breedingpopulation size, and generated a weakly densitydependent model: RPtt=+−071660108300373... and a strongly densitydependent model: RPtt= +−113900137601131... Predicted recruitment was then rescaled to reflect the current definition of midcontinent mallards which now excludes birds from Alaska but includes mallards observed in the Great Lakes region. 34 Pond Dynamics We modeled annual variation in Canadian pond numbers as a firstorder autoregressive process. The estimated model was: 35 t PPtt+=++12212703420..ε where ponds are in millions and εtis normally distributed with mean = 0 and variance = 1.2567. Variance of Prediction Errors Using the balance equation and submodels described above, predictions of breedingpopulation size in year t+1 depend only on specification of population size, pond numbers, and harvest rate in year t. For the period in which comparisons were possible, we compared these predictions with observed population sizes. We estimated the predictionerror variance by setting: ()()()()()[]()eNNeNNNnttobstprettobstpret=−=−Σlnln~,$lnlnthen assumingand estimating012σσ22 where Nobs and Npre are observed and predicted population sizes (in millions), respectively, and n = the number of years being compared. We were concerned about a variance estimate that was too small, either by chance or because the number of years in which comparisons were possible was small. Therefore, we calculated the upper 80% confidence limit for σ2 based on a Chisquared distribution for each combination of the alternative survival and reproductive submodels, and then averaged them. The final estimate of σ2 was 0.0280, equivalent to a coefficient of variation of about 18%. Model Implications The population model with additive hunting mortality and weakly densitydependent recruitment (SaRw) leads to the most conservative harvest strategy, whereas the model with compensatory hunting mortality and strongly densitydependent recruitment (ScRs) leads to the most liberal strategy. The other two models (SaRs and ScRw) lead to strategies that are intermediate between these extremes. Under the models with compensatory hunting mortality (ScRs and ScRw), the optimal strategy is to have a liberal regulation regardless of population size or number of ponds because at harvest rates achieved under the liberal alternative, harvest has no effect on population size. Under the strongly densitydependent model (ScRs), the density dependence regulates the population and keeps it within narrow bounds. Under the weakly density dependent model (ScRw), the densitydependence does not exert as strong a regulatory effect, and the population size fluctuates more. Model Weights Model weights are calculated as Bayesian probabilities, reflecting the relative ability of the individual alternative models to predict observed changes in population size. The Bayesian probability for each model is a function of the model’s previous (or prior) weight and the likelihood of the observed population size under that model. We used Bayes’ theorem to calculate model weights from a comparison of predicted and observed population sizes for the years 1996–2009, starting with equal model weights in 1995. APPENDIX 3: Eastern Mallard Models Model Structure We also revised the population models for eastern mallards in 2002 (Johnson et al. 2002a, USFWS 2002). The current set of six models: (1) relies solely on federal and state waterfowl surveys (rather than the Breeding Bird Survey) to estimate abundance; (2) allows for the possibility of a positive bias in estimates of survival or reproductive rates; (3) incorporates competing hypotheses of strongly and weakly densitydependent reproduction; and (4) assumes that hunting mortality is additive to other sources of mortality. As with midcontinent mallards, all population models for eastern mallards share a common balance equation to predict changes in breedingpopulation size as a function of annual survival and reproductive rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅+−⋅+⋅⋅+⋅⋅⋅11ψ where: N = breedingpopulation size, p = proportion of males in the breeding population, Sam, Saf, Sym, and Syf = survival rates of adult males, adult females, young males, and young females, respectively, Am = ratio of young males to adult males in the harvest, d = ratio of young male to adult male direct recovery rates, ψ = the ratio of male to female summer survival, and t = year. In this balance equation, we assume that p, d, and ψ are fixed and known. The parameter ψ is necessary to account for the difference in anniversary date between the breedingpopulation survey (May) and the survival and reproductive rate estimates (August). This model also assumes that the sex ratio of fledged young is 1:1; hence Am/d appears twice in the balance equation. We estimated d = 1.043 as the median ratio of young:adult male bandrecovery rates in those states from which wing receipts were obtained. We estimated ψ = 1.216 by regressing through the origin estimates of male survival against female survival in the absence of harvest, assuming that differences in natural mortality between males and females occur principally in summer. To estimate p, we used a population projection matrix of the form: ()()⎥⎦⎤⎢⎣⎡⋅⎥⎦⎤⎢⎣⎡⋅⋅⋅+=⎥⎦⎤⎢⎣⎡++ttafyfmymmamttFMSSdASdASFMψ011 where M and F are the relative number of males and females in the breeding populations, respectively. To parameterize the projection matrix we used average annual survival rate and age ratio estimates, and the estimates of d and ψ provided above. The right eigenvector of the projection matrix is the stable proportion of males and females the breeding population eventually would attain in the face of constant demographic rates. This eigenvector yielded an estimate of p = 0.544. We also attempted to determine whether estimates of survival and reproductive rates were unbiased. We relied on the balance equation provided above, except that we included additional parameters to correct for any bias that might exist. Because we were unsure of the source(s) of potential bias, we alternatively assumed that any bias resided solely in survival rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅⋅+−⋅+⋅⋅+⋅⋅⋅11Ωψ 36 (where Ω is the biascorrection factor for survival rates), or solely in reproductive rates: ()()()()()()()()NNpSpSpAdSpAdStttamtaftmtymtmtyf+=⋅⋅+−⋅+⋅⋅⋅+⋅⋅⋅⋅11αα (where α is the biascorrection factor for reproductive rates). We estimated Ω and α by determining the values of these parameters that minimized the sum of squared differences between observed and predicted population sizes. Based on this analysis, Ω = 0.836 and α = 0.701, suggesting a positive bias in survival or reproductive rates. However, because of the limited number of years available for comparing observed and predicted population sizes, we also retained the balance equation that assumes estimates of survival and reproductive rates are unbiased. Survival Process For purposes of AHM, annual survival rates must be predicted based on the specification of regulationspecific harvest rates (and perhaps on other uncontrolled factors). Annual survival for each age (i) and sex (j) class under a given regulatory alternative is: ()()Shvctijjtamij,,=⋅−⋅−⎛⎝⎜⎜⎞⎠⎟⎟θ11 where: S = annual survival, jθ= mean survival from natural causes, ham = harvest rate of adult males, v = harvest vulnerability relative to adult males, and c = rate of crippling (unretrieved harvest). This model assumes that annual variation in survival is due solely to variation in harvest rates, that relative harvest vulnerability of the different agesex classes is fixed and known, and that survival from natural causes is fixed at its sample mean. We estimated jθ= 0.7307 and 0.5950 for males and females, respectively. Reproductive process As with survival, annual reproductive rates must be predicted in advance of setting regulations. We relied on the apparent relationship between breedingpopulation size and reproductive rates: ()RabNtt=⋅⋅exp where Rt is the reproductive rate (i.e., dAmt), Nt is breedingpopulation size in millions, and a and b are model parameters. The leastsquares parameter estimates were a = 2.508 and b = 0.875. Because of both the importance and uncertainty of the relationship between population size and reproduction, we specified two alternative models in which the slope (b) was fixed at the leastsquares estimate ± one standard error, and in which the intercepts (a) were subsequently reestimated. This provided alternative hypotheses of strongly densitydependent (a = 4.154, b = 1.377) and weakly densitydependent reproduction (a = 1.518, b = 0.373). 37 Variance of Prediction Errors Using the balance equations and submodels provided above, predictions of breedingpopulation size in year t+1 depend only on the specification of a regulatory alternative and on an estimate of population size in year t. For the period in which comparisons were possible (1991–96), we were interested in how well these predictions corresponded with observed population sizes. In making these comparisons, we were primarily concerned with how well the biascorrected balance equations and reproductive and survival submodels performed. Therefore, we relied on estimates of harvest rates rather than regulations as model inputs. We estimated the predictionerror variance by setting: ()()()()()[]eNNeNNNttobstprettobstpret=−=−Σlnln~,$lnlnthen assumingand estimating02σσ22 where Nobs and Npre are observed and predicted population sizes (in millions), respectively, and n = 6. Variance estimates were similar regardless of whether we assumed that the bias was in reproductive rates or in survival, or whether we assumed that reproduction was strongly or weakly densitydependent. Thus, we averaged variance estimates to provide a final estimate of σ2 = 0.006, which is equivalent to a coefficient of variation (CV) of 8.0%. We were concerned, however, about the small number of years available for estimating this variance. Therefore, we estimated an 80% confidence interval for σ2 based on a Chisquared distribution and used the upper limit for σ2 = 0.018 (i.e., CV = 14.5%) to express the additional uncertainty about the magnitude of prediction errors attributable to potentially important environmental effects not expressed by the models. Model Implications Modelspecific regulatory strategies based on the hypothesis of weakly densitydependent reproduction are considerably more conservative than those based on the hypothesis of strongly densitydependent reproduction. The three models with weakly densitydependent reproduction suggest a carrying capacity (i.e., average population size in the absence of harvest) >2.0 million mallards, and prescribe extremely restrictive regulations for population size <1.0 million. The three models with strongly densitydependent reproduction suggest a carrying capacity of about 1.5 million mallards, and prescribe liberal regulations for population sizes >300 thousand. Optimal regulatory strategies are relatively insensitive to whether models include a bias correction or not. All modelspecific regulatory strategies are “knifeedged,” meaning that large differences in the optimal regulatory choice can be precipitated by only small changes in breedingpopulation size. This result is at least partially due to the small differences in predicted harvest rates among the current regulatory alternatives (see the section on Regulatory Alternatives later in this report). Model Weights We used Bayes’ theorem to calculate model weights from a comparison of predicted and observed population sizes for the years 1996–2009. We calculated weights for the alternative models based on an assumption of equal model weights in 1996 (the last year data was used to develop most model components) and on estimates of yearspecific harvest rates (Appendix 5). 38 APPENDIX 4: Western Mallard Models In contrast to midcontinent and eastern mallards, we did not model changes in population size of western mallards as an explicit function of survival and reproductive rate estimates (which in turn may be functions of harvest and environmental covariates). We believed this socalled “balanceequation approach” was not viable for western mallards because of insufficient banding in Alaska to estimate survival rates, and because of the difficulty in estimating stockspecific fall age ratios from a sample of wings derived from a mix of breeding stocks. We therefore relied on a discrete logistic model (Schaefer 1954), which combines reproduction and natural mortality into a single parameter r, the intrinsic rate of growth. The model assumes densitydependent growth, which is regulated by the ratio of population size, N, to the carrying capacity of the environment, K (i.e., equilibrium population size in the absence of harvest). In the traditional formulation, harvest mortality is additive to other sources of mortality, but compensation for hunting losses can occur through subsequent increases in production. However, we parameterized the model in a way that also allows for compensation of harvest mortality between the hunting and breeding seasons. It is important to note that compensation modeled in this way is purely phenomenological, in the sense that there is no explicit ecological mechanism for compensation (e.g., densitydependent mortality after the hunting season). The basic model for both the Alaska and California–Oregon stocks had the form: ()AMttttttthdwhereKNrNNN⋅=−⎥⎦⎤⎢⎣⎡⎟⎠⎞⎜⎝⎛−+=+αα111 and where t = year, hAM = the harvest rate of adult males, and d = a scaling factor. The scaling factor is used to account for a combination of unobservable effects, including unretrieved harvest (i.e., crippling loss), differential harvest mortality of cohorts other than adult males, and for the possibility that some harvest mortality may not affect subsequent breedingpopulation size (i.e., the compensatory mortality hypothesis). Estimation Framework We used Bayesian estimation methods in combination with a statespace model that accounts explicitly for both process and observation error in breeding population size. This combination of methods is becoming widely used in natural resource modeling, in part because it facilitates the fitting of nonlinear models that may have nonnormal errors (Meyer and Millar 1999). The Bayesian approach also provides a natural and intuitive way to portray uncertainty, allows one to incorporate prior information about model parameters, and permits the updating of parameter estimates as further information becomes available. We first scaled N by K as recommended by Meyer and Millar (1999), and assumed that process errors et were lognormally distributed with mean 0 and variance σ2. Thus, the process model had the form: ()()[](){}()21111,0~11loglogσNewhereehdPrPPPKNPttAMtttttttt+⋅−−+==−−−− The observation model related the unknown population sizes (PtK) to the population sizes (Nt) estimated from the breedingpopulation surveys in Alaska and CaliforniaOregon. We assumed that the observation process yielded additive, normally distributed errors, which were represented by: 39 BPOPtttKPNε+=, ),0(~2BPOPBPOPtNwhereσε. Use of the observation model allowed us to account for the sampling error in population estimates, while permitting us to estimate the process error, which reflects the inability of the model to completely describe changes in population size. The process error reflects the combined effect of misspecification of an appropriate model form, as well as any unmodeled environmental drivers. We initially examined a number of possible environmental covariates, including the Palmer Drought Index in California and Oregon, spring temperature in Alaska, and the El Niño Southern Oscillation Index (http://www.cdc.noaa.gov/people/klaus.wolter/MEI/mei.html). While the estimated effects of these covariates on r or K were generally what one would expect, they were never of sufficient magnitude to have a meaningful effect on optimal harvest strategies. We therefore chose not to further pursue an investigation of environmental covariates, and posited that the process error was a sufficient surrogate for these unmodeled effects. Parameterization of the models also required measures of harvest rate. Beginning in 2002, harvest rates of adult males were estimated directly from the recovery of reward bands. Prior to 1993, we used direct recoveries of standard bands, corrected for bandreporting rates provided by Nichols et al. (1995b). We also used the bandreporting rates provided by Nichols et al. (1995b) for estimating harvest rates in 1994 and 1995, except that we inflated the reporting rates of fulladdress and tollfree bands based on an unpublished analysis by Clint Moore and Jim Nichols (Patuxent Wildlife Research Center). We were unwilling to estimate harvest rates for the years 1996–2001 because of suspected, but unknown, increases in the reporting rates of all bands. For simplicity, harvest rate estimates were treated as known values in our analysis, although future analyses might benefit from an appropriate observation model for these data. In a Bayesian analysis, one is interested in making probabilistic statements about the model parameters (θ), conditioned on the observed data. Thus, we are interested in evaluating P(θ data), which requires the specification of prior distributions for all model parameters and unobserved system states (θ) and the sampling distribution (likelihood) of the observed data P(data θ). Using Bayes theorem, we can represent the posterior probability distribution of model parameters, conditioned on the data, as: ) ()() (θθθdataPPdataP×∝. Accordingly, we specified prior distributions for model parameters r, K, d, and P0, which is the initial population size relative to carrying capacity. For both stocks, we specified the following prior distributions for r, d, and σ2: ()()()001.0,001.0~2,0~4427.1,0397.1~2gammaInverseUniformdnormalLogr−−−σ The prior distribution for r is centered at 0.35, which we believe to be a reasonable value for mallards based on lifehistory characteristics and estimates for other avian species. Yet the distribution also admits considerable uncertainty as to the value of r within what we believe to be realistic biological bounds. As for the harvestrate scalar, we would expect d ≥ 1 under the additive hypothesis and d < 1 under the compensatory hypothesis. As we had no data to specify an informative prior distribution, we specified a vague prior in which d could take on a wide range of values with equal probability. We used a traditional, uninformative prior distribution for σ2. Prior distributions for K and P0 were stockspecific and are described in the following sections. 40 We used the publicdomain software WinBUGS (http://www.mrcbsu.cam.ac.uk/bugs/) to derive samples from the joint posterior distribution of model parameters via MarkovChain Monte Carlo (MCMC) simulations. We obtained 510,000 samples from the joint posterior distribution, discarded the first 10,000, and then thinned the remainder by 50, resulting in a final sample of 10,000. Alaska mallards Data selection.Breeding population estimates of mallards in Alaska (and the Old Crow Flats in Yukon) are available since 1955 in WBPHS strata 1–12 (Smith 1995). However, a change in survey aircraft in 1977 instantaneously increased the detectability of waterfowl, and thus population estimates (Hodges et al. 1996). Moreover, there was a rapid increase in average annual temperature in Alaska at the same time, apparently tied to changes in the frequency and intensity of El Niño events (http://www.cdc.noaa.gov/people/klaus.wolter/MEI/mei.html). This confounding of changes in climate and survey methods led us to truncate the years 1955–1977 from the time series of population estimates. Modeling of the Alaska stock also depended on the availability of harvestrate estimates derived from bandrecovery data. Unfortunately, sufficient numbers of mallards were not banded in Alaska prior to 1990. A search for covariates that would have allowed us to make harvestrate predictions for years in which bandrecovery data were not available was not fruitful, and we were thus forced to further restrict the time series to 1990–2005. Even so, harvest rate estimates were not available for the years 1996–2001 because of unknown changes in bandreporting rates. Because available estimates of harvest rate showed no apparent variation over time, we simply used the mean and standard deviation of the available estimates and generated independent samples of predictions for the missing years based on a logit transformation and an assumption of normality: ()200119960830.0,3265.2~1ln−=−⎟⎟⎠⎞⎜⎜⎝⎛−tforNormalhhtt Prior distributions for K and P0.—We believed that sufficient information was available to use mildly informative priors for K and P0. In recent years the Alaska stock has contained approximately 0.8 million mallards. If harvest rates have been comparable to that necessary to achieve maximum sustained yield (MSY) under the logistic model (i.e., r/2), then we would expect K ≈ 1.6 million. On the other hand, if harvest rates have been less than those associated with MSY, then we would expect K < 1.6 million. Because we believed it was not likely that harvest rates were >r/2, we believed the likely range of K to be 0.8–1.6 million. We therefore specified a prior distribution that had a mean of 1.4 million, but had a sufficiently large variance to admit a wide range of possible values: ()41224.0,13035.0~LognormalK Extending this line of reasoning, we specified a prior distribution that assumed the estimated population size of approximately 0.4 million at the start of the timeseries (i.e., 1990) was 20–60% of K. Thus on a log scale: ()5108.0,6094.1~−−UniformPo Parameter estimates.—The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed timeseries of population estimates. The posterior means of K and r were similar to their priors, although their variances were considerably smaller (Table 1). However, the posterior distribution of d was essentially the same as its prior, reflecting the absence of information in the data necessary to reliably estimate this parameter. Table 1. Estimates of model parameters resulting from fitting a discrete logistic model with MCMC to a time 41 series of estimated population sizes and harvest rates of mallards breeding in Alaska, 1990–2008. Parameter Mean SD 95% credibility interval K 1.125 0.323 0.659–1.890 P0 0.342 0.097 0.207–0.560 d 1.066 0.541 0.091–1.946 r 0.310 0.131 0.095–0.594 σ2 0.024 0.013 0.008–0.057 California–Oregon mallards Data selection.—Breedingpopulation estimates of mallards in California are available starting in 1992, but not until 1994 in Oregon. Also, Oregon did not conduct a survey in 2001. To avoid truncating the time series, we used the admittedly weak relationship (P = 0.11) between California and Oregon population estimates to predict population sizes in Oregon in 1992, 1993, and 2001. The fitted linear model was: ()CAtORtNN0848.074486+= To derive realistic standard errors, we assumed that the predictions had the same mean coefficient of variation as the years when surveys were conducted (n = 14, CV = 0.082). The estimated sizes and variances of the California–Oregon stock were calculated by simply summing the statespecific estimates. We pooled banding and recovery data for California and Oregon and estimated harvest rates in the same manner as that for Alaska mallards. Although banded sample sizes were sufficient in all years, harvest rates could not be estimated for the years 1996–2001 because of unknown changes in bandreporting rates. As with Alaska, available estimates of harvest rate showed no apparent trend over time, and we simply used the mean and standard deviation of the available estimates and generated independent samples of predictions for the missing years based on a logit transformation and an assumption of normality: ()200119960335.0,9633.1~1ln−=−⎟⎟⎠⎞⎜⎜⎝⎛−tforNormalhhtt Prior distributions for K and P0.— Unlike the Alaska stock, the California–Oregon population has been relatively stable with a mean of 0.48 million mallards. We believed K should be in the range 0.48 – 0.96 million, assuming the logistic model and that harvest rates were ≤ r/2. We therefore specified a prior distribution on K that had a mean of 0.7 million, but with a variance sufficiently large to admit a wide range of possible values: ()41224.0,5628.0~−LognormalK The estimated size of the California–Oregon stock was 0.48 million at the start of the timeseries (i.e., 1992). We used a similar line of reasoning as that for Alaska for specifying a prior distribution P0, positing that initial population size was 40–100% of K. Thus on a log scale: ()0.0,9163.0~−UniformPo Parameter estimates.—The logistic model and associated posterior parameter estimates provided a reasonable fit to the observed time series of population estimates. The posterior means of K and r were similar to their priors, although the variances were considerably smaller (Table 2). Interestingly, the posterior mean of d was <1, 42 43 suggestive of a compensatory response to harvest; however the standard deviation of the estimate was large, with the upper 95% credibility limit >1. Table 2. Estimates of model parameters resulting from fitting a discrete logistic model with MCMC to a timeseries of estimated population sizes and harvest rates of mallards breeding in California and Oregon, 1992–2008. Parameter Mean SD 95% credibility interval K 0.655 0.178 0.445–1.117 P0 0.741 0.161 0.429–0.986 d 0.599 0.410 0.038–1.588 r 0.335 0.220 0.067–0.883 σ2 0.013 0.012 0.001–0.044 APPENDIX 5: Modeling Mallard Harvest Rates Midcontinent We modeled harvest rates of midcontinent mallards within a Bayesian hierarchical framework. We developed a set of models to predict harvest rates under each regulatory alternative as a function of the harvest rates observed under the liberal alternative, using historical information. We modeled the probability of regulationspecific harvest rates (h) based on normal distributions with the following parameterizations: ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LfLLMfLMMRLRRCCCNhpNhpNhpNhpυδμυδμγυμγυμ++ For the restrictive and moderate alternatives we introduced the parameter γ to represent the relative difference between the harvest rate observed under the liberal alternative and the moderate or restrictive alternatives. Based on this parameterization, we are making use of the information that has been gained (under the liberal alternative) and are modeling harvest rates for the restrictive and moderate alternatives as a function of the mean harvest rate observed under the liberal alternative. For the harvestrate distributions assumed under the restrictive and moderate regulatory alternatives, we specified that γR and γM are equal to the prior estimates of the predicted mean harvest rates under the restrictive and moderate alternatives divided by the prior estimates of the predicted mean harvest rates observed under the liberal alternative. Thus, these parameters act to scale the mean of the restrictive and moderate distributions in relation to the mean harvest rate observed under the liberal regulatory alternative. We also considered the marginal effect of frameworkdate extensions under the moderate and liberal alternatives by including the parameter δf. To update the probability distributions of harvest rates realized under each regulatory alternative, we first needed to specify a prior probability distribution for each of the model parameters. These distributions represent prior beliefs regarding the relationship between each regulatory alternative and the expected harvest rates. We used a normal distribution to represent the mean and a scaled inversechisquare distribution to represent the variance of the normal distribution of the likelihood. For the mean (μ) of each harvestrate distribution associated with each regulatory alternative, we use the predicted mean harvest rates provided in USFWS (2000:13–14), assuming uniformity of regulatory prescriptions across flyways. We set prior values of each standard deviation (ν) equal to 20% of the mean (CV = 0.2) based on an analysis by Johnson et al. (1997). We then specified the following prior distributions and parameter values under each regulatory package: Closed (in U.S. only): )0018.0,6(~)()60018.0,0088.0(~)(2222χνμ−InvScaledpNpCC These closedseason parameter values are based on observed harvest rates in Canada during the 1988–93 seasons, which was a period of restrictive regulations in both Canada and the United States. For the restrictive and moderate alternatives, we specified that the standard error of the normal distribution of the scaling parameter is based on a coefficient of variation for the mean equal to 0.3. The scale parameter of the inversechisquare distribution was set equal to the standard deviation of the harvest rate mean under the restrictive and moderate regulation alternatives (i.e., CV = 0.2). 44 Restrictive: )0133.0,6(~)()615.0,51.0(~)(2222χνγ−InvScaledpNpRR Moderate: )0222.0,6(~)()626.0,85.0(~)(2222χνγ−InvScaledpNpRM Liberal: )0261.0,6(~)()60261.0,1305.0(~)(2222χνμ−InvScaledpNpRL The prior distribution for the marginal effect of the frameworkdate extension was specified as: )01.0,02.0(~)(2Npfδ The prior distributions were multiplied by the likelihood functions based on the last seven years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided in Table 1. Table 1. Parameter estimates for predicting midcontinent mallard harvest rates resulting from a hierarchical, Bayesian analysis of midcontinent mallard banding and recovery information from 1998 to 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0088 0.0021 h1998 0.1021 0.0069 Cν 0.0019 0.0005 h1999 0.0981 0.0071 Rγ 0.5107 0.0622 h2000 0.1249 0.0083 Rν 0.0129 0.0033 h2001 0.0923 0.0088 Mγ 0.8525 0.1066 h2002 0.1128 0.0062 Mν 0.0215 0.0055 h2003 0.1131 0.0069 μL 0.1103 0.0067 h2004 0.1200 0.0109 Lν 0.0198 0.0035 h2005 0.1157 0.0085 fδ 0.0085 0.0075 h2006 0.1094 0.0078 h2007 0.1034 0.0066 h2008 0.1112 0.0065 45 Eastern We modeled harvest rates of eastern mallards using the same parameterizations as those for midcontinent mallards: ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LfLLMfLMMRLRRCCCNhpNhpNhpNhpυδμυδμγυμγυμ++ We set prior values of each standard deviation (ν) equal to 30% of the mean (CV = 0.3) to account for additional variation due to changes in regulations in the other Flyways and their unpredictable effects on the harvest rates of eastern mallards. We then specified the following prior distribution and parameter values for the liberal regulatory alternative: Closed (in U.S. only): )024.0,6(~)()6024.0,08.0(~)(2222χνμ−InvScaledpNpCC Restrictive: )0404.0,6(~)()6228.0,76.0(~)(2222χνγ−InvScaledpNpRR Moderate: )0488.0,6(~)()628.0,92.0(~)(2222χνγ−InvScaledpNpRM Liberal: )0531.0,6(~)()60531.0,1771.0(~)(2222χνμ−InvScaledpNpRL A previous analysis suggested that the effect of the frameworkdate extension on eastern mallards would be of lower magnitude and more variable than on midcontinent mallards (USFWS 2000). Therefore, we specified the following prior distribution for the marginal effect of the frameworkdate extension for eastern mallards as: )01.0,01.0(~)(2Npfδ 46 The prior distributions were multiplied by the likelihood functions based on the last four years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided in Table 2. Table 2. Parameter estimates for predicting eastern mallard harvest rates resulting from a hierarchical, Bayesian analysis of eastern mallard banding and recovery information from 2002 to 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0797 0.0257 h2002 0.1617 0.0125 Cν 0.0233 0.0061 h2003 0.1456 0.0105 Rγ 0.7619 0.0919 h2004 0.1360 0.0114 Rν 0.0393 0.0100 h2005 0.1308 0.0119 Mγ 0.9217 0.1145 h2006 0.1034 0.0132 Mν 0.0473 0.0118 h2007 0.1209 0.0134 μL 0.1459 0.0153 h2008 0.1206 0.0120 Lν 0.0433 0.0087 fδ 0.0034 0.0096 Western We modeled harvest rates of western mallards using a similar parameterization as that used for midcontinent and eastern mallards. However, we did not explicitly model the effect of the framework date extension because we did not use data observed prior to when framework date extensions were available. In the western mallard parameterization, the effect of the framework date extensions are implicit in the expected mean harvest rate expected under the liberal regulatory option. ),(~)(:Liberal),(~)(:Moderate),(~)(:eRestrictiv),(~)(:Closed2222LLLMLMMRLRRCCCNhpNhpNhpNhpυμυμγυμγυμ We set prior values of each standard deviation (ν) equal to 30% of the mean (CV = 0.3) to account for additional variation due to changes in regulations in the other Flyways and their unpredictable effects on the harvest rates of western mallards. We then specified the following prior distribution and parameter values for the liberal regulatory alternative: Closed (in U.S. only): )00264.0,6(~)()600264.0,01.0(~)(2222χνμ−InvScaledpNpCC Restrictive: 47 )01683.0,6(~)()6153.0,51.0(~)(2222χνγ−InvScaledpNpRR Moderate: )02805.0,6(~)()6255.0,85.0(~)(2222χνγ−InvScaledpNpRM Liberal: )033.0,6(~)()6033.0,11.0(~)(2222χνμ−InvScaledpNpRL The prior distributions were multiplied by the likelihood functions based on the last four years of data under liberal regulations, and the resulting posterior distributions were evaluated with Markov Chain Monte Carlo simulation. Posterior estimates of model parameters and of annual harvest rates are provided Table 3. Table 3. Parameter estimates for predicting western mallard harvest rates resulting from a hierarchical, Bayesian analysis of western mallard banding and recovery information from 2008. Parameter Estimate SD Parameter Estimate SD Cμ 0.0102 0.0187 h2008 0.1435 0.0243 Cν 0.0181 0.0045 Rγ 0.5103 0.0636 Rν 0.0164 0.0042 Mγ 0.8524 0.1025 Mν 0.0272 0.0069 μL 0.1155 0.0124 Lν 0.0318 0.0076 48 APPENDIX 6: Scaup Model We use a statespace formulation of scaup population and harvest dynamics within a Bayesian estimation framework (Meyer and Millar 1999, Millar and Meyer 2000). This analytical framework allows us to represent uncertainty associated with the monitoring programs (observation error) and the ability of our model formulation to predict actual changes in the system (process error). 8.1 Process Model Given a logistic growth population model that includes harvest (Schaefer 1954), scaup population and harvest dynamics are calculated as a function of the intrinsic rate of increase (r), carrying capacity (K), and harvest (Ht). Following Meyer and Millar (1999), we scaled population sizes by K (i.e., Pt = Nt/K) and assumed that process errors (εt) are lognormally distributed with a mean of 0 and variance. The state dynamics can be expressed as 2Processσ 197401974εePP= ()008,1975,...,2t,/)1(1111t=−−+=−−−−teKHPrPPPttttε where P0 is the initial ratio of population size to carrying capacity. To predict total scaup harvest levels, we modeled scaup harvest rates (ht) as a function of the pooled direct recovery rate (ft) observed each year with ./tttfhλ= We specified reporting rate (λt) distributions based on estimates for mallards (Anas platyrhynchos) from large scale historical and existing reward banding studies (Henny and Burham 1976, Nichols et. al. 1995b, P. Garrettson unpublished data). We accounted for increases in reporting rate believed to be associated with changes in band type (e.g., from AVISE and new address bands to 1800 toll free bands) by specifying year specific reporting rates according to 9961974,...,1 )04.0,38.0(~=tNormaltλ . 0081997,...,2 )04.0,70.0(~=tNormaltλ We then predicted total scaup harvest (Ht) with ()[].0081974,...,2t,1t=−+=KPrPPhHttt 8.2 Observation Model We compared our predictions of population and harvest numbers from our process model to the observations collected by the Waterfowl and Breeding Habitat Survey (WBPHS) and the Harvest Survey programs with the following relationships, assuming that the population and harvest observation errors were additive and normally distributed. May breeding population estimates were related to model predictions by ,2008,...,1974),,0(~where,2,==−tNKPNBPOPtBPOPtBPOPttObservedtσεε where is specified for each year with the variance estimates resulting from the WBPHS. 2,BPOPtσ 49 We adjusted our harvest predictions to the observed harvest data estimates with a scaling parameter (q) according to ()[]()).,0(~ where008,1974,...,2t,12,HarvesttHtHttttObservedNε/qKPrPPhHtσε==−+− We assumed that appropriate measures of the harvest observation error could be approximated by assuming a coefficient of variation for each annual harvest estimate equal to 0.15 (Paul Padding pers. comm.). The final component of the likelihood included the year specific direct recovery rates that were represented by the rate parameter (ft) of a Binomial distribution indexed by the total number of birds banded preseason and estimated with. 2,Harvesttσ ),(~,/ttttttfMBinomialmMmf= where mt is the total number of scaup banded preseason in year t and recovered during the hunting season in year t and Mt is the total number of scaup banded preseason in year t. 8.3 Bayesian Analysis Following Meyer and Millar (1999), we developed a fully conditional joint probability model, by first proposing prior distributions for all model parameters and unobserved system states and secondly by developing a fully conditional likelihood for each sampling distribution. Prior Distributions For this analysis, a joint prior distribution is required because the unknown system states P are assumed to be conditionally independent (Meyer and Millar 1999). This leads to the following joint prior distribution for the model parameters and unobserved system states ).,,,,, (), ()()()()()()()(),,,,,,,(2Pr21112Pr0102Pr,...,102ProcessntttttocessocessttTocessttfKrPPpPPpPpppfpqpKprpPPfqKrPσλσσλσλΠ=−−−×= In general, we chose noninformative priors to represent the uncertainty we have in specifying the value of the parameters used in our assessment. However, we were required to use existing information to specify informative priors for the initial ratio of population size to carrying capacity (P0) as well as the reporting rate values (λt) specified above that were used to adjust the direct recovery rate estimates to harvest rates. We specified that the value of P0, ranged from the population size at maximum sustained yield (P0 = NMSY/K = (K/2)/K = 0.5) to the carrying capacity (P0= N/K= 1), using a uniform distribution on the log scale to represent this range of values. We assumed that the exploitation experienced at this population state was somewhere on the righthand shoulder of a sustained yield curve (i.e., between MSY and K). Given that we have very little evidence to suggest that historical scaup harvest levels were limiting scaup population growth, this seems like a reasonable prior distribution. We used noninformative prior distributions to represent the variance and scaling terms, while the priors for the 50 population parameters r and K were chosen to be vague but within biological bounds. These distributions were specified according to P0 ~ Uniform(ln(0.5),0), K ~ Lognormal (2.17, 0.667), r ~ Uniform (0.00001, 2), ft ~ Beta(0.5,0.5), q~ Uniform(0.0, 2), ~ Inverse Gamma (0.001, 0.001). 2Processσ Likelihood We related the observed population, total harvest estimates, and observed direct recoveries to the model parameters and unobserved system states with the following likelihood function: )., (),,,,,, (),, (),,,,,,, ,,,(12121,...,122Pr,...,1,...,1,...,1,...,1ΠΠΠ===××=TtttHarvestTttttBPOPTtttTHarvestocessttTTTTfMmpqfKrPHpKPNpPqfKrMmHNPttσλσσσλ Posterior Evaluation Using Bayes theorem we then specified a posterior distribution for the fully conditional joint probability distribution of the parameters given the observed information according to ).,, (),,,,,, (),, (),,,,, (), ()()()()()()()(),,, ,,,,,,,(121212Pr21112Pr0102Pr,...,1,...,1,...,1,...,1,...,102PrtTtttHarvestTttttBPOPTtttocessntttttocessocessttTTTTTocess 
Tag  LibrarySourceBird_publications 
Date created  20130123 

