Original scientific article UDK 594.1: 591.134.2(262.3) Received: 2013-04-03 DEFINITION OF A NEW FORMULA FOR THE CALCULATION OF THE TOTAL HEIGHT OF THE FAN SHELL PINNA NOBILIS IN THE MIRAMARE MARINE PROTECTED AREA (TRIESTE, ITALY) Milena TEMPESTA Miramare Marine Protected Area, I-34151 Trieste, Viale Miramare 349, Italy E-mail: milena@riservamarinamiramare.it Donatella DEL PIERO University of Trieste, Department of Life Science, I-34127 Trieste, via L. Giorgieri 10, Italy Saul CIRIACO Shoreline coop., Area Science Park, I-34149 Trieste, Padriciano 99, Italy ABSTRACT Four years of monitoring of Pinna nobilis population in Miramare Marine Protected Area (Trieste, Italy) have led to the development of a new equation for indirect calculation of total shell height (H ) for this bivalve, endemic to the Mediterranean Sea. As the formulas listed in the cited literature are not corroborated by in situ measured growth data, it was necessary to describe a specific formula for Miramare population. This exponential equation was called MirExp. A specific growth curve has also been created with this formula and the data collected every year from 2008 to 2011. This curve shows that in Miramare P. nobilis grows faster in its early years than what is reported in other studies carried out in Croatia and Spain. Population density is also the highest among those reported in consulted papers, which make reference to studies conducted in various Mediterranean protected and non-protected areas. Key words: Pinna nobilis, population dynamics, growth, marine protected area, Adriatic Sea DEFINIZIONE DI UNA NUOVA FORMULA PER IL CALCOLO DELLA LUNGHEZZA TOTALE DEGLI ESEMPLARI DI PINNA NOBILIS PRESENTI NELLA AMP DI MIRAMARE (TRIESTE, ITALIA) SINTESI Dopo 4 anni di monitoraggio del popolamento di Pinna nobilis all'interno dell'Area Marina Protetta di Miramare (Trieste) si è giunti alla definizione di una nuova equazione per il calcolo indiretto della lunghezza totale (H) della conchiglia di questo bivalve endemico del Mediterraneo. Le formule riportate in bibliografia non trova-no riscontro nei dati di accrescimento misurati in situ e quindi è stato necessario descrivere una formula specifica per il popolamento di Miramare, denominata MirExp trattandosi di una equazione esponenziale. Grazie a questa formula e ai dati raccolti annualmente dal 2008 al 2011, si è definita anche una curva di accrescimento specifica che ha evidenziato come questo bivalve si accresca più velocemente nei primi anni a Miramare rispetto a quanto riportato in altri studi svolti in Croazia e Spagna. La densità di popolamento, inoltre, è la più alta tra quelle ripor-tate nei lavori consultati effettuati in varie aree del Mediterraneo, protette e non. Parole chiave: Pinna nobilis, dinamica di popolazione, accrescimento, area marina protetta, Adriatico INTRODUCTION Pinna nobilis L., 1758 is an organism of the phylum Mollusca, class Bivalvia and subclass Pteriomorphia, endemic to the Mediterranean Sea, where it is present since the end of the Miocene (Gómez-Alba, 1988). It is one of the largest living bivalves and among these, like other Pinnids, shows the fastest growth of the shell (Richardson et al., 2004): specimens with a length of 120 cm have been found (Zavodnik et al., 1991; Richardson et al, 1999; Garcia-March et al, 2007). The fan shell live partially buried upright in the sandy and pebbly substrata particularly in association with marine phanerogams such as Posidonia oceanica, Cymodocea nodosa, Zo-stera marina e Nanozostera noltii (Zavodnik et al, 1991). Despite the interest of the scientific community for this species and its conservation, knowledge about the ecology and population dynamics of P. nobilis are still fragmentary. Although in France and in Spain there are several studies on various aspects of this bivalve, such as its ecology (Combelles et al., 1986; De Gaulejac & Vicente, 1990; Butler et al., 1993), morphology, growth and age assessment (Moreteau & Vicente, 1980; García--March & Ferrer, 1995, Richardson et al., 1999; Hendriks et al., 2012), gametogenesis and reproduction (De Gaulejac, 1995; De Gaulejac et al., 1995a, b), there are only sporadic publications on the ecology and morphology of this species in the Adriatic Sea (Mihailinovic, 1955; Zavodnik, 1967; Zavodnik et al, 1991; Peharda et al, 2002; Siletic & Peharda, 2003; Prestinenzi, 2009). The importance of studying this bivalve mollusc is shown, inter alia, by its protection status under official documents such as the Barcelona Convention (1995), ratified by the Italian government with Act no. 175 of 25/05/1999, and the EU Habitats Directive (43/92). In particular, P. nobilis is listed under Annex II of the Barcelona Convention including endangered species requiring a specific conservation strategy, and under Annex IV of the Habitats Directive on the conservation of natural and semi-natural habitats and of wild flora and fauna, where species requiring strict protection measures are listed. Although protective measures have been introduced already in the early nineties, there is still uncertainty on the conservation status of P. nobilis populations at Mediterranean scale (Addis et al., 2009). The amateur collection of shells, the destruction due to anchoring and fishing both commercial and recreational, but especially the regression of the seagrass meadows as its elective habitat, have contributed to the decline of the population of P. nobilis (Zavodnik et al., 1991; Centoducati et al., 2007; Katsanevakis, 2007b). On the contrary, an increase of the number of individuals of this species has been reported in the Gulf of Trieste (Northern Adriatic Sea) in recent years (Prestinenzi, 2009), particularly in Mirama-re Marine Protected Area as for the juveniles is concerned. Therefore marine reserves - which, regardless of their size, are associated to high values of density, biomass, individual size and specific diversity (Halpern, 2003) - can be considered essential to protect a breeding stock that can increase marine living resources through spill-over (Tewfik & Béné, 2003). Estimates of P. nobilis recruitment and spill-over rates on populations in small-size MPAs (Mljet National Park, Croatia) have already been made by using models (Peharda & Vilibic, 2008), which however need further refinement through in situ experiments and systematic monitoring. To better understand this new positive trend in the presence of P. nobilis, a monitoring was carried out on a non-stop basis in Miramare MPA for 4 years, from 2008 to 2011. Its aim was to assess the population structure and dynamics inside the Reserve by counting and measuring the size of individuals, estimating their density and studying their spatial distribution. Data collected have been compared over four years thanks to the use of tagging techniques, through which monitored individuals were identified. Measurements were made in situ by SCUBA divers as this is a method of non-destructive sampling and therefore an excellent tool for the analysis of species or habitats protected or rare (Katsanevakis, 2007a). MATERIALS AND METHODS The sampling methodology follows the indications provided by the Protocol to study and monitor Pinna no-bilis populations within marine protected areas written by MEPA (Malta Environment and Planning Authority) in the framework of the MedPAN project (García-March & Vicente, 2006). Data were collected through measurements made by SCUBA divers. In Miramare, a round area of nearly 5.6 m radius (approximately 100 m2) was defined: the area was sampled by two SCUBA divers swimming from its central point - whose GPS coordinates were recorded - outwards and stopping whenever they saw a P. nobilis individual. Its distance from the centre was marked by using a tape, degrees to the north were noted down and a plastic numerical tag was put at the base of the shell bonded with a nylon cable tie. The height of the unburied shell (H), maximum width at the point of maximum dorso-ventral length of the shell (L ) and minimum width at the bottom (L . ) were then max min measured. Morphometric measurements were performed by using a 0.1 cm tape, whereas the total anteroposterior shell height (H ) was estimated by using equations to relate the measurements of unburied shell height and total height. During the dives unbroken shells of dead individuals were collected to validate total height data resulting from equations. Double sampling In 2011, two monitoring actions were carried out with the help of the Coast Guard Divers Team: P. nobilis specimens were surveyed both in May by coast guard divers and in August by Miramare biologists. Each individual was therefore counted and measured twice in a period of three months, and the same morphometric measurements were performed. The comparison between the two sampling periods showed a difference in the measured minimum width at the bottom (L . ), that min ' leads to a statistically significant difference (Mann-Whitney test U = 964, z = -1.974, p(same)< 0,05). Statistical analysis was performed using the open source data analysis package PAST ver. 2.17 (Hammer et al., 2001). Lmin is necessary to calculate total shell height (Htot) with García-March & Ferrer equation (García-March & Ferrer, 1995), which was chosen in the past as its results were the closest to the measurements made on dead individuals with unbroken shells collected during the dives. Miramare's new formula for H , calculation tot Hence a different formula was created to calculate H, which could better adapt to Miramare's measurement method. So, the new formula derived directly from the measurements on Miramare's specimens. In order to avoid any mistakes due to the difficulty of measuring minimum width at the bottom, the new formula considered parameters that can be better measured such as the height of the unburied shell (H) and maximum width at the point of maximum dorso-ventral length of the shell (L ), thus neglecting L . . The new formula was max ' o o min calculated following the indications given by Katsane-vakis (2006), who obtained the formula to convert the size-at-age and growth-rate results into total shell height from the relationship between H , and L of some tot max dead individuals with unbroken shell. The formula was based on the data of 10 individuals with unbroken shells and no concretions that were col- lected in Miramare during the samplings in 2007, 2008 and 2010. Total shell height and maximum width were measured by using a 0.1 cm tape and data were put in a graph with H , values on the y-axis and L values o r tot / max on the x-axis. Probability measures (r) showed that the equation best describing the relationship between H -Lmax is the exponential equation (MirExp) given below: (MirExp) y = 10.259 e0 R2 = 0.9118 where y is H value in cm and x (to be inserted) is L value in cm, both measured in situ. max From the specimen measured during the four year sampling period, a growth curve was elaborated according to the Siletic & Peharda (2003) methodology applied to a study of P. nobilis in the Mljet National Park (Croatia). Total shell height data used were so derived from MirExp resulting values. Siletic & Peharda (2003) estimated P. nobilis growth parameters by using von Bertalanffy equation: L = L -k( t-tl 0)) where Lt is H at a given time, L^ is the x-intercept on Gulland-Holt plot, k is the regression line slope on Gulland-Holt plot and t-t0 is the time since the first measurement made in t0 (Gulland & Holt, 1959). The Gulland-Holt plot was constructed by using the data resulting from the measurements performed on the 24 tagged shells over two years (2008 and 2011). In the plot, the growth rate expressed in cm/year (i.e., the difference between the values measured over the two years) is shown in the y axis, whereas the average total shell height obtained from the mean value of H of each individual is shown in the x axis. from 0 to 20 cm from 20 to 40 cm from 40 to 60 cm i □ 2008 □ 2009 □ 2010 □ 2011 Fig. 1: Size frequency distribution of Pinna nobilis over four years of sampling (2008,2009,2010 and 2011) Sl. 1: Velikostna porazdelitev leščurjev v obdobju 2008-2011 Tab. 1: Total height (Htot) in cm measured in 2008, 2009,2010 and 2011 on 24 individuals tagged in 2008 and re-measured every year Tab. 1: Celotna višina (Htot) v cm, izmerjena leta 2008, 2009, 2010 in 2011 pri 24 primerkih, označenih leta 2008 in na novo izmerjenih vsako leto No. Tag No. H. „2008 H. „2009 H. „ 2010 H. „ 2011 1 17 29.4 35.9 37.4 43.3 2 15 34.5 40.6 44.0 46.9 3 14 34.5 39.0 37.4 39.3 4 2 37.4 39.0 44.0 51.7 5 4 19.6 25.0 34.5 42.3 6 21 25.0 29.4 34.5 37.4 7 5 15.4 24.0 30.6 35.9 8 7 15.4 25.0 31.8 34.5 9 29 15.4 24.0 31.8 35.9 10 20 14.2 23.0 29.4 37.4 11 32 27.1 34.5 37.4 37.4 12 18 29.4 34.5 34.5 39.0 13 67 29.4 33.2 37.4 38.7 14 74 29.4 31.8 34.5 34.0 15 58 27.1 31.8 34.5 36.8 16 65 29.4 35.9 37.4 40.6 17 39 18.1 28.2 34.5 42.6 18 230 16.0 23.0 30.6 34.5 19 11 16.7 27.1 31.8 40.6 20 100 18.1 29.4 34.5 41.9 21 61 19.6 37.4 39.0 43.3 22 97 40.6 47.7 37.4 44.0 23 175 21.2 31.8 40.6 42.3 24 189 19.6 25.0 36.2 44.0 RESULTS Over four years, a total of 323 individuals were tagged counted and measured (91 in 2008, 79 in 2009, 77 in 2010 and 76 in 2011). The 26 % of the 2008 tagged pool, 24 specimens, maintain the numerical tag along the four years and so were re-measured year after year, thus providing an indication of their in situ growth. Figure 1 shows the values of average total shell height calculated with the new MirExp equation, with shells divided into 4 size classes of 20 cm each. According to literature, 20 cm of H is the threshold under which ' tot individuals should be considered juveniles (Combelles et al., 1986; Butler et al., 1993; Richardson et al, 1999). Classes are divided into 0-20 cm (juveniles), 20-40 cm (adult juveniles), 40-60 cm (adults) and more than 60 cm (old individuals). It is evident that in 2008 most individuals were juveniles, whereas in 2009 and 2010 they were adult juveniles and in 2011 adults. There are no individuals longer than 60 cm, which at present is the maximum height for specimens found in Miramare. The comparison between four sampling years shows a general increase in H values of individuals from one year to the next, as well as a different size distribution as juveniles were only present in 2008 and 2011, whereas recruitment and the birth of new shells almost did not occur in 2009 and 2010. As a whole, a shift can be seen if size class frequencies for all sampled specimens are considered: in 2008, 59 % of the population consisted of individuals smaller than 20 cm in total height; in 2009 and 2010, 94 % and 91 % respectively consisted of shells ranging from 20 to 40 cm, which became only 58 % in 2011. Population density distribution was marked by fluctuations over the years. The highest average density was recorded in 2008 with 20.84 individuals/100 m2, followed by 2009 with 18.88 ind./100 m2, 2010 with 15.69 ind./100 m2 and 2011 with 13.57 ind./100 m2. Growth As the sampling lasted 4 years, the growth curve was calculated by using only the specimens that were found and re-measured all four years, from 2008 to 2011. In particular, 2008 saw the presence of many juveniles that presumably settled in that year and, consequently, were likely to be at least age four in 2011. The morphometric data from the pool of the same 24 specimens were collected in situ in 2008, 2009, 2010 and 2011 (Tab. 1) and were then inserted in Miramare's exponential formula (MirExp) to calculate H The resulting total shell height values for 2008 and 2011 were used to obtain a growth curve, according to the method used by Siletic & Peharda (2003) in their population study of P. nobilis in Mljet National Park (Croatia). The Gul-land-Holt plot created with the Miramare's shells data on 2008 and 2011 (Fig. 2), shows that the point where the line crosses the x axis corresponds to the L^ value (52.6 cm) and the line slope corresponds to a k value of 0.3408. Von Bertalanffy equation - adjusted to the measurements performed on P. nobilis individuals found in Miramare MPA - was therefore as follows: Lt = 52.6 (1-e-°-3408