I IId‘ uuloo u... .....L 1... .-_... . . - .. .. w.....d.1...«t.... i». v... a L... . I . .. 1 ..3. I}... ... t I . urteaI .. 2: .1} .I I..?) . . ‘30:..I1 ’2. 7.0%..5 g‘OII‘ . O o ...-903‘. o..._n‘.o“. . a‘h‘qu‘fiv .QWanug‘.” . K. 1”,y if . I I.. {a 1‘; ......a :97; ....I... 1‘. . It ...... to o..l'QtD1¢-&.It£lat. II<-.¢.O.PV ..‘O k‘fl8ils .... .. I . .. ... . .... . of ...-...... , . .3... .3 ..3I.. ~t:I.?z..l-.:..O. .... ...). .t.:l...:.7_. 1b3_d8‘.9’.~3 I. . .20.!‘739 . . to... .. .. . 043.19... $.12: 9. .0 .0 ... . 1...... .r . .0... .......l .63.. ...l . 10"... .O..£u.ll~n,.h.!. . .-J\c.. . . . .. a} a... .... . I! OI...‘ 500:... O...3‘.J‘. .3 u .0. . . in..." us}! 9.» a:.0...9~0-.“.‘l. .00 v.‘380.‘kooo~o'\~‘uwc .... 9 0.. . .52... I 8.013.. ‘o 5 u 5' 0.1.2.: .0. . ... I! .. ton... .. at! .C ..t ..‘S. -‘ Int... 11. J In... -. (I..-606100. 0.": .80! 1...! o‘o .021; .. .zclafioli" .0 .s,.‘d.lsl d 3‘.“ II.I-.-. I p.. a... ‘0! I: O. uPJ.. .. I.. .o I» .0 ’ol‘.‘ .I ... .....o. ..I.0I.n'.‘ofi.0§v. . a... 00!. . Q3I_I.I.! . l‘ . - ‘.O.Iurf... Ch...v',~...‘nfh0.l on... 0 Au. ‘ o. .a..-!-I..l... "L. .‘ .... .... .. .3r....o . . . 98.9.3203‘.:-.t£:..IS . ..‘Iél. lat?! lel.._. ,...c‘:.\cc...?..rv0r 0.. 1......373. . .. . . 25‘9013.’..la.19l‘- . o 7.70 «.... .......,..a"..ot\. . . I.. . ..... . . ... ... . II 3......00-‘30393 ..sC. .... I o .. .. .. . . . . . , .lv ....3Inv.‘~on. ...-’u-ce... .17. . s .. . ., . . . I. . ... _ . . .pocub“ .1. 36,3! I H“.uut’.t.£ . a... (31‘. Iisuhn-‘Wmogfi .u JO; A? b. . 1.0...- o D. I. . . Co. {...-3“ .Im.‘ .. I" o . .0. I .. .....I .12....”u...) . 3L? 03.9.5.» k‘.cdk$\3\..zoli b). . ... :4 o . , .....11... ~ ...-I. u. I“... .01! . 9.0L! .35... ... A Voolt.’ . u. .2!!! .. . .6 .(I .1 ‘2’. .0 lav}. .2...“ vi . fl.” . . V’I I ' ‘ . . ~ . . . 0“" 206......‘3‘ a. . lu..'9.o .... . . . ‘0:- IO. .. I.....Os.... 1".6ol. , ... . . ‘oé . I ......G . I . . . ... . .. . ... isles...) 40‘ O....‘...‘CQ§.I!‘ . Q. .... I.....i .9! u . .. .. . . I..... . I. ((8... .Uw.l‘ 5.03%! ..III.‘ . . . u. o 6...? ...... . C . . . 9 \‘33299ll . . . . _ .I’I.‘l ...-CI o\co.o0. l .12.. .33... 9.. . . .. . . . . ..IQ.I‘ .6. ...I! Q. I...) .33 0A,.- .) .....Iérou.‘ . . . . , .. .- .... .. . . . 8r . . I . o0. , . .. ...: ...-l :4. ...... 43.. o .l .0 3.4.; ... 0.3.184 - . Ik‘. I . . . . r . .0. v. I. , ...-n. .. ... 4.61.30... .0 ‘OI‘ 3""l Ol'lh. ' _ . . . . . .1. . . . . . . . - .1. .. . I... b... ...I. .. . . , ..‘Ior- .... 4:53.. ....‘I U I 0.. . . .‘ll- '3’;.II.\:.1 . O. . .. . . , . .. . . _ ..V... . . I ...: ...-9.9 13.3. I... ... . ...... . . ... . ... ...‘II 6:... 2"» . A ... .3'l\.. . . 5...! .-.. .0010?!) 4.... 93“! ”Z I . .u o. 9- ..~.. ~...- 09..../. 90.1. .. I. . I.....l. . 0'9. 0. b. I: . . . . . ....‘0. Im.‘coI-O.$l..| ¢ .IQV’IK..I.‘113:-.Ib v... IL 0 . ... 1" 9...: v. .... 30.. . ‘. ~ I . . . . . o... ..v I ... ..v..u.1... V . . .. . . . . .V...‘ ‘ I...0o.l... . 3.1. o, I. 01:? ow... .... 5“ o I . 6/7'0 I. . v ‘ fl . : 0.! .23... n .o.' . . . . . . . 1 . . . ... u at I of. .‘.-.. «.... . . . . I .1 2:33 I... is... O I It . I ‘0‘! “II. .I... it! .0.- ‘I W . ,I‘c.. .— to. c . fit... 4 .o .I. . . . .10 ‘u {I I... of: of ,0. I . .‘d~. l. .. I'...‘Q ...: I.. 0.50,...03 3‘ ‘c v , ... .. . . . .. . . . 9... HI. I...» 2 ..l‘ .- cos..- .3! . o a. £139 . .l .00 .IWolb t a... ...-.... . . I! I .0. n I? ... .I I I o .. on. Al I ...}...ot’o:“‘ IIZC‘. 1!. in ll 9L. 14""... n . f. i . . S .26.. .3. n .V .88... r 5.0! 3‘ . .7 . '3... I.&.~.o§‘ u. I .... ......I.t. ... o. .. ..{1 .... x ...-... .A . . LC‘ to!” 3-,}! .‘13 u .o a .. lo... .4 1. . i. I.. .. r}. .. I .; . .3... o 11. u o . ..I. l ...z‘ .204 . .lA-Qtoc ... I 3' .3: o‘...» . y'.- 0. ... ...... .7 .._.otlvg Z-I . I... ..I.' J 01.-.}.-.210013“ 4:... 0‘ ... I. ~01...’ :. I. . .w... 321 ... I... I 0 I ‘ .0... . .I .D o 19...... .39.... .... ‘I‘ .. , qudfll7 v :I.I..I _ .0 12" .0. 10.28.... .C o 3 .0. u». I. 3c. an... uflf . . . I. . 0", ul .. I.. .I‘u; . I . “I A. . O‘o.oOI..."... .3"... 1. .10 n Tit. ..1 av. . .3...“ .‘13!’. . .0. .PID.IDI~..3I31%.I~.1 ... . . 2 .IA‘ of. 0. 97.1....81‘ 3......itv o.l..' :3.- It.‘. at»? . ' r: .. ...o. .... .. .. . . . . . . ...!1 .... ...-1.9. ot‘fit‘: .‘I.....N..Il'0...‘. 3.4“. .I ..Iu...’ lili1.lup‘ .1!“ ‘ 0 o'. 79.3.... .... . . , . .. . I. . . . .. , . c 3. 17.. .0. .91.! ...a Q . o‘414I. a.."lOII.IiI..l...kl:fiu. d ..u... ivl“o‘.“. I , . . .v . . ... . I". A. . . . . . . . . . I. .o u. .A.. ... C‘ . .l... '3... . .xc‘ICx‘... 9:. o...- .‘ ....Lth-‘I ... ‘5».‘pn73‘u r I“ . .5. 1 o .. ... ... . . .. . ..l . . . . . . . . . . . no... 5...! .l.- ...-I.. A. I .3'k.01..:, ...! no. . ....o . ... 'hfl} {.l!x.!‘2iz‘a Oiéutill‘b .1 . .I .. .. to. .. .. ... I. .. . . . . . «9.9.5.... O. .... lo I..-...... .- .‘ ...! L20... ... Ill .0 4,}..ct3b ...-...cd-“Q‘.“$uu3.‘..l . .. I 2.. I. .. h. .... v . . . . . . . .. I. t ... 0‘3 . ...JXéO‘G 9.301.. .100... I 3:3... 5:». J .‘g..-. .... 1i .olnl...2.c ILtII!’. .h. v. «o .12 P.. . a.‘ .. wl. . . .. . . .. . . ...-I....l l '5 I. .90... at... . .k.‘.“.ha§ ‘a.ooll'\v.o."..o. o . .v .. ... O . . o . ..(o o . {I o. . .. . . . . . .. . . . . OI. )....$9. 4 :- “.€. 1'03. . u“. .. ..r‘. I O . I . I . . . . .. . i .... ...530100 0"“. _ . r a . . 9 , . : . . .. . .. .. , .. . . 2.5.5.1! .....;.o.l.3511I2. 3. co ....2 . r . . . o . . , I r _ ....L. . ., . .o . o v .7 . .. . o . . . v. . . . .I. .. I u . 9 I301 . ‘01- ...‘.v!.... . . . . . . 5.1.1 . LOO-I.....‘Ioi I.. .. ... . . . . - . n . . . . . . . . . u . .13'. . . ... . . , . . . . .. . . ‘c‘n... v0. . . . I 3...? o‘ I... . 17-310.. oo.‘ 10...! 3 . I.. 32235.! ‘- 10....3..t .’ .Z.¢.fr..a.:§§k. . .. ...... a. on-..) . ‘0; . .. f3|t ~Q’O ole ‘1‘ up . . II. . 39.1.3... . ...-3.... ...; . . . .1 211’s...“qu J\.:na§ ....9 .... :3»... l.‘.:o?.o I 7/ .9 . .. a:: .PI ’5 900.45 .....£, ...l’ 70.... ... . I.. .. .‘IJv... . . . . . ‘2‘)? .UQIIIA4’1‘IOI Q‘fiini. .- .. 3v .. . . . . . .. .... ..‘r . . .00... I‘C.....o€25. 9‘19 In... 30.0.12 . . I I I.. . . .. I. I A I. 1. .....‘Ir . .. o ..o 2... .83 . . . . . so. oki‘—'l.1 tIis.‘.fil\v-. ‘5.1 p '6- _. .I a r cl-.5....I. .0. .I o... ...! .53 ..I 10.. u. .. a. 7.1.1... I. . o 01.. DJ 1 1:. ...-Q. . :9. an. I}... .iIt: .‘l I. I . o . .. .. . . .... lo I. I. t... .I.‘ .. .ab v. :0. . I! I... ...Al ....4. . 3.‘Iltco“:v.:'~5 .I-. ‘3. ‘ 3. ..a I . r , . o. .33.. . ... ....o . ...: . .1: .. :1....3..!. . o. :7 P. I . . . . o... . ..I$§-...2.I.C’....l7\, .Péfii.i’1‘§3ku3'! .. .... .I . . .. . . _ . , . . .....a. . .. ...: All . I Is. 3. I .o \I..II....c - . . . .. ...II‘ I... o no. 20th $91. 5‘03“ 5-0.3..oél. 3 ‘70.). . ..a ._ ..c . .. . . . . . o . I. ...... o. o 2 . 3.'3.-..I.o..._. Io ... o .. -. I.....v . . . . ..I. \n-I‘Uuqo? If. coo. no.9.n - Q.o$'.‘0o..£ ustzc‘l‘u q 0‘ . .. . .. . , , . . . . .. .. . . .. .. ... ... 3.0... to. .o .... 2..., o. . .1 . ... . o It ...... I... .1. Did .... . . .. 1.. :(_’..’I.'F\|i&9.udoi 0.4:. .3..I~.‘3. 03.- .12 .I . .. . . . . .. . .. a . ......u .. ... . ... .. 3 ., .. I 0.34.. o..- ...::I. .. . '1! I... ...!‘vl I'anto.’_..l!...232..ilnnup it; h I .luv.. . _ . . . . . . . . .J .l.’,; In... 13lu. ...: ...); .53.; ..II. .II... 2...? . .32. .. . . . 3..-}: ..llnl‘n.....‘.l..I..\S1$l .. r . II I . . . . . . {.I .. . .0 . ..o. . I f . SI}... .82. o J. :60 .- £20 5}... 012.010 .. . 5...“! 3'0 o.€!$i"':’§.$t . .. 0.. ,‘II. . .,.. . 0.. ....40 :o .. . v ‘Q-‘. ‘.' . . .... . 3'... .o" . ..b‘? . 0.1. . «a..- ...-.4" ".OO.-LY n1..o..o . ... I. ... -. ... I.. . .9. II... -¢o- 13...... ..2. .33. . 2.1.. ...... I... . . . 1~ On I. I u 0. .. a. . ... .A. ... . 3|... 9. II. ... .- .. I o A I. . I I It _ 6 . . I U - ... . . I o . .9 I. I . . . I .6 I I o I _ s . o.\ o o o. ._ I o . u . . :5 n n I. I ~ I o I - t I I . o . . . D .... I... ‘ ~ .. g ‘I IA. 0 . . o o. a. o . .a . ... . v . . I o . . I u I o .I. .... . o . .. ... .I .9 no.0... sI. . I I . a, . O I . . .... I ...... ‘I. .II. 0 I. u . OI - o o I O . . c on. u v . . .. s c K Na .1 0 ‘ c v . I . . . O. I .. . . n . .... . . . . I.. . I ... .. o. -v 5. o. o . o . I .0 I I . I T O 5! I to . I .o I..... ......O 3.. .- I. . .. . . v . I ..I l . o . . I o . . a . .. . . . . I 6 I; 3‘" I . o . . .I. . . I I. .x . . 1 . ... . [Iii-all I...-Yl v a I I O I! ea; 0' . I I I . a I ....... . . ton”: O . . . | OI O o . . . _ I . I o I .o a. I. o. l. to. . o ...... . .‘I. IcVO . . ..v . 3...: 4 . .v . . . . . . . I. g I o a O I 0 I J . I . . n.0.’!ll' I . . . :34 ."IOOI...'OI"I'I.ID 3'0; 0 I I . . . not. 2.4.. ’aififieoz.‘z.o.0.$l§ . . . . .. . . . . ... . ..‘COEAI Iii... .I‘ngizé- . . ... . I . .. . . . . . £”15:!:‘ I 01.00.... I I. I . . . . . C . . , . . u . .. . . . {OIQtzlf ' I 32.0“?l" ": . . . . _ .. . ., . . . . . . . ..‘n \I! 03.01.100.00 iII§’¥.I..l-Irl . I .. . . . . . . .. It!!! 11.?0‘d Qlefit’I‘t’,‘.Llli: . _ . . . . . . . .oIa‘f1I-Rt‘cgf ‘i‘:“.‘ 7 I .. . . . .. . _ . . .. . .. A. .oI’i‘l.9!,.t9 .siovll.§b.llt I 0. a . 0 .I . ..zo-i,lafl”‘v‘d¢:" 1.0.x: . . . . 16.0 IQIIOI'IQQO‘..§’IOOI. :03... al' . 0...: 4 $600). ' Ilv‘tiov‘i..‘“'lo i‘foaz’ I .0 hi“’.0(":!".l‘" QiIIQ.v‘ .. .. . o . .- (6:.000'ic 31"}:i’VlI-‘Ooloill I 7 I 0"! S..‘I"§’.’"’O ’YI‘fIII' . I... :Q?:009nzla.if 2.3;.‘0’ I‘ll! I . I ’I...‘IO“."( I'i‘f‘gfiz' . 'r . .I .00 0“ wool-D'I‘IQGI’IItic I _ O. o 0“; . t . o 9 I . . . . ~ g‘:fl 5., II. I . 30‘ o S'ggof‘l.!.u£ tl‘IoO‘CII. 6.. ’O‘at; It... I 3329.};9‘. ‘0‘ a .( 1’.I‘$’al& -..olt‘nl‘a" 4,059.}: If? II.OI.CIIO‘IOII '3‘,‘|V‘1’£ ’ O.- ‘Oio n‘vO’I-.IPO:O‘ ‘40; $’. I .3 310.03! tall?!§‘3i’!.*fi ‘l. I . 4‘ 0 III 0.2096110“ ... .f a II:?IIQ.}I: no}. 0.0: $11!!! loot! «fl. (Iii..- if!o.‘l. - Ill-lootinl ....)‘QI'IoQIi .‘IMIA‘I. I‘D-PI\‘IQ - . C o (......Io' loo-IIIOI.I'O‘ICIO.I"‘ {f‘o O} I 04 I}: .21....0 .’O':‘O-|‘I‘I’"l613-.n 'OI‘ vs. .0 0‘ ‘r I10 0.. II I.|.o.!§'o‘.' .‘o 0‘: ‘IOoLI‘OC I I6 of.f.’II-.S$"I l g I..! 1......4 ’1]! . I . ,f".§!.0 :6 ...... 5.0"": II. ‘0 I....fngolII‘ 31!. one (I..-053......I' Ifi’x- lif,’..¢(1’4v0‘o‘:‘ I I x - I . . . '1‘..- . .. . I ... I. u I o 0:00... ..u 0.3: . no! I ... . , I. ...-0.0.0.0 I‘V‘ .II u . . I. OO.v.O.OoQ.I.O ... .....-.n.ac .90.... II’ . .0 f-..“tf.‘(l¢c"b‘ I-‘I0.IQ.‘3“IO ...! . .... .06)- ..u“ 0.0.0! .000 ll 0-1-06». 0.. I.. . o .o I I I I. O u v . . . I I . . I . v I. I o.’;. O . I p o o. . . . .. . . . . . I . .- .n.~ . .. . . ......o. ... I '(‘-.I L .....IIOthfi.‘ .6 C. I. ".16 .2 :1 .. b N . . ... .t¢‘ . . I ... I . .. I.f 3 ... 15$ _ 4.. 0%.. A... 0». {SI}. .... ...? .vLIw...’ .l...do$,..§.:t ...fl 816(23ttoihwcflzooxfo‘s go... It... .. 53:51.82 :2. ..n’lobkmg.&fiku 036...! 6’JVIYMT.5. l a. 0 Q:— .I . . ... .II. .0 . t l a O .C .c ..J .0 Jo... .35.;N‘» IC. In... t a4£D3 Ln. ...-.60....b a.) 5533‘... I.. ’0 4:»..Io-‘u. Q 5. ‘1‘:6"C.I 0.0.0.11. O ...l .. 3“33’ 0 a: I..... o. . 35:53 2 .,.o.II.Ib-..D I I :2. 6T 0 .I 6 ,ID . 1.1.0:. :0 b . cl... ..v»... 62...? o...o...!...tb:.b..o ..Ia. .0 25‘ .IleJ I .3 it 1032'... .8. po. 7 0.7..- ...lt..l...vt1.. I. . .. , . . . .L. {i . .. ’2. . I“ . . v 05.!" . 80...!.Gn: 4.... O. . I o ' a . . I . 0.. l . I .0 - ~ O ..o. C . 6.3:? v i ' I H H l —-\ l 3009 I LIBRARIES MICHIGAN STATE UNIVERSITY EAST LANSING, MI 48824-1048 This is to certify that the thesis entitled HOW MANY FISH ARE THERE AND HOW MANY CAN WE KILL? IMPROVING CATCH PER EFFORT INDICES OF ABUNDANCE AND EVALUATING HARVEST CONTROL RULES FOR LAKE WHITEFISH IN THE GREAT LAKES presented by Jonathan J. Deroba has been accepted towards fulfillment of the requirements for the PhD. degree in Fisheries and Wildlife, and Ecology, Evolutionary Biology, and Behavior rx )‘CNW’D— R (/éflflkfi...l (Major Professor’s Signature 11L“, VLL ( 7 ! (9‘ (Vb? Date MSU is an Affirmative Action/Equal Opportunity Employer PLACE IN RETURN BOX to remove this checkout from your record. To AVOID FINES return on or before date due. MAY BE RECALLED with earlier due date if requested. DATE DUE DATE DUE DATE DUE 5/08 KuProi/Achros/ClRC/DateDqundd HOW MANY FISH ARE THERE AND HOW MANY CAN WE KILL? IMPROVING- CATCH PER EFFORT INDICES OF ABUNDANCE AND EVALUATING HARVEST CONTROL RULES FOR LAKE WHITEFISH IN THE GREAT LAKES By Jonathan J. Deroba A DISSERTATION ‘ Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Fisheries and Wildlife, and Ecology, Evolutionary Biology, and Behavior 2009 ABSTRACT HOW MANY FISH ARE THERE AND HOW MANY CAN WE KILL? IMPROVING CATCH PER EFFORT INDICES OF ABUNDANCE AND EVALUATING HARVEST CONTROL RULES FOR LAKE WHITEFISH IN THE GREAT LAKES By Jonathan J. Deroba My dissertation has two main objectives: 1) to explore alternative ways to use commercial lake Whitefish fishery catch per effort (CPE) data as an index of abundance in 1836 Treaty-ceded waters of the Great Lakes, and 2) to evaluate alternative harvest control rules for lake Whitefish. Chapter 1 was directed at exploring alternative ways to use commercial lake Whitefish fishery CPE data, while Chapters 2 and 3 covered topics related to harvest control rules. Fishery CPE data is often used to assess relative fish abundance, and assessments used in 1836 Treaty-ceded waters of the Great Lakes assume that commercial CPE (i.e., ratio of aggregate catch to aggregate effort in each year) from gill-net and trap-net fisheries is proportional to abundance. However, CPE may change due to factors other than abundance. In Chapter 1, I developed general linear mixed models (GLMMS) to account for sources of variation in CPE unrelated to abundance, and used the least- squares means (LSMs) for each year as an alternative to the current index of abundance. Effects such as license holder, boat Size, and month accounted for much of the variation in CPE. LSMs and the current CPE index displayed different temporal trends among years in some areas, suggesting the importance of adjusting fishery CPE for effects like boat Size, season, and license holder. Harvest policies use control rules to dictate how fishing mortality or catch and yield levels are determined. Common control rules include constant catch, constant fishing mortality rate, and constant escapement. The “best” control rules for meeting common fishery objectives (e.g., maximizing yield) is a source of controversy in the literature, and results are seemingly contradictory. In Chapter 2, I conducted a detailed review of the relevant harvest control rule literature to compare control rules for their ability to meet widely used fishery objectives and identify potential causes for contradictory results. The relative performance of control rules at meeting common fishery objectives was affected by: fishery objectives, whether uncertainty in estimated stock sizes was included in analyses, whether the maximum recruitment level was varied in an autocorrelated fashion over time, how policy parameters were chosen, and the amount of compensation in the stock—recruit relationship. More research is needed to compare control rules while considering these and related factors. In Chapter 3, I used an age-structured simulation model that incorporated stochasticity in life history traits and multiple uncertainties to compare the current harvest control rule for lake Whitefish (constant fishing rate; CF) with a range of alternative control rules, including conditional constant catch (CCC), biomass-based (BB), and CF and BB rules with a 15% limit on the interannual change in the target catch. The CF and BB rules Simultaneously attained higher average yield and spawning stock biomass than other control rules, while the CCC rule and limiting the target catch changes by 15% had the lowest yearly variability in yield. The low yearly variability in yield provided by limiting target catch changes to 15% comes at the cost of frequently reducing biomass to low levels, so that in many situations other control rules would be preferred. ACKNOWLEDGEMENTS I am grateful to many people for helping me throughout my doctoral program. I especially want to thank my advisor, Dr. James Bence, for not only teaching me about quantitative fisheries analyses, but for also making the process fun. Jim was always available and helpful whenever I had questions, and I am honored to have studied with him and call him my fi'iend. I thank the other members of my graduate committee, Dr. William Taylor, Dr. Jean Tsao, and Dr. Michael Jones for greatly improving my research. The Modeling Sub—Committee of the 1836 Treaty Waters partially guided my research and answered many questions about the Great Lakes. In particular, Mark Ebener provided data, discussion, and entertainment. The Michigan Department of Natural Resources and US. Fish and Wildlife Service provided funding. A scholarship from the International Association of Great Lakes Research and a fellowship from the Department of Fisheries and Wildlife also made my life much easier. Several staff members of the Department of Fisheries and Wildlife also ensured that my degree progressed smoothly. Lastly, I would like to thank my family and friends. My parents, Bill and Judy Deroba, have supported me in countless ways and no amount of thanks could adequately express my gratitude for their hard work and sacrifices. My brother, Aaron, provided me with countless hours of fun and laughter. My grandparents, Bill and Judy Hoyer and Irene Deroba, were always available to listen and provide anything I needed. Julie Nieland helped maintain my sanity and made me a better person. My many friends from the QFC, Bence/Jones lab, M.I.R.T.H. lab, F W Department, related Departments, and MSU provided many necessary distractions. Thank you all. iv TABLE OF CONTENTS LIST OF TABLES .......................................................................................................... vii LIST OF FIGURES ......................................................................................................... ix INTRODUCTION AND SUMMARY ............................................................................. 1 Indices of Abundance .......................................................................................... 1 Chapter 1: Improving indices of abundance for lake Whitefish ................ 2 Harvest Control Rules ........................................................................................... 4 Chapter 2: A review of harvest control rules ............................................ 4 Chapter 3: Evaluating harvest control rules for lake Whitefish ................. 5 Overall Conclusions and Future Directions .......................................................... 8 References ........................................................................................................... 10 CHAPTER 1 DEVELOPING MODEL-BASED INDICES OF LAKE WHITEFISH ABUNDANCE USING COMMERCIAL FISHERY CATCH AND EFFORT DATA IN LAKES HURON, MICHIGAN AND SUPERIOR ...................................................................... 12 Abstract ............................................................................................................... 13 Introduction ......................................................................................................... 14 Methods ............................................................................................................... 17 Study area ................................................................................................ 17 Data and analyses .................................................................................... 18 Results ................................................................................................................. 24 Gill-Net Fishery ...................................................................................... 24 Trap-Net Fishery ..................................................................................... 25 Discussion ........................................................................................................... 26 Acknowledgements ............................................................................................. 34 References ........................................................................................................... 36 CHAPTER 2 A REVIEW OF HARVEST POLICIES: UNDERSTANDING RELATIVE PERFORMANCE OF CONTROL RULES ................................................................... 5 3 Abstract ............................................................................................................... 54 Introduction ......................................................................................................... 55 Common control rules ......................................................................................... 59 Common fishery objectives ................................................................................ 61 Relative performance with “perfect” information .............................................. 64 Comparing control rules ......................................................................... 64 Effect of the stock-recruit relationship ................................................... 70 Relative performance with “imperfect” information .......................................... 71 Policy parameters unadjusted for uncertainty ......................................... 72 V Uncertainty adjusted policy parameters .................................................. 74 Selecting catch, fishing mortality, and threshold levels ..................................... 76 Available options —Simulations or biological reference points ............... 76 Constant catch levels ............................................................................... 77 Constant fishing mortality rate F levels .................................................. 78 Threshold levels ...................................................................................... 80 Summary and conclusions .................................................................................. 82 Acknowledgements ............................................................................................. 87 References ........................................................................................................... 89 CHAPTER 3 EVALUATING HARVEST CONTROL RULES FOR LAKE WHITEFISH IN THE GREAT LAKES: ACCOUNTING FOR VARIABLE LIFE-HISTORY TRAITS ...... 108 Abstract ............................................................................................................. 108 Introduction ....................................................................................................... 109 Harvest control rules ............................................................................. 109 Lake Whitefish Fishery History and Management ............................... 109 Methods ............................................................................................................. 112 Stocks in 1836 Treaty waters ................................................................ 112 Overview of simulations ....................................................................... 112 Growth sub-models ............................................................................... 114 Stock-recruitment sub-model ................................................................ 116 Maturity sub-model ............................................................................... 118 Population abundance sub-model ......................................................... 119 Assessment and implementation error .................................................. 120 Control rules .......................................................................................... 121 Performance metrics ............................................................................. 124 Results ............................................................................................................... 126 Overview of results ............................................................................... 126 Rank-order performance of control rules .............................................. 127 Sensitivity to future growth uncertainty ............................................... 128 Sensitivity to assessment error parameters ........................................... 131 Discussion...............; ......................................................................................... 132 References ......................................................................................................... 1 3 7 APPENDIX A ............................................................................................................... 157 Classifying management units as fast or slow .................................................. 157 Growth sub-models ........................................................................................... 158 Stock-recruitment sub-model ............................................................................ 159 Maturity sub-model ........................................................................................... 161 Maximum sustainable yield and unfished SSB ................................................ 162 APPENDIX B ............................................................................................................... 166 COMPREHENSIVE LIST OF REFERENCES ........................................................... 183 Vi LIST OF TABLES CHAPTER 1 Table 1. Years of lake Whitefish catch and effort data analyzed from gill-net and trap-net fisheries conducted in 1836 Treaty of Washington-ceded waters of Lakes Superior (management units designated as WFS), Huron (WFH), and Michigan (WF M). .......... 49 Table 2. Difference (A) in the corrected Akaike’s information criterion (AICc) between the final lake Whitefish abundance model and a means model (i.e., a model with only a year effect; AAICc = means model AICc - final model AICc) for gill-net and trap-net fisheries conducted in 1836 Treaty-ceded waters of Lakes Superior (management units designated as WFS), Huron (WFH), and Michigan (WF M). ......................................... 50 Table 3. Average variance component (0' 2) estimates for a model Of lake Whitefish catch per effort (CPE) in the gill-net fishery (residual error = 0'5”" b g, ; license holder = of ; month x year = 030,) and for a model of trap-net fishery CPE (residual error = Grim]; year x license holder = ail ; month x year = 0,30,; month x year x license holder = 0' in”) for 1836 Treaty-ceded waters of Lakes Superior, Huron, and Michigan. Variance component estimates were averaged across management units. ..................... 51 Table 4. Average coefficient estimates for three boat size categories (small: <20 ft; medium: 20—30 R; large: >30 ft) in a model of lake Whitefish catch per effort for the gill- net fishery in 1836 Treaty-ceded waters of Lakes Superior, Huron, and Michigan. Coefficients were averaged across management units within each lake. ....................... 52 CHAPTER 2 Table 1. The rank order performance of control rules for meeting each of several different fishery objectives. Results given in columns of the table correspond to cases assuming no error in estimates of stock size, with the inclusion of error in estimates of stock size, with and without policy parameters adjusted for uncertainty, and with and without autocorrelation in the maximum level of recruitment (see text). When errors in stock size estimates were incorporated, studies that compared performance for control rules using the policy parameters that were optimal without these errors are “unadjusted”; studies that sought policy parameters that were optimal over the uncertainty are “uncertainty adjusted.” When uncertainty in stock assessments was incorporated, rank order reflects finding for the highest levels of assessment error that were evaluated. vii When for a given table column there are no studies that evaluated relative performance of a control rule, these policies are missing (-). ................................................................ 101 Table 2. Papers that compared harvest policies for meeting common fishery objectives assuming no error in estimates of stock Size, with the inclusion of error in estimates of stock size, and with or without autocorrelation in the maximum level of recruitment (i.e., asymptote of a Beverton-Holt stock-recruit function). Specific control rules and characteristics included in each paper are indicated with a X. ..................................... 104 Table 3. Studies that evaluated various biological reference points (BRP) ................. 105 CHAPTER 3 Table 1. The different assessment and implementation error parameter values evaluated. ...................................................................................................................... 155 Table 2. The performance of optimal policy parameters (see text) chosen assuming the incorrect future about lake Whitefish growth (autocorrelated versus recent) for each of four trade-offs in performance metrics and for fast and slow growing stocks in 1836 Treaty-ceded waters. The pairs of values under each grth scenario sub-heading are the performance of the policy parameters chosen assuming the wrong future lake Whitefish growth, and are in the same respective order as the performance metrics defined by the trade-off plot for each row. The values in parentheses are the percent change from the optimal set of policy parameters chosen assuming the correct future lake Whitefish growth. Spawning stock biomass values are reported as a fraction of the unfished level. .............................................................................................................................. 156 APPENDIX A Table A1. Growth category and the range of years for which a random sample of commercial trap-net data were available from lake Whitefish management units in the 1836 Treaty-ceded waters of Michigan used in pararneterizing growth and maturity sub- models. .......................................................................................................................... l 63 Table A2. Expected mean length (mm) at age of lake Whitefish for four growth models (see text for details). ...................................................................................................... 164 Table A3. Upper and lower caps placed on simulated length (mm) at age-3 lake Whitefish for fast and slow growth category Simulations. ............................................ 165 viii LIST OF FIGURES CHAPTER 1 Figure 1. Map of 1836 Treaty—ceded waters and lake Whitefish management units in Lakes Superior (units designated as WF S), Huron (WFH), and Michigan (WF M; Ebener et al. 2005). ..................................................................................................................... 39 Figure 2. Average coefficient estimates (iSE representing uncertainty resulting from variability among management units shown in Figure 1) for the month effect from a general linear mixed model standardizing lake Whitefish catch per effort (catch = aggregate round mass of fish) for the gill-net fishery (top panel; effort = aggregate net length in thousands of feet) and trap-net fishery (bottom panel; effort = number of lifts) in 1836 Treaty-ceded waters of Lakes Superior, Huron, and Michigan. Coefficient estimates were averaged across various years (generally 1981—2001) and management units ...... 40 Figure 3. Annual (generally 1981—2001) proportional difference between an index of lake Whitefish abundance from a general linear mixed mode] (i.e., least-squares mean) and catch per effort (ratio of fish round mass to aggregate net length) from the gill-net fishery in 1836 Treaty-ceded waters of Lakes Superior (management units designated as WFS), Huron (WFH), and Michigan (WFM). ................................................................ 44 Figure 4. Annual (generally 1981—2001) proportional difference between an index of lake Whitefish abundance from a general linear mixed model (i.e., least-squares mean) and catch per effort (ratio of fish round mass to aggregate number of lifts) from the trap- net fishery in 1836 Treaty—ceded waters of Lakes Superior (management units designated as WFS), Huron (WFH), and Michigan (WFM) ............................................................. 48 CHAPTER 2 Figure 1. Basic control rules and how fishing mortality generally changes with biomass or abundance for each type. .......................................................................................... 106 Figure 2. VaIiants of basic control rules and how fishing mortality generally changes with biomass or abundance for each type. .................................................................... 107 CHAPTER 3 Figure l. 1836 Treaty-ceded waters and lake Whitefish management units in Lakes Superior, Huron, and Michigan (Ebener et a1. 2005) .................................................... 141 ix Figure 2. Example of the biomass based control rule used in this analysis (solid line). Dashed lines are provided as a reference for defining the policy parameters. ............. 142 Figure 3. Median spawning stock biomass versus yield (left column) and interannual variability in yield versus the proportion of years with spawning stock biomass less than 20% of the unfished level (right column) for baseline levels of assessment and implementation error parameters and fast growth similar to recent levels (see text for details), for the conditional constant catch control rule (top row), constant-F (black dots, middle row), constant-F with a 15% limit on the interannual change in target catch (grey dots, middle row), biomass based (black dots, bottom row), and biomass based with a 15% limit on the interannual change in target catch (grey dots, bottom row). Each dot corresponds to 1000 simulations for one combination of policy parameters for the given control rule. ................................................................................................................... 144 Figure 4. AS in figure 3, except for median yield versus interannual variability in yield (left column) and yield versus the proportion of years with spawning stock biomass less than 20% of the unfished level (right column). ............................................................ 146 Figure 5. Median Interannual variability in yield versus the proportion of years with spawning stock biomass less than 20% of the unfished level for baseline assessment and implementation error variance, the extent of autocorrelation in assessment error set equal to 0.0 or 0.9, and fast grth similar to recent levels (see text for details). Control rules are displayed as in Figure 3. ......................................................................................... 148 Figure 6. As in Figure 5, except with the extent of autocorrelation in assessment error set equal to the baseline level (0.7), and for assessment error variance equal to 0.01 or 0.20. .......................................................................................................................... 150 Figure 7. As in Figure 5, except with yield versus the interannual variability in yield ............................................................................................................................... 152 Figure 8. As in Figure 5, except with the extent of autocorrelation in assessment error equal to the baseline level (0.7), for assessment error variance equal to 0.01 or 0.20, and for yield versus the interannual variability in yield. ..................................................... 154 APPENDIX E Figure Bl. Median spawning stock biomass versus yield for baseline assessment and implementation error variance, the extent of autocorrelation in assessment error set equal to 0.0 or 0.9, and fast growth similar to recent levels (see text for details), for the conditional constant catch control rule (top row), constant-F (black dots, middle row), constant-F with a 15% limit on the interannual change in target catch (grey dots, middle row), biomass based (black dots, bottom row), and biomass based with a 15% limit on the interannual change in target catch (grey dots, bottom row) .......................................... 168 X Figure B2. As in Figure 1B, except with the extent of autocorrelation in assessment error set equal to the baseline level (0.7), and for assessment error variance equal to 0.01 or 0.20 ................................................................................................................................ 170 Figure B3. As in Figure 1B, except for median yield versus the proportion of years with spawning stock biomass less than 20% of the unfished level. ...................................... 172 Figure B4. Median yield versus the pIOportion of years with spawning stock biomass less than 20% of the unfished level with the extent of autocorrelation in assessment error set equal to the baseline level (0.7), for assessment error variance equal to 0.01 or 0.20, and fast growth similar to recent levels (see text for details). Control rules are displayed as in Figure 1B. ............................................................................................................. 174 Figure B5. Median yield versus spawning stock biomass for baseline levels of assessment and implementation error parameters and for autocorrelated fast growth, autocorrelated slow growth, and slow growth similar to more recent levels. Control rules are displayed as in Figure 1B. ....................................................................................... 176 Figure B6. As in Figure B5, except for variability in yield versus risk. ...................... 178 Figure B7. As in Figure B5, except for yield versus variability in yield. .................... 180 Figure B8. As in Figure BS, except for yield versus the proportion of years with spawning stock biomass less than 20% of the unfished level ....................................... 182 xi INTRODUCTION AND SUMMARY Many fisheries are managed by using estimates of abundance and other parameters from model-based stock assessments (e.g., fitted statistical catch at age models) for setting annual fishery harvest quotas. Stock assessments are often fit to an index of abundance, and so the estimates from the stock assessments can critically rely on the accuracy of both the index and a measure of uncertainty for the index (Maunder and Starr, 2003). Harvest control rules are often used to set a quota as a function of the current estimate of the system state (e. g., an abundance estimate from an assessment). These topics, indices of abundance and harvest control rules, were the main foci of my research. 1. Indices of Abundance Catch per effort (CPE) is usually used as the index of abundance for most fisheries, and the common assumption is that CPE changes in proportion to abundance, which is also referred to as “constant catchability” (Quinn and Deriso, 1999). Violations of this assumption can lead to inaccurate estimates of abundance from stock assessments, and consequently ineffective management, which sometimes results in fishery collapse (Rose and Kulka, 1999; Harley et al., 2001). To avoid violations of this assumption, CPE indices of abundance are ideally based on fishery independent survey data (e.g., Helser et al., 2004). Such surveys are not available for many fisheries and so many indices of abundance used in assessments are based on fishery dependent data. Fishery dependent data is more likely to violate the constant catchability assumption due to things such as systematic changes in characteristics of the fishing fleet (e. g., technological 1 advancements, entrance and exit of individual vessels), non-random search effort, and the spatial distribution of the fish stock (Rose and Kulka, 1999; Harley et al., 2001; Maynou et al., 2003; Battaile and Quinn, 2004; Bishop et al., 2004; Campbell, 2004). Even stock assessment models that allow for some temporal changes in catchability will tend to work better when such temporal variation is lower (Wilberg and Bence, 2006; Wilberg et al., 2008) To account for some of the variation in CPE not attributable to changes in abundance, and provide a more accurate index, CPE data can be “standardized” by fitting statistical models to the catch and effort data, and then using “year-effect” estimates as the index of abundance (Maunder and Punt, 2004; Venables and Dichmont, 2004). Year- effect estimates are commonly used because detecting trends in abundance over time is usually the objective (Maunder and Punt, 2004). Frequently, some form of general or generalized linear model is used to standardize the CPE data (Maunder and Punt 2004). 1.1. Chapter I : Improving indices of abundance for lake Whitefish My main objective in Chapter 1 was to produce standardized indices of abundance for lake Whitefish in 1836 Treaty-ceded waters of Lakes Huron, Michigan, and Superior, but this work also allowed me to develop expertise in statistical techniques (e.g., mixed models) that I used to parameterize the simulation model of chapter 3. Currently, statistical catch at age assessments are fit in each of 18 management units, and a quota is also set for each unit. The assessments are fit using two separate CPE indices of abundance from gill-nets and trap-nets, with CPE estimated as the ratio of sum of aggregate catch to sum of aggregate effort in each year. I developed general linear mixed models (GLMM) for each gear type to standardize the fishery CPE data. Factors 2 included in the GLMMS were fixed effects of year, month, and boatsize (gill-net fishery only), and random effects of license holder (i.e., analogous to boat captain), grid (i.e., location), and all two and three way interactions. The effect of the standardization by using the GLMM method was evaluated by examining the temporal trends in the proportional difference (PD) between the least squares means for each year (LSM) and CPE (i.e., aggregate catch divided by aggregate effort for each year). Since both the LSMS and CPE are relative indices, changes in PD over time were of interest and not whether average PD differed from 1.0. Factors that were particularly influential in the GLMM models were month, boat size, and license holder, which was similar to factors important for marine commercial fisheries where standardization is more widely applied than in freshwater systems. The proportional difference between the LSMS and CPE trended through time in some management units, suggesting that adjusting fishery CPE for effects such as boat size, season, and license holder was important. So, I concluded that model-based indices of abundance should replace non-standardized CPE in some lake Whitefish stock assessment models, especially those management units where the proportional difference trended through time. In management units where the proportional difference did not trend through time, using a model-based index of abundance may still be beneficial. Accounting for variability due to random effects led to year specific estimates of uncertainty (e. g., the standard errors for the LSMS) that were not available when using non-standardized CPE. Using improved years-specific estimates of uncertainty to weight the influence of indices of abundance can increase the accuracy of stock assessment estimates (Helser et al., 2004; Maunder and Starr, 2003). 2. Harvest Control Rules Harvest control rules are guidelines that specify an amount of catch, fishing effort, or fishing mortality as a specific, and usually simple, function of a current estimate of the system state (e.g, spawning biomass; Deroba and Bence, 2008). Common control rules include constant catch, constant fishing mortality rate, constant escapement, or a few variations of these. Each control rule is also defined by a number of policy parameters. For example, the constant fishing mortality rate control rule is defined by one policy parameter, the target level of fishing mortality. Ideally, a harvest control rule is chosen because it meets fishery objectives (e. g., maximize yield, minimize interannual variability in yield). However, which rules are best at meeting certain fishery objectives is a source of controversy in the literature. Furthermore, the relative performance of control rules depends on Specific characteristics of the fishery and underlying population dynamics that are incorporated into an evaluation. Consequently, selecting a harvest control rule and policy parameters can be a difficult task. 2.]. Chapter 2: A review of harvest control rules In Chapter 2 I reviewed the harvest control rule literature with two Objectives: 1) to compare and contrast the relative performance of various control rules at meeting common fishery objectives, and 2) to identify reasons for what seem to be contradictory results. The findings were also relevant for designing the harvest control rule evaluation of Chapter 3 (see below). I found that the relative performance of control rules at meeting common fishery objectives was affected by: the given fishery objective, whether uncertainty in estimated stock sizes was included in analyses (i.e., assessment error), whether the maximum recruitment level (e.g., the asymptote of a Beverton—Holt stock— 4 recruit fiInction) varied in an autocorrelated fashion over time, and the amount of compensation in the stock—recruit relationship. Also, few studies have compared control rules using optimal parameters (e.g., those that maximize some objective function) that were found while including assessment error. More commonly, parameters that are optimal without assessment error are used in a comparison of control rules that includes assessment error. This approach can produce misleading results. Lastly, more research is needed to compare control rules when accounting for uncertainty in key population parameters, when stock—recruitment or other population dynamic parameters vary over time, and for fisheries with non-yield-based or competing objectives. 2.2. Chapter 3: Evaluating harvest control rules for lake Whitefish Chapter 3 addressed some of the harvest control rule research needs identified in Chapter 2, and was based on a simulation analysis with the objective of evaluating the ability of alternative control rules to meet fishery objectives for lake Whitefish in 1836 Treaty-ceded waters. Currently, a quota is set for each management unit so that total annual mortality rate equals 65% for ages experiencing the highest levels of fishing mortality. Because assessments in these waters assume a constant natural mortality rate across ages and time (Ebener et al., 2005), this is equivalent to a constant fishing mortality rate (constant-F) control rule. The constant-F control rule and the parameter for the control rule (i.e., 65% total annual mortality rate) are based on analyses conducted over 30 years ago (Healey, 1975), and so may not be optimal for meeting fishery objectives. Lake Whitefish stocks in 1836 Treaty-ceded waters are characterized by temporal and Spatial variation in various population parameters. For example, lake Whitefish 5 growth in some areas of the Great Lakes declined during the 19905 and 20003, coincident with declines in an important prey source, Diporeia (Hoyle et al., 1999; Pothoven et al., 2001; Mohr and Nalepa, 2005), but similar declines have not occurred everywhere despite Similar ecosystem changes (e.g., Cook et al., 2005; Lumb et al., 2007). Growth rates, maturity ogives, natural mortality, and stock-recruit relationships also likely differ spatially among some of the management units (e.g., Wang et al., 2008). Drawing from my experiences with GLMMS from Chapter 1 and partially based on the results of Chapter 2, I developed a stochastic age-structured simulation model that incorporated stochasticity in life history traits, uncertainty in future lake Whitefish growth, and other sources of uncertainty to compare the current harvest control rule with a range of alternative control rules, including conditional constant catch (CCC), constant- F, biomass-based (BB), and constant-F and BB rules with a 15% limit on the interannual change in the target catch. Separate sets of growth parameters were estimated for fast and slow growth stocks, and separate sets of simulations were done for these two categories of individual stocks. Furthermore, I developed two variants of a growth model to represent alternative hypotheses about future lake Whitefish growth; one with temporally autocorrelated changes in growth and another where growth remained similar to more recent patterns. Uncertainty in the stock-recruitment relationship was incorporated by drawing stock-recruit parameters for each simulation from a set of possible values, which were based on data from each management unit and estimated using a GLMM (i.e., similar statistical model used in Chapter 1). The simulations also included assessment and implementation error. Some of the model features mentioned above were included because the results of Chapter 2 indicated that these can affect 6 relative control rule performance, in particular, accounting for uncertainty in the stock- recruit relationship and assessment error. Each control rule was evaluated over a range of the policy parameters that define the control rules. The performance of the control rules was evaluated by examining trade-off plots of spawning stock biomass (SSB) versus yield (Y), interannual variability in yield (Yvar) versus the proportion of years that SSB fell below 20% of the unfished level (SSBF=0), Y versus Yvar, and Y versus the proportion of years that SSB fell below 20% of SSBF=0. While treating future growth as known, the rank order performance of the control rules for each of the performance metrics was generally robust to sources of uncertainty. For example, the constant-F and BB rules simultaneously attained higher average yield and spawning stock biomass than all other control rules. The CCC rule and limiting the constant-F or BB mles to a 15% change in target catch had the lowest yearly variability in yield. The low yearly variability in yield provided by limiting target catch changes to 15%, however, came at the cost of frequently reducing biomass to low levels, so that in many situations other control rules would be preferred. The sensitivity of results to uncertainty about future lake Whitefish growth was control rule specific and depended on whether stock growth was fast or slow. For fast growth stocks, selecting control rules and policy parameters by incorrectly assuming that future growth will be autocorrelated resulted in little cost from the optimum levels relative to the alternative of incorrectly assuming future growth will be similar to recent levels. For slow growth stocks, however, the robustness to choosing policy parameters based on an erroneous assumption about future lake Whitefish growth depended on the control rule and trade-off plot. The decision about how best to select control rules and 7 policy parameters will ultimately depend on how competing fishery objectives are weighted relative to each other. Generally, however, control rules and policy parameters for fast growth stocks Should likely be selected assuming future growth will be autocorrelated, but a universal recommendation for slow growth stocks is less clear (i.e., depends on the control rule and fishery objectives). Depending on how important different fishery objectives are, a control rule and policy parameters other than the one currently in use (i.e., constant-F based on a total annual mortality rate of 65%) may be worth considering. For example, a BB control rule with appropriately selected policy parameters could likely produce nearly the same or more yield, spawning stock biomass, and less risk with little cost in variability in yield relative to the currently used policy. Similarly, the CCC control rule can likely provide less variability in yield, but at the cost of yield. So, if maintaining low variability in yield is more desirable than maximizing yield, a CCC control rule may want to be considered. 3. Overall Conclusions and Future Directions The results of this dissertation have implications for the improved management of lake Whitefish in the Great Lakes, but the results are also more generally applicable. In Chapter 1, I found that model-based indices of abundance should likely replace non- standardized indices in fitting stock assessment models. The factors important to the standardization process also seem to be consistent among systems, and so should be considered when standardizing CPE data for most fisheries. Likewise, updating stock assessments for most fisheries to include standardized indices of abundance and associated measures of uncertainty would likely produce more accurate estimates of abundance and other population parameters, and so reduce assessment error, which in 8 Chapter 2 was shown to affect relative control performance. In addition to assessment error, Chapter 2 highlighted several other characteristics and uncertainties of harvest policy evaluations that have affected control rule performance, and so should be considered when developing harvest policy analyses for any fishery. The results of Chapter 2, however, also revealed that little research has historically considered these characteristics. Chapter 3 added to the body of research that has considered factors important to control rule performance. The CCC control rule, which was first published in an analysis of Pacific halibut Hippoglossus stenolepis, had never been evaluated while considering assessment error (Clark and Hare, 2004). Similarly, few published analyses have considered control rules with limits on the interannual change in target catch. Lake trout Salvelinus namaycush in 1836 Treaty-ceded waters are managed with such a restraint, but given the generally poor performance of these control rules another option may be warranted. Chapter 3 also evaluated the sensitivity of relative control rule performance to one form of time-varying growth that had never been considered before, and time-varying population parameters have been Shown to affect control rule performance (Chapter 2). The results in regards to the rank order and sensitivity of the control rules to this source of uncertainty are likely generally applicable to any fishery experiencing similar conditions. References Battaile, BC. and T.J. Quinn 11. 2004. Catch per unit effort standardization of the eastern Bering Sea walleye Pollock fleet. Fisheries Research 70 (2004): 161-177. Bishop, 1., W.N. Venables, and Y-G. Wang. 2004. Analyzing commercial catch and effort data from a Penaeid trawl fishery: A comparison of linear models, mixed models, and generalized estimating equations approaches. Fisheries Research 70(2004): 179-193. Campbell, RA. 2004. CPUE standardisation and the construction of indices of stock abundance in a spatially varying fishery using general linear models. Fisheries Research 70(2004): 209-227. Clark, W.G., and SR. Hare. 2004. A conditional constant catch policy for managing the Pacific halibut fishery. North American Journal of Fisheries Management 24: 106-113. Cook, H.A., T.B. Johnson, B. Locke, and BI. Morrison. 2005. Status of lake Whitefish in Lake Erie. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 87-104. Ebener, M.P., J.R. Bence, K. Newman, and P. Schneeberger. 2005. Application of statistical catch-at-age models to assess lake Whitefish stocks in the 1836 treaty- ceded waters of the upper Great Lakes. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 271-309. Harley, S.J., R.A. Myers, and A. Dunn. 2001. Is catch-per-unit-effort proportional to abundance? Canadian Journal of Fisheries and Aquatic Sciences 58: 1760-1772. Healey, MC. 1975. Dynamics of exploited Whitefish populations and their management with special reference to the Northwest Territories. Journal of the Fisheries Research Board of Canada. 32: 427-448. Helser, T.E., A.E. Punt, and RD. Methot. 2004. A generalized linear mixed model analysis of a multi-vessel fishery resource survey. Fisheries Research 70(2004): 251-264. Hoyle, J .A., T. Schaner, J .M. Casselman, and R. Derrnott. 1999. Changes in lake Whitefish stocks in eastern Lake Ontario following Dreissena mussel invasion. Great Lakes Research Review 4: 5-10. 10 Lumb, C.E., T.B. Johnson, H.A. Cook, and J .A. Hoyle. 2007. Comparison of lake Whitefish growth, condition, and energy density between Lakes Erie and Ontario. Journal of Great Lakes Research 33: 314-325. Maunder, MN. and AB. Punt. 2004. Standardizing catch and effort data: a review of recent approaches. Fisheries Research 70(2004): 141-159. Maunder, MN. and P.J. Starr. 2003. Fitting fisheries models to standardized CPUE abundance indices. Fisheries Research 63(2003): 43-50. Maynou, F., M. Demestre, and P. Sanchez. 2003. Analysis of catch per unit effort by multivariate analysis and generalised linear models for deep-water crustacean fisheries off Barcelona. Fisheries Research 65(2003): 257-269. Mohr, LG, and Nalepa, T.F. (Editors). 2005. Proceedings of a workshop on the dynamics of lake Whitefish (Coregonus clupeaformis) and the amphipod Diporeia spp. in the Great Lakes. Great Lakes Fishery Commission Technical Repport 66. Pothoven, S.A., T.F. Nalepa, P.J. Schneeberger, and SB. Brandt. 2001. Changes in diet and body condition of lake Whitefish in southern Lake Michigan associated with changes in benthos. North American Journal of Fisheries Management: 21: 876- 883. Quinn, T.J., II, and RB. Deriso. 1999. Quantitative Fish Dynamics. Oxford University Press Inc. New York, New York. Rose, GA. and D.W. Kulka. 1999. Hyperaggregation of fish and fisheries: how catch- per-unit-effort increased as the northern cod declined. Canadian Journal of Fisheries and Aquatic Sciences 56(supplement 1): 118-127. Venables, W.N., and CM. Dichmont. 2004. GLMS, GAMS, and GLMMS: an overview of theory for applications in fisheries research. Fisheries Research 70(2004): 319- 337. Wang, H-Y, T.O. HOOk, M.P. Ebener, L.C. Mohr, and P.J. Schneeberger. 2008. Spatial and temporal variation of maturation schedules of lake Whitefish in the Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 65: 2157-2169. Wilberg, M.J., and J .R. Bence. 2006. Performance of time-varying catchability estimators in statistical catch-at-age analysis. Canadian Journal of Fisheries and Aquatic Sciences 63: 2275-2285. Wilberg, M.J., B.J. Irwin, M.L. Jones, and LR. Bence. 2008. Effects of source-Sink dynamics on harvest policy performance for Yellow Perch in southern Lake Michigan. Fisheries Research 94: 282-289. 11 CHAPTER 1 Deroba, J .J . and J .R. Bence. 2009. Developing model-based indices of lake Whitefish abundance using commercial fishery catch and effort data in Lakes Huron, Michigan, and Superior. North American Journal of Fisheries Management 29: 50-63. The content of this chapter is intended to be identical to the cited publication and is based on the accepted manuscript with changes that reflect corrections made during copy editing. Any differences Should be minor and are unintended. 12 Abstract Fishery catch per effort (CPE) is often used to assess relative fish abundance, and in many Great Lakes and other freshwater applications this is based on either an average or the ratio of sum of aggregate catch to sum of aggregate effort. In particular, assessments used to estimate the abundance of lake Whitefish and recommend harvest quotas in the 1836 Treaty-Ceded waters of Lakes Huron, Michigan, and Superior assume that commercial CPE from gill-net and trap-net fisheries is proportional to abundance, but CPE may change due to factors other than abundance, leading to violations of this assumption. To account for sources of variation in CPE not attributable to abundance, general linear mixed models (GLMMS) were developed for each management unit, and least squares means (LSMS) for each year were used as the index of abundance. The effect of the standardization by using the GLMM method was evaluated by examining the temporal trends in the proportional difference between the LSMS and CPE (i.e., aggregate catch divided by aggregate effort for each year). Of the random effects included in the final GLMM for the gill-net fishery, license holder accounted for the most variation. The fixed effect of boat size category on CPE depended on lake, where on average in Lake Superior there was little difference, but in Lakes Michigan and Huron large boats had lower CPE than medium and small boats. CPE was on average higher from October to December than in other months. The proportional difference between the LSMS and CPE trended through time in some management units, suggesting that adjusting fishery CPE for effects such as boat size, season, and license holder is important. Factors influential to lake Whitefish commercial fishery CPE are similar to factors that have been shown to be important in marine commercial fisheries. 13 Introduction Lake Whitefish, C oregonus clupeaformis, has supported a historically important fishery for Native American bands and a highly valued commercial fishery in the upper Great Lakes (Lakes Huron, Michigan, and Superior). In the late 18003 and early 19003, lake Whitefish were often the most highly valued commercial species and usually comprised the greatest proportion of total yield from each of the upper Great Lakes (Koelz 1926; Brown et a1. 1999). Lake Whitefish stocks collapsed in each of these lakes in the 19303 and 403 due to overexploitation, sea lamprey, Petromyzon marinus, predation, and pollution (Smiley 1882; Koelz 1926; Jensen 1976; Brown et al. 1999; Ebener and Reid 2005). From the 19603 through the 19803, lake Whitefish stocks rebounded in each of the lakes largely due to improved management of commercial harvest, sea lamprey control, pollution remediation, and the introduction of salmonines that reduced the abundance of the invasive alewife, Alosa pseudoharengus, and rainbow smelt, Osmerus mordax (Ebener 1997; Mohr and Ebener 2005a). In the 19903, lake Whitefish once again became the main commercial species, particularly in Lake Huron where the species comprised over 80% of the total commercial yield (Mohr and Ebener 2005b) In 1979, the rights of Native American bands to fish in the Michigan waters of the upper Great Lakes, as reserved in a treaty signed in 183 6, were reaffirmed by US. federal courts. Since the reaffirmation of treaty fishing rights, periodic stock assessments have been conducted for stocks within Spatially defined management units, with the fishery data and harvest from within each management unit treated as applying to a reproductively isolated stock (Figure 1; Ebener et al. 2005). Stock assessments are 14 conducted and harvest recommendations based on the assessments are made annually for each individual management unit. Within each management unit commercial fishery catch and effort data are reported on a 10-minute by 10-minute statistical grid basis, which allows for some spatial resolution within management units. Since 2000, guidelines for the management of lake Whitefish have been set according to a Consent Decree. The 2000 Consent Decree created a Technical Fisheries Committee (TFC) and its Modeling Subcommittee (MSC) to conduct stock assessments and Specify total allowable catches (TACs) and harvest regulating guidelines (HRGS, see below). TACS are limits to catch, and are used in management units where some yield is allocated to the state licensed fishery and some to the tribal fishery. HRGS are targets for yield used to guide regulations for lake Whitefish in units where all yield is allocated to the tribal fishery. The MSC fits statistical catch-at-age (CAA) models to commercial fishery data to estimate population numbers, mortality rates, fishery harvest, and other population parameters of interest. The estimates of the population parameters are then used to project each stock’s abundance into the future, and then a TAC or HRG is calculated by applying a reference mortality rate to the estimate of the next year’s abundance. The CAA models use fishery effort data and an assumed relationship between fishing mortality and fishery effort. Age (a) and year (y) specific fishing mortality rates (F) are estimated as the product of age specific selectivity (S) and year Specific “fishing intensity” (f) for each of two fishery gears, gill-nets and trap-nets: Fi,a,y = Si,afi,y; (1) where i denotes gear type and, 15 fi, y = Er, yqi, yer, y ; (2) where E is fishery effort specific to each gear type, q is catchability, and e is multiplicative observation error. The details of the CAA models have been described in Ebener et a1. (2005). Equation 2 is equivalent to assuming that the commercial fishery catch per effort (CPE), estimated as the ratio of sum of aggregate catch to sum of aggregate effort in each year, is on average proportional to average abundance over the fishing year, and that deviations from this average relationship are independent variations from year to year. Violations of the assumption that CPE is proportional to average abundance can occur due to changes in fishing power of gear, or if the Spatial and temporal distribution of fishery effort is non-random (Quinn and Deriso 1999). Violations of this assumption are called hyperdepletion when CPE declines faster than abundance at high stock sizes, and hyperstability when CPE does not decline as drastically as abundance at high stock sizes (Quinn and Deriso 1999). For example, an increase in the number of fishing operations could cause some fishermen to operate in lower quality habitat. Thus, CPE could decline even if fish abundance did not, resulting in hyperdepletion. Hyperstability is the more common occurrence and leads to overestimation of biomass and underestimation of fishing mortality, which has too ofien gone unrecognized and led to fishery collapses (Rose and Kulka 1999; Harley et a1. 2001). To account for some of the variation in CPE not attributable to changes in abundance, and improve assessments and associated fishery management, CPE can be “Standardized” by fitting statistical models to the catch and effort data, and then using “year-effect” estimates as the index of abundance (Maunder and Punt 2004; Venables and 16 Dichmont 2004). Commonly, some form of general or generalized linear model is used to standardize the CPE data (Maunder and Punt 2004). Year is usually included as one of the explanatory variables because detecting trends in abundance over time is usually the objective (Maunder and Punt 2004). Other explanatory variables often include a spatial element or some measure of individual vessel fishing power (e.g., boat size) (Battaile and Quinn 2004; Bishop et al. 2004). Our objectives were (1) to standardize lake Whitefish CPE data in the upper Great Lakes to attain an index of abundance that more accurately reflected changes in lake whitefish biomass than CPE; (2) gain an improved understanding of factors that influence commercial fishery CPE for lake Whitefish; and (3) compare the factors that are important for this fishery with those found to influence CPE in other fisheries of the world. Currently for lake trout, Salvelinus namaycush, in these waters, indices of abundance are based on the least squares means (LSMS) for each year from a general linear mixed model (GLMM; Deroba and Bence in press). Consequently, we explored the use of a similar GLMM for lake Whitefish, and compared the temporal trends in the LSMS for each year to that of the CPE. Our concern here is that the LSMS account for sources of variation in CPE not considered when CPE is estimated as a ratio of sum of aggregate catch to sum aggregate effort in each year, and might reveal substantially different interannual trends in apparent relative abundance. Methods Study Area Our study area was the waters relevant to the 1836 Treaty, which encompassed the majority of Michigan waters of Lakes Superior, Huron, and Michigan (Figure l). 17 These waters were stratified into 18 management units with individual surface areas ranging from 69,000 to 733,000 ha, and a total surface area of 5.8 million ha (Figure l; Ebener et a1. 2005). Analyses were done separately for each management unit because these are treated as reproductively isolated stocks and define the resolution of spatial stratification used to manage lake Whitefish (see introduction; Ebener et a1. 2005). Data and Analyses Data were collected from commercial fishing operations as part of a requirement for all licensed vessels to submit monthly reports that describe for each day of the month the weight of fish landed, the amount of gear lifted, the 10-minute by 10-minute statistical grid where the catch and effort occurred, and other auxiliary information (Ebener et a1. 2005). Monofilament large-mesh gill-nets with 2 114-mm stretched mesh and 6-14 m tall trap-nets accounted for nearly 100% of the lake Whitefish commercial harvest, and analyses were only conducted on these two gear types. The range of years included in this study differed by management unit and gear type, and some years are missing because no catch or effort was reported (Table 1). Analyses were only conducted on 12 of the 18 management units for the gill-net fishery, and 10 of the 18 management units for the trap-net fishery because few or no observations were recorded within most years for some management units and gears. CPE was estimated separately for gill-nets and trap-nets as the ratio of sum of aggregate catch to sum of aggregate effort in each year, as is currently used in the CAA models. Catch was measured as the round mass of Whitefish for both gears, while effort was measured in 10003 of feet of net for gill-nets, and number of lifts for trap-nets. 18 GLMMS were fit separately for gill-nets and trap-nets, with log.= (CPE+1) as the dependent variable. We applied a log, transformation because examination of the distribution of the data showed that this was necessary to meet the assumption of normality for general linear models (McCulloch and Searle 2001; Gelman and Hill 2007). We added 1.0 to all CPE observations prior to transformation to address the (infrequent, ~0.001% for both gear types) occurrence of zero CPE observations. This added constant represents a low CPE for gill nets and the lowest possible CPE for trap nets, and more than 99% of CPE values exceeded 1.0 (the constant) for both gear types. Our initial full model for gill-nets included fixed effects of year, month, and boat size, and random effects of license holder, grid, and all possible two and three way interactions. In preliminary analyses, interactions of a higher order than three ways were not estimable for any management units, and so were excluded from firrther consideration. Because not enough individual license holders fished with multiple boat sizes, license holder and boat size were confounded when two and three way interactions with license holder and two and three way interactions with boat size were included in the same model. Furthermore, in preliminary analyses interactions with license holder were only estimable for two management units, while interactions with boat size were estimable in all management units. Consequently, all interactions with license holder were also excluded from further consideration. Thus, the new “full” model included fixed effects of year ((1),), month (13,"), boat size (71,), and random effects of license holder (6;), grid (kg), and all two and three way interactions except those with license holder: loge(CPE+l)=,u+ay +,Bm +yb+c1+kg +oym +pyb +qyg +rmb +smg +tbg + umbg + dgmy + hgyb + fymb + 5iymbgl§ l9 where u is the overall mean, 0y", is the interaction of year and month, pyb is the interaction of year and boat size, qyg is the interaction of year and grid, rmb is the interaction of month and boat size, Smg is the interaction of month and grid, t bg is the interaction Of boatsize and grid, umbg is the interaction of month and boat size and grid, dgmy is the interaction of grid and month and year, hgyb is the interaction of grid and year and boat size, jymb is the interaction of year and month and boat size, and 8iymbgl is residual error for each observation, 1'. This model assumes that the random effects and residual error are all independent and identically distributed as normal with a mean of zero. Boat size was a categorical effect and sizes were defined as: small (3 20 ft), medium (20-30 It), and large (2 30 ft). The full model for trap—nets included fixed effects of year and month, and random effects of license holder, grid, and all two and three way interactions: loge(CPE+l)= p+ay +,6m +c1+kg +0ym +vyl+wm1+smg +ng +qyg +Zyml +dgmy +aygl +emgl +5iymgl§ where vyl is the interaction of year and license holder, Wm] is the interaction of month and license holder, xg’ is the interaction of grid and license holder, Zyml is the interaction of year and month and license holder, ayg, is the interaction of year and grid and license holder, 8mg] is the interaction of month and grid and license holder, and all other terms are defined as for gill-nets. In four of the 10 management units analyzed for the trap-net fishery, all of the observations came from one boat size category, and so this effect was not evaluated. 20 Final models for both gear types were determined by evaluating which effects could be removed using corrected Akaike’s information criterion (AICc) (Bumham and Anderson 2002). Our model selection approach was to first consider which random effects would be removed from the final model while keeping all fixed effects in the model (N go and Brand 1997). Random effects were selected prior to fixed effects so that the final models had the Simplest error structure possible (i.e., a random effect would be eliminated rather than a fixed effect that explained similar sources of variation). Our approach to selecting random effects was to drop each random effect one at a time, while keeping all other effects in the model. Once a random effect was removed, AAlCc was then calculated by subtracting AICc for the reduced model from AICc for the full model. If AAICc was greater than 2.0 (Bumarn and Anderson 2002), the factor not present in the reduced model was eliminated from the final model, otherwise the factor was retained. We followed this approach because with 22 management unit and gear combinations and 12 potential random effects to consider for each, fitting and comparing all possible models was not practical. A random effect was also dropped from the final model if the variance estimate for that factor was zero. Restricted maximum likelihood (REML) was used for model fitting when comparing models with different random effects, given its superior performance in estimating random effects (McCulloch and Searle 2001). Once the best set of random effects was selected, the best set of fixed effects was selected by comparing AICc values for all possible combinations of fixed effects. Models were fit using maximum likelihood (ML) instead of REML because comparisons with AICc based on REML are not valid when comparing models with different fixed effects (SAS 2003). During this process the previously determined best random effects 21 portion of the model was used. Year ((1,) was not evaluated during model selection because the objective is to estimate a yearly index of abundance, and 30 year must be retained in the final model. The AAICc values are not reported in the results because this would require reporting a value for each factor that was included in the firll models for each management unit and gear type (i.e., 298 values). Rather, we report the AAICc values between a means model (i.e., a model with only a year effect) and the final model (AAICc = AICc means model — AICc final model) to quantify the likely improvement that the final models offer over the current indices of abundance that do not account for factors other than year. Generally, the same effects were included in the final model for each management unit, but the models for some management units could be improved by the elimination of an effect that improved model fit for the majority of the management units, or inclusion of an effect that did not improve model fit for the majority of the management units. For the Simplicity of reporting results in these analyses, we eliminated an effect in all management units if it only improved model fit in a minority of management units. LSMS for each year were calculated by summing the overall mean (u), the coefficient estimate for each year ((1,), and the average of the coefficient estimates over all levels of fixed effects other than year in the final models (SAS 2003). The LSMS for each year from the final model, as determined by the majority, were nearly identical to the LSMS from other models that improved model fit for a minority of management units. Consequently, we believe that the conclusions of these analyses are robust to this approach. However, if the estimated uncertainty (e.g., standard errors) associated with 22 LSMS (or alternatively year effects or other functions of model parameters) is important, as in fitting stock assessment models to indices of abundance where the standard errors are used to weight the indices of abundance relative to other data (e. g., Maunder 2001; Maunder and Starr 2003), a different model than that reported as the final model here may be warranted for some management units. Differences in the back-transformed LSMS for each year and CPE+1 were qualitatively examined by plotting the proportional difference (PD) between the two measures across years for each management unit included in this analysis. PD was calculated as: D _ (CPE +1) exp(LSM) ' The PD is a measure of how much larger or smaller CPE is than the LSMS. For example, if PD=2 then the CPE is two times larger than the index of abundance based on the mixed mode]. Since both the LSMS and CPE are relative indices, changes in PD over time are of interest and not whether average PD differs from 1.0. Consequently, if PD varied without trend we concluded that the two approaches generally suggested similar trends in abundance through time, although differences may have existed for a given year. Conversely, if PD trended through time we concluded that the index of abundance provided by the two approaches suggested different temporal trends. The relative effect of factors included in the final model on CPE was determined by averaging coefficient estimates across management units and comparing the average values. For random effects, the variance component estimates for each effect were used in estimating the average; while for fixed effects, the coefficient estimates for each level 23 of a factor were used. For boat Size, the averages were estimated separately for each lake because different boat sizes may perform differently in each lake. Results Gill-net Fishery The final model for the gill-net fishery included fixed effects of year, month, and boat Size, and random effects of license holder and the interaction of year and month: loge(CPE+1) = ,u + ay + ,8," +7b +c1+oym +5iymbgl- The final model improved model fit over a means model in all but one management unit, with an average AAICc value of 362.5 and values ranging from -10 to 2514 (Table 2). The final model may not have improved fit over a means model in WFM-06 because this management unit had the smallest sample size (N=308; mean N=1452), which may not provide enough data to adequately capture the variability in CPE caused by the various factors. Of the random effects, the license holder effect accounted for the most variation in CPE (Table 3). The effect of boat size depended on lake (Table 4). In Lake Superior, CPE did not vary much among boat size classes. On Lake Huron, small and medium boats had similar CPE, which was less than that for large boats. On Lake Michigan, CPE ordered as medium > small > large boats. CPE was generally low during January through September, highest in October and November, and intermediate between these levels in December (Figure 2). The index of abundance provided by the GLMMS suggested different temporal patterns than CPE (i.e., PD trended through time) over some or all of the time series in some management units for the gill-net fishery (Figure 3). In Lake Huron, the PD for management units WFH-Ol and WFH-04 generally varied without trend, while in WFH- 24 02 PD declined during 1982-1983, but varied without trend for the remainder of the time series. In Lake Michigan, PD in WFM-02 increased during 1987-1988 and then decreased. In WFM-03, PD increased in variability over the time series and increased during 1999-2001. PD in WFM-04 generally declined through time. In WFM-OS, PD generally varied without trend, but declined during 1997-1999 and then increased. In WFM-06, PD declined during 1993-1997. In Lake Superior, the PD in WFS-05, WFS- 06, WFS-07, and WFS-08 generally varied without trend, except during 1999-2001 in WFS-05 when PD declined. Trap-net Fishery . The final model for the trap-net fishery included fixed effects of year and month, and random effects of the interactions of month and year, year and license, and month and year and license: loge(CPE +1) = ,u + or), + ,6," + kmy + vy1+ pmy1+ aiml. The final model improved model fit over a means model in all management units by an average AAICc value of 170.1, with values ranging from 2.2 to 478.2 (Table 2). Of the random effects, the interaction of year and license holder accounted for the most variation in logc(CPE+l), even more than residual error (Table 3). CPE was generally low during January through September, with the exception of May, highest in October and November, and intermediate between these levels in December (Figure 2). The index of abundance provided by the GLMMS showed different temporal trends than CPE (i.e., PD trended through time) over all or some of the time series in some management units for the trap-net fishery (Figure 4). In Lake Huron, the PD in WFH-Ol and WFH-02 generally varied without trend, while the PD in WFH-04 varied 25 without trend until 1998 when PD increased to 2000 and then decreased. In Lake Michigan, the PD in WFM-Ol, WFM-02, and WFM-03 generally varied without trend, except during 2000-2001 in WFM-Ol when PD increased. In WFM-04 and WFM-OS, PD varied cyclically with a period of approximately two years in WFM-04 and six years in WFM-OS. In Lake Superior, the PD in WFS-07 generally varied without trend, while the PD in WFS-08 increased during 1984-1986, but varied without trend during the few other years of data. Discussion CPE is often assumed to be proportional to abundance, but CPE can change due to factors other than abundance that cause violations of this assumption (Quinn and Deriso 1,999; Battaile and Quinn 2004). Violations of the assumption of proportionality can lead to inaccurate estimates of abundance from stock assessments, and in particular hyperstability can increase the risk for fishery collapse (Rose and Kulka 1999; Harley et a1. 2001). Indices of abundance based on commercial fishery catch and effort data are at an especially high risk of violating the assumption of proportionality due to things such as systematic changes in characteristics of the fishing fleet (e.g., technological advancements, entrance and exit of individual vessels), non-random search effort, and the spatial distribution of the fish stock (Rose and Kulka 1999; Harley et a1. 2001; Maynou et a1. 2003; Battaile and Quinn 2004; Bishop et al. 2004; Campbell 2004). For these reasons, fishery CPE data from many major marine fisheries are now often standardized using various statistical models (e.g., general linear mixed models, generalized linear models) that account for some of the variation in CPE not attributable to abundance, 30 that the “year-effect” becomes a more accurate index of abundance (Maunder and Punt 26 2004; Venables and Dichmont 2004). Factors commonly included in models used to standardize CPE data include factors for time (usually year), location (e.g., grid in this study), individual vessels, characteristics of vessels that affect catchability (e. g., vessel size, horsepower, GPS), among other factors (Maunder and Punt 2004). The temporal trends exhibited by standardized CPE data (e.g., LSMS) have differed from that of non-standardized CPE data (e.g., ratio of aggregate catch to aggregate effort in each year) in other studies (Maynou et a1. 2003; Battaile and Quinn 2004), as was true for some management units in our evaluation of Great Lakes Whitefish fisheries. Thus, we believe that model-based indices of abundance should replace non- standardized CPE in some lake Whitefish stock assessment models, especially those management units where PD was shown to trend through time. Converting to the use of model-based indices of abundance in the stock assessment models for these management units would likely produce more accurate estimates (e.g., abundance estimates) than the current approach of treating raw effort as an index of fishing mortality (equivalent to using CPE as an abundance index). This outcome would also likely hold true for other freshwater systems, where model based methods for standardizing CPE data have not been used as frequently as in marine systems. The reason for the changes in PD in this study can be partially explained by when most fishing occurred and who fished in each year. For example, in 1988 in the WFM-02 gill-net fishery, fewer observations were made in the spring (i.e., when CPE is lower relative to other times of year) and more observations were taken from license holders with relatively high CPE than in other years, which may explain the spike in PD. Similarly, in the WFM-04 trap-net fishery, peaks in PD occurred in years when more 27 observations came from license holders who did well in that year relative to other license holders. Consequently, indices of abundance based on CPE in these and other areas would most likely be driven by differences in the number of observations taken among seasons or from difference license holders, and not due to changes in abundance as is being assumed in stock assessments. In addition to providing a more accurate index of abundance, the use of mixed effects models also allows the uncertainty around the indices of abundance to be more accurately quantified for each year, and this can be especially important if these estimates of uncertainty are used to weight the importance of the yearly CPE indices in stock assessment models (Helser et a1. 2004, Maunder and Starr 2003). Maunder and Starr (2003) describe methods for how yearly indices of abundance can be weighted by their coefficient of variation in fitting stock assessment models, and also found that stock assessment estimates (e.g., abundance estimates) can be less accurate when each yearly index of abundance is weighted equally, instead of using a year specific weight. Furthermore, Helser et a1. (2004) found that ignoring the variability due to random effects, including vessel and the interaction of vessel and year, similar to the effects of license and the interaction of license and year in this study, may lead to an underestimation of uncertainty in indices of abundance. Thus, if the CPE data used in fitting lake Whitefish stock assessment models were replaced with model-based standardized CPE indices and an associated estimate of uncertainty for each year (e.g., the standard errors around the LSMS), uncertainty in the indices of abundance would be more accurately quantified and CAA stock assessment estimates would also likely be more accurate. This benefit would accrue even in areas where CPE and model-based 28 indices Showed Similar temporal patterns (i.e., PD did not show any trends or systematic temporal patterns). We do not believe that calculating a fishery CPE index, by combining CPE each year over strata defined based on statistical modeling, provides a viable alternative to the use of indices directly derived from model-based methods. This conclusion applies especially in the presence of the types of random effects we saw for Great Lakes lake Whitefish data and that appear to be common to fishery CPE data from marine systems. A large advantage of a model-based approach is that the complex correlated error structure resulting from such random effects can be parsimoniously accounted for. The studies cited above suggest that a stratification approach would either underestimate uncertainty in the indices of abundance and lead to inaccurate stock assessment results by ignoring variability attributable to random effects, or would require so many strata with so few observations per stratum that the resulting indices would be poorly estimated. For example, our model for the gill-net fishery would suggest strata need to account for seasonality, boat size, and individual license, but available data only consist of monthly summaries by license. Even if data were combined over similar months, few observations would be available per stratum. Perhaps in some situations (e.g., if random effects were less important), data from each year could be post-stratified into relatively few strata. In such a situation, calculating indices based on combining raw results over strata might be a viable approach, with the advantage of not requiring refitting of statistical models each time a new year of data is collected. An alternative approach to using model-based output as an index of abundance in stock assessments is to integrate the standardization process into the estimation procedure 29 of the stock assessment models (Maunder 2001; Maunder and Langley 2004). Such an approach still models CPE data in the same way as in our analysis here, but integrates the CPE model as a sub-model of the overall assessment. Maunder (2001) found that integrating the CPE standardization into the estimation procedure of the stock assessment model provided a more accurate representation of the uncertainty in stock assessment parameter estimates. The reason for this result, however, was unclear, and so more research is needed in this area, especially given the programming and data management challenges associated with integrating complex GLMM and related models for fishery CPE into assessment models. Standardization techniques used for fishery CPE data cannot ensure that all sources of variation in CPE not attributable to changes in abundance have been considered. For example, changes that are confounded with year and universally affect the fishing fleet, or density dependent changes in catchability, cannot be accounted for using model based standardization methods. Factors left untreated by standardization methods should be addressed in the stock assessments where the CPE indices of abundance are used, for example by allowing for time-varying catchability (Wilberg and Bence 2006). The factors in the final models for both the gill-net and trap-net fishery were similar to models developed for other fisheries (Maynou et a1. 2003; Battaile and Quinn 2004; Bishop et a1. 2004; Helser et a1. 2004). This commonality suggests that similar factors are likely to be important and necessary for consideration when standardizing CPE data for most fisheries. Year is usually included as one of the explanatory variables because detecting trends through time is often the objective for developing indices of 30 abundance, as in this study (Maunder and Punt 2004). Temporal factors on a finer scale than year have also been included in statistical models used for CPE standardization in order to account for systematic temporal patterns in fish abundance or catchability (Battaile and Quinn 2004). Battaile and Quinn (2004) used a fixed effects analysis of variance to standardize CPE data for the eastern Bering Sea walleye pollock, T heragra chalcogramma, trawl fishery, and found a significant effect of time of day (i.e., categorical variable for daylight versus nighttime hours), with higher catch rates dming the daylight hours. They suggested that catch rates were higher during daylight hours because walleye pollock school during those times, but spread out to feed during nighttime, which reduces catchability. In this study, month was included in the final model for the gill-net and trap-net fisheries, with higher catch rates from October to December. The higher catch rates in those months were likely caused by an increase in the catchability of lake Whitefish facilitated by spawning aggregations, which usually occurs during those times in most areas of the Great Lakes (Becker 1983). The results of these studies suggest that temporal factors that account for systematic changes in fish aggregating behaviors should be considered in models used to standardize CPE data whenever possible Various measures of vessel “power” have also been included in models used for standardizing CPE data. Vessel “power” is any measure of the boat or crew that likely affects catchability, and so affects the indices of abundance that result from CPE data taken from those vessels. In the eastern Bering Sea walleye pollock trawl fishery, longer vessels tended to have higher catch rates than shorter vessels as indicated by the coefficient estimates for each vessel participating in the fishery (Battaile and Quinn 31 2004). For the trawl fishery directed at Norway lobster, Nephrops vorvegicus, and deep- water red shrimp, Aristeus antennatus, in the northwestern Mediterranean Sea, generalized linear models used for CPE standardization included measures of the gross tonnage of vessels, engine horsepower, and total length (Maynou et a1. 2003). Generally, longer more powerful vessels had higher catch rates. In the absence of direct measures of vessel power, some surrogate could also be used. For example, Punt et a1. (1996) included the number of crew on the vessel as a surrogate for vessel length in generalized linear models used to standardize albacore, T hunnus alalunga, longline CPE data. For the lake Whitefish fishery in this study, a categorical effect of vessel length was used for the gill-net fishery as a measure of vessel power, but the affects on CPE were inconsistent across lakes. This inconsistency makes broad conclusions about the relative success of various vessel Sizes difficult, but the explanation may be in the characteristics of the lakes themselves. The depth gradient of Lake Superior is relatively steep and permits access to fishing grounds by all boat sizes, and so all boat sizes performed similarly. Conversely, Lake Michigan offers more shallow fishing grounds that are more accessible to small and medium Sized boats, and this may have resulted in higher catch rates than longer boats in that lake. The reason for the relative performance of each boat size in Lake Huron, however, is not clear. A factor for individual vessel, such as license holder in this study, is also commonly included in models for CPE standardization (Maynou et a1. 2003; Battaile and Quinn 2004; Bishop et a1. 2004; Cooper et a1. 2004; Helser et a1. 2004). Similar to results here, an individual vessel factor explained the most variability in CPE in the eastern Bering Sea walleye pollock trawl fishery (Battaile and Quinn). Generalized linear 32 models that included vessel also explained the most variation in CPE for the deep-water red shrimp trawl fishery in the Mediterranean (Maynou et a1. 2003). Cooper et a1. (2004) and Helser et a1. (2004) also found that individual vessel and interactions with vessel should be included in the final models used to standardize U.S. west coast groundfish bottom trawl surveys. The results of Cooper et a1. (2004) and Helser et al. (2004) suggest that even with survey data, standardizing CPE may be necessary, and the availability of model-based indices should not replace the use of consistent survey sampling. The consistent inclusion of an individual vessel effect indicates that individual vessel may serve as a “catch all” for characteristics of boats not included in models (Battaile and Quinn 2004). For example, Maynou et a1. (2003) suggested that the inclusion of individual vessel likely accounts for the expertise of individual fishers or unmeasured technical characteristics, such as investment in technology. The large amount of variation explained by the random effect of license holder and interactions with license holder in this study for both fishery gears also suggests that this factor is accounting for the effects of some unmeasured characteristics, such as those suggested by Maynou et al. (2003). Making inference about the causal or biological mechanisms for some of the two- and three-way interactions included in the final models in this study is not straightforward. However, as Battaile and Quinn (2004) note, identifying causal mechanisms is not required when standardizing CPE data, because the purpose is to account for effects coincident with the variables included in the model. So, the specific higher order interactions may not be indicative of anything biologically meaningful, only 33 that CPE varies coincident with combinations of those factors, either due to those factors themselves or other variables that co-vary with them. The random effect of grid was not included in the final models for either the gill- net or trap-net fisheries, which is surprising considering that typically there is spatial variation in fish density or fishing success. Campbell (2004) found that non-randomly sampled locations led to biased indices of abundance, unless the total habitat area of the stock was spatially stratified and each CPE observation was weighted by the relative amount of sampling effort in the strata from where the observation was taken. This result suggests that not accounting for spatial variation in sampling effort can lead to biased indices of abundance. The effect of grid in this study may have not been included in final models because the analyses were already run on spatially stratified stocks delineated by management unit. However, the results of Campbell (2004) and the spatial variability that likely exists in fish density and fishing success for most fisheries suggests that spatial effects should always be considered when standardizing CPE data. Acknowledgements We are grateful to Michigan Department of Natural Resources (MDNR) and Chippewa-Ottawa Resource Authority personnel who collected and provided data for this project. We thank the 1836 treaty waters modeling subcommittee and personnel at the Quantitative Fisheries Center at Michigan State University for commenting on presentations based on this research. We would especially like to thank Mark Ebener and Phil Schneeberger for their insights on characteristics of the lake Whitefish fishery and these data. Discussions with Ty Wagner, Gretchen Anderson, and Melissa Mata on mixed models and model selection aided this research. Funding for this project was 34 provided by the MDNR and US. Fish and Wildlife Service - Federal Aid for Sportfish Restoration (Project F -80-R, Study 230713). This manuscript is publication 2009-06 of the Quantitative Fisheries Center at Michigan State University. 35 References Battaile, BC. and TJ. Quinn 11. 2004. Catch per unit effort standardization of the eastern Bering Sea walleye Pollock fleet. Fisheries Research 70(2004): 161-177. Becker, G.C. 1983. Fishes of Wisconsin. The University of Wisconsin Press. Madison, Wisconsin. Bishop, J ., W.N. Venables, and Y-G. Wang. 2004. Analyzing commercial catch and effort data from a Penaeid trawl fishery: A comparison of linear models, mixed models, and generalized estimating equations approaches. Fisheries Research 70(2004): 179-193. Brown, R.W., M. Ebener, and T. Gorenflo. 1999. Great Lakes commercial fisheries: historical overview and prognosis for the future. Pages 307-354 in W. W. Taylor and C. P. Ferreri, editors. Great Lakes fishery policy and management: a binational perspective. Michigan State University Press, East Lansing. Burnham, K.P., and DR. Anderson. 2004. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, Second Edition. Springer-Verlag New York, Inc. New York. Campbell, RA. 2004. CPUE standardisation and the construction of indices of stock abundance in a spatially varying fishery using general linear models. Fisheries Research 70(2004): 209-227. Cooper, A.B., A.A. Rosenberg, G. Stefansson, and M. Mangel. 2004. Examining the importance of consistency in multi-vessel trawl survey design based on the US. west coast groundfish bottom trawl survey. Fisheries Research 70(2004): 23 9- 250. Deroba, J.J., and J .R. Bence. In press. Assessing model-based indices of lake trout abundance in 1836 Treaty waters of Lakes Huron, Michigan, and Superior. Michigan Department of Natural Resources: Fisheries Research Report, Ann Arbor, MI. Ebener, MP. 1997. Recovery of lake Whitefish populations in the Great Lakes: a story of successful management and just plain luck. Fisheries 22: 18-20. Ebener, M.P., J .R. Bence, K. Newman, and P. Schneeberger. 2005. Application of statistical catch-at-age models to assess lake Whitefish stocks in the 1836 treaty- ceded waters of the upper Great Lakes. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 271-309. 36 Ebener, MP, and D.M. Reid. 2005. Historical context. In The state of Lake Huron 1999. Edited by MP Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 9-18. Gelman, A., and J. Hill. 2007. Data Analysis Using Regression and Multilevel Hierarchical Models. Cambridge University Press. New York, New York. Harley, S.J., R.A. Myers, and A. Dunn. 2001. Is catch-per-unit-effort proportional to abundance? Canadian Journal of Fisheries and Aquatic Sciences 58: 1760-1772. Helser, T.E., A.E. Punt, and RD. Methot. 2004. A generalized linear mixed model analysis of a multi-vessel fishery resource survey. Fisheries Research 70(2004): 25 1-264. Jensen, AL. 1976. Assessment of the United States lake Whitefish fisheries of Lake Superior, Lake Michigan, and Lake Huron. Journal of the Fisheries Research Board of Canada 33: 747-759. Koelz, W. 1926. Fishing industry of the Great Lakes. Pages 554-617, In Report of the US. Commissioner of Fisheries for 1925. Maunder, MN. 2001 . A general framework for integrating the standardization of catch per unit of effort into stock assessment models. Canadian Journal of Fisheries and Aquatic Sciences 58: 795-803. Maunder, MN. and AD. Langley. 2004. Integrating the standardization of catch-per- unit-of-effort into stock assessment models: testing a population dynamics model and using multiple data types. Fisheries Research 70: 389-395. Maunder, MN. and A.E. Punt. 2004. Standardizing catch and effort data: a review of recent approaches. Fisheries Research 70(2004): 141-159. Maunder, MN. and P.J. Starr. 2003. Fitting fisheries models to standardized CPUE abundance indices. Fisheries Research 63(2003): 43-50. Maynou, F., M. Demestre, and P. Sanchez. 2003. Analysis of catch per unit effort by multivariate analysis and generalised linear models for deep-water crustacean fisheries off Barcelona. Fisheries Research 65(2003): 257-269. McCulloch, CE. and SR. Searle. 2001. Generalized, Linear, and Mixed Models. New York: John Wiley and Sons, Inc. Mohr, LC, and MP. Ebener. 2005a. The coregonine community. In The state of Lake Huron 1999. Edited by MP. Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 69-76. 37 Mohr, LC, and MP. Ebener. 2005b. Description of the fisheries. In The state of Lake Huron 1999. Edited by MP. Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 19-26. Ngo, L., and R. Brand. 1997. Model selection in linear mixed effects models using SAS proc mixed. SAS Institute Inc., Proceedings of the 22nd Annual SAS Users Group International Conference: 1335-1340. Punt, A.E., A.J. Penney, and R.W. Leslie. 1996. Abundance indices and stock assessment of south Atlantic albacore. Collective Volume of Scientific Papers of the International Commission for the Conservation of Atlantic Tunas 43: 225-245. Quinn, T.J., II, and RB. Deriso. 1999. Quantitative Fish Dynamics. Oxford University Press Inc. New York, New York. Rose, GA. and D.W. Kulka. 1999. Hyperaggregation of fish and fisheries: how catch- per-unit-effort increased as the northern cod declined. Canadian Journal of Fisheries and Aquatic Sciences 56(supplement 1): 118-127. SAS. 2003. SAS version 9.1 help and documentation. Cary, North Carolina: SAS Institute, Inc. Smiley, CW. 1882. Changes in the fisheries of the Great Lakes during the decade, 1870-1880. Transactions of the American Fish-Cultural Association 11: 28-37. Venables, W.N., and CM. Dichmont. 2004. GLMS, GAMS, and GLMMS; an overview of theory for applications in fisheries research. Fisheries Research 70(2004): 319- 337. Wilberg, M.J. and J .R. Bence. 2006. Performance of time-varying catchability estimators in statistical catch-at-age analysis. Canadian Journal of Fisheries and Aquatic Sciences 63: 2275-2285. 38 g. 0000“ O .. N 3 ”‘ ° {>- Z‘ ' w B s S W PS -04 wrs-os M8416 -. Ca nada t Michigan WES-08 WFH-OI WFM-Ol WFM-03 {S ‘1“ WFH-03 DD I . 9 WFM-O-I 'e, :0. WFM-OZ . wru -04 ”04% (‘4 9&0» ‘ 1696 ’\ reaty Bounds ‘ . %, Mthban WFM-07 Treaty Boundary LWFM-OS \ Treaty Boundary ? Figure l.— 1836 Treaty-ceded waters and lake Whitefish management units in Lakes Superior, Huron, and Michigan (Ebener et al. 2005). 39 1.0 7 0.8 . Gill-net Fishery Average Coefficient Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 2): 1 Trap-net Fishery 0.4 — I” 0.2 — 0.0 T I I . -0.2 - -0.4 - Average Coefficient -08 - r- -1.0 - Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Figure 2.—-Average coefficient estimates (i SE representing uncertainty resulting from variability among management units) for the effect of month from a general linear mixed effects model standardizing catch per effort (catch = aggregate round mass of lake Whitefish) for the lake Whitefish gill-net fishery (top panel; effort = aggregate length of net in 10003 of feet) and trap-net fisheries (bottom panel; effort = number of lifts) in the 1836 treaty-ceded waters of Lakes Superior, Huron, and Michigan. Coefficient estimates were averaged across various years (generally 1981-2001) and lake Whitefish management units included in this analysis. 40 EmEo w.l_u~ouo&o=m_ @5238 c0250: :8 :33 o». ”655880 W08 m manna. :58: :1on Song ch; Sea 3:88 3025 mg «an: mass can 88: van ammo: Ans—mo 0». SEE 53m om E8 «$183: 8 mmmnomfia man» 0». 35m": 0». no” men own: «85 ~33 m mEéfi E8 Scream: 3303. men «5105 v.33 Amazons—J. $398: mam E8 2385: BgmmoBaB 5:8 5 90 5.3 333,638 (<80; om buxom manger ISO? mam Ewe—ammo: 50:58 E "Em gamma. 1 4 OJ 1 Proportional Difference —l N Ewan». w AAHIbA A03 38 Ammm Awmu Ammo AooA Ammw Amom Aomw. Aomm NooA Proportional Difference A A A A A A A A A A 52.2 A A A Al 4 A H A A A J 0 A03 Ammw Ammm Ammw Aomm AomA Amow Aomm Aomu Ammo NooA A A1 4 J A A H A a A A A A A A 1 A H4 A A A I aba <&.:m80 00:2 10:3 €388.32 883083 0083088 twang-888 2:580 roar 028:8 350on 3053.8 <85 on .0338 003083 30003.08. 00:2 mmomumama 003053..” 1.8805 00:083-.” 0030833u - 380305563000 00883 003062 438030556380 00308.: F 0800 8800303 mmomumama 0800 80000303 - 0085:: 00:083-." - 00:082-.” 80000302 - - 2080 003063 008: - - - - - 05.038 1.80:0? 3.3.3.8 Em: 0m 084.86.08.83 2 2.058.: 0.03000 >88 0 4:80:05 003063 003083 0032 888303 - mmomumsma 00:083-." 00:083..” - ...—380303563000 380305563000 00:08:» h 8000 - 00:05:...” 08000 8800303 - 00:08:. 000 - - 8800302 - - 00:083-.” - - - - - 5030 003053 008: - - - - - 100 H020 _ A8000. 0032 5.030 00.80 .2080 20 mung. .0 080: 0.00 mama—080 5.3.. ... 030.. 0.00 803080 0020.080 c8020.080 2.0x 30x :000=0_080 30.. 008.080 30x .8033303 3003.583 ”003.58.: $003.58.: 9.00.0080 €808.32 980.880 00:3 >0.:080 v0.24 0.0.8.. €808.52 00003083 0083088 10.03083 >0..=080 00.0008 35000:" 2.8.3.00 080: ”08.5.5 ....30 4200305563000 00000 - - - - - 00:08:...n - - - - - 020008 35000:” 3.3.3.00 8002.3. .0 50.0 0.. .3008 3083.030. 00:08:. 008: - $808.0 - 00:08:...u - 00:08:. 00:08:” 00:08:...” - 80000303. - 80000303 - 380305563000 00000 - 00:08:.-.” - - - 00:08:. 80000302. 101 HmEn N. $353 an: 8:633 320$ 3:98 won 5895 8:55: maroQ ogaomgm ”Ewan—Sm so 9.3:: 8:338 ow floor £39 2:: 9a 5288: 2.9.3:: 83828 om man—A an? wan can. an 25.2: 3883330: 5 :8 88355: :22 ow 82.5883 ab; ”55688 cm a ago—427:0: mfioowémoam 3:983. mvoommo 838— 2.3 25 05388133 523g 5 80: 333 80 59083 23.. w X. 38:35: ”32:32: «30: min 55. mag-33m 02:3. ”53 3- «2.3.3 Canofiamnwn 02.3.32. 20 E42 mason on 0." am ...—:6ng mm 000 .3 03.835 «55:02 33.3.3 <35 2. 13:8 30.22 8mm 5%.: man £26. 5?. ...mcn m» m_. 88 ammo man mama. 3.5 mmma Bum . ..mooomo: man ._.m<_o.. 8mm xooaom man mzcfimq 33 12.32. 83 . Own man 129: Sam wczmgonz man away 88 005m man mnmzomq Emu momma 88 mnmiwzmas 8mm ..maam m." 2. 88 $268 mag nmasm 80m mam: m" m_. 33 ..mzam m" m_. 83 95:09 83 0.22% mag Ema—imam: 8mm . 3.83% m» m... .68 Emam man 1330: Non: 230x man mo=m§maa moon <><><><><><>< ><><><><><><><><><><><><><>< ><><>< ><><><><>< ><><><><><><><><>< ><><>< >< .><><><><><><><><><>< ><><><>< ><><><>< ><><><><>< ><><><><>< X X ><><><><><><>< 102 132a A8313 382.32: ”33.5.2: 900.» «.3 .55. mania-8m 02:3. wig .8- 9:23 Cane—3.22. 02.3.32. 20 m22. m32. no 0.“ 0m 4:325... mm 000 ..c 0632.3 .6323... 3.3.3.3 2.»: 2 0572222820: 2. $9.39.: 20393 >303 m .3325... X x x x 40.qu man Imam ~09. xmficxmém moon .....mmma m. m_. woom $3.88: m2. 2.3. 83 ..moocmo: 2:. 43.2 .8... .8038 m2. 9.22 .03 02... m2. Ian... .80 0:69: .00» 9.6263... man. wma... 88 Amcoma Bow . 992 m2. 3.8.8 Bow mfismzmaz Bow N326 m. m.. 88m ..msam m. m_. .58 ..maam m. m... .03 0.22% 2:. 52.63: 88 3.258.» m. m.. 88 m_mamx m2. mo=m§m3= moon <><><><>< ><><>< ><><><><><><><><><><><><><><><><><>< ><><><><><><>< X X X x X x x x 0362.5 352.2: 3.3.3.3 m8...» ”2.55% ......8 ><><><><><><><><><><><><><><><><><*><><><>< x X XX 7.6.22.2 m2. 03333 .03 x x x - x . 0:5: 2 m... .80 x . x w «650.60.. 2 m_. .80 x xx— X x x " 103 Ada—n N. A8513 33.35: ”33:32: ra<><><><><><>< ><><><><><><>< x X xx xXxx ><><><><>< 033:5 mcaozoan 3.33:3 5:32.? ... <3... on 3.01.8 303.353 ><><><><>< >< X><><><>< XX X ><><>< ><><><>< ><><>< X x XXX x ><><><><><><>< ><><><><><>< >< ><><><><><><><>< '><>< 104 Hugo w. mammam ...m. 23:32. £595 30.8.8. Bmfiosoo .555 9%. amt 2038308 Oman: r03.» Zm< 0.2x .39 83.88: mi. 2.9. .83. 5%.: .03. mammaiam .39 933.53 man msmcsma .omfl 1.603 mi. 2268 .8». Oman... m2. .<.m:o: .8... .s>< 9332.3 .39 Own 0. m.. .03. 2.59.6: m2. 50.3 .89 >255 .8w .<_m< 291mm manage: mag 000.8 .89 >::m_m .88 23.3 32.8.3. ..o228: man «.238 .88. 32.05 8.3. Dmamo .8» 33me m o.. \HXo\o 30 man. 333. 3.6... 30:. .um. 9233 5890.3 3.0.8. .39 933.53 m2. mzmnzma .omfl 09.8 .8». Omnmo .83. 0.3. .8.“ ..mxoomm: .8». 008.32 .88. ..mmBm: .8w. Imam. m2. manumx .89 00:6 man Gimme: moo.“ 02:83 noon. mm:.xm.nm3 man 93:98: noon 983.5% man 2626:. .83. 0320. m. m.. .88. 0.9.x .8.. 0003.3. .88. ..mxoamm: .89 ....mom mag mammaism .89 .uc.._oxm m. m.. .8N maammx mag 2.18:. .89 0.2.» .89 0.qu .88. 005m m2. 0.5.80: woo.“ 0.qu noon. <<.__.m3m woo». woo... woo... mm...xm.:m: m2. mficsmamoa moi mammaéim man mzmuzma .83 $33.25... m2. mzmuzma .83. .5533 .89 ....mom man 933.53 .89 .5858 man 2.25 .88. 00:6 man Gimme: moo.“ 0:5: man wNmfl. .89 :33 ~08 ...—.3255 .93.» 9...... m. m_. .89 9m} .8.“ ....mom m2. m.mmm:<<.:m .89 O~E=a~=8 wchumm S. >c:=&=.8 52:93.... >w==as=8 Ema?" rimmmwo 833— ::8 mag :05 mmasm Bog—5 masaazw 0:858 2:: 92.5% 2. 3.5958 m3 man: 38. 106 3 drama—SE 5 0c=a3o=£ 0258.: 038.. — 0v 02:—ES.»— 025359 0:8: N Fishing Mortality 5 wmcawmm $62.. 5 ESE—mm wumaa- UV vac—523:». diam—SE m. 5:8... was—=68 0&8: 5:8... Hus—3.8 Fishing Mortality wmcauuu 2. 2:59:88 ESE—mm S. >U==aa=oa ESE-mm S. >r==aa=nm 3min N.l150‘y y 1l; (19) Zy>150Yy where 124 —z Fa,yNa,y(1— e a’y] W (20) a,y ' Za’y _ a=12+ Yy — Za=3 The trade-off plots described above were then constructed for each of three percentiles calculated among simulations: the median, the percentile value less than 75% of the values among simulations (i.e., 25th percentile for SSB and Y, 75th percentile for risk and Yvar), and the percentile value less than 90% of the values among simulations (i.e., 10‘h percentile for SSB and Y, 90th percentile for risk and Yvar). The rank order performance of the control rules was evaluated for each growth i scenario and trade-off plot by examining the combination of policy parameters that resulted in the best performance for each control rule. This method equates to determining the optimal control rule as if future growth (i.e., autocorrelated versus recent growth patterns) was a known certainty. To evaluate the sensitivity of control rule performance to uncertainties about future growth, an optimal set of policy parameters (see below) was chosen for each control rule, trade-off plot, and growth scenario. The optimal set of policy parameters for each future growth scenario was then applied to the alternative fiiture growth scenario (i.e., autocorrelated versus recent growth patterns), but for the same type of stock growth (i.e., fast or slow growing). The difference in performance between the optimal set of policy parameters and the set of policy parameters being applied under the wrong future growth scenario was used as the measure of sensitivity. This method equates to evaluating how well a control rule would perform if policy parameters were chosen assuming one type of future growth was true, when in fact the alternative growth future 125 was true. The optimal set of policy parameters for each trade-off plot was defined as the set that: maximized yield and maintained SSB above 50% of SSBF=0, minimized variability in yield and produced risk less than 0.40, maximized yield and produced variability in yield less than 0.4, and maximized yield and produced risk less than 0.40. This method of selecting optimal policy parameters was only used to illustrate the sensitivity of control rule performance to uncertain future growth, and because this does not necessarily reflect desired tradeoffs, should not be used for management without careful consideration of fishery objectives. 3. Results 3.1. Overview of results Varying the level of implementation error had little effect on the relative or absolute performance of the control rules, and consequently, all the results below are for the baseline level. Varying the parameters related to assessment error affected the relative performance of some control rules, but only for trade-off plots that included variability in yield (Section 3.4). For the other performance metrics, the relative and absolute performance of the control rules varied little among the different levels of assessment error parameters (see Appendix B), and this was consistent among growth scenarios. Thus, the choice of optimal parameters would also generally be robust to the level of assessment error. As a consequence, results for the sensitivity of the control rules to varying assessment error parameters (Section 3.4) include example graphs based on simulations of fast growth similar to more recent levels, and the results in all other sections are for baseline levels of assessment error. Results for the different percentiles showed similar trade-offs among performance metrics and relative differences among 126 control rules; so only results for the median are considered further. The rank order performance of the control rules for each trade-off plot was also the same for all levels of each source of uncertainty (see Appendix B), and these rank orders are presented in Section 3.2 with example graphs based on simulations of fast growth similar to more recent levels. 3. 2. Rank-order performance of control rules For the plot of SSB versus Y, the BB control rule performed best, providing more Y at a given level of SSB and higher SSB at a given level of Y than other control rules (Figure 3). The BB control rule was followed in performance by constant-F, CCC, and the BB-lim and constant-Flim control rules (Figure 3). For the plot of Yvar versus risk, the CCC, BB-lim, and constant-Flim control rules provided less Yvar at a given level of risk than other control rules (Figure 3). These same control rules, however, were also more risky at a given level of Yvar than other control rules (Figure 3). For the plot of Y versus Yvar, the CCC, BB-lim, and constant-Flim control rules provided less or similar Yvar at a given level of Y than other control rules (Figure 4). The BB control rule, however, attained more Y at a given level of Yvar, and was followed in performance by the constant-F, CCC, and the BB-lim and constant-Flim control rules (Figure 4). For the plot of Y versus risk, the BB control rule performed best, providing more Y at a given level of risk and less risk at a given level of Y (Figure 4). The BB control rule was followed in performance by constant-F, CCC, and the BB-lim and constant-Flim control rules (Figure 4). 127 3.3. Sensitivity to future growth uncertainty 3.3.1. Fast growth stocks For SSB versus Y, yield was generally insensitive to future growth, and the BB and constant-F control rules were less sensitive than other control rules. All of the percent changes in yield from the optimum were less than 6%, as were the changes in SSB for the BB and constant-F control rules (Table 2). For the CCC, BB-lim, and constant-Flim control rules, the percent changes in SSB were at least 10%, and SSB decreased from the optimum and below 50% of SSBF=0 (i.e., the level used to define optimal) when policy parameters were chosen as though grth would be similar to recent levels and autocorrelated growth was the true future, but increased from the optimum for the opposite situation (Table 2). So, control rules and policy parameters chosen by incorrectly assuming future growth will be autocorrelated cost little in yield and produced more SSB than the optimum, relative to incorrectly assuming future growth will be similar to recent levels, which produced less SSB than the optimum. For Yvar versus risk, results were insensitive to the future growth scenario. For all control rules, the percent change from the optimal levels was 0.00 (Table 2). For Y versus Y var, results were generally insensitive to the future growth scenario. All of the percent changes in yield were less than 5%, and the changes in variability in yield were all less than 8% (Table 2). Variability in yield increased from the optimum and, for the BB rule, above 0.40, when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true future, but decreased from the optimum for the opposite situation (Table 2). So, control rules and policy parameters chosen by incorreCtly assuming fiiture growth will be autocorrelated 128 cost little in yield and produced less variability in yield than the optimum, relative to incorrectly assuming future growth will be similar to recent levels, which produced more variability in yield than the optimum. For Y versus risk, yield was generally insensitive to the future growth scenario, but risk was more sensitive. All of the percent changes in yield were less than 5%, but the changes in risk were more variable (Table 2). Risk increased from the optimum and, for the BB-lim and constant-Flim rules, to 0.40, when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true future, but decreased from the optimum for the opposite situation (Table 2). So, control rules and policy parameters chosen by incorrectly assuming future growth will be autocorrelated cost little in yield and produced less risk than the optimum, relative to incorrectly assuming future growth will be similar to recent levels, which produced more risk than the optimum. 3.3.2. Slow growth stocks For SSB versus Y, the CCC, BB, and constant-F control rules were less sensitive to the future growth scenario than the BB-lim and constant-Flim rules. All of the percent changes in SSB and yield were less than 5% for the CCC, BB, and constant-F control rules (Table 2). For the BB-lim and constant-Flim rules, however, yield decreased and SSB increased from the optimum when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true future (Table 2). Results for these control rules were less consistent when policy parameters were chosen as though growth would be autocorrelated, but growth similar to recent levels was the true future. For the BB-lim rule, both yield and SSB decreased from the optimum, and 129 SSB was less than 50% of SSBF=0 (Table 2). For the constant-Flim rule, yield increased but SSB decreased below 50% of SSBF=0 (Table 2). So, the costs and benefits of choosing policy parameters based on assuming an incorrect future for lake Whitefish growth depended on the control rule. _ For Y var versus risk, results were insensitive to the future growth scenario, except for the CCC control rule. For the CCC control rule, variability in yield and risk increased _ L from the optimum when policy parameters were chosen as though growth would be 5 similar to recent levels and autocorrelated growth was the true fiiture (Table 2). When H policy parameters were chosen as though growth would be autocorrelated and growth similar to recent levels was the true future, the percent changes were less than 3%, and variability in yield increased while risk decreased from the optimum levels (Table 2). So for the CCC control rule, policy parameters chosen by incorrectly assuming future growth will be autocorrelated cost little in variability in yield and produced slightly less risk than the optimum, relative to incorrectly assuming future growth will be similar to recent levels, which produced more variability in yield and risk than the optimum. For Y versus Y var, sensitivity to the future growth scenario depended on the control rule, but the BB control rule was the most sensitive. For the BB and constant-F control rules, yield and variability in yield increased from the optimum (above 0.40 for variability in yield) when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true future, but decreased from the optimum for the opposite situation (Table 2). For the other control rules, yield and variability in yield decreased from the optimum when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true 130 future (Table 2). When policy parameters were chosen as though growth would be autocorrelated and growth similar to recent levels was the true future, yield decreased and variability in yield increased from the optimum (Table 2). So, the costs and benefits of choosing policy parameters based on assuming an incorrect future for lake Whitefish growth will depend on the control rule. For yield versus risk, results were generally more sensitive to the future growth scenario than other trade-off plots. For all control rules, yield and risk decreased from the optimum when policy parameters were chosen as though growth would be similar to recent levels and autocorrelated growth was the true future, but in the opposite situation, yield decreased or was the same and risk increased from the optimum, and above 0.40 for the BB-lim rule (Table 2). So, control rules and policy parameters chosen by incorrectly assuming future growth will be similar to recent levels produced costs in yield at the benefit of less risk than the optimum, relative to incorrectly assuming future growth will be autocorrelated, which produced costs in yield and risk from the optimum. 3. 4. Sensitivity to assessment error parameters Varying assessment error parameters affected the performance of the control rules for trade-off plots that included Yvar. For the plot of Yvar versus risk, control rules that performed well at each performance metric increased in superiority over other control rules as pn was decreased (Figure 5) or 0'3 was increased (Figure 6). Specifically, the CCC, BB-lim, and constant-Flim control rules increased in superiority in terms of Yvar at a given level of risk, and the BB and constant-F rules increased in superiority in terms of attaining less risk at a given level of Yvar. For the plot of Y versus Yvar, the CCC, BB- lim, and constant-Flim control rules increased in superiority in terms of Yvar at a given 131 level of Y as Pn was decreased (Figure 7) or 0'3 was increased (Figure 8), but an increase in the superiority of the BB and constant-F control rules in terms of Y at a given level of Yvar did not materialize. 4. Discussion In this study, the BB control rule provided more yield and less risk than other control rules with all else being equal, and was followed in rank-order by constant-F, CCC, and the BB-lim and constant-Flim control rules, which is consistent with similar research (Irwin et al., 2008; Punt et al., 2008). Irwin et al. (2008) used a similar BB control rule for a recreational yellow perch Percaflavescens fishery in southern Lake Michigan and found that the BB control rules produced higher yields and less risk than a constant-F control rule at given levels of SSB. Punt et al. (2008) used a BB control rule that reduced catch (instead of F) linearly with SSB for groundfish off the US. west coast and found that the BB control rules produced higher yield and less risk than a constant-F control rule at low levels of productivity, but performance was nearly equal for high productivity. The CCC, BB-lim, and constant-Flim control rules provided less interannual variation in yield than other control rules, but at the cost of yield and risk. Clark and Hare (2004) found that the CCC control rule could produce similar yields and SSB than a constant-F control rule, but with less variability in yield for Pacific halibut Hippoglossus stenolepis. Their simulations, however, did not include parameter uncertainty in the stock-recruit relationship or assessment error. The results of this study suggest that including these sources of uncertainty affects the relative performance of the CCC control rule, as has been shown for other control rules (Deroba and Bence, 2008). For roundfish 132 stocks managed by the International Council for the Exploration of the Sea, limits on the interannual change in target catch that prevented quotas from being decreased often resulted in less yield and greater frequency of low SSB than a constant-F control rule (Kell et al., 2006), which is consistent with the findings of this study where the BB-lim and constant-Flim control rules performed worst in terms of yield and risk of low spawning stock with all else being equal (e. g., constrained to provide some specified level of spawning stock or yield). The relative performance of limiting the interannual change in target catch, however, can depend on how tight the restraint is on the interannual change, productivity, variability in recruitment or growth, and current status of the stock (Punt et al., 2002; Kell et al., 2006). The rank order performance of the control rules, while treating future growth as known, was robust to the type of autocorrelated growth evaluated in this study, and the growth scenarios in general, but time varying population dynamics have been shown to affect relative performance (Walters and Parma, 1996; Deroba and Bence, 2008). Walters and Parma (1996) showed that a constant fishing mortality rate control rule performed better in terms of yield when the asymptote parameter of a Beverton—Holt stock—recruit model was autocorrelated, and this was counter to when this parameter of the stock-recruit model was constant. Whether catchability is treated as time-varying in stock assessments also has an effect on control rule performance (e. g., Dichmont et al., 2006). Few studies have evaluated the effect of time-varying parameters, but more research is warranted in this area (Deroba and Bence, 2008). The robustness of control rule performance to uncertainty about future lake Whitefish growth, and the most robust future lake Whitefish growth to assume for 133 selecting policy parameters, depended on whether stock growth was fast or slow. For fast growth stocks, selecting control rules and policy parameters by incorrectly assuming that future growth will be autocorrelated resulted in little cost from the optimum levels relative to the alternative of incorrectly assuming future growth will be similar to recent levels. For slow growth stocks, however, the robustness of which future lake Whitefish growth to assume for selecting policy parameters depended on the control rule and trade- off plot. The decision to select control rules and policy parameters based on some assumption about future lake Whitefish growth will ultimately depend on how competing fishery objectives are weighted relative to each other. Generally, however, control rules and policy parameters for fast growth stocks should likely be selected assuming future grth will be autocorrelated, but a universal recommendation for slow growth stocks is less clear (i.e., depends on the control rule and fishery objectives). The results were generally robust to the level of assessment error. As reported here, Punt et al. (2008) found that assessment error affected results for Yvar, but other performance metrics similar to those included in this study were robust to this source of uncertainty. Irwin et al. (2008) did not include a measure of yield variability, but also found that other performance metrics similar to those included in this study were insensitive to varying assessment error. Conversely, Katsukawa (2004) reported that the superiority of BB control rules over constant-F control rules in terms of yield, diminished with increasing assessment error variance. This contradiction may have occurred because Katsukawa (2004) did not include other sources of uncertainty (e. g., shape of the stock- recruitrnent curve) that may outweigh the effect of assessment error on the relative performance of control rules. Alternatively, Katsukawa (2004) summarized tradeoffs in 134 terms of yield versus the minimum biomass over a time-horizon, and stocks managed in the face of greater assessment uncertainty might suffer more extremes in stock size. In this study, limiting the interannual change in target catch by 15% was inferior or at best similar to other control rules for all performance metrics, and in some cases was more sensitive to uncertainty in future lake Whitefish growth. So, using alternative control rules would likely cost little relative to limiting the interannual change in target catch by 15%, and would produce benefits in most situations. This result may not be general, however, as benefits associated with limiting the interannual variability in target catch depend on fishery objectives, the degree to which the interannual change is restrained, stock status, and other population parameters such as growth variability (Kell etaL,2006) Depending on the relative weight of fishery objectives, a control rule and policy parameters other than the one currently in use (i.e., constant-F based on a total annual mortality rate of 65%) may want to be considered for lake Whitefish populations in Lakes Huron, Michigan, and Superior. For example, a BB control rule with appropriately selected policy parameters could likely produce nearly the same or more yield, spawning stock biomass, and less risk with little cost in variability in yield relative to the currently used policy. Similarly, the CCC control rule can likely provide less variability in yield, but at the cost of yield. So, if maintaining low variability in yield is more desirable than maximizing yield, a CCC control rule may want to be considered. Not all dynamics or uncertainties about lake Whitefish were included in this study, and unanticipated changes in the future may require periodic reviews of this evaluation (Butterworth, 2008). For example, this study did not include density dependent growth, 135 which 1115 i this 4 L1‘ 4‘ 0 cause a like coal it on But: Mu ("014115 Iltl? which has been shown to occur for lake Whitefish (Henderson et al., 1983; Kratzer et al., 2005). Density dependence was not included because the recent declines in lake Whitefish growth that have occurred in some areas generally happened over a time period when abundance has declined or remained relatively stable, and so could not have been caused by density dependence (e. g., Lumb et al., 2007). Density dependence is, however, a likely compensatory response in lake Whitefish and so should be considered if conditions arose to make such dependence important. To address such changes that may be outside the realm of uncertainties included in a management strategy evaluation, Butterworth (2008) recommended scheduling periodic reviews to consider whether evaluations should be updated. If radical unanticipated changes occur, Butterworth (2008) recommended making an ad hoc adjustment to the pre-agreed control rule until the management strategy evaluation can be updated and tested for robustness to the recent changes. Alternatively, stakeholders could agree on a pre-determined default management plan that would be applied temporarily until the management strategy evaluation is reviewed (Butterworth, 2008). Such scheduled maintenance of this evaluation would also be prudent given the uncertainties and variability in lake Whitefish population dynamics. 136 Bear: 3901’: Bro“ Bun: C00 Der References Beauchamp, K.C., N.C. Collins, and BA. Henderson. 2004. Covariation of growth and maturation of lake Whitefish. Journal of Great Lakes Research 30(3): 451-460. Booth, A.J. 2004. Determination of cichlid specific biological reference points. Fisheries Research 67: 307-316. Brown, R.W., M. Ebener, and T. Gorenflo. 1999. Great Lakes commercial fisheries: historical overview and prognosis for the future. Pages 307-354 in W. W. Taylor and C. P. Ferreri, editors. Great Lakes fishery policy and management: a binational perspective. Michigan State University Press, East Lansing. Butterworth, D.S. 2008. Some lessons from implementing management procedures. Pages 381-397 in K. Tsukamoto, T. Kawamura, T. Takeuchi, T.D. Beard, Jr., and M.J. Kaiser, editors. Fisheries for global welfare and environment, 5th world fisheries congress 2008. TERRAPUB, Tokyo. Clark, W.G. 1991. Groundfish exploitation rates based on life history parameters. Canadian Journal of Fisheries and Aquatic Sciences 48: 734-750. Clark, W.G., and SR. Hare. 2004. A conditional constant catch policy for managing the Pacific halibut fishery. North American Journal of Fisheries Management 24: 106-113. Cook, H.A., T.B. Johnson, B. Locke, and B.J. Morrison. 2005. Status of lake Whitefish in Lake Erie. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 87-104. Deroba, J .J . and J .R. Bence. 2008. A review of harvest policies: understanding relative performance of control rules. Fisheries Research 94: 210-223. Dichmont C.M., A. Deng, A.E. Punt, W. Venables, and M. Haddon. 2006. Management strategies of short-lived species: The case of Australia’s Northern Prawn Fishery 3. Factors affecting management and estimation performance. Fisheries Research 82: 23 5-245. Ebener, M.P., J .R. Bence, K. Newman, and P. Schneeberger. 2005. Application of statistical catch-at-age models to assess lake Whitefish stocks in the 1836 treaty- ceded waters of the upper Great Lakes. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 271-309. 137 Ebener '1. ltleh Heels Hand H031 lm‘ir Ebener, MP, and D.M. Reid. 2005. Historical context. In The state of Lake Huron 1999. Edited by MP. Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 9-18. Fujioka, J .T., J. Heifetz, and M.F. Sigler. 1997. Choosing a harvest strategy for sablefish based on uncertain life-history parameters. Pages 247-251 in NCAA Technical Report NMFS 130 Biology and Management of Sablefish; Papers from the International Symposium on the Biology and Management of Sablefish, Seattle. Healey, MC. 1975. Dynamics of exploited Whitefish populations and their management with special reference to the Northwest Territories. J. Fish. Res. Board Can. 32: 427-448. Henderson, B.A., J .J . Collins, and J .A. Reckahn. 1983. Dynamics of an exploited population of lake Whitefish in Lake Huron. Canadian Journal of Fisheries and Aquatic Sciences 40: 1556-1567. Hoyle, J .A., T. Schaner, J .M. Casselman, and R. Dermott. 1999. Changes in lake Whitefish stocks in eastern Lake Ontario following Dreissena mussel invasion. Great Lakes Research Review 4: 5-10. Irwin, B.J., M.J. Wilberg, J .R. Bence, and ML. Jones. 2008. Evaluating alternative harvest policies for yellow perch in southern Lake Michigan. Fisheries Research 94: 267-281. Jacobson, RC. and W.W. Taylor. 1985. Simulation of harvest strategies for a fluctuating population of lake Whitefish. North American Journal of Fisheries Management 5: 537-546. Jensen, AL. 1976. Assessment of the United States lake Whitefish fisheries of Lake Superior, Lake Michigan, and Lake Huron. Journal of the Fisheries Research Board of Canada 33: 747-759. Katsukawa, T. 2004. Numerical investigation of the optimal control rule for decision- making in fisheries management. Fisheries Science 70: 123-131. Kell, L.T., G.M. Pilling, G.P. Kirkwood, M.A. Pastoors, B. Mesnil, K. Korsbrekke, P. Abaunza, R. Aps, A. Biseau, P. Kunzlik, C.L. Needle, B.A. Roel, and C. Ulrich. 2006. An evaluation of multi-annual management strategies for ICES roundfish stocks. ICES Journal of Marine Science 63: 12-24. Koelz, W. 1926. Fishing industry of the Great Lakes. Pages 554-617, In Report of the US. Commissioner of Fisheries for 1925. 138 Kratzer, J .F ., W.W. Taylor, C.P. Ferreri, and M.F. Ebener. 2005. Factors affecting growth of lake Whitefish in the upper Laurentian Great Lakes. Advances in Limnology 60: 459-470. Lumb, C.E., T.B. Johnson, H.A. Cook, and J .A. Hoyle. 2007. Comparison of lake Whitefish growth, condition, and energy density between Lakes Erie and Ontario. Journal of Great Lakes Research 33: 314-325. Mohr, L.C., and MP. Ebener. 2005. Description of the fisheries. In The state of Lake Huron 1999. Edited by MP. Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 19-26. Mohr, L.C., and Nalepa, T.F. (Editors). 2005. Proceedings of a workshop on the dynamics of lake Whitefish (Coregonus clupeaformis) and the amphipod Diporeia spp. in the Great Lakes. Great Lakes Fishery Commission Technical Repport 66. Myers RA, Bowen KG, and Barrowman NJ. 1999. Maximum reproductive rate of fish at low population sizes. Can J Fish Aquat Sci 56: 2404-2419. Myers, R.A., G. Mertz, and J. Bridson. 1997. Spatial scales of interannual recruitment variations of marine, anadromous, and freshwater fish. Can J Fish Aquat Sci 54: 1400-1407. Nieland, J .L., M.J. Hansen, M.J. Seider, and J .J . Deroba. 2008. Modeling the sustainability of lake trout fisheries in eastern Wisconsin waters of Lake Superior. Fisheries Research 94: 304-314. Pothoven, S.A., T.F. Nalepa, P.J. Schneeberger, and SB. Brandt. 2001. Changes in diet and body condition of lake Whitefish in southern Lake Michigan associated with changes in benthos. North American Journal of Fisheries Management: 21: 876- 883. Punt, A.E., M.W. Dom, and M.A. Haltuch. 2008. Evaluation of threshold management strategies for groundfish off the US west coast. Fisheries Research 94: 251-266. Punt, A.E., A.D.M. Smith, and G. Cui. 2002. Evaluation of management tools for Australia’s South East Fishery 3. Towards selecting appropriate harvest strategies. Marine and Freshwater Research 53: 645-660. Quinn, T.J., II, and RB. Deriso. 1999. Quantitative Fish Dynamics. Oxford University Press Inc. New York, New York. Quinn, T.J., II, R. Fagen, and J. Zheng. 1990. Threshold management policies for exploited populations. Canadian Journal of Fisheries and Aquatic Sciences 47: 2016-2029. 139 Smiley, CW. 1882. Changes in the fisheries of the Great Lakes during the decade, 1870-1880. Transactions of the American Fish-Cultural Association 11: 28-37. Smith, A.D.M., K.J. Sainsbury, and RA. Stevens. 1999. Implementing effective fisheries management systems — management strategy evaluation and the Australian partnership approach. ICES Journal of Marine Science 56: 967-979. Thompson, G.G. 1993. A proposal for a threshold stock size and maximum fishing mortality rate. Pages 303-320 in S.J. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Walters, C.J., and A.M. Parma. 1996. Fixed exploitation rate strategies for coping with effects of climate change. Canadian Journal of Fisheries and Aquatic Sciences 53: 148-158. Wang, H-Y, T.O. Hook, M.P. Ebener, L.C. Mohr, and P.J. Schneeberger. 2008. Spatial and temporal variation of maturation schedules of lake Whitefish in the Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 65: 2157-2169. 140 lrcan ll urndan O o as W FS-ll-l w rs-os w rsoo . C a n a d a % «\u—g/ kw ’t. Michigan WFS-08 ‘l. WFH-OI \ A WFMM “FM-03 " - . M51) WFH-03 3 . J ' b a " V o' -‘ . 0 wrist—04 ”a. 0; \H‘M-oz , , ‘00,, g . “HI-04 0.4,“. E- o'- g Q :n E. 8 if? WFM-OS Mbhban 2"“ wrM-ov / D Treaty Buundan W F.“ -08 l‘reatv Buu Man of Figure l.— 1836 Treaty-ceded waters and lake Whitefish management units in Lakes Superior, Huron, and Michigan (Ebener et al. 2005). 141 Fsat Fishing Mortality 0.05 SSBLT ssiaHT Spawning biomass relative to the unfished level Figure 2.—Example of the biomass based control rule used in this analysis (solid line). Dashed lines are provided as a reference for defining the policy parameters. 142 EmEo 91780.0: mnmésmsm 0.00.. 0.0.080 88.5 800.8. 0.85 an» 003...... ...... 82...... 0.83.500. $500.:Q ... 0.85 588 33.0.. 300030.. 0m 0503 2.... magnum £00.. 30803 8mm ...»: 83 0.. .30 00.8.50 _m<0_ Em... 0075.... m0. 0.50:8 .95; 0.. 0088.00... 0.... 0010300830.. 0.40. 0080.083 0.... ....m. was... $320.. .0 Ron... .085 $00 8... m0. 00.0.5. m0. .30 000030..»— m 00.5.»... 88.. 000.8. .08 ..00 3.5. 8380.-.“. ...—mow 008. 8.008 3.3. 008.03%. 23. 0 3.x. :0... 0.. ...0 0.83.500. orgma ... SE... 008: $80. 008. 3.008 «0.5. 0.030% 00000 3.00.. 008. 00.8... 82.. ...... 0.080% 00000 2.... 0 3.x. :0... 0.. ...a 5.30.0000— nrwsmo ... Sam... 38: £30. 008. 00.8... «0.5. mun: 00. 00.48.0003 .0 .30 0.8.0.5 ..8... Eco «mac—0.80m m0. 000 80.00.0030: 0.. 00:8. 0030.083 ..0.. 90 m.0 ... .1. 9.8 0. 0.800. ..0. 300.00 .00... 0.25 0.00.00 5.0.0.500. 0100...... ... 0.0... :00. 00.00.... 0.... 300.00 «.0... 0300 0.00.00 0300....0: 0.. 000.0 2.... 0002.....m 0.00.. 0.0.0000 .000 ...00 3.x. 0.. ...0 00.15.00 .00. am... 00.08.... 145 0 r O l 0 (6)1) mam 0 (6):) PlalA 000007 000003 000009 '17 omB! :1 Variability in yield Variability in yield Variability in yield 0 o _. _. N O O -* -* N .0 .O 7‘ .-‘ N 'o '01 'o '01 'o 'o '01 b '01 'o o 01 o 01 o l J l J 000002 000002 000002 000007 00 o o ’0" '0 0000017 0000017 000009 000009 000009 Proportion of years with Proportion of years with Proportion of years with SSB < 20% unfished SSB SSB < 20% unfished SSB SSB < 20% unfished SSB 9.99.0.0? omxsmooo J l l 1 J A o o N N O O O O O O O O O 0 -§ h C O O O O O O O O O O) O) O O O O O O O O O O 146 .u.m0.0 m.l.<.00.0.. 5.0.0.500. 00.00:...N ... 0.0.0 .0000 800.00 0.000....00 0., 000.0 2.... 0002.....m 0.00.. 0.0.0000 .000 0.0.. 8.x. 0.. ...0 0.....0..00 .0<0. ..0. 0000....0 000000.00... 000 ....0.0...0...0..0.. 00.0.. <0..0..00. ...0 9:0... 0.. 00.000.40.000: ... 000000.00... 00.0. 00. 0000. .0 ob 0. 0.0. 000 .,00. E02... 0......0. .0 .000... .0<0.0 .000 .0... .0. 00.0..0v. 000.8. 00.00 0.0 0.00.0000 00... 050.0 m. 147 91915 05111qule maul ua mumemA '9 omfiyd Proportion of years with Proportion 0' years With Proportion of years with 383 < 20% unfished sse 833 < 20% unfished SSB sse < 20% unfished ssa 9.0.0.090 -°-°-°-°-°—“ 9.09.090 ONAmmOOON‘bmmo 010001000 p;1111J'J1¥1119L11141 CI 0 o \ ~. 9‘ .o 0‘ 01 3 .1 -1 ll 0 0 o O -* "« e- 01 0‘ 01 N N~ .N- o O o Proportion of years with Proportion of years with Proportion of years with 883 < 20% unfished SSB SSB < 20% unfished SSB SSB < 20% unfished 333 o o o o o _. o o o o o _. . . . . . . . . . . . . o o o o .. o N :s o: oo o o N .0 o: oo o '0 i0 '0 ‘o: g 'o L 1 1 1 1 4 P 1 1 ¢ L 1 1 L 1 1 1 1 J o .0 0 ° I Q. / .o- o 01 '01" 3 || A s- a T‘. - 01 '0,“ NJ N o 'o‘ 148 050.0 0..I>0 .: 050.0 0. 0.800. 2.... .00 0.30:. 0. 00.000.410.33: 000000.003 0.8. 00. 0000. .0 ...0 0000.50 .0<0. 8.... 0:0 .0. 000000.003 0.3. <0...0:00 0000. .0 ob. 0. one. 149 916M U! mem Q'l new U! 01101211211 S'l '9 amfii :1 Proportion of years with Proportion of years with Proportion of years with 883 < 20% unfished SSB SSB < 20% unfished SSB SSB < 20% unfished SSB .0 .0 .0 .0 .0 r" o N «5 O) on o 9 1 1 1 1 1 1 0 90 UL 0'2 Proportion of years with SSB < 20% unfished SSB 9.0.0990 010001030 | i k 1 4L J 00 i 0'2 O O O O O . '0 '0 '0 '0 0 'o P 4 i L4 0 .0- 01 .O- .01- .N‘ 0 Proportion of years with 888 < 20% unfished SSB OOOOO-t 'oiubbabo'o 9111111 0 .04 UI dd 0 A 01 01 o 150 9.9.0.09.“ 003001000 1 l l L 1 3L7"... 90 — .0 0|. u ”-2 [0'0 Q'l 0'2 Proportion of years with SSB < 20% unfished SSB PPPPP.‘ 01010010110 1 J L 1 1 J 1...... J .9 o .0 OI .0 3N '0“ II p N _. o '0‘ .N. O 050.0 0.110 .0 050.0 0. 0.800. 2.... 800.00 0.0.0 <0.000 800.0: 80.00000. <0...00....0 .0 0.0.0. 151 050.0 0. In yield Variabili In yield ield Variabili ' Iny Variabili Na.1 ..0 1 00 1 00.1 1.05.11» 1. . 1 .00” 0.0 O ““ 0.. N... 00 1 90 1 00w1 . . ~00000 000000 J'- . . .. .ut\v . .. _ 000000 .vo N001 00 1 90 1 00w1 . . 000000 ~00000 0 000000 0.. N00000 000000 50.0 :6. . 000000 N001 00.. ..0 1 .0: 0.0 0 00's N001 ..0 1 00.1 00.1 J1 . ~00000 000000 I o. o 00 0.x. 0003‘ 0 000000 0.. 0 N001 00 1 00.1 .90- 0.. 0 . 000000 N00000 000000 888 388 <00 :6. 000000 152 30:8 m.|1>m 3 Emcaa 0. 98000 $09 05 9:9: om $883330: 5 ammommBoa 039 35— 8 90 ammo—mam _m080. 0:88 035: 0003010 0.5088885 .8003 >Sooo:0_080 30003 00302830. 02883 008: <85 <0. 000 3003 0-3 0.03 :3 000000 T3 0.00 013 30000 3 0.00 3 004000 3 can 3 <<0<0m u a m a N a o .3 ._ ._ ..N Awm #3 wow mam mam mom 2 V mum mmm mum 3 u REV A: Ag 3 o man man— mmo mam mm.“ A3 ANA m._ ._ gm mac 9 w 9.5 a: mom Nam mam wow kmu 3w 5.3 #8 m: A mmm m3. mg 164 Table A3—Upper and lower caps placed on simulated length (mm) at age-3 lake Whitefish for fast and slow growth category simulations. Growth Category Cap Fast Slow Upper 539 601 Lower 228 21 0 165 APPENDIX B This appendix presents additional trade-off plots beyond those presented in the main text. Displaying all the results for four performance metrics, five control rule variants, three levels of productivity, seven levels for the parameters of assessment and implementation error, and four growth scenarios would require 1,512 panels. Displaying all these results in this appendix was not feasible. The subset of trade-off plots shown here support statements in the main text. Figures B 1 -B4 illustrate that varying the parameters related to assessment error had little effect on the relative or absolute performance of control rules for trade-off plots not involving variability in yield. Figures BS-B8 illustrate that the rank order performance of the control rules is robust to the growth scenarios. 166 Emfio 21:32:»: mvmgmsm £0wa 30353 <22; 365 man 38:30 mamammaoa Ba mac—«338:0: 63.2. 553.8. 25 0x82 0». $883388: 3 $8382: 3.8 mg 35; 8 ob 9. ob. Ba mum” mag: £328. 8 883 $85 Ammo 8x" m3 mag—my m3 Em 83525— 8833 88: 838— ::a Qou 32v. 8385:“. Aime—n no? Baa—o 32v. 838:7»... <55 m 3.x. :8: 0: En 583555— 25me 5 Samoa 88: Amaow mow. Baa—o 82V. Eon—gm 633 95% no? 6283 8,5. 25 325mm 3mg 23. m 3.x. :3: o: 90 W 58555— orgmn m: Eamon 88: 38% no? cone:— 32v. Ewan" 2. Yield (kg) Yield (kg) Yield (kg) mooooo wooooo Aooooo wooooo Nooooo ‘ooooo o mooooo mooooo Aooooo wooooo Nooooo aooooo o mooooo mooooo Aooooo wooooo mooooo dooooo o 1 1 1 1 1 1 o J L I. 1 l 1 1 q q d a mooo .5000 .5000 moooo mmooo o I. 1 1 1 1 1 a q _ 4 . wooo aoooo .5600 N025 Mmooo 5 o q — ~ q _ 88 :88 :88 ~88 ~88 mum 3:82.: Yield (kg) Yield (kg) Yield (kg) mooooo mooooo #ooooo wooooo mooooo Aooooo o 4 a _ a _ o mooo A 0000 A mooo Noooo mmooo #00000 3.... uooooo .r:::xra..e NOOOOO . o nooooo d e o . _ . . _ o mooo ..oooo Amoco Noooo Nmooo mooooo-n mooooo 500000 wooooo Nooooo Aooooo.r u o _ _ . _ _ o mooo 40000 Amoco Moooo mmooo mmm 315033 1 J2; 1 1 1 1 mooooo mooooo J l 1 1 1 . 168 1 J 1 L Ewe—8 w~.l >m m: Emfio :w. 0x83 2:: :8 8x83 came—883138: 5 88889: 0:3 mm" BE: 8 :8 @9858 _om 5 Ems?“ 5. 38¢? m3 885: $05 58% :8 32838: o». «8% 3:: mumgmam mfioow 323% 8mm :5: N03 om Em 5:15er _oSooo:m_m8¢ 888 1 888 1 888 A 888 1 M. 888 1 888 1 W 888 1 ~88 1 m ~88 1 888 1 88 1 o _ _ 4 _ . o o 88 88 38 ~88 ~88 8 000000 1 000000 1 000000 1 8008 1 888 1 W. 888 1 000000 1 ll: W... 000000 1 888 1 l....... m ~88 1 1.00000 1 . . . o . 000000 1 0 A d . _ 0 _ 0 o 88 88 a8 88 ~88 8 000000 1 000000 1 888 1 888 0 888 1 m .88 1 888 1 W 888 1 888 1 V m 888 1 388 1 . 388 1 o _ 1 . _ 1 1 o o 88 88 88 ~88 ~88 o 000 8.525. 1V . mooo 808 $08 moooo $08 000 8.50:3 q . 176 30:8 00.186 8 30:8 00» 8838 8a 8103:? 8 $05 888 amw. 177 30:8 00. Proportion of years with Proportion of years with SSB < 20% unfished SSB SSB < 20% unfished SSB SSB < 20% unfished SSB Proportion of years with mmm. >Soooz1m_m8a 0.92 >Sooo:m_m88 ..o 1 0.92 0.3.8.18 .883 0.0 4 ._.o 1. 0.0 1 0.0 8 0.0 1 8 1 8 . pm 1 s 0.# 1 0.# 4 0.# 1 0.0 11 _ 4 . _ 0.0 4 _ _ .0 0.0 A; . . . . ob oh 8 ..m ~.o 8 9m 8 ..m ~. 0... 3 ..o ..m ~.o 0.0 1 0.0 1 0.0 L Pm 1 oh 4 0.0 1 oh 1 0.# 0 0.# 1 oo 7 o.» 1 o.» 1 on 1 0.0 (I 2 4 q . 0.0 1 l q _ 8 . o C 8 cm 3 ..m ~o oo 3 ..o ..m ~.o . -1 1 1 1 1 0.0 0.0 ..o ..m No ...0 1 1_.0 4 ._.o 1 0.0 1 0.0 0 ob 1 01m 11 pm 1. O.w .1 C15. 11 0.# 1 0.5. 11 o.~ 1 o.~ 1 o.~ 1 o... 0 0.0 1 0.0 1 1 1 ob on ..o ..m ~.o ob 9m 8 ..m ~.o ob 9m .8 ..m No <8... <8 <15... <3 <15... <3 Emfio wu.|>m m: Ewan mm. 983: won 305 <22; <35:sz m: 305. 179 Ewan mu. mmmfi >Sooo=m_m8a No 1 M e w 8 1 .m i so 1 .m .m 3 1 W I CI‘ 90 1 1 1 o N808 888 888 m No 1 .w Y 3 1 .m u 1.0 1 m .m .l O.m 11 a j v e. o ‘00». O O. 0.0 4 a A o 888 388 888 No 1 M .m Y 2. 1 .m i 2 1 m .n 3 1 a v ob 1 1 1 o 888 .888 888 51.2.3: N o 1 m_0<< >Sooo:m_m8a No 1 was 93:2 8 mmoma ._.w 1 9w 1 ...o 1 4.0 1 om 1 pm 1 e 00‘ oo o o 0“ ob 1 1 _ . 1 1 1 o 888 388 888 o 888 888 888 No 1 N; S 1 2. 1 S 1 S 1 00° 3 1 om 1 000 “oo‘OL .0 e coo-‘04 0.0 q 1 . 0.0 _ _ 1 o 888 $88 888 o 888 888 888 No 1 No 1 3m 1 ._.m 1 4.0 1 ._.o 1 9m 1 oo 1 ob 1 1 1 ob 1 1 1 o 888 .888 888 o 188° 888 888 <3... :6. <55 3: Ewan mm.|>m 5 3mg» wme 982: m3 <35 <22; 90 @8830: 03133 <39 mumgmsm floor 3925 5mm 9»: Neck. cm :5 cammroa _oc882m_m88 2 1 Ab .1 pm 1 pm 1 9m 1 Pm 1 v o;.1 o; 1 ohm- o.» 1 .. u pp. ....uinfl ob 1 oav1 1 1 o wooooo gooooo mooooo 0 9b 1 4.04 gm 1 F’ 0) 1L okw1 OAW1 oxu1 F3 15 1 on~1 F3 n) l s . oku o 888 888 888 o Po 1 F) c: .A O 4 0gw1 Avm1 Ohw1 Ohw1 o;a1 Ayn egg 888 o Proportion 0" years With Proportion of years with Proportion of years with 888 < 20% unfished SSB 333 < 20% unfished SSB 533 < 20% unfished SSB o 888 .888 so: :6. m_0<< >c88:m_m88 #o.4 ohw1 chw1 oxa1 ¢N 1 nooooo 1‘ d mooooo nooooo — 1. Aooooo 9o 1 mooooo o S 1 Avm1 chw1 o;a1 ON 1 9o a — C #00009 mooooo A.o 1 Cd 1 Ca 1 9#.1 oh~1 9b 1 o # mooooo a Aooooo <88 :6. 992 93:2 8 «883 mooooo mooooo A booooo ‘ nooooo A mooooo 182 . mooooo mooooo <88 :6. _ Aooooo . mooooo COMPREHENSIVE LIST OF REFERENCES Alverson, D.L., and W.T. Pereyra. 1969. Demersal fish exploration in the northeastern Pacific Ocean — an evaluation of exploratory fishing methods and analytical approaches to stock size and yield forecasts. Journal of the Fisheries Research Board of Canada 26: 1985-2001. Annala, J .H. 1993. Fishery assessment approaches in New Zealand’s ITQ system. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 791-805. Barange, M., M. Bemal, M.C. Cercole, L.A. Cubillos, C.L. Cunningham, G.M. Daskalov, J .A.A. De Oliveira, M. Dickey-Collas, K. Hill, L.D. Jacobson, F.W. deter, J. Masse, H. Nishida, M. Niquen, Y. Oozeki, I. Palomera, S.A. Saccardo, A. Santojanni, R. Serra, S. Somarakis, Y. Stratoudakis, C.D. van der Lingen, A. Uriarte, and A. Yatsu. In press. Current trends in the assessment and management of small pelagic fish stocks. Chapter 10 in D. Checkley, D.M. Jr, C. Roy, Y. Oozeki and J. Albeit, editors. Climate Change and Small Pelagic Fish. Cambridge University Press. Battaile, BC. and T.J. Quinn 11. 2004. Catch per unit effort standardization of the eastern Bering Sea walleye Pollock fleet. Fisheries Research 70 (2004): 161-177. Beauchamp, K.C., N.C. Collins, and BA. Henderson. 2004. Covariation of growth and maturation of lake Whitefish. Journal of Great Lakes Research 30(3): 451-460. Becker, G.C. 1983. Fishes of Wisconsin. The University of Wisconsin Press. Madison, Wisconsin. Beddington, J. R., and J. G. Cooke. 1983. The potential yield of fish stocks. FAO Fisheries Technical Paper 242. Beddington, J .R., and RM. May. 1977. Harvesting natural populations in a randomly fluctuating environment. Science 197: 463-465. Bishop, J ., W.N. Venables, and Y-G. Wang. 2004. Analyzing commercial catch and effort data from a Penaeid trawl fishery: A comparison of linear models, mixed models, and generalized estimating equations approaches. Fisheries Research 70(2004): 179-193. Booth, A.J. 2004. Determination of cichlid specific biological reference points. Fisheries Research 67: 307-316. 183 Brodziak, J. 2002. In search of optimal harvest rates for West Coast groundfish. North American Journal of Fisheries Management 22: 25 8-271. Brodziak, J ., and C .M. Legault. 2005. Model averaging to estimate rebuilding targets for overfished stocks. Canadian Journal of Fisheries and Aquatic Sciences 62: 544- 562. Brown, R.W., M. Ebener, and T. Gorenflo. 1999. Great Lakes commercial fisheries: historical overview and prognosis for the future. Pages 307-354 in W. W. Taylor and C. P. F erreri, editors. Great Lakes fishery policy and management: a binational perspective. Michigan State University Press, East Lansing. Burnham, K.P., and DR. Anderson. 2004. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, Second Edition. Springer-Verlag New York, Inc. New York. Butterworth, D.S. 2008. Some lessons from implementing management procedures. Pages 381-397 in K. Tsukamoto, T. Kawamura, T. Takeuchi, T.D. Beard, Jr., and M.J. Kaiser, editors. Fisheries for global welfare and environment, 5th world fisheries congress 2008. TERRAPUB, Tokyo. Butterworth, D.S., and M0. Bergh. 1993. The development of a management procedure for the South African anchovy resource. Pages 83-99 in S.J. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Caddy, J .F., and R. Mahon. 1995. Reference points for fisheries management. FAO Fisheries Technical Paper Number 347. Caddy, J. F ., and R. McGarvey. 1996. Targets or limits for management of fisheries? North American Journal of Fisheries Management 16: 479-487. Campana, S.E., W. Joyce, L. Marks, L.J. Natanson, N.E. Kohler, C.F. Jensen, J .J . Mello, and H.L. Pratt, Jr. 2002. Population dynamics of the porbeagle in the northwest Atlantic ocean. North American Journal of Fisheries Management 22: 106-121. Campbell, RA. 2004. CPUE standardisation and the construction of indices of stock abundance in a spatially varying fishery using general linear models. Fisheries Research 70(2004): 209-227. Clark, C. W. 1973. The economics of overexploitation. Science 181: 630-634. 184 Clark, C.W., and GP. Kirkwood. 1986. On uncertain renewable resource stocks: optimal harvest policies and the value of stock surveys. Journal of Environmental Economics and Management 13: 23 5-244. Clark, W.G. 1991. Groundfish exploitation rates based on life history parameters. Canadian Journal of Fisheries and Aquatic Sciences 48: 734-750. Clark, W.G. 1993. The effect of recruitment variability on the choice of a target level of spawning biomass per recruit. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 233-246. Clark, W.G. 1999. Effects of an erroneous natural mortality rate on a simple age- structured stock assessment. Canadian Journal of Fisheries and Aquatic Sciences 56:1721-1731. Clark, W.G. 2002. F35% revisited ten years later. North American Journal of Fisheries Management 22: 251-257. Clark, W.G., and SR. Hare. 2004. A conditional constant catch policy for managing the Pacific halibut fishery. North American Journal of Fisheries Management 24: 106-113. Cochrane, K.L., D.S. Butterworth, J .A.A. De Oliveira, and BA. Reel. 1998. Management procedures in a fishery based on highly variable stocks and with conflicting objectives: experiences in the South Afi‘ican pelagic fishery. Reviews in Fish Biology and Fisheries 8: 177-214. Collie, J .S., and H. Gislason. 2001. Biological reference points for fish stocks in a multispecies context. Canadian Journal of Fisheries and Aquatic Sciences 58: 2167-2176. Collie, J .S., and PD. Spencer. 1993. Management strategies for fish populations subject to long-term environmental variability and depensatory predation. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93- 02: 629-650. Conover, DO, and SB. Munch. 2002. Sustaining fisheries yields over evolutionary time scales. Science 297(5): 94-96. Cook, H.A., T.B. Johnson, B. Locke, and B.J. Morrison. 2005. Status of lake Whitefish in Lake Erie. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 87-104. 185 Cooper, A.B., A.A. Rosenberg, G. Stefansson, and M. Mangel. 2004. Examining the importance of consistency in multi-vessel trawl survey design based on the US. west coast groundfish bottom trawl survey. Fisheries Research 70(2004): 239- 250. Deriso, RB. 1982. Relationship of fishing mortality to natural mortality and growth at the level of maximum sustainable yield. Canadian Journal of Fisheries and Aquatic Sciences 39: 1054-1058. Deriso, RB. 1987. Optimal F 0,. criteria and their relationship to maximum sustainable yield. Canadian Journal of Fisheries and Aquatic Sciences 44(Supplement 2): 339-348. Deroba, J .J ., and J .R. Bence. In press. Assessing model-based indices of lake trout abundance in 1836 Treaty waters of Lakes Huron, Michigan, and Superior. Michigan Department of Natural Resources: Fisheries Research Report, Ann Arbor, MI. Deroba, J .J ., and J .R. Bence. 2008. A review of harvest policies: understanding relative performance of control rules. Fisheries Research 94: 210-223. Dichmont C.M., A. Deng, A.E. Punt, W. Venables, and M. Haddon. 2006a. Management strategies of short-lived species: The case of Australia’s Northern Prawn Fishery 1. Accounting for multiple species, spatial structure and implementation uncertainty when evaluation risk. Fisheries Research 82: 204- 220. Dichmont C.M., A. Deng, A.E. Punt, W. Venables, and M. Haddon. 2006b. Management strategies of short-lived species: The case of Australia’s Northern Prawn Fishery 2. Choosing appropriate management strategies using input controls. Fisheries Research 82: 221-234. Dichmont C.M., A. Deng, A.E. Punt, W. Venables, and M. Haddon. 2006c. Management strategies of short-lived species: The case of Australia’s Northern Prawn Fishery 3. Factors affecting management and estimation performance. Fisheries Research 82: 235-245. DiNardo, G.T., and J .A. Wetherall. 1999. Accounting for uncertainty in the development of harvest strategies for the Northwestern Hawaiian Islands lobster trap fishery. ICES Journal of Marine Science 56: 943-951. Dorn, M.W. 2002. Advice on west coast rockfish harvest rates form Bayesian meta- analysis of stock-recruit relationships. North American Journal of Fisheries Management 22: 280-300. 186 Ebener, M .P. 1997. Recovery of lake Whitefish populations in the Great Lakes: a story of successful management and just plain luck. Fisheries 22: 18-20. Ebener, MP, and D.M. Reid. 2005. Historical context. In The state of Lake Huron 1999. Edited by MP. Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 9-18. Ebener, M.F., J .R. Bence, K. Newman, and P. Schneeberger. 2005. Application of statistical catch-at-age models to assess lake Whitefish stocks in the 1836 treaty- ceded waters of the upper Great Lakes. In Proceedings of a workshop on the dynamics of lake Whitefish and the amphipod Diporeia spp. in the Great Lakes. Edited by LC. Mohr and T.F. Nalepa. Great Lakes Fishery Commission Technical Report 66. pp. 271-309. Eggers, D.M. 1993. Robust harvest policies for Pacific salmon fisheries. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93- 02: 85-106. Enberg, K. 2005. Benefits of threshold strategies and age-selective harvesting in a fluctuating fish stock of Norwegian spring spawning herring. Marine Ecology Progress Series 298: 277-286. Engen, S., R. Lande, and B-E. Saether. 1997. Harvesting strategies for fluctuating populations based on uncertain population estimates. Journal of Theoretical Biology 186: 201-212. F ieberg, J. 2004. Role of parameter uncertainty in assessing harvest strategies. North American Journal of Fisheries Management 24: 459-474. Francis, RC. 1974. Relationship of fishing mortality to natural mortality at the level of maximum sustainable yield under the logistic stock production model. Journal of the Fisheries Research Board of Canada 31: 1539-1542. Francis, RC. 1993. Monte Carlo evaluation of risks for biological reference points used in New Zealand fishery assessments. Pages 221-230 in S.J. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Frederick, S.W., and RM. Peterman. 1995. Choosing fisheries harvest policies: when does uncertainty matter? Canadian Journal of Fisheries and Aquatic Sciences 52: 291-306. 187 Fujioka, J .T., J. Heifetz, and M.F. Sigler. 1997. Choosing a harvest strategy for sablefish based on uncertain life-history parameters. Pages 247-251 in NOAA Technical Report NMFS 130 Biology and Management of Sablefish; Papers from the International Symposium on the Biology and Management of Sablefish, Seattle. Gabriel, W.L., M.P. Sissenwine, and W.J. Overholtz. 1989. Analysis of spawning stock biomass per recruit: an example for Georges Bank haddock. North American Journal of Fisheries Management 9: 383-391. Gatto, M., and S. Rinaldi. 1976. Mean value and variability of fish catches in fluctuating environments. Journal of the Fisheries Research Board of Canada 33: 189-193. Gelman, A., and J. Hill. 2007. Data Analysis Using Regression and Multilevel Hierarchical Models. Cambridge University Press. New York, New York. Getz, W. M., R. C. Francis, and G. L. Swartzman. ,1987. On managing variable marine fisheries. Canadian Journal of Fisheries and Marine Sciences 44: 1370-1375. Getz, W.M., and R.G. Haight. 1989. Population Harvesting: Demographic Models of Fish, Forest, and Animal Resources. Princeton University Press, Princeton, New Jersey. Gibson, A.J.F., and RA. Myers. 2004. Estimating reference fishing mortality rates from noisy spawner-recruit data. Canadian Journal of Fisheries and Aquatic Sciences 61:1771-1783. Goodyear, GP. 1993. Spawning stock biomass per recruit in fisheries management: foundation and current use. Pages 67-81 in S]. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Hall, D. L., R. Hilborn, M. Stocker, and C. J. Walters. 1988. Alternative harvest strategies for Pacific herring. Canadian Journal of Fisheries and Aquatic Sciences 45: 888-897. Harley, S.J., R.A. Myers, and A. Dunn. 2001. Is catch-per-unit-effort proportional to abundance? Canadian Journal of Fisheries and Aquatic Sciences 58: 1760-1772. Hayes, DB. 2000. A biological reference point based on the Leslie matrix. Fisheries Bulletin 98: 75-85. Healey, MC. 1975. Dynamics of exploited Whitefish populations and their management with special reference to the Northwest Territories. Journal of the Fisheries Research Board of Canada. 32: 427-448. 188 Heino, M. 1998. Management of evolving fish stocks. Canadian Journal of Fisheries and Aquatic Sciences 55: 1971-1982. Helser, T.E., and J .K.T. Brodziak. 1998. Impacts of density-dependent growth and maturation on assessment advice to rebuild depleted U.S. silver hake stocks. Canadian Journal of Fisheries and Aquatic Sciences 55: 882-892. Helser, T.E., A.E. Punt, and RD. Methot. 2004. A generalized linear mixed model analysis of a multi-vessel fishery resource survey. Fisheries Research 70(2004): 251-264. Henderson, B.A., J .J . Collins, and J .A. Reckahn. 1983. Dynamics of an exploited population of lake Whitefish in Lake Huron. Canadian Journal of Fisheries and Aquatic Sciences 40: 1556-1567. Hightower, J. E., and G. D. Grossman. 1987. Optimal policies for rehabilitation of overexploited fish stocks using a deterministic model. Canadian Journal of Fisheries and Aquatic Sciences 44: 803-810. Hightower, J .E., and W.H. Lenarz. 1989. Optimal harvesting policies for the widow rockfish fishery. American Fisheries Society Symposium 6: 83-91. Hilbom, R., and C.J. Walters. 1992. Quantitative Fisheries Stock Assessment: Choice, Dynamics, and Uncertainty. Chapman and Hall, New York. Hilbom, R., A. Parma, and M. Maunder. 2002. Exploitation rate reference points for west coast rockfish: are they robust and are there better alternatives? North American Journal of Fisheries Management 22: 365-375. Hj erne, O., and S. Hansson. 2001. Constant catch or constant harvest rate? The Baltic Sea cod fishery as a modelling example. Fisheries Research 53: 57-70. Hoyle, J .A., T. Schaner, J .M. Casselman, and R. Dermott. 1999. Changes in lake Whitefish stocks in eastern Lake Ontario following Dreissena mussel invasion. Great Lakes Research Review 4: 5-10. Irwin, B.J., M.J. Wilberg, J .R. Bence, and ML. Jones. 2008. Evaluating Alternative Harvest Policies for Yellow Perch in Lake Michigan. Fisheries Research 94:267- 281. Ishimura, G., A.E. Punt, and DD. Huppert. 2005. Management of fluctuating fish stocks: the case of Pacific whiting. Fisheries Research 73: 201-216. 189 Jacobson, PC, and W.W. Taylor. 1985. Simulation of harvest strategies for a fluctuating population of lake Whitefish. North American Journal of Fisheries Management 5: 537-546. J akobsen, T. 1992. Biological reference points for northeast Arctic cod and haddock. ICES Journal of Marine Science 49: 155-166. J akobsen, T. 1993. The behavior of F 10.1,, Fmed, and F high in response to variation in parameters used for their estimation. Pages 119-125 in S.J. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Jensen, AL. 1976. Assessment of the United States lake Whitefish fisheries of Lake Superior, Lake Michigan, and Lake Huron. Journal of the Fisheries Research Board of Canada 33: 747-759. Katsukawa, T. 2004. Numerical investigation of the optimal control rule for decision- making in fisheries management. Fisheries Science 70: 123-131. Kell, L.T., G.M. Pilling, G.P. Kirkwood, M.A. Pastoors, B. Mesnil, K. Korsbrekke, P. Abaunza, R. Aps, A. Biseau, P. Kunzlik, C.L. Needle, B.A. Roel, and C. Ulrich. 2006. An evaluation of multi-annual management strategies for ICES roundfish stocks. ICES Journal of Marine Science 63: 12-24. Koelz, W. 1926. Fishing industry of the Great Lakes. Pages 554-617, In Report of the US. Commissioner of Fisheries for 1925. Koonce, J .F., and B.J. Shuter. 1987. Influence of various sources of error and community interactions on quota management of fish stocks. Canadian Journal of Fisheries and Aquatic Sciences 44(Supplement 2): 61-67. Kratzer, J .F ., W.W. Taylor, C.P. Ferreri, and MP. Ebener. 2005. Factors affecting growth of lake Whitefish in the upper Laurentian Great Lakes. Advances in Limnology 60: 459-470. Lande, R., B-E. Saether, and S. Engen. 1997. Threshold harvesting for sustainability of fluctuating resources. Ecology 78(5): 1341-1350. Lande, R., S. Engen, and B-E. Saether. 1995. Optimal harvesting of fluctuating populations with a risk of extinction. The American Naturalist 145: 728-745. Larkin, P. A. 1977. An epitaph for the concept of maximum sustainable yield. Transactions of the American Fisheries Society 106: 1-11. 190 Larkin, P.A., and WE. Ricker. 1964. Further information on sustained yields from fluctuating environments. Journal of the Fisheries Research Board of Canada 21(1):1-7. Leaman, BM. 1993. Reference points for fisheries management: the western Canadian experience. Pages 15-30 in S.J. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Lillegard, M., S. Engen, B-E Saether, and R. Toresen. 2005. Harvesting strategies for Norwegian spring-spawning herring. Oikos 110: 567-577. Lowe, SA. and G.G. Thompson. 1993. Accounting for uncertainty in the development of exploitation strategies for the atka mackerel resource of the Aleutian Islands. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 203-231. Lumb, C.E., T.B. Johnson, H.A. Cook, and J .A. Hoyle. 2007. Comparison of lake Whitefish growth, condition, and energy density between Lakes Erie and Ontario. Journal of Great Lakes Research 33: 314-325. Mace, RM. 2001. A new role for MSY in single-species and ecosystem approaches to fisheries stock assessment and management. Fish and Fisheries 2: 2-32. Mace, P.M., and MP. Sissenwine. 1993. How much spawning per recruit is enough? Pages 101-118 in 8.1. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Maguire, J .J ., and PM. Mace. 1993. Biological reference points for Canadian Atlantic gadoid stocks. Pages 321-331 in SI Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Matsuda, H., and PA. Abrams. 2006. Maximal yields from multispecies fisheries systems: rules for systems with multiple trophic levels. Ecological Applications 16: 225-237. Maunder, MN. 2001. A general framework for integrating the standardization of catch per unit of effort into stock assessment models. Canadian Journal of Fisheries and Aquatic Sciences 58: 795-803. 191 Maunder, MN. and AD. Langley. 2004. Integrating the standardization of catch-per- unit-of-effort into stock assessment models: testing a population dynamics model and using multiple data types. Fisheries Research 70: 389-395. Maunder, MN. and A.E. Punt. 2004. Standardizing catch and effort data: a review of recent approaches. Fisheries Research 70(2004): 141-159. Maunder, MN. and P.J. Starr. 2003. Fitting fisheries models to standardized CPUE abundance indices. Fisheries Research 63(2003): 43-50. Maynou, F., M. Demestre, and P. Sanchez. 2003. Analysis of catch per unit effort by multivariate analysis and generalised linear models for deep-water crustacean fisheries off Barcelona. Fisheries Research 65(2003): 257-269. McCulloch, CE. and SR. Searle. 2001. Generalized, Linear, and Mixed Models. New York: John Wiley and Sons, Inc. McGlade, J .M. 1989. Integrated fisheries management models: understanding the limits to marine resource exploitation. American Fisheries Society Symposium 6: 139- 165. Milner-Gulland, E.J., K. Shea, H. Possingham, T. Coulson, and C. Wilcox. 2001. Competing harvesting strategies in a simulated population under uncertainty. Animal Conservation 4: 157-167. Mohr, L.C., and MP. Ebener. 2005a. The coregonine community. In The state of Lake Huron 1999. Edited by MP Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 69-76. Mohr, L.C., and MP Ebener. 2005b. Description of the fisheries. In The state of Lake Huron 1999. Edited by MP Ebener. Great Lakes Fishery Commission Special Publication 05-02, pages 19-26. Mohr, L.C., and Nalepa, T.F. (Editors). 2005. Proceedings of a workshop on the dynamics of lake Whitefish (Coregonus clupeaformis) and the amphipod Diporeia spp. in the Great Lakes. Great Lakes Fishery Commission Technical Repport 66. Murawski, S.A., and J .S. Idoine. 1989. Yield sustainability under constant-catch policy and stochastic recruitment. Transactions of the American Fisheries Society 118: 349-367. Myers, R.A., K.G. Bowen, and NJ. Barrowman. 1999. Maximum reproductive rate of fish at low population sizes. Canadian Journal of Fisheries and Aquatic Sciences 56: 2404-2419. 192 Myers, R.A., A.A. Rosenberg, P.M. Mace, N. Barrowman, and V.R. Restrepo. 1994. In search of thresholds for recruitment overfishing. ICES Journal of Marine Science 51: 191-205. Myers, R.A., and B. Worm. 2005. Extinction, survival or recovery of large predatory fishes. Philosophical Transactions of the Royal Society B 360(2005): 13-20. Myers, R.A., G. Mertz, and J. Bridson. 1997. Spatial scales of interannual recruitment variations of marine, anadromous, and freshwater fish. Can J Fish Aquat Sci 54: 1400-1407. Nieland, J .L., M.J. Hansen, M.J. Seider, and J .J . Deroba. 2008. Modeling the sustainability of lake trout fisheries in eastern Wisconsin waters of Lake Superior. Fisheries Research 94: 304-314. Ngo, L., and R. Brand. 1997. Model selection in linear mixed effects models using SAS proc mixed. SAS Institute Inc., Proceedings of the 22nd Annual SAS Users Group International Conference: 1 3 35-1 340. Nostbakken, L. 2006. Regime Switching in a fishery with stochastic stock and price. Journal of Environmental Economics and Management 51: 231-241. NRC (National Research Council). 1994. Improving the Management of US. Marine Fisheries. National Academy Press, Washington, DC. Overholtz, W.J. 1999. Precision and uses of biological reference points calculated from stock recruitment data. North American Journal of Fisheries Management 19: 643-657. Overholtz, W.J., S.F. Edwards, and J .K.T. Brodziak. 1993. Strategies for rebuilding and harvesting New England groundfish resources. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 507-527. Pacific Fishery Management Council. 1998. Options and analyses for the coastal pelagic species fishery management plan: appendix B to amendment 8. 134 pages. http://www.pcouncil.org/cps/cpsfmp/a8apdxb.pdf Parma, A. 1993. Retrospective catch-at-age analysis of Pacific halibut: implications on assessment of harvesting policies. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 247-265. 193 Peterman, R.M., and Anderson, J.L. 1999. Decision analysis: a method for taking uncertainties into account in risk-based decision making. Human and Ecological Risk Assessment 5: 231-244. Polacheck, T., N.L. Klaer, C. Millar, and AL. Preece. 1999. An initial evaluation of management strategies for the southern bluefin tuna fishery. ICES Journal of Marine Science 56: 811-826. Pothoven, S.A., T.F. Nalepa, P.J. Schneeberger, and SB. Brandt. 2001. Changes in diet and body condition of lake Whitefish in southern Lake Michigan associated with changes in benthos. North American Journal of Fisheries Management: 21: 876- 883. Potter, E.C.E., J .C. MacLean, R.J. Wyatt, and R.N.B. Campbell. 2003. Managing the exploitation of migratory salmonids. Fisheries Research 62: 127-142. Punt, A.E. 1997. The performance of VPA based management. Fisheries Research 29: 217-243. Punt, A.E. 2003. Evaluating the efficacy of managing west coast groundfish resources through simulations. Fisheries Bulletin 101: 860-873. Punt, A.E., M.W. Dom, and M.A. Haltuch. 2008. Evaluation of threshold management strategies for groundfish off the US. West Coast. Fisheries Research 94:251-266. Punt, A.E., A.D.M. Smith, and G. Cui. 2002a. Evaluation of management tools for Australia’s South East Fishery 2. How well can management quantities be estimated? Marine and Freshwater Research 53: 631-644. Punt, A.E., A.D.M. Smith, and G. Cui. 2002b. Evaluation of management tools for Australia’s South East Fishery 3. Towards selecting appropriate harvest strategies. Marine and Freshwater Research 53: 645-660. Punt, A.E., A.D.M. Smith, and G. Cui. 2002c. Evaluation of management tools for Australia’s South East Fishery 1. Modelling the South East Fishery taking account of technical interactions. Marine and Freshwater Research 53: 615-629. Punt, A.E., A.J. Penney, and R.W. Leslie. 1996. Abundance indices and stock assessment of south Atlantic albacore. Collective Volume of Scientific Papers of the International Commission for the Conservation of Atlantic Tunas 43: 225-245. Quiggin, J. 1992. How to set catch quotas: a note on the superiority of constant effort rules. Journal of Environmental Economics and Management 22: 199-203. 194 Quinn, T.J., II, and J .S. Collie. 2005. Sustainability in single-species population models. Philosophical Transactions of the Royal Society B 360: 147-162. Quinn, T.J., II, R. Fagen, and J. Zheng. 1990. Threshold management policies for exploited populations. Canadian Journal of Fisheries and Aquatic Sciences 47: 2016-2029. Quinn, T. J. 11, and N. J. Szarzi. 1993. Determination of sustained yield in Alaska's recreational fisheries. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 61-84. Quinn, T.J., II, and RB. Deriso. 1999. Quantitative Fish Dynamics. Oxford University Press Inc. New York, New York. Rahikainen, M., and R.L. Stephenson. 2004. Consequences of growth variation in northern Baltic herring for assessment and management. ICES Journal of Marine Science 61: 338-350. Ralston, S., J .R. Bence, W.G. Clark, R.J. Conser, T. Jagielo, and T.J. Quinn 11. 2000. West Coast groundfish harvest rate policy workshop. Panel Report, Seattle, Washington. Reed, W.J. 1979. Optimal escapement levels in stochastic and deterministic harvesting models. Journal of Environmental Economics and Management 6: 350-363. Ricker, W.E. 1958. Maximum sustained yields from fluctuating environments and mixed stocks. Journal of the Fisheries Research Board of Canada 15(5): 991- 1006. Ricker, W.E. 1975. Computation and interpretation of biological statistics of fish populations. Bulletin of the Fisheries Research Board of Canada 191. Rose, GA. and D.W. Kulka. 1999. Hyperaggregation of fish and fisheries: how catch- per-unit-effort increased as the northern cod declined. Canadian Journal of Fisheries and Aquatic Sciences 56(supplement 1): 118-127. Ruppert, D., R. L. Reish, R. B. Deriso, and R. J. Carroll. 1985. A stochastic population model for managing the Atlantic menhaden fishery and assessing managerial risks. Canadian Journal of Fisheries and Aquatic Sciences 42: 1371-1379. SAS. 2003. SAS version 9.1 help and documentation. Cary, North Carolina: SAS Institute, Inc. 195 Sethi, G., C. Costello, A. Fisher, M. Hanemann, and L. Karp. 2005. Fishery management under multiple uncertainty. Journal of Environmental Economics and Management 50: 300-318. Siddeek, M.S.M., and A.H.S. Al-Hosni. 1998. Biological reference points for managing kingfish in Oman waters. Naga: the ICLARM Quarterly 32-36. Sigler, M.F., and J .T. F ujioka. 1993. A comparison of policies for harvesting sablefish in the Gulf of Alaska. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93-02: 7-19. Sissenwine, MP. 1978. Is MSY an adequate foundation for optimum yield? Fisheries 3: 22-24, 37-42. Sissenwine, MP, and J .G. Shepherd. 1987. An alternative perspective on recruitment overfishing and biological reference points. Canadian Journal of Fisheries and Aquatic Sciences 44: 91 3-91 8. Sladek Nowlis, J ., and B. Bollermann. 2002. Methods for increasing the likelihood of restoring and maintaining productive fisheries. Bulletin of Marine Science 70: 715-731. Smiley, CW. 1882. Changes in the fisheries of the Great Lakes during the decade, 1870-1880. Transactions of the American F ish-Cultural Association 11: 28-37. Smith, A.D.M., K.J. Sainsbury, and RA. Stevens. 1999. Implementing effective fisheries management systems — management strategy evaluation and the Australian partnership approach. ICES Journal of Marine Science 56: 967-979. Spencer, P. D. 1997. Optimal harvesting of fish populations with nonlinear rates of predation and autocorrelated environmental variability. Canadian Journal of Fisheries and Aquatic Science 54: 59-74. Steinshamn, SI. 1993. Management strategies: fixed or variable catch quotas. Pages 373-385 in SI. Smith, J .J. Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Steinshamn, 8.1. 1998. Implications of harvesting strategies on population and profitability in fisheries- Marine Resource Economics 13: 23-36. Swain. D.P., A.F. Sinclair, and J .M. Hanson. Evolutionary response to size-selective mortality in an exploited fish population. Proceedings of the Royal Society B 274: 1015-1022. 196 Tautz, A., P.A. Larkin, and WE. Ricker. 1969. Some effects of simulated long-term environmental fluctuations on maximum sustained yield. Journal of the Fisheries Research Board of Canada 26: 2715-2726. Thompson, G.G. 1993. A proposal for a threshold stock size and maximum fishing mortality rate. Pages 303-320 in S]. Smith, J .J . Hunt, and D. Rivard, editors. Risk Evaluation and Biological Reference Points for Fisheries Management. Canadian Special Publication of Fisheries and Aquatic Sciences 120. Vasconcellos, M. 2003. An analysis of harvest strategies and information needs in the purse seine fishery for the Brazilian sardine. Fisheries Research 59: 363-378. Venables, W.N., and CM. Dichmont. 2004. GLMS, GAMS, and GLMMs: an overview of theory for applications in fisheries research. Fisheries Research 70(2004): 319- 337. Walters, C.J. 1986. Adaptive management of renewable resources. MacMillan, New York, New York, USA. Walters, C .J ., and A.M. Parma. 1996. Fixed exploitation rate strategies for coping with effects of climate change. Canadian Journal of Fisheries and Aquatic Sciences 53: 148-158. Walters, C.J., and PH. Pearse. 1996. Stock information requirements for quota management systems in commercial fisheries. Reviews in Fish Biology and Fisheries 6: 21-42. Walters, C.J., V. Christensen, 8.]. Martel], and J .F. Kitchell. 2005. Possible ecosystem impacts of applying MSY policies from single-species assessment. ICES Journal of Marine Science 62: 558-568. Wang, H-Y, T.O. Hédk, M.P. Ebener, L.C. Mohr, and P.J. Schneeberger. 2008. Spatial and temporal variation of maturation schedules of lake Whitefish in the Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences 65: 2157-2169. Wilberg, M.J., and J .R. Bence. 2006. Performance of time-varying catchability estimators in statistical catch-at-age analysis. Canadian Journal of Fisheries and Aquatic Sciences 63: 2275-2285. Wilberg, M.J., B.J. Irwin, M.L. Jones, and J .R. Bence. 2008. Effects of source-sink dynamics on harvest policy performance for Yellow Perch in southern Lake Michigan. Fisheries Research 94: 282-289. 197 Williams, EH. 2002. The effects of unaccounted discards and misspecified natural mortality on harvest policies based on estimates of spawners per recruit. North American Journal of Fisheries Management 22: 311-325. Zheng, J., F.C. Funk, G.H. Kruse, and R. Fagen. 1993a. Evaluation of threshold management strategies for Pacific herring in Alaska. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93- 02: 141-165. Zheng, J ., T.J. Quinn 11, G.H. Kruse. 1993b. Comparison and evaluation of threshold estimation methods for exploited fish populations. Proceedings of the International Symposium on Management Strategies for Exploited Fish Populations, University of Alaska Sea Grant College Program Report Number 93- 02: 267-289. 198 5 l l TY LIBRARIES H 255 601 1293 03 ll 3 l m R E w N U E T A T S N A m H m M . . . . o. . s o u . . . .. . o . . . . . . . . .. . . . u . . . . . . . - . .1 .. .. . . .... ..... 1 . 1. .1 . 1 . .. . n . . . 1 . _ 1 . .1. 1 . . _ .... . . . ... .. . . .1 v u . 4 ' a v e V n.1 . I . . . . .1 ~ . . . a .1 . . . . . . . . . 1 . 1 1 . 1 .. . . . .1 . . . .. . . . o . 1 n r r e . . l e . I v. .u 1 o 1 . . r. . . . 1 . 1. . .. u. a ul ~ 1 r . n. v . . . . . .. . . 1 . 1 . 1 . . . 1 . . . . . . . 1 1... . .. . . . 1 . 1 1 .1 . . . . n c .. 1 . . . . .. . . . U . ~ 0 1 1 . 1 . .. .. . . . . .. . . . .. . . u . . . 1 . . . . 1 . .. . .. . . . . .. .. 1. .. . 1. . u a u e v . r .. . u 1 . . 1 . ... .. . 1 . . . . .. . . .. . 1 . . . . s . . . . . .. . . . n . u a n . u ‘ . v s . u . . . . 1 1 . J . 1 . . .. . . . . . . . . . . . . 1 . . . . . . . . . . .. . . ‘ ~ - n .n I. n . a. .0 u I . .1 . 1 . . 1 . . . u u t . . . 1 . 1 . . . .. n J. . . . 1 . .v . .... .. .. 1 . . . . .1... .. . . . . . . . .. . 1 . .1 . 1 .. 1 .. . . . 1 . . . . . . .. . . o . .. . . . . . . . . . . . . . . 1 . .1 . 1.. ... . . . . . s . .. . . . . . . .. . 1.1.1 . . . . . .. . . . . 1 . . . . 1 . . . . . . .. . . .. . . . . . . . . .u .n . ... o . 1- . . .. a la n n . . I v I . . . . . . . . . . . . . . . u I u . . . . . r .1 . . .. . . .. .1 . . . . . . .. .. . . . . . .1 . .. .. . . .. . n . . .. 1 . ... . . s . . . ... .. . . . . 1 . . . .1... . . .1 1 . . . . ... .. 1 .. . 1 n I . a o e. 1- . I I. l . . o. . v . . .. . .. .... . . . . .. . . . . 1 .1 . . . . . 1... . . . . . 1 . _ 1 . ... .. . . .1. 1 . . . . . . . . . . .. . . 1 .. . .. .1 . . . . . . . - . . 1 ... . 1 . . 1 . .1 . . .1 . . .. . .. . . . . .. 1. . . . . . n .. 1‘. . . . .1 4 ...... 1 ‘ .... . . . ... .. . 1. . . .. . . 1. . 1 . .... . . . . 1. 1. . r . .. 1 1 1. 1... . . . . . .. . 1. 1 . . 1. .. . 1 . . . . . ... . .. . ... 1 . ... . . . . ... . . ... . . . . . ... . . .. .1 . . . . . . . . .1 . a . ~ u e u e . . . . . . 1. . . . . . . 1 . .. . . . .. .1 .... 1 . . . . .... . .1 1 u. . c. .1 . . .. .. . . . .1... . .. . 1 . . . . 1 . . . . .. .. .. .. . . . .. . . . . . ... . . ~ . ... . .. .. . . . . . vu. .. ...1 . .o .1. .. .. . . . . . 1. . . 1 . . ......o .u 4 c l.o.I-.ro\o~1. ... .