This article was downloaded by: [Department Of Fisheries] On: 20 July 2015, At: 19:14 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: 5 Howick Place, London, SW1P 1WG
North American Journal of Fisheries Management Publication details, including instructions for authors and subscription information: http://www.tandfonline.com/loi/ujfm20 Raft and Floating Radio Frequency Identification (RFID) Antenna Systems for Detecting and Estimating Abundance of PIT-tagged Fish in Rivers Eric R. Fethermana, Brian W. Avilab & Dana L. Winkelmanc a Colorado Parks and Wildlife, 317 West Prospect Road, Fort Collins, Colorado 80526, USA b Colorado State University, Colorado Cooperative Fish and Wildlife Research Unit, Department of Fish, Wildlife and Conservation Biology, Colorado State University, Fort Collins, 80523, Colorado, USA c U.S. Geological Survey, Colorado Cooperative Fish and Wildlife Research Unit, Department of Fish, Wildlife and Conservation Biology, Colorado State University, Fort Collins, 80523, Click for updates Colorado, USA Published online: 10 Oct 2014.
To cite this article: Eric R. Fetherman, Brian W. Avila & Dana L. Winkelman (2014) Raft and Floating Radio Frequency Identification (RFID) Antenna Systems for Detecting and Estimating Abundance of PIT-tagged Fish in Rivers, North American Journal of Fisheries Management, 34:6, 1065-1077, DOI: 10.1080/02755947.2014.943859 To link to this article: http://dx.doi.org/10.1080/02755947.2014.943859
PLEASE SCROLL DOWN FOR ARTICLE
Taylor & Francis makes every effort to ensure the accuracy of all the information (the “Content”) contained in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http:// www.tandfonline.com/page/terms-and-conditions North American Journal of Fisheries Management 34:1065–1077, 2014 C American Fisheries Society 2014 ISSN: 0275-5947 print / 1548-8675 online DOI: 10.1080/02755947.2014.943859
ARTICLE
Raft and Floating Radio Frequency Identification (RFID) Antenna Systems for Detecting and Estimating Abundance of PIT-tagged Fish in Rivers
Eric R. Fetherman* Colorado Parks and Wildlife, 317 West Prospect Road, Fort Collins, Colorado 80526, USA Brian W. Avila Colorado State University, Colorado Cooperative Fish and Wildlife Research Unit, Department of Fish, Wildlife and Conservation Biology, Colorado State University, Fort Collins, Colorado 80523, USA Dana L. Winkelman U.S. Geological Survey, Colorado Cooperative Fish and Wildlife Research Unit, Department of Fish, Wildlife and Conservation Biology, Colorado State University, Fort Collins, Colorado 80523, USA
Abstract Portable radio frequency identification (RFID) PIT tag antenna systems are increasingly being used in studies examining aquatic animal movement, survival, and habitat use, and their design flexibility permits application in a wide variety of settings. We describe the construction, use, and performance of two portable floating RFID PIT tag antenna systems designed to detect fish that were unavailable for recapture using stationary antennas or electrofishing. A raft antenna system was designed to detect and locate PIT-tagged fish in relatively long (i.e., ≥10 km) river reaches, and consisted of two antennas: (1) a horizontal antenna (4 × 1.2 m) installed on the bottom of the raft and used to detect fish in shallower river reaches (<1 m), and (2) a vertical antenna (2.7 × 1.2 m) for detecting fish in deeper pools (≥1 m). Detection distances of the horizontal antenna were between 0.7 and 1.0 m, and detection probability was 0.32 ± 0.02 (mean ± SE) in a field test using rocks marked with 32-mm PIT tags. Detection probability of PIT-tagged fish in the Cache la Poudre River, Colorado, using the raft antenna system, which covered 21% of the wetted area, was 0.14 ± 0.14. A shore-deployed floating antenna (14.6 × 0.6 m), which covered 100% of the wetted area, was designed for use by two operators for detecting and locating PIT-tagged fish in shorter (i.e., <2 km) river reaches. Detection distances of the shore-deployed floating antenna were between 0.7 and 0.8 m, and detection probabilities during field deployment in the St. Vrain River exceeded 0.52. The shore-deployed floating antenna was also used to estimate abundance of PIT-tagged fish. Results suggest that the shore-deployed floating antenna could be used as an
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 alternative to estimating abundance using traditional sampling methods such as electrofishing.
Passive integrated transponder tag technology has many ad- behavior studies, especially those examining habitat selection vantages over traditional marking techniques. They allow indi- or movement (Nunnallee et al. 1998; Zydlewski et al. 2006; vidual identification, have an infinite life as long as the tag is not Bond et al. 2007; Compton et al. 2008; Connolly et al. 2008; damaged, are easily applied and well retained, and have minimal Aymes and Rives 2009). effects on growth and survival of fishes (Gries and Letcher 2002; Stationary antennas are typically used to detect PIT-tagged Zydlewski et al. 2006; Ficke et al. 2012). Traditionally, the utility fish, but the use of portable antennas is becoming more com- of PIT tagging has been limited to physical recapture events mon. Portable antennas have been used in studies examining using methods such as electrofishing (Zydlewski et al. 2006). aquatic animal movement, survival, and habitat use, and their Stationary antennas allow passive recaptures of PIT-tagged design flexibility permits application in a wide variety of set- fish and have recently been used to detect PIT-tagged fish in tings. Initial technological advances in portable PIT tag antenna
*Corresponding author: [email protected] Received December 19, 2013; accepted July 1, 2014 1065 1066 FETHERMAN ET AL.
systems enabled effective detection in shallow, wadable rivers reach of the Cache la Poudre River, Colorado. The array was of salmonid fishes such as Atlantic Salmon Salmo salar (Rous- constructed using a 4.9-m, self-bailing, inflatable river raft. The sel et al. 2000; Zydlewski et al. 2001), Brown Trout S. trutta first antenna was installed in the bottom of the raft (horizontal) (Cucherousset et al. 2005), and steelhead Oncorhynchus mykiss and was designed for continuous deployment to detect fish in (anadromous Rainbow Trout) (Hill et al. 2006). Additionally, shallower (<1 m) sections of the river. The second antenna was a portable antennas have been developed to detect and locate dropper antenna (vertical) designed for intermittent deployment age-0 Northern Pike Esox lucius (Cucherousset et al. 2007) to detect fish in deeper (≥1 m) pools. and various age-classes of European Eel Anguilla anguilla and Both antennas consisted of two, loosely bound, continu- Common Dace Leuciscus leuciscus (Cucherousset et al. 2010). ous loops of 12-gauge thermoplastic, high-heat-resistant, nylon- Previous designs of portable antenna systems have limited their coated (THHN) wire. Binding of the wire occurred at specific effectiveness to shallow, wadable floodplains or small streams; attachment points with the raft (horizontal) or support beams however, a boat-mounted antenna system was recently devel- used to maintain antenna shape (vertical). Between these points, oped for monitoring mussel populations in larger, nonwadable the wire formed loose gaps of varying distances. The horizon- rivers (Fischer et al. 2012). tal antenna was an elliptical antenna 4 × 1.2 m in size and Portable antenna systems are limited by factors affecting de- located in the self-bailing channel of the raft (Figure 1). An- tection efficiency including tag size, power source, tag orienta- tenna shape and loop proximity were maintained by threading tion, antenna proximity if multiple antennas are used (Zydlewski the wire through sections of flexible plastic tubing secured to et al. 2006), tag collision (Axel et al. 2005; O’Donnell et al. the self-bailing holes in the floor of the raft with soft nylon 2010), and disruption of the magnetic field by the presence of cord. The vertical antenna was a 2.7- × 1.2-m rectangle, main- metal (Greenberg and Giller 2000; Bond et al. 2007). For exam- tained by four, 19-mm PVC crossbeams secured to the antenna ple, tag orientation relative to the antenna field affects detection wire with expandable spray-foam insulation (Figure 1). Holes and is higher when the tag is oriented perpendicular rather than were drilled in each crossbeam to allow water entry and 51-mm parallel to the antenna (Nunnallee et al. 1998; Morhardt et al. PVC caps filled with cement were attached to the lowermost 2000; Zydlewski et al. 2006; Compton et al. 2008; Aymes and crossbeam for ballast. Foam pipe insulation placed on the first Rives 2009). Disruption due to antenna proximity can be re- crossbeam allowed the antenna to hang vertically in the water duced through the use of multiplexers (Aymes and Rives 2009), column while remaining at the water surface without becom- and disruption from metal can be reduced by utilizing nonin- ing completely submerged. Connectors produced for welding ductive materials such as epoxy coil encasements or nylon nuts applications were used to allow for a quick disconnection of and bolts (Fischer et al. 2012). Potential limitations should be the vertical antenna in the event that the antenna got entangled accounted for during the design process, and we elaborate on and caught on submerged rocks or vegetation while deployed. the ways in which these limitations were accounted for in our The vertical antenna design facilitated easy deployment at the portable antenna designs. head of a pool, retrieval at the tail end of the pool, and onboard We describe the design and construction of two, portable, storage in an accordion-like fashion for swift deployment on floating, radio frequency identification (RFID) PIT tag antenna subsequent occasions. systems: a raft antenna system and a shore-deployed floating The horizontal and vertical antennas were both connected antenna system. To assess the performance of each antenna sys- to an Oregon RFID half-duplex (HDX) multiplex reader (al- tem, we estimated detection distance and detection probability ternating read cycle of six times per second for each antenna), using both experimental and field data. Our research objective which helped prevent proximity detection errors (Aymes and
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 for the raft antenna system was to determine the location and Rives 2009). The HDX reader stored detections for the array fate of PIT-tagged fish in relatively long (i.e., ≥10 km) river along with date and time of detection. Two, 12-V, marine, deep- reaches. For the shore-deployed floating antenna system, our cycle batteries, connected in parallel, powered the raft antenna research objectives were twofold. First, we wanted to compare system. Batteries, tuner boxes, and the reader were placed in abundance estimates to those obtained via electrofishing to de- plastic, top-locking containers and strapped to a rigid plastic termine whether the antenna could be used as an alternative for deck located on the floor of the raft, which prevented equipment estimating abundance of PIT-tagged fish. Second, we wanted to shifts and submersion during deployment (Figure 1). determine the location and fate of fish that were not available Detection distance.—To measure the detection distance of for recapture via electrofishing due to movement from release the horizontal antenna, the raft was elevated on stands, and locations in shorter (i.e., <2 km) river reaches. detection distances were measured by running a 32-mm PIT tag in a perpendicular and parallel orientation past the antenna on ◦ METHODS horizontal, vertical, and 45 detection planes. The horizontal detection plane was defined as the plane extending 180◦ to the Raft Antenna System sides of the antenna, the same plane in which the antenna existed Design and construction.—We designed a two-antenna ar- lying flat on the bottom of the raft, and was used to simulate ray for detecting PIT-tagged fish in a relatively long (≥10 km) detection of PIT-tagged fish near the water surface or in shallow PORTABLE RFID ANTENNA SYSTEMS 1067 Downloaded by [Department Of Fisheries] at 19:14 20 July 2015
FIGURE 1. (A) Schematic representation of the raft antenna system (not to scale). Two continuous loops of 12-gauge THHN wire were used to create the horizontal (heavy dotted line) and vertical (heavy solid line) antennas. Both antennas were connected to tuning boxes (T), which were in turn connected to a multiplex reader (MR) and batteries (B) housed inside plastic, top-locking containers (light gray) and strapped to a rigid plastic deck (dark gray). (B) Diagram of the shore-deployed floating antenna system (not to scale). The antenna consisted of a single loop of 8-gauge, multistrand, speaker wire connected to a tuner box (T), which interfaced with the reader (R) and battery (B) enclosed in the sling-load backpack.
water on the edges of the river. The vertical detection plane was of the antenna detection field. When the tag was detected, a defined as the plane directly under the antenna, 90◦ to the bottom piezoelectric buzzer connected to the reader produced an audible of the raft, and was used to simulate detection when passing beep. Maximum continuous detection distance was defined as directly over a PIT-tagged fish. The 45◦ detection plane was the distance between when a beep was heard for every movement defined as the plane extending at a 45◦ angle to the bottom of the of the tag past the antenna (100% detection rate) and when a raft and was used to simulate detection of fish on the periphery lack of beep indicated that detections were being missed. The 1068 FETHERMAN ET AL.
presence of a diamond-plated, floor-covering, aluminum oar TABLE 1. Placement of rocks marked with PIT tags, with regard to transect frame was thought to affect the horizontal array as metal within number (T), depth (cm), and distance from center (DFC; cm), of the 50 rocks (two on each of 25 transects) used in the raft antenna system detection probability the detection field had been shown to disrupt both detection experiment conducted in the Parvin Lake inlet stream. occurrence and detection distance (Greenberg and Giller 2000; Bond et al. 2007). To test this, maximum detection distances T Depth DFC T Depth DFC T Depth DFC were measured with and without the frame and compared. To evaluate the effects of the oar frame, tag orientation, and 1 38.1 53.3 10 34.3 94.0 19 5.1 346.7 detection plane on maximum detection distance, we used a gen- 1 38.1 190.5 10 20.3 82.6 19 30.5 19.1 eral linear model (GLM) as implemented in SAS ProcGLM 2 64.8 22.9 11 16.5 330.2 20 25.4 171.5 (SAS Institute 2010). We considered an intercept-only model 2 25.4 257.8 11 34.3 25.4 20 22.9 224.8 as well as models that included effects of oar frame only, tag 3 53.3 20.3 12 45.7 180.3 21 25.4 36.8 orientation only, detection plane only, additive effects between 3 35.6 147.3 12 30.5 53.3 21 12.7 450.9 oar frame, orientation, and detection plane, and an interaction 4 44.5 121.9 13 33.0 104.1 22 15.2 237.5 model. Models were ranked using Akaike’s information crite- 4 48.3 138.4 13 30.5 195.6 22 53.3 616.0 rion corrected for small sample sizes (AIC ), compared using 5 52.1 85.1 14 33.0 190.5 23 38.1 2.5 c 5 48.3 82.6 14 39.4 144.8 23 27.9 271.8 AICc differences ( AICc) and weights (wi), and we report pa- rameter estimates and associated SEs from the top-supported 6 33.0 269.2 15 41.9 180.3 24 47.0 22.9 models (Burnham and Anderson 2002). 6 39.4 55.9 15 35.6 94.0 24 26.7 205.7 Detection probability.—To determine the detection proba- 7 44.5 90.2 16 22.9 342.9 25 33.0 381.0 bility of the horizontal antenna under riverine conditions, we 7 25.4 214.6 16 30.5 358.1 25 17.8 15.2 epoxied fifty 32-mm HDX PIT tags to rocks and placed them 8 17.8 299.7 17 35.6 108.0 in a 95.5-m (5.6 m average width) reach of an inlet stream lo- 8 25.4 35.6 17 30.5 257.8 cated at Parvin Lake, Red Feather Lakes, Colorado. We divided 9 17.8 194.3 18 22.9 298.5 the reach into transects (N = 25) and two rocks were placed 9 44.5 293.3 18 29.2 158.8 on each transect. The first transect was located 6.8 m down- stream from the raft put-in site, which allowed the raft to be underway before detection occurred. Subsequent transects were in program MARK (White and Burnham 1999). The Huggins located at random distances from the first by using a random closed capture–recapture estimator differs from a traditional number generator. Depths and locations of the rocks were cho- closed capture–recapture estimator in that only two parameters, sen deliberately to provide a variety of distances and depths for capture or detection probability (p) and recapture probability analysis. Rocks were placed such that PIT tags were oriented (c) are included in the likelihood; N is conditioned out of the parallel to the banks, similar to the orientation of fish facing likelihood and estimated as a derived parameter using estimates into the current. Distance from the south bank and water depth of p (Huggins 1989). This quality allows individual covariates was recorded for each tag, and the metric distance from center affecting p to be included in Huggins estimator (Huggins 1991). (DFC) was calculated by dividing the transect length in half and Primary assumptions are that tags are not lost, are correctly subtracting the distance from south bank (Table 1). identified, and that the system is closed. Ten passes were conducted to estimate detection probability. Encounter histories were constructed for each rock such that A crew of three was used to operate the raft, and the raft was if it were detected on a pass, it would be given a value of 1, and
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 maneuvered down the center of the inlet stream on each pass to if it were not detected, it would be given a value of 0. Depth reduce bias because the raft operators were the same people that and DFC were included as individual covariates for each rock, placed the tagged rocks. Additionally, approximately the same and models in which p was constant, varied by depth, DFC, course was followed on all 10 passes to reduce bias. The raft was or the additive combination of the two, were included in the maneuvered at a 45◦ angle to the banks, providing a detection model set; p equaled c in all models, as c was not expected to be field roughly 4.2 m wide, including both the 1.2- and 4-m axes affected by previous detection. Models were ranked using AICc, of the antenna, and allowing raft operators to maneuver through compared using AICc and wi (Burnham and Anderson 2002), meanders and avoid obstacles. Unfortunately, the shallow nature and a model-averaged parameter estimate and unconditional SE of the inlet stream precluded the use of the vertical antenna in were reported (Anderson 2008). Information from all models the detection probability field test. Detection probability for the with wi > 0 were included in the model-averaged parameter raft antenna system, including continuous deployment of the estimate. In addition, cumulative AICc weights were used to horizontal antenna and intermittent deployment of the vertical assess the relative importance of each covariate. antenna, was later estimated following deployment and PIT- Location and fate of PIT-tagged fish.—Our objective was to tagged fish detection in the Cache la Poudre River. detect and locate PIT-tagged fish that had been marked and re- Detection probability (p) for the horizontal antenna was es- leased in the Cache la Poudre River near Glen Echo, Colorado, timated using the Huggins closed capture–recapture estimator during a previous study (Fetherman 2013). We were primarily PORTABLE RFID ANTENNA SYSTEMS 1069
FIGURE 2. Study area in the Cache la Poudre River, Colorado. The map shows the reach in which the raft antenna system was deployed to determine PIT-tagged fish location and fate, the locations where Brown Trout were and were not removed from the river as part of the study conducted by Fetherman (2013), the stationary antenna locations, and the four study segments (CPR1, CPR2, CPR3, and CPR4) in which abundance was estimated using both the shore-deployed floating antenna system and electrofishing.
interested in fish that were not available for detection by station- also documented a decrease in detection distance associated ary antennas deployed in the river because they had migrated with the metal oar frame (see below). The raft antenna system from the initial study sections, or they remained in the study was deployed in low-water conditions, which were conducive sections and could not be detected by the stationary antennas. to higher detection probabilities but made maneuvering the raft During the initial tagging event (Fetherman 2013), we PIT- difficult; however, we attempted to maneuver the raft within the tagged 5,271 Rainbow Trout and Brown Trout with 32-mm HDX river’s thalweg. The horizontal antenna was deployed contin- tags, inserted posterior to the pectoral fin through the midventral uously to detect fish in shallower reaches (<1 m), while the body wall into the peritoneal cavity using a hypodermic needle vertical antenna was intermittently deployed in pools ≥ 1m Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 (Prentice et al. 1990; Acolas et al. 2007). At the time of tagging, deep. Rainbow Trout averaged 177 mm TL, ranging from 116 to Two passes were made through the 11.3-km reach on sub- 236 mm TL, and Brown Trout averaged 274 mm TL, ranging sequent days. Raft course, restricted within the deeper water from 107 to 500 mm TL. Fish were marked 1 year prior to raft of the thalweg, and vertical antenna deployment within deep antenna deployment. pools were similar among the passes. Before the start of each The raft antenna system was deployed in an 11.3-km section pass, operators’ watches were synchronized with the PIT tag of the Cache la Poudre River (19.6 m average width; Figure 2). reader clock. We recorded start and stop times, as well as times Prior to deployment, the horizontal and vertical antennas were at which recognizable landmark features adjacent to the river tuned and tested with a PIT tag to ensure proper operation. A were passed, allowing us to pair PIT tag detection times and crew of six was used to maneuver the raft: a captain, four pad- locations for analysis. dlers, and a person to operate the antenna equipment and deploy Following deployment in the Cache la Poudre River, two- the vertical antenna in pools. Four paddlers were needed because pass detection data from the horizontal and vertical antennas previous research indicated that a metal oar frame would inter- were pooled to obtain an overall estimate of p and c for the fere with the operation of the antennas (Greenberg and Giller raft antenna system using the Huggins closed capture–recapture 2000; Bond et al. 2007). In assessing our antenna systems, we estimator in program MARK. To meet the closure assumption, 1070 FETHERMAN ET AL.
passes were made on consecutive days over a relatively long maneuvered downstream, allowing the river current to carry the stretch of river, making it unlikely that fish died or moved out of antenna over the majority of obstacles, primarily large boulders. the study section between passes. Although fish were previously Detection distance.—Detection distance was tested by run- marked with PIT tags and released, the raft antenna deployment ning a 32-mm PIT tag over the antenna in the horizontal, vertical, was treated as a traditional mark–recapture study in which the and 45◦ detection planes, holding the tag at both a perpendicular first pass was the “mark” pass and the second pass was the and parallel orientation to the antenna. Both sides of the antenna “recapture” pass. Encounter histories were constructed for each were tested to determine whether there were differences in de- fish such that fish encountered only on the first pass had a history tection symmetry. When the tag was detected by the reader, a of “10,” fish encountered only on the second pass had a history piezoelectric buzzer attached to the reader produced an audible of “01,” and fish encountered on both passes had a history of beep. Maximum continuous detection distance was determined “11.” To estimate p and c, the model set included models in as the distance between when a beep was heard for every move- which p and c were constant and equal, constant but not equal, ment of the tag past the antenna (100% detection rate); when or varied by species. Models were ranked using AICc, compared an audible beep was not produced, this indicated that detections using AICc and wi, and we report model-averaged parameter where being missed. estimates and associated SEs (models with wi > 0). To evaluate antenna symmetry and influence of detection Locations of PIT-tagged fish detected by the raft antenna plane on maximum detection distance, we used a general linear system were determined by comparing the times at which the model (GLM) as implemented in SAS ProcGLM. We considered fish were detected to the times at which recognizable landmark an intercept-only model, as well as models that included effects features adjacent to the river were passed. To determine fate, of side only and detection plane only and models with additive tag numbers detected by the raft antenna system were compared and interactive effects between side and detection plane. Models with the release information associated with the initial PIT tag- were ranked using AICc and compared using AICc and wi, and ging event and the tag numbers obtained from the stationary we report parameter estimates and associated SEs from the top antennas deployed in the Cache la Poudre River. Numbers of supported models. fish detected by the raft antenna system, fish location, and fish Abundance estimation.—One of our primary objectives was fate are reported. to determine whether the shore-deployed floating antenna could be used to estimate abundance of PIT-tagged fish. To accom- Shore-deployed Floating Antenna System plish this, we compared estimates obtained from antenna de- Design and construction.—We designed a river-spanning, tection data with traditional electrofishing estimates. We de- shore-deployed, floating antenna system for detecting PIT- ployed the antenna in 10 river segments, four segments in the tagged fish in shorter (i.e., <2 km) river reaches. The antenna Cache la Poudre River (Figure 2) and six segments in the St. was designed to not only determine the location and fate of fish Vrain River, Lyons, Colorado (Figure 3). Passive integrated that were not available for detection via electrofishing, but also transponder-tagged fish within the Cache la Poudre River were to test its utility for use in place of traditional methods, such as the same as those previously described for the raft antenna sys- electrofishing, for estimating the abundance of PIT-tagged fish. tem (N = 5,271). Rainbow Trout and Brown Trout in the St. The shore-deployed floating antenna was rectangular in Vrain River were PIT-tagged and released 1 year prior to abun- shape (14.6 × 0.6 m) and consisted of a single loop of in- dance estimation (N = 1,569). At the time of tagging, Rainbow sulated, 8-gauge, multistrand, copper speaker wire. Antenna Trout averaged 192 mm TL, ranging from 113 to 272 mm TL, shape was maintained by threading the wire through foam pipe and Brown Trout averaged 190 mm TL, ranging from 120 to
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 insulation (used for flotation) and 13-mm PVC crossbeams, lo- 480 mm TL. Fish that were PIT-tagged within the St. Vrain cated every 1.8 m along the length of the antenna (Figure 1). River were being used to evaluate and compare movement rates Floating nylon rope was threaded through the upstream side of through a whitewater park reach (WWP) containing artificially the foam pipe insulation, allowing operators to maintain tension constructed flow-control structures and through a natural reach and antenna shape during deployment. An Oregon RFID HDX (natural) where no such flow-control structures existed. single reader and tuner box, located in the top compartment of a Sampling segments varied in length, width, and habitat char- plastic-framed, sling-load pack and a 12-V, marine, deep-cycle acteristics. In the Cache la Poudre River, the two most upstream battery, secured to the pack via the sling, were used to power segments (CPR1: 114 m long × 17 m wide; CPR2: 165 × the antenna (Figure 1). 12 m) were characterized by slower-velocity pool habitat with Antenna design facilitated two-person deployment, with both higher-velocity riffles on the upstream end. The two most operators walking along the banks of the river to avoid scaring downstream segments (CPR3: 91 × 21 m; CPR4: 124 × 19 m) fish out of the study segments by wading during deployment. were both characterized by moderate riffle habitat with higher- One person carried the sling-load pack and was the primary velocity riffle habitat on the upstream end (Figure 2). The Cache guide for the antenna. The second person retained tension on la Poudre River was sampled at an average discharge of 3.4 m3/s, the nylon rope in order to maintain antenna shape and guide and average segment depth was 0.6 m. Segments within the the antenna over obstacles. The fully extended antenna was St. Vrain River consisted of naturally occurring and artificially PORTABLE RFID ANTENNA SYSTEMS 1071
FIGURE 3. Study area in the St. Vrain River, Colorado. The map shows the two 0.8-km reaches in which the shore-deployed floating antenna system was deployed to determine PIT-tagged fish location and fate, and the six study segments in which abundance was estimated using both the shore-deployed floating antenna system and electrofishing.
constructed pools of varying depth. The three artificial pools Abundance estimates of PIT-tagged fish per segment were (lower, middle, and upper WWP) were 1.6, 1.4, and 2.1 m obtained using the Huggins closed capture–recapture estimator deep, respectively, and the three natural pools (lower, middle, in program MARK. Encounter histories and model set were and upper natural) were 0.5, 0.4, and 1.0 m deep, respectively constructed similar to those described for estimating detection (Figure 3). Average width of segments was 7.6 m, and the probability of the raft antenna system. Data for each segment segments were sampled at an average discharge of 0.7 m3/s. was analyzed separately (10 analyses in total). Models were Fish abundance was estimated by making two passes with ranked using AICc and compared using AICc and wi, and we the fully extended, shore-deployed, floating antenna through a report model-averaged parameter estimates and associated SEs > Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 study segment; the antenna was folded up and returned to the (models with wi 0). top of the segment between passes. The first pass with the an- Estimates of PIT-tagged fish abundance obtained from the tenna was considered the “mark” pass, while the second pass shore-deployed floating antenna were compared with those ob- acted as the “recapture” pass. In the St. Vrain River, closure tained from the same study segment using electrofishing. Im- for the estimates was achieved by using block nets at both the mediately following antenna deployment, two-pass (Cache la upstream and downstream ends of the study segments. In the Poudre River) or three-pass (St. Vrain River) removal abun- Cache la Poudre River, we used stream features, such as high- dance estimates were conducted in the same segment using a velocity riffles on the upstream end of the segments, to restrict bank electrofishing unit with four electrodes. All fish captured fish movement during the estimates. Use of natural barriers is were weighed, measured, and scanned for PIT tags using a hand- the protocol currently used by Colorado Parks and Wildlife to held reader. Tag numbers were compared with those recorded conduct electrofishing population estimates, and we wanted to by the shore-deployed floating antenna to determine whether use a similar protocol when deploying the antenna. However, the same fish were detected by both gears. we realized that we were assuming closure and it was not guar- Abundance estimates were similarly obtained from the elec- anteed. As in the St. Vrain River, passes were conducted one trofishing data using the Huggins closed capture–recapture es- immediately after the other, and wading was avoided to prevent timator in program MARK, with the exception that c was fixed scaring fish from the segment. to zero in all models because fish were removed after capture 1072 FETHERMAN ET AL.
TABLE 2. Model selection results for factors influencing maximum detection distance for the horizontal antenna of the raft antenna system. The maximized log-likelihood [log(L)], the number of parameters (K) in each model, and AICc values are shown. Models are ranked by their AICc differences ( i) relative to the best model in the set and Akaike weights (wi) quantify the probability that a particular model is the best model in the set given the data and the model set.
2 Model R log(L) K AICc i wi Orientation × Frame × Plane 0.90 −207.46 12 441.30 0.00 1.00 Orientation + Plane 0.87 −225.11 5 460.65 19.35 0.00 Orientation + Frame + Plane 0.87 −224.22 7 463.25 21.96 0.00 Orientation 0.84 −237.02 2 478.13 36.83 0.00 Orientation + Frame 0.84 −236.27 4 480.82 39.52 0.00 Intercept only 0.00 −370.24 1 742.52 301.22 0.00 Plane 0.02 −368.50 3 743.16 301.86 0.00 Frame 0.00 −370.13 2 744.34 303.04 0.00 Frame + Plane 0.03 −368.38 5 747.19 305.89 0.00
(Hense et al. 2010; Saunders et al. 2011). Encounter histories tagging event and tag numbers of fish captured via electrofish- were constructed for only PIT-tagged fish and were constructed ing. Numbers of fish detected by the shore-deployed floating such that fish captured on the first pass had a history of “10” and antenna system, fish location, and fish fate are reported. fish caught on the second pass had a history of “01.” Rainbow Detection probability for the shore-deployed floating antenna Trout and Brown Trout were included as groups in the same was estimated for each reach following deployment. Encounter analysis and p was modeled as constant or varied by length, histories and model sets were constructed similar to those de- species, or both. Models were ranked using AICc and compared scribed for estimating detection probability of the raft antenna using AICc and wi, and we report model-averaged parameter system. Models were ranked using AICc and compared using estimates and associated 95% CIs (models with wi > 0). Abun- AICc and wi, and we report model-averaged parameter esti- dance estimates obtained with the two gear types were compared mates and associated SEs (models with wi > 0). within each study segment, and differences between the gears were determined by a lack of overlap in 95% CIs. RESULTS Location and fate of PIT-tagged fish.—Another primary ob- jective of the shore-deployed floating antenna was to determine Raft Antenna System the location and fate of PIT-tagged fish that were not available Detection distance.—Tag orientation, oar frame presence, for recapture via electrofishing in the St. Vrain River study seg- and detection plane had the largest influence on hand-held PIT ments. To address this objective we deployed the antenna in tag maximum detection distance of the horizontal antenna of two 0.8-km reaches of the St. Vrain River (Figure 3). In the the raft antenna system (AICc weight = 1.00; Table 2). Average natural reach, the antenna was deployed just upstream from the detection distance was greater for tags oriented perpendicular upper natural study segment to just downstream from the lower than for those oriented parallel to the antenna. The presence natural study segment. In the WWP reach, the antenna was de- of the oar frame affected detection distance along the vertical, ployed just upstream from the upper WWP study segment to ◦
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 horizontal, and 45 detection planes for tags oriented perpendic- just downstream from the lower WWP study segment. All six of ular to the antenna. Detection distance in the vertical detection the previously described study segments were contained within plane was affected by the presence of the oar frame for tags ori- the two reaches. Shallow riffle habitat constituted the majority ented parallel to the antenna, but were not affected by oar frame of the habitat between study segments in both reaches. presence in the horizontal or 45◦ detection planes (Figure 4). Two passes with the antenna were made through each reach. Detection probability.—The use of PIT-tagged rocks in the Due to the narrow width of the reaches, the array was deployed Parvin Lake inlet stream showed that both DFC and depth af- ◦ at a 45 angle to the stream banks, allowing full extension of fected p of the horizontal antenna of the raft antenna system the array for proper tuning. Before the start of each pass, op- (model = DFC + Depth; AICc weight = 0.99). A negative re- erators’ watches were synchronized with the reader clock. We lationship was observed between DFC and detection probability recorded start and stop times, as well as times at which recog- (βˆ =−0.011 ± 0.001) indicating that the farther a tag was from nizable landmark features adjacent to or within the river were the center of the inlet stream, and consequently the center of the passed, allowing us to pair PIT tag detection times and locations antenna’s detection field, the less likely it was to be detected. for analysis; average deployment time was roughly 45 min per Interestingly, there was a positive relationship between p and pass. To determine PIT-tagged fish fate, tag numbers were com- depth (βˆ = 0.02 ± 0.005). This relationship likely occurred pared with the release information associated with the initial PIT because tags located near the center of the inlet stream were PORTABLE RFID ANTENNA SYSTEMS 1073
Shore-deployed Floating Antenna System Detection distance.—Detection plane had the largest influ- ence on hand-held PIT tag maximum detection distance for the shore-deployed floating antenna (model 1 = Plane, AICc weight = 0.48; model 2 = Side + Plane, AICc weight = 0.44; Table 3). On average, tags in the vertical detection plane were detected at a greater distance (79.9 ± 0.8 cm) than tags in the horizontal or 45◦ detection planes. Horizontal (71.6 ± 0.7 cm) and 45◦ (72.3 ± 0.8 cm) detection planes exhibited av- erage maximum detection distances. There was some evidence that detection distances were not symmetrical about the antenna (model 2 = Side + Plane, AICc weight 0.44; Table 3). However, average maximum detection distances did not appear to differ between the right (75.3 ± 1.2 cm) and left (74.0 ± 0.9 cm) FIGURE 4. Maximum detection distances for the horizontal antenna of the sides of the antenna. raft antenna system with and without the aluminum oar frame and with tags oriented perpendicular and parallel to the antenna in the vertical (V), horizontal Abundance estimation.— Abundances of PIT-tagged fish ob- (H), and 45◦ detection planes. Error bars represent SE. tained from the shore-deployed floating antenna were similar to those obtained via electrofishing, as indicated by overlapping 95% CIs, in two of the four study segments in the Cache la deeper than those placed closer to the banks and, thus, more Poudre River (CPR2 and CPR3; Figure 5). However, it is im- likely to be detected regardless of depth. The model-averaged portant to note the 95% CIs overlapped in these two segments estimate of p for the horizontal antenna was 0.32 ± 0.02 (mean because of large 95% CIs for the electrofishing estimates (CPR2) ± SE). or the antenna estimates (CPR3), suggesting that the estimates Location and fate of PIT-tagged fish.—A total of 44 PIT- were fairly imprecise, likely due to the low numbers of PIT- tagged fish (32 Rainbow Trout and 12 Brown Trout) were de- tagged fish detected in these segments. In CPR2, 19 PIT tags tected by the 4.2-m-wide antenna detection field (21% of the were detected by the antenna, whereas eight PIT-tagged fish wetted area) over the 11.3-km reach of the Cache la Poudre were detected via electrofishing; four PIT-tagged fish were de- River. Two unique fish were detected by the intermittently de- tected by both gear types. In CPR3, seven PIT tags were detected ployed vertical antenna, whereas 42 unique fish were detected by the antenna, whereas five PIT-tagged fish were detected via by the continuously deployed horizontal antenna. Three fish electrofishing; two PIT-tagged fish were detected by both gear were detected on both passes, which allowed both p and c to be types. estimated for the raft antenna system. Detection probability of Abundance estimates differed in two of the four study seg- the raft antenna system, including continuous deployment of the ments in the Cache la Poudre River (CPR1 and CPR4; Figure 5). horizontal antenna and intermittent deployment of the vertical The estimate of PIT-tagged fish abundance in CPR4 was higher antenna, was 0.14 ± 0.14, and for c was 0.13 ± 0.07. with electrofishing than with the antenna, which was potentially Twenty-seven of the fish detected had never been detected a function of differences in maneuverability of the gear types by the stationary antennas after their release 1 year prior to de- around the large boulders located within this study segment. In ployment of the raft antenna system. Of these, 15 were detected CPR4, six PIT tags were detected by the antenna, whereas 11 Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 in their release location and had never moved past the station- PIT-tagged fish were detected via electrofishing; two PIT-tagged ary antennas. The remaining 12 fish had moved upstream or fish were detected by both gear types. Abundance estimation downstream, but were not detected by the stationary antennas. via electrofishing was not possible in CPR1 due to depletion The complementary use of the raft antenna system enabled us failure (i.e., three PIT-tagged fish caught on both passes); how- to locate and determine the fate of those fish not detected by the ever, abundance estimates were obtainable using the antenna stationary antennas. (Figure 5). In CPR1, 11 PIT tags were detected by the antenna, The other 17 fish detected by the raft antenna system had whereas six PIT-tagged fish were detected via electrofishing; been previously detected by the stationary antennas. Detections five PIT-tagged fish were detected by both gear types. at stationary antennas occurred between 8 and 12 months prior The shore-deployed floating antenna failed to obtain compa- to detection by the raft antenna system, and use of the system rable abundance estimates to those obtained by electrofishing allowed us to confirm the new locations of these fish. Fish moved in the WWP study segments in the St. Vrain River (Figure 5), both upstream and downstream from their release locations. which was likely a function of segment depth. The shallowest On average, fish were located within 1.2 km of the locations study segment (middle WWP, 1.4 m) exceeded the maximum from which they had been released; however, several fish were read range of the array by 0.6 m. All three PIT tags detected by detected farther from their release location including one fish the antenna in the lower and middle WWP segments were not located over 4 km upstream from its release site. detected via electrofishing. 1074 FETHERMAN ET AL.
TABLE 3. Model selection results for factors influencing maximum detection distance for the shore-deployed floating antenna system. The maximized log- likelihood [log(L)], the number of parameters (K) in each model, and AICc values are shown. Models are ranked by their AICc differences ( i) relative to the best model in the set and Akaike weights (wi) quantify the probability that a particular model is the best model in the set given the data and the model set.
2 Model R log(L) K AICc i wi Plane 0.64 −43.43 3 93.49 0.00 0.48 Side + Plane 0.66 −42.27 4 93.63 0.14 0.44 Side × Plane 0.68 −41.33 6 97.06 3.57 0.08 Intercept only 0.00 −64.98 1 132.05 38.56 0.00 Side 0.02 −64.57 2 133.45 39.96 0.00
Similar abundance estimates were obtained in the lower nat- of the antenna by 0.2 m, the abundance estimate obtained us- ural and middle natural study segments in the St. Vrain River ing the antenna was lower than that obtained by electrofishing (Figure 5). Maximum pool depth in these study segments did (Figure 5). Twelve PIT tags were detected by the antenna in not exceed 0.5 m. In the lower natural segment, 20 PIT tags the upper natural segment, whereas 21 PIT-tagged fish were de- were detected by the antenna, whereas 18 PIT-tagged fish were tected via electrofishing; only two PIT-tagged fish were detected detected via electrofishing; six PIT-tagged fish were detected by both gear types. by both gear types. In the middle natural segment, 11 PIT tags Location and fate of PIT-tagged fish.—Thirty-two PIT- were detected by the antenna, whereas 10 PIT-tagged fish were tagged fish—16 Rainbow Trout and 16 Brown Trout—were detected via electrofishing; four PIT-tagged fish were detected detected by the shore-deployed floating antenna, which covered by both gear types. In the upper natural segment, which had a 100% of the wetted area within the 0.8-km natural reach of maximum depth of 1 m and exceeded the maximum read range the St. Vrain River. In the 0.8-km WWP reach, 49 PIT-tagged fish were detected by the shore-deployed floating antenna: 18 Rainbow Trout and 31 Brown Trout. Estimated p for the shore- deployed floating antenna did not differ between the reaches: p = 0.52 ± 0.15 in the natural pool reach and p = 0.60 ± 0.10 in the WWP reach. Of the 32 fish detected by the shore-deployed floating antenna in the natural reach of the St. Vrain River, 12 were detected in the location to which they had been released 6 to 12 months prior to antenna deployment. Eighteen of the 32 fish that had been released within the reach exhibited upstream and down- stream movements, and were detected in locations that differed from their release location. Six of these fish were located within the lower, middle, or upper natural segments and were available for detection, whereas the other 12 fish were located between study segments and were not available for detection via elec-
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 trofishing. The new locations of these fish would not have been known had they not been detected by the shore-deployed float- ing antenna. Finally, two of the fish that had been released in the WWP reach 12 months earlier were located between study segments and were not available for detection via electrofish- ing. The detection of these fish confirmed that fish were able to navigate the artificial structures of the WWP and move up- stream, an observation that was only made possible because the locations of these fish were confirmed by the shore-deployed floating antenna. Of the 49 fish detected by the shore-deployed floating antenna in the WWP reach of the St. Vrain River, 31 were detected in FIGURE 5. Estimated number of PIT-tagged salmonids per section estimated the location in which they had been released. Twelve of the via electrofishing (Cache la Poudre River: two passes; St. Vrain River: three passes) and the shore-deployed floating antenna (two passes) within the study 49 fish were detected downstream from where they had been segments in (A) the Cache la Poudre River and (B) the St. Vrain River. Error released within the WWP reach. Two fish had been released bars represent 95% CI. in the natural reach upstream from the WWP reach and had PORTABLE RFID ANTENNA SYSTEMS 1075
exhibited downstream movement into the WWP reach. Four distance from the detection field of the raft, which suggests that fish were detected upstream from the location in which they had if coverage had been wider, p would have been higher. Coverage been released, again confirming that fish were able to navigate was less in the Cache la Poudre River where the raft antenna the artificial structures of the WWP. All four were located be- system covered only 21% of the wetted area. This resulted in tween the WWP study sections and as such were not available alowerp (0.14) than that of the Parvin Lake field test. Rafting for detection via electrofishing. Therefore, the shore-deployed technique, e.g., pointing the nose of the raft straight downriver floating antenna helped confirm upstream movement through versus maintaining the raft at a 45◦ angle to the banks, could the WWP structures, which would not have been confirmed also affect coverage and therefore p of the raft antenna system. using traditional sampling methods. In contrast, the shore-deployed floating antenna system, which covered 100% of the wetted area, produced higher es- timates of p during deployment, and p was similar between the DISCUSSION natural and WWP reaches in the St. Vrain River. Similarities in p Both our raft and shore-deployed floating antenna systems were likely a result of where fish were detected in these reaches. allowed us to determine the location and fate of PIT-tagged fish In the natural reach, depths rarely exceeded the read range of that had not moved, moved but had not been detected by station- the antenna, which allowed nearly all fish in the reach to be ary antennas or electrofishing, or had migrated from the release available for capture on one or both passes. In the WWP reach, location entirely. Our raft antenna system was deployed over fish were most commonly detected in the shallow riffle habitat an 11.3-km section of the Cache la Poudre River, detecting fish between the deeper study segments, the depth of which did not in locations that would not have otherwise been sampled using exceed the read range of the antenna. Estimation of p requires traditional sampling methods in shorter study segments. Our that at least one fish needs to be detected on both passes (White shore-deployed floating antenna allowed us to obtain estimates et al. 1982). In this study repeat detections on subsequent passes of PIT-tagged fish abundance in shorter study segments in both occurred in the shallower sections of the WWP reach. Fish in the rivers and determine the location and fate of fish in two 0.8- deeper study segments were generally unavailable for detection km reaches of the St. Vrain River. Overall, the raft and shore- on either pass, which likely artificially inflated the estimates of deployed floating antenna systems have overcome limitations p in the WWP reach. Overall, p could be increased by increasing recognized with other portable systems, namely antenna size, the number of passes made through a study section. However, a stream distance surveyed, coverage, and detection distances. balance is needed between the number of passes made, the time Detection distances for both antenna systems were between it takes to make a pass, and the information gained by adding 0.7 and 1.0 m and can be partially attributed to using 32-mm more passes to the study. tags (Zydlewski et al. 2006). Other studies have used 12- Only two fish were detected by the vertical antenna of the or 23-mm tags, which resulted in lower detection distances raft antenna system, potentially contributing to the low p.Lin- (Roussel et al. 2000; Zydlewski et al. 2001; Cucherousset et al. nansaari and Cunjak (2007) suggested that there may be a 2005; Hill et al. 2006). However, one possible disadvantage of fright bias associated with larger submerged antennas if fish greater detection distance is an increased chance of multiple are tracked in their active state. Therefore, our vertical antenna tags being present in the detection field of the array, resulting may have caused some behavioral avoidance when deployed in in no tags being detected (tag collision; Axel et al. 2005; pools. In addition, the vertical antenna was tuned fully extended O’Donnell et al. 2010). In addition, greater detection distance prior to deployment. Antenna shape may have differed when could result in the detection of ghost tags, i.e., tags lodged in deployed in pools, or intermittent deployment and storage may
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 the substrate through a combination of tag loss, predation, and have caused the antenna to become detuned, thereby reducing natural mortality (O’Donnell et al. 2010). Finally, although detection distance of the vertical antenna upon deployment. The detection distances were greater than other portable antenna small number of fish detected by the vertical antenna suggests designs, both the raft and the shore-deployed floating antenna that antenna contributed little to the overall detection of fish systems were still limited to use in shallow (<1m)river by the raft antenna system. However, the utility of the vertical reaches. The utility of both antenna systems for use in deeper antenna could be increased by correcting issues regarding de- river reaches could be increased by manipulating voltage or tuning of the antenna by constructing a more rigid frame that adjusting wire arrangement (e.g., multiple loops or adjusting maintains antenna shape during deployment. spacing between loops) to increase their detection distance. Despite potential issues with detection probability, the ability Despite greater detection distances, detection probability for to determine the location and fate of PIT-tagged fish over short the raft antenna system was relatively low. Low p could be a (0.8 km) and long (11.3 km) distances is a major advantage of result of the coverage of the raft antenna system relative to these antenna systems relative to other portable antenna designs. the width of the river. In the Parvin Lake field test, the raft Most portable antenna designs, with the exception of the boat- antenna system covered an average of 32% of the wetted area mounted antenna for monitoring mussels (Fischer et al. 2012), of the inlet stream. Consequently, p for the horizontal antenna have been constrained to use in shallow, wadable streams; as was 0.32, and modeling indicated that p was affected by tag a result, survey length was limited by the length of river an 1076 FETHERMAN ET AL.
operator could walk. Here, we demonstrated that both antenna also illicit an avoidance response, especially when conditions systems could detect and locate fish that had migrated from are such that shadows are cast by the antenna (e.g., sunny days; smaller study segments and were no longer available for detec- Ellis et al. 2013), fish are less likely, relative to smaller designs, tion by stationary antennas or electrofishing. In several cases, to move completely out of the detection field due to antenna we were able to document upstream movement of fish through coverage (i.e., 100% of the wetted area). The effect of operator structures thought to be a barrier to movement (i.e., WWP struc- experience is also reduced due to antenna coverage as the op- tures in the St. Vrain River), and over long distances (up to 4 km erator is not required to identify specific locations or habitats in the Cache la Poudre River). However, little inference can be to sample. However, fish located directly behind large boulders made regarding fish fate from tags that were detected in the or other obstacles may not be detected by the floating antenna same location or downstream from where they had been re- when passing over these obstacles. leased. Ghost tags (O’Donnell et al. 2010) could potentially be The design flexibility of these antenna systems provides an detected by these antenna systems, leading to incorrect inter- opportunity to potentially combine designs and create a larger pretations of the data regarding fish location and fate. Ghost tag detection field for greater river coverage. For example, wing- detections cannot be removed from the data without locating like floating antennas could be combined with the raft antenna the tags using a smaller wand-type antenna or electrofishing to system to create a larger array that could cover more of the confirm whether tags were retained or lost by the fish to which wetted area over long distances. Multiplexers, or a well-designed they are associated. master–slave set-up, could be used to power the system and Portable PIT-tag antenna systems are fairly accurate in prevent proximity detection errors (Aymes and Rives 2009). estimating abundance of PIT-tagged fish in small streams However, the larger size of the system could potentially result in (O’Donnell et al. 2010; Sloat et al. 2011). Portable antennas a greater chance of entanglement with obstacles such as boulders have the advantage of allowing frequent sampling for abun- or submerged trees, and this would need to be considered during dance estimation and fish location without subjecting the fish to the design of these larger systems. excessive handling stress or mortality, and by minimizing dis- Our portable antenna systems provide noninvasive methods turbance to individuals (Sloat et al. 2011). However, this feature that minimize disturbance to individual fish for determining the also excludes the ability to examine fish for growth or physi- fate of PIT-tagged fish in both small (hundreds of meters) and ological parameters (Zydlewski et al. 2001) or to estimate the large (kilometers) river reaches. Through the use of marker tags, overall abundance (marked and unmarked) of fish within a des- accurate timing devices, and submeter GPS, the location of fish ignated area. Our results suggest that if estimates of tagged fish can be determined fairly accurately using these systems. In addi- are desired, and handling fish (beyond tagging) is not neces- tion, the shore-deployed floating antenna provides an alternative sary to collect individual information (e.g., fish size or signs of to traditional sampling methods for estimating PIT-tagged fish disease), portable antennas such as the shore-deployed floating abundance. More research is needed to examine the effects of antenna present an alternative to traditional sampling methods ghost tags on inferring the fate of fish without physical recap- such as electrofishing. However, estimates may be imprecise. In tures, to assess fish behavioral responses to antenna systems, addition, detection distance limitations must be considered, as and to determine ways to increase p and reduce variability in abundance may be greatly underestimated in deeper study seg- abundance estimates. ments, e.g., those of the WWP reach. Finally, our data regarding the small number of tags detected by both gears suggest that ghost tags (O’Donnell et al. 2010) could influence abundance ACKNOWLEDGMENTS
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 estimates obtained with the shore-deployed floating antenna This work was sponsored in part by Colorado Parks and system. Wildlife and the Colorado Cooperative Fish and Wildlife Re- The shore-deployed floating antenna system overcomes some search Unit at Colorado State University, and funding was pro- of the limitations observed with other antenna systems, such as vided in part by the Federal Aid in Sport Fish Restoration pro- antenna size, or those caused by lack of operator experience gram, project F-394R. We thank G. Schisler, K. Davies, M. and fish behavior (O’Donnell et al. 2010). Many of the previ- Kondratieff, and C. Praamsma of Colorado Parks and Wildlife, ously described portable antenna systems were small, designed and L. Bailey, G. Fraser, and N. Shannon for their help with to be operated by one person in a small stream (Roussel et al. antenna design and deployment. Any use of trade, firm, or prod- 2000; Cucherousset et al. 2005; Hill et al. 2006), and as such, uct names is for descriptive purposes only and does not imply antenna coverage was small relative to the width of the river. endorsement by the U.S. Government. This study was approved Our shore-deployed floating antenna is the largest two-person under Colorado State University Institutional Animal Care and portable antenna described to date, as previous two-person an- Use Committee 10-1956A. tennas did not exceed 5 m in length (Linnansaari and Cunjak 2007). Submersion of the antenna is not required, theoretically REFERENCES reducing the chance of a behavioral response to an antenna lo- Acolas, M. L., J. -M. Roussel, J. M. Lebel, and J. L. Bagliniere.` 2007. Laboratory cated within the water column. Although overhead stimuli may experiment on survival, growth and tag retention following PIT injection into PORTABLE RFID ANTENNA SYSTEMS 1077
the body cavity of juvenile Brown Trout (Salmo trutta). Fisheries Research Hense, Z., R. W. Martin, and J. T. Petty. 2010. Electrofishing capture efficien- 86:280–284. cies for common stream fish species to support watershed-scale studies in Anderson, D. R. 2008. Model based inference in the life sciences: a primer on the central Appalachians. North American Journal of Fisheries Management evidence. Springer, New York. 30:1041–1050. Axel, G. A., E. F. Prentice, and B. P. Sandford. 2005. PIT-tag detection system Hill, S. A., G. B. Zydlewski, J. D. Zydlewski, and J. M. Gasvoda. 2006. Develop- for larger-diameter juvenile fish bypass pipes at Columbia River basin hydro- ment and evaluation of portable PIT tag detection units: PITpacks. Fisheries electric dams. North American Journal of Fisheries Management 25:646–651. Research 77:102–109. Aymes, J. C., and J. Rives. 2009. Detection efficiency of multiplexed passive Huggins, R. M. 1989. On the statistical analysis of capture experiments. integrated transponder antennas is influenced by environmental conditions Biometrika 76:133–140. and fish swimming behavior. Ecology of Freshwater Fish 18:507–513. Huggins, R. M. 1991. Some practical aspects of a conditional likelihood ap- Bond, M. H., C. V. Hanson, R. Baertsch, S. A. Hayes, and R. B. MacFarlane. proach to capture experiments. Biometrics 47:725–732. 2007. A new low cost instream antenna system for tracking passive integrated Linnansaari, T. P., and R. A. Cunjak. 2007. The performance and efficacy of a transponder (PIT)-tagged fish in small streams. Transactions of the American two-person operated portable PIT-antenna for monitoring spatial distribution Fisheries Society 136:562–566. of stream fish populations. River Research and Applications 23:559–564. Burnham, K. P., and D. R. Anderson. 2002. Model selection and multimodel Morhardt, J. E., D. Bishir, C. I. Handlin, and S. D. Mulder. 2000. A portable inference: a practical information-theoretic approach, 2nd edition. Springer- system for reading large passive integrated transponder tags from wild trout. Verlag, New York. North American Journal of Fisheries Management 20:276–283. Compton, R. I., W. A. Hubert, F. J. Rahel, M. C. Quist, and M. R. Bower. 2008. Nunnallee, E. P., E. F. Prentice, B. F. Jonasson, and W. Patten. 1998. Evaluation Influences of fragmentation on three species of native warmwater fish in a of a flat-plate PIT tag interrogation system at Bonneville Dam. Aquaculture Colorado River basin headwater stream system, Wyoming. North American Engineering 17:261–272. Journal of Fisheries Management 28:1733–1743. O’Donnell, M. J., G. E. Horton, and B. H. Letcher. 2010. Use of portable an- Connolly, P. J., I. G. Jezorek, K. D. Martens, and E. F. Prentice. 2008. Measuring tennas to estimate abundance of PIT-tagged fish in small streams: factors the performance of two stationary interrogation systems for detecting down- affecting detection probability. North American Journal of Fisheries Man- stream and upstream movement of PIT-tagged salmonids. North American agement 30:323–336. Journal of Fisheries Management 28:402–417. Prentice, E. F., T. A. Flagg, C. S. McCutcheon, D. F. Brastow, and D. C. Cross. Cucherousset, J., J. R. Britton, W. R. C. Beaumont, M. Nyqvist, K. Sievers, 1990. Equipment, methods, and an automated data-entry station for PIT and R. E. Gozlan. 2010. Determining the effect of species, environmental tagging. Pages 335–340 in N. C. Parker, A. E. Giorgi, R. C. Heidinger, D. B. conditions and tracking method on the detection efficiency of portable PIT Jester Jr., E. D. Prince, and G. A. Winans, editors. Fish marking techniques. telemetry. Journal of Fish Biology 76:1035–1045. American Fisheries Society, Symposium 7, Bethesda, Maryland. Cucherousset, J., J.-M. Paillisson, and J.-M. Roussel. 2007. Using PIT technol- Roussel, J.-M., A. Haro, and R. A. Cunjak. 2000. Field test of a new method ogy to study the fate of hatchery-reared YOY Northern Pike released into for tracking small fishes in shallow rivers using passive integrated transpon- shallow vegetated areas. Fisheries Research 85:159–164. der (PIT) technology. Canadian Journal of Fisheries and Aquatic Sciences Cucherousset, J., J.-M. Roussel, R. Keeler, R. A. Cunjak, and R. Stump. 2005. 57:1326–1329. The use of two new portable 12-mm PIT tag detectors to track small fish SAS Institute. 2010. SAS system software, release 9.2. SAS Institute, Cary, in shallow streams. North American Journal of Fisheries Management 25: North Carolina. 270–274. Saunders, W. C., K. D. Fausch, and G. C. White. 2011. Accurate estimation of Ellis, T. R., T. Linnansaari, and R. A. Cunjak. 2013. Passive integrated transpon- salmonid abundance in small streams using nighttime removal electrofish- der (PIT) tracking versus snorkeling: quantification of fright bias and com- ing: an evaluation using marked fish. North American Journal of Fisheries parison of techniques in habitat use studies. Transactions of the American Management 31:403–415. Fisheries Society 142:660–670. Sloat, M. R., P.F. Baker, and F. K. Ligon. 2011. Estimating habitat-specific abun- Fetherman, E. R. 2013. Introduction and management of Myxobolus cerebralis- dances of PIT-tagged juvenile salmonids using mobile antennas: a compari- resistant Rainbow Trout in Colorado. Doctoral dissertation. Colorado State son with standard electrofishing techniques in small streams. North American University, Fort Collins. Journal of Fisheries Management 31:986–993. Ficke, A. D., M. C. Kondratieff, and C. A. Myrick. 2012. The effects of PIT White, G. C., D. R. Anderson, K. P. Burnham, and D. L. Otis. 1982. Capture- tagging on the swimming performance and survival of three nonsalmonid recapture and removal methods for sampling closed populations. Los Alamos freshwater fishes. Ecological Engineering 48:86–91. National Laboratory, LA-8787-NERP, Los Alamos, New Mexico. Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 Fischer, J. R., T. E. Neebling, and M. C. Quist. 2012. Development and eval- White, G. C., and K. P. Burnham. 1999. Program MARK: survival esti- uation of a boat-mounted RFID antenna for monitoring freshwater mussels. mation from populations of marked animals. Bird Study 46(Supplement): Freshwater Science 31:148–153. 120–138. Greenberg, L. A., and P. S. Giller. 2000. The potential of flat-bed passive Zydlewski, G. B., A. Haro, K. G. Whalen, and S. D. McCormick. 2001. Perfor- integrated transponder antennae for studying habitat use by stream fishes. mance of stationary and portable passive transponder detection systems for Ecology of Freshwater Fish 9:74–80. monitoring of fish movements. Journal of Fish Biology 58:1471–1475. Gries, G., and B. H. Letcher. 2002. Tag retention and survival of age-0 Atlantic Zydlewski, G. B., G. Horton, T. Dubreuil, B. Letcher, S. Casey, and J. Zydlewski. Salmon following surgical implantation with passive integrated transponder 2006. Remote monitoring of fish in small streams: a unified approach using tags. North American Journal of Fisheries Management 22:219–222. PIT tags. Fisheries 31:492–502. This article was downloaded by: [Department Of Fisheries] On: 20 July 2015, At: 19:14 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered office: 5 Howick Place, London, SW1P 1WG
North American Journal of Fisheries Management Publication details, including instructions for authors and subscription information: http://www.tandfonline.com/loi/ujfm20 A Market Model of Eastern Bering Sea Alaska Pollock: Sensitivity to Fluctuations in Catch and Some Consequences of the American Fisheries Act James W. Stronga & Keith R. Criddlea a School of Fisheries and Ocean Sciences, University of Alaska Fairbanks, 17101 Point Lena Loop Road, Juneau, Alaska 99801, USA Published online: 16 Oct 2014.
Click for updates
To cite this article: James W. Strong & Keith R. Criddle (2014) A Market Model of Eastern Bering Sea Alaska Pollock: Sensitivity to Fluctuations in Catch and Some Consequences of the American Fisheries Act, North American Journal of Fisheries Management, 34:6, 1078-1094, DOI: 10.1080/02755947.2014.944678 To link to this article: http://dx.doi.org/10.1080/02755947.2014.944678
PLEASE SCROLL DOWN FOR ARTICLE
Taylor & Francis makes every effort to ensure the accuracy of all the information (the “Content”) contained in the publications on our platform. However, Taylor & Francis, our agents, and our licensors make no representations or warranties whatsoever as to the accuracy, completeness, or suitability for any purpose of the Content. Any opinions and views expressed in this publication are the opinions and views of the authors, and are not the views of or endorsed by Taylor & Francis. The accuracy of the Content should not be relied upon and should be independently verified with primary sources of information. Taylor and Francis shall not be liable for any losses, actions, claims, proceedings, demands, costs, expenses, damages, and other liabilities whatsoever or howsoever caused arising directly or indirectly in connection with, in relation to or arising out of the use of the Content. This article may be used for research, teaching, and private study purposes. Any substantial or systematic reproduction, redistribution, reselling, loan, sub-licensing, systematic supply, or distribution in any form to anyone is expressly forbidden. Terms & Conditions of access and use can be found at http:// www.tandfonline.com/page/terms-and-conditions North American Journal of Fisheries Management 34:1078–1094, 2014 C American Fisheries Society 2014 ISSN: 0275-5947 print / 1548-8675 online DOI: 10.1080/02755947.2014.944678
ARTICLE
A Market Model of Eastern Bering Sea Alaska Pollock: Sensitivity to Fluctuations in Catch and Some Consequences of the American Fisheries Act
James W. Strong and Keith R. Criddle* School of Fisheries and Ocean Sciences, University of Alaska Fairbanks, 17101 Point Lena Loop Road, Juneau, Alaska 99801, USA
Abstract A simultaneous-equation equilibrium model of domestic and international markets for fillets, surimi, and roe of Alaska pollock (i.e., Walleye Pollock Gadus chalcogrammus) is used to examine the combined effect of variability in landings and changes in product mix on first wholesale prices and revenues. From the commencement of American- ization in the mid-1970s through the joint venture era of the 1980s and following full Americanization in early 1990s, U.S. producers of Alaska pollock focused mainly on selling surimi and roe to Japan. Passage of the American Fisheries Act of 1998 (AFA) unleashed an economic transformation that more than doubled fishery revenues without increasing exploitation rates or catches. The increased value came about because freedom from the race for fish allowed fishers and processors to increase product quality, improve product recovery rates, diversify their mix of product forms, and develop new markets in the USA and Europe. Despite the size and importance of the Alaska pollock fishery—among U.S. fisheries, it consistently ranks number one by volume and in the top five by value—little attention has been given to the development of a formal model of Alaska pollock markets or the market impacts of the AFA. The results of the model indicate that first wholesale revenues are maximized when harvests are maximized, but the value of the harvest is particularly sensitive to changes in the demand for Alaska pollock harvested in Russian waters, price differentials between fillet and surimi, and Japanese inventories of Alaska pollock products. Comparative static simulations com- paring pre- and post-AFA fishery product recovery rates and product mix indicate revenues would have increased by 62% from 1998 to 2007.
The eastern Bering Sea fishery for Alaska pollock (i.e., Wall- the fishery. This led to rapid increases in the number of U.S.- Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 eye Pollock Gadus chalcogrammus) yields gross ex-vessel rev- flagged harvesters and catcher–processors and the development enues of over US$300 million and a first wholesale value of over of extensive shore-based processing capacity in the 1980s and $1.2 billion. Over the last 30 years, annual total allowable catch 1990s. In the mid-1980s, the total allowance level of foreign (TAC) has ranged from 813,000 metric tons to 1,494,900 metric fishing was zeroed-out in favor of allocations to joint ventures. tons, making the Alaska pollock fishery the largest U.S. fishery By 1991, the joint-venture allocations had also been eliminated in tonnage. During the 1960s and 1970s, Alaska pollock was to make room for harvesting by U.S.-flagged catcher vessels that largely harvested by a Japanese mothership fleet that extracted delivered to U.S.-flagged motherships or processing facilities in roe and processed the meat into surimi, a fish-based protein the USA and for harvesting by U.S.-flagged catcher–processors paste. A mostly Russian mothership fleet produced a frozen and U.S.-flagged catcher vessels that delivered to them. Con- headed and gutted (H&G) product destined for consumption tinued increases in harvesting and processing capacity during in the Eastern Bloc. Passage of the Fishery Conservation and the 1990s precipitated a race for fish and engendered heated po- Management Act in 1976 stimulated an “Americanization” of litical battles over allocations between the at-sea (motherships,
*Corresponding author: [email protected] Received June 21, 2013; accepted July 1, 2014
1078 ALASKA POLLOCK MARKETS 1079
catcher–processors, and associated catcher vessels) and inshore range of catches, implying that ex-vessel revenues could be max- sectors (Strong and Criddle 2013). While a few processors in- imized by harvesting the maximum permissible level of catch. vested in fillet technology, roe and surimi remained the primary However, because the AFA led to substantial reorganization of products. During the race for fish, it was not advantageous for production, that model is unlikely to provide a good representa- processors to produce time-intensive product such as fillets. In- tion of the effect of variation in catch limits under current market stead, they focused on high-throughput production methods to conditions. Fell (2008) used cointegration methods to examine accommodate the time-compressed volumes of catch that re- the responsiveness of Alaska pollock processors to changes in sulted from the race for fish. surimi and fillet prices before and after implementation of the These conditions changed following passage of the Amer- AFA and found that price responsiveness increased in both the ican Fisheries Act (AFA) in 1998. The AFA established per- inshore and at-sea sectors after the AFA was implemented. Mor- manent sector allocation shares for each fleet sector and autho- rison Paul et al. (2009) used a revenue function model to examine rized the formation of harvesting cooperatives in which fishing which product form was chosen by catcher–processors pre- and operations privately negotiated contracts that guaranteed each post-AFA and found that the mix of regular and deep-skin fillets, vessel a portion of the TAC allocated to that sector (Criddle surimi, and roe was responsive to variations in product prices and Macinko 2000; Matulich et al. 2001; Anderson 2002). The and increased operational flexibility that followed implementa- catcher–processor sector and catcher vessels that delivered to tion of the AFA. These choices of product form, in turn, affected them formed cooperatives in 1999. Catcher vessels that deliv- operation choices related to tow duration, crewing, and vessel ered to inshore processors formed cooperatives in 2000, as did power. the catcher vessels that delivered to motherships. Of the Alaska This study uses monthly time series data from 2000 to 2008 pollock allocated to these groups, the catcher–processor sector to extend the single product–single market model of Herrmann received 40%, the inshore sector 50%, and the mothership sector et al. (1996) to a represent post-AFA demand and allocation of 10%. The AFA also expanded the Western Alaska Community fillets, roe, and surimi to markets in Japan, Europe, and the USA. Development Quota (CDQ) program from an allocation of 7.5% In contrast to the cointegration approach applied in Fell (2008) of the eastern Bering Sea pollock TAC to 10% of the TAC for to identify structural breaks before and after the implementation each eastern Bering Sea groundfish species. The six CDQ enti- of AFA, and in contrast to the revenue function approach used ties represent 65 western Alaska communities and were created by Morrison Paul et al. (2009) to explicate the joint effects of to provide opportunities for community residents to benefit from product mix and productivity that arose from implementation of the wealth being derived from exploitation of eastern Bering the AFA, this study uses a structural model of post-AFA markets Sea fishery resources (Ginter 1995). These groups have, until to account for the effects of variability in landings and changes recently, leased their quota to AFA-qualified fishing vessels op- product demand on the mix of products produced, patterns of erating in the at-sea and inshore sectors. Under the AFA Alaska trade, and first wholesale prices and revenues. Simulations with pollock operations have enjoyed increased operational stability pre- and post-AFA product mixes, recovery rates, and markets and increased revenues (Felthoven 2002; Fell 2008; Wilen and based on the estimated model are used to explore the effects Richardson 2008; Morrison Paul et al. 2009). The AFA enabled on revenue of moving from a race-for-fish fishery to operating processors to look beyond simply processing the largest number under a cooperative structure. of fish before the end of the fishing season and to instead focus on maximizing revenues through increased capital expenditures into technology, expanded global markets, and optimized prod- METHODS Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 uct mix. The Alaska pollock fishery is the largest U.S. fishery in vol- Model Specification ume and among the top five in value. Despite its relative impor- The five most important markets for Alaska pollock prod- tance, markets for Alaska pollock products have received very ucts were modeled simultaneously: fillet sales to the U.S. and little attention in formal models. While Hiatt et al. (2010) pro- Europe, surimi shipped to Japan and U.S., and roe wholesaled vides an excellent descriptive overview of Alaska pollock prod- to Japan. The first four of these markets were represented by ucts and their markets, it does not model supply-and-demand inverse demand equations representing amount of the product relationships that drive observed patterns of production or trade and its respective price in that market and allocation functions flows and is thus ill suited for representing future production or that distributed quantities of product among competing markets. value. Using data from 1986 to 1993, Herrmann et al. (1996) de- The roe market was also represented by an inverse demand func- veloped a simultaneous equation model of the Japanese demand tion. Although U.S. export statistics suggest that 34% of Alaska for U.S. surimi, the U.S. supply of surimi to Japan, inventory pollock roe is exported to Korea, Japanese import statistics in- holdings of surimi in Japan and the United States, and U.S. ex- dicate that 85% of that roe is re-exported to Japan. Thus, over vessel price for Alaska pollock. They found that although surimi 95% of U.S. Alaska pollock roe production (2000–2008) was exports to Japan were unresponsive to changes in export price, ultimately sold into the Japanese market. The demand equations ex-vessel demand was inelastic over the biologically plausible were specified as inverse (price dependent) demands because, 1080 STRONG AND CRIDDLE
in a typical year, fishers harvest the entire TAC, which is ex- Market Clearing Identities ogenously determined, and global prices adjust to the quantity Market clearing identities are equations that transform input produced. Each demand equation includes per capita imports, data series into appropriate forms for inclusion in the behavioral substitute prices, and per capita income. The producer price equations. These adjustments account for changes in price due index for intermediate foods and feeds was used to normalize to inflation and fluctuating exchange rates. Quantities supplied prices to 2008 levels. The equation system was modeled using were not adjusted to account for changes in human populations monthly data for 2000–2008. due to the relatively short time span modeled. The first 12 market The supply of Alaska pollock of U.S. origin was treated as clearing identities included in the model are used to obtain time a stochastic process moderated by variations in survival and series for the endogenous variables (Table 3). The 13th identity growth of Alaska pollock, various environmental factors, and adjusts prices to 2008 price levels to eliminate the influence of changes in harvests. The upper bound on U.S. harvests (the inflation. The final six identities describe the time series used for TAC) is set by the National Marine Fisheries Service (NMFS) Alaska pollock fillets and surimi allocated to the United States. on behalf of the U.S. Secretary of Commerce and based on recommendations of the North Pacific Fishery Management Allocation Equations Council (NPFMC). The NPFMC, in turn, is advised by a Sci- Allocation of Alaska pollock fillets to the United States.—The entific and Statistical Committee invested with authority (under United States and Europe are the two largest markets for fillets the Magnuson–Stevens Fishery Conservation and Management produced from U.S. harvests of Alaska pollock. Nearly 44% of Act) to set binding upper limits for the acceptable biological all fillets produced from 2000 to 2008 were retained for sale in catch and the overfishing level, both of which serve as upper the U.S. market. Equation (1) represents the quantity of these fillet US bounds on the TAC. In practice, the TAC is often well below fillets, Qt , as a linear function of monthly dummy vari- the overfishing level because the Bering Sea and Aleutian Is- ables, mi ; the U.S. real price for single frozen Alaska pollock US land groundfish fisheries are subject to an aggregate cap of 2 fillet fillet blocks, Pt ; the real price of U.S. Alaska pollock fillet million metric tons. The Scientific and Statistical Committee US→EU fillet acts as a peer-review panel for stock assessments prepared by exports to Europe, Pt ; the quantity of Alaska pollock fillet SS NMFS staff who develop models that incorporate information fillets produced by the inshore processing sector, Qt ;the from the fishery and fishery-independent surveys to estimate quantity of Alaska pollock fillets produced by the at-sea pro- fillet AS and project recruitment. While the fishery occasionally harvests cessing sector, Qt ; the lagged real U.S. import price for IM→US less than the TAC, landings are generally within a few percent- fillet Alaska pollock fillets, Pt−1 ; and the lagged quantity of ages of the TAC; thus, the supply of Alaska pollock is not well fillet US Alaska pollock fillets allocated to the USA, Q − : represented by a traditional supply function. Instead, the dispo- t 1 sition of surimi to the USA and Japan and the disposition of 11 US fillets to the USA and Europe were represented by allocation fillet US = β + β + β fillet Qt 10 1i mi 112 Pt equations that each considered own price, prices in compet- i=1 ing markets, levels of production in the at-sea processing and US→EU + β fillet + β fillet SS shore-based processing sectors, and monthly seasonal variables. 113 Pt 114 Qt Additional variables, such as lagged quantities and competing IM→US product prices, were included in some of the allocation equa- + β fillet AS + β fillet 115 Qt 116 Pt−1 tions to account for variability due to other factors. Linear and fillet US Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 + β . nonlinear model forms were explored for each of the behavioral 117 Qt−1 (1) equations. The ultimate choice of functional form was based on goodness of fit for the individual equations and on performance The coefficients βni represent unknown parameters to be esti- of the overall equation system. Data sources used in the model mated, where the subscript associates the coefficient with the are listed in Table 1. Import and export statistics for the USA ith explanatory variable in the nth equation. Because there are were obtained from the U.S. International Trade Commission, no public data on the monthly quantities of U.S.-origin Alaska which maintains records of export and import quantities and pollock fillets sold within the United States, a proxy time series prices. Japanese imports and exports of roe and surimi, as well was constructed based on publically available data and advice as Japanese consumption statistics and welfare measures, were from industry representatives. These representatives indicated obtained from Trade Statistics of Japan. European import and that most pin-bone-out and pin-bone-in fillets are sold within export statistics were obtained from the Eurostat database. Pro- 90 d and that deep-skin fillets are nearly all sold within 180 d. duction data for the at-sea and inshore sectors of the U.S. Alaska A complete description of the identity (number 19) is listed in pollock fishery were derived from the NMFS catch account- Table 3. ing system (P. Britza, NMFS/Alaska Regional Office, personal While some U.S.-origin Alaska pollock are headed and gut- communication). Variable definitions and sources are listed in ted, reprocessed in China, and then reimported into the USA, Table 2. a majority of Alaska pollock fillets imported into the USA are ALASKA POLLOCK MARKETS 1081
TABLE 1. Data sources, frequencies, and ranges for information used in the market model of eastern Bering Sea Alaska pollock fisheries.
Number Type Description Source Frequency and range 1 Japanese roe consumption Total value and quantity Trade statistics of Japan Monthly, 1999–2008 purchased per household. 2 Roe auction prices Seattle roe auction prices. BANR reports Intermittent spring and fall, 1999–2009 3 Japanese surimi and roe Amount in inventory. Trade statistics of Japan Monthly, 1992–2009 cold-storage holdings 4 Roe, fillet, and surimi Quantity produced by. NOAA, personal Monthly, 1998–2008 production inshore and at-sea sector. communication 5 Japanese roe import quantity Average price and quantity of Trade statistics of Japan Monthly, 1999–2008 and price imports from Russia, USA, and Korea. 6 Japanese surimi import Average price and quantity of Trade statistics of Japan Monthly, 2000–2008 quantity and price imports from all countries. 7 U.S. surimi import and export Average price and quantity U.S. International Trade Monthly, 1992–2009 quantity and price for individual and total Commission countries. 8 U.S. fillet price 66-lb block, US$/kg Urner Barry Monthly, 1998–2008 9 U.S. fillet export and import Average price and quantity U.S. International Trade Monthly, 1992–2009 quantity and price for individual and total Commission countries. 10 European fillet and surimi Average price and quantity of Eurostat Monthly, 1999–2009 import quantity and price imports from USA, Russia, China, and Japan. 11 U.S. per capita income Real terms U.S. Bureau of Labor Monthly, 1999–2008 statistics 12 Japanese index of real Real terms Trade statistics of Japan Monthly, 1999–2008 consumption expenditures 13 German per capita Income Real terms Eurostat Monthly, 1999–2008 14 U.S. halibut fillet prices 10–20 lb, US$/kg Urner Barry Monthly, 1998–2008 15 U.S. producer price index Intermediate foods and feeds U.S. Bureau of Labor Monthly, 1998–2008 statistics 16 Exchange rates U.S. dollar, yen, and euro Board of Governors of the Monthly, 1998–2008 U.S. Federal Reserve
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 fillet US→EU twice-frozen fillets of Russian origin processed in China. There- Alaska pollock fillet exports to Europe, Qt , as a linear fore the lagged real price of Alaska pollock imports to the function of monthly dummy variables, mi ; the real price of U.S.- USA was included in equation (1) to represent the impact of US→EU origin Alaska pollock fillet exports to Europe, fillet P ;the Russian-origin Alaska pollock fillet imports on sales of U.S.- t real U.S. price (from Urner Barry reports) for single, frozen, origin Alaska pollock fillets in the United States. The lagged US fillet quantity of U.S.-origin Alaska pollock fillets allocated to the Alaska pollock fillet blocks, Pt ; the lagged quantity of fillet SS USA represents the effect that past fillet sales to the United fillets produced by the inshore sector, Qt−1; and the lagged fillet AS States have on current sales through, for example, establish- quantity of fillets produced by the at-sea sector, Qt−1: ment of long-term contracts with major purchasers. 11 Allocation of Alaska pollock fillets to Europe.—The sales US→EU fillet US→EU = β + β + β fillet volume of U.S.-origin Alaska pollock fillets exported to Europe Qt 20 2i mi 212 Pt exceeded the volume of Alaska pollock fillets retained in the i=1 United States for the first time in 2002 and has frequently done US +β fillet + β fillet SS so since 2004. Overall, from 2000 to 2008, more than 49% 213 Pt 214 Qt−1 of U.S.-origin Alaska pollock fillets were sold into European + β fillet AS . markets. Equation (2) represents the quantity of U.S.-origin 215 Qt−1 (2) 1082 STRONG AND CRIDDLE
TABLE 2. Variable definitions and sources of information used in the market model of eastern Bering Sea Alaska pollock fisheries. Source numbers refer to those listed in Table 1.
Variable Definition Source
fillet US Pt U.S. Urner Barry price of single frozen Alaska pollock fillet block. 8 fillet US→EU Pt Price of U.S.-origin Alaska pollock fillet exports to Europe. 9 fillet SS Qt Alaska pollock fillets produced by the inshore sector. 4 fillet AS Qt Alaska pollock fillets produced by the at-sea sector. 4 fillet IM→US Pt U.S. import price of Alaska pollock fillets. 9 fillet US→EU Qt U.S.-origin Alaska pollock fillets exported to Europe. 9 fillet totalEx Qt Total exports of U.S.-origin Alaska pollock fillets. 7 fillet total Qt Total quantity of Alaska pollock fillets produced. 4 surimi US→JN Qt U.S.-origin Alaska pollock surimi imported to Japan. 6 surimi US→JN P(JP )t Price of U.S.-origin Alaska pollock surimi imported to Japan. 6 surimi SS Qt Alaska pollock surimi produced in the inshore sector. 4 surimi AS Qt Alaska pollock surimi produced in the at-sea sector. 4 surimi US Pt U.S. surimi import price. 7 surimi totalEx Qt Total exports of U.S.-origin Alaska pollock surimi. 7 surimi total Qt Total quantity of U.S.-origin Alaska pollock surimi. 4 I1 and I2 Indicator variables for missing observations of U.S.-origin Alaska pollock surimi prices. US Inct Real U.S. per capita disposable personal income. 11 fillet € CR→EU P(EU )t European import price of Alaska pollock fillets from China and Russia. 10 € GE Inc( EU )t Index of disposable income of Germany in euros. 13 €→$ Ext Europe–U.S. exchange rate: one euro to U.S. dollars. 16 →$ Ext Japan–U.S. exchange rate: one U.S. dollar to Japanese yen. 16 roe US→JN P( JP )t Japanese import price of U.S.-origin Alaska pollock roe. 5 roe SS Qt Alaska pollock roe produced in the inshore sector. 4 roe AS Qt Alaska pollock roe produced in the at-sea sector. 4 roe CR→JN P ( JP )t Japanese import price of Alaska pollock roe from China and Russia. 5 JN Inct Japanese index of real consumption expenditures. 12 roe JN St Japanese cold-storage holdings of Alaska pollock roe. 3 I3 Indicator variable for missing observations of Alaska pollock roe prices in Japan. surimi JN St Japanese cold-storage holdings of Alaska pollock surimi. 3 surimi other→JN P( JP )t Japanese import price of surimi other than U.S.-origin Alaska pollock surimi. 6 surimi € US→EU P ( EU )t European import price of U.S.-origin Alaska pollock surimi. 10 halibut US Pt Price of halibut in the USA. 14 I5and I6 Indicator variable for data point with no observations on U.S.-origin Alaska
Downloaded by [Department Of Fisheries] at 19:14 20 July 2015 pollock surimi prices. roe US→JN Qt Japan imports of Alaska pollock roe. 5 US = PPIt U.S. producer price index for intermediate foods and feeds (base 1982). 15
Lagged quantities of fillets produced by the inshore and at-sea surimi SS surimi produced in the inshore sector, Qt−2; the two-period sectors were used to represent the difference in time between lagged quantity of Alaska pollock surimi produced in the at-sea fillet production and their distribution to Europe. surimi AS sector, Qt−2; and the real U.S. price (from Urner Barry Allocation of Alaska pollock surimi to Japan.—The primary US fillet market for Alaska pollock surimi is Japan, which purchased reports) for single, frozen, Alaska pollock fillet blocks, Pt : over 63% of total U.S.-origin Alaska pollock surimi from 2000 11 US→JN to 2008. Equation (3) models the quantity of Alaska pollock → surimi QUS JN = β + β m + β surimi P surimi imported from the USA by Japan, surimi QUS→JN,as t 30 3i i 312 t t i=1 a linear function of monthly dummy variables, mi ; the real +β surimi Q SS + β surimi Q AS price of Alaska pollock surimi imports by Japan from the USA, 313 t−2 314 t−2 US→JN US surimi + β fillet . Pt ; the two-period lagged quantity of Alaska pollock 315 Pt (3) ALASKA POLLOCK MARKETS 1083
TABLE 3. Market-clearing identities in the market model of eastern Bering Sea Alaska pollock fisheries.
Identity Equation