Brachiopods , sipunculans , enteropneusts and metals from two estuarine tidal flats , Pacific , Costa Rica

Reports on the abundances and on metal concentrations in intertidal estuarine invertebrates from the Eastern Tropical Pacific are rare. Thus, the objectives of this report are to make accessible data on the abundances (1984-1987, 49 dates; 2013, 12 dates) of sipunculans, brachiopods and hemichordates from a sand-mud flat; and on trace metals (1996, 2000) and abundances (2015, 3 dates) of sipunculans and brachiopods at a sand flat in the Gulf of Nicoya estuary (10o N-85o W). Cores (17.7 cm2) were collected at the sand-mud flat, and quadrats (0.2 m2) at the sand flat. The flats contrasted in their sand (65 % vs 90 %) and silt+clay (31.5 % vs 5.6 %) contents. At the sand-mud flat (1984-87: 1.83 m2) the sipunculans were represented by 13 individuals, the brachiopods by 129 and the acorn worms by 185, with estimated maximum densities of: 5.7, 29, and 40 ind./m2, respectively. Trace metal (Fe, Mn, Ni, Cr, Cd, Zn, and Pb) analysis (Atomic Absorption Spectrometry) were conducted in specimens of Sipunculus nudus (Sipuncula) and Glottidia audebarti (Brachiopoda). Maximum mean concentrations in S. nudus were: For non-depurated worms, Fe (16.0 mg/g dw) > Mn (165 μg/g dw) > Zn (81 μg/g dw) > Cu (26 μg/g dw) > Cr (11 μg/g dw) > Ni (10.4 μg/g dw) > Pb (9.3 μg/g dw) > Cd (1.2 μg/g dw). For 72 hour depurated worms: Fe (5.0 mg/g dw) > Mn (61 μg/g dw) > Zn (39 μg/g dw) > Cu (24 μg/g dw) > Ni (8.4 μg/g dw) > Pb (2.7 μg/g dw) > Cd (0.62 μg/g dw). For G. audebarti: Fe (1.6 mg/g dw-soft parts) > Zn (123.5 μg/g dw-soft parts) > Cu (31.4 μg/g dw-pedicles) > Pb (21.0 μg/g dw-shells) > Cd (5.2 μg/g dw-soft parts) > Cr (4.7 μg/g dw-shells). For sediments; Fe (46 mg/g dw) > Mn (41.3 μg/g dw) > Zn (63 μg/g dw) > Cu (36.2 μg/g dw) > Cr (31.5 μg/g dw) > Pb (21.1 μg/g dw) > Ni (16.1 μg/g dw) > Cd (1.1 μg/g dw). These concentrations were expected for a nonindustrialized estuary. At the sand flat (Area sampled: 10.6 m2 ) 76 individuals of G. audebarti, 112 of G. albida, and 366 of S. nudus were collected in 2015, with estimated maximum densities of: 7.1, 10.5, and 31 ind./m2, respectively. Densities of G. audebarti and G. albida were relatively low, while those of S. nudus were relatively high when compared with other reports. The shell lenght of G. audebarti ranged from 9.0 mm to 38.0 mm and from 6.0 mm to 29.0 mm for G. albida. These ranges were within those found for these lingulides elsewhere. The mean length of S. nudus was 41 mm and the maximum weight was 1.6 g, which are small. No brachiopods were found at the sand-mud flat in 2013, nor enteropneusts at the sand flat in 2015. G. audebarti had a relatively stable presence, while G. albida almost vanished from the samples at the end of 2015. The spatial distributions of the three invertebrates were found aggregated at both intertidal flats. Strong ENSO warming events during 1983 and 2015, and red tides in 1985 may have influenced the abundances. Rev. Biol. Trop. 64 (3): 1311-1331. Epub 2016 September 01.

At the nearby Cocorocas sand flat, a preliminary evaluation of the presence of metals in S. nudus and G. audebarti was presented by Vargas and Abdullah (1997) and the final results are included herein.A comparison of the Cocorocas macrofauna with that of similar flats in other latitudes was made by Dittmann and Vargas (2001) and by Dittmann (2002).At Cocorocas the two species of brachiopods and S. nudus coexist and are important components of the macrofauna.
The presence of imposex in Acanthais brevidentata snails from the Gulf was informed by Gravel, Johanning, McLachlan, Vargas, and Oberdörster (2006) and concentrations of PCBs in S. nudus and other sipunculans were reported by Spongberg (2006).The abundance and spatial distribution of the giant onuphid polychaete worm Americonuphis reesei in sand flats from the upper Gulf of Nicoya was described by Rojas and Vargas (2008).The presence of Pharmaceuticals and Personal Care Products (PCCP) in rivers draining into the Gulf of Nicoya and other coastal areas was evaluated by Spongberg et al. (2011).More recently, Vargas, Acuña-González, Gómez and Molina (2015) summarized data on trace metals from snails and clams collected at the four coastal sites, including the razor clam Tagelus affinis from the Cocorocas sand flat.
No recent published data is available on the abundances of brachiopods and sipunculans from the Gulf of Nicoya estuary.Reports are also rare on the abundances of these groups and acorn worms for the Eastern Tropical Pacific region.Data on trace metals in sipunculans and brachiopods are scarce worldwide.Therefore, concentrations of trace metals in S. nudus and G. audebarti from the 2000 survey are also included herein.This information may be crucial in studies assessing the potential recovery of a disturbed system like the Gulf of Nicoya estuary (Vargas, 2016).This assessment needs information on the structure (diversity and abundance) and function (energy flow) of its benthic fauna, and the presence-absence of certain groups of invertebrates may be key diagnostic factors as the review by Borja, Dauer and Elliott (2010) indicates.
Thus, the objectives of this study were to make accessible data on the abundances of sipunculans, brachiopods and hemichordates (1984)(1985)(1986)(1987)2013) at a sand-mud flat; and on trace metals (1996,2000) and abundances (2015) of sipunculans and brachiopods at a sand flat in the upper Gulf of Nicoya estuary.

MATERIAL AND METHODS
Collection sites: Samples were collected at two intertidal sedimentary environments located on the Southern shore of the Punta Morales Peninsula, upper Gulf of Nicoya estuary (10 o N -85 o W), Pacific coast of Costa Rica (Fig. 1A) during 1984-1987and 2013(Punta Morales sand-mud flat) and 1996, 2000, and 2015 (Cocorocas sand flat).
The Punta Morales flat is characterized by soft sediments.A person walking on them quickly sinks to knee high (Fig. 1B).According to Vargas (1987Vargas ( , 1988Vargas ( , 1996) ) conspicuous biogenic structures are tubes of the onuphid polychaetes Diopatra ornata and Americonuphis reesei (Fig. 1B) and fecal mounds of hemichordates.Snails like the scavenger Nassarius luteostoma and the predator Natica unifasciata crawl on the sediment surface year around, while the sand dollar Encope stokessi is more frequent during the dry seasons.At low tide the sediments remain wet (Fig. 1B) and small tide pools serve as refuge for crabs, shrimps, and gobiid fishes.Several species of shore birds feed at low tide on this fauna.
The Cocorocas sand flat is characterized by dark-gray sandy sediments allowing a person to walk and stand on them without sinking (Fig. 1C).The flat is under the influence of freshwater discharges from the nearby Lagarto River that carries heavy loads of sediments during the rainy season, and coarse sediments are deposited on the flat.Ripple marks are observed at low tide as well as occasional ray feeding pits.Tubes of onuphid worms are scarce and are mainly of D. ornata.Hemichordate fecal mounds are absent.Snails and sand dollars are also rare on the sediment surface.At low tide, the sediment surface gets almost dry (Fig. 1C) and ray feeding pits act as refuges for small fish and crustaceans.Shore birds were rarely seen on this flat.At both the Punta Morales and Cocorocas flats the infauna includes the brittle star Ophiopholis geminata.At Punta Morales the small razor clam Tagelus bourgeoisae is frequent, while at Cocorocas the larger  In this upper region of the Gulf (Fig. 1A) the average tidal range is 2.3 m and tides are semi-diurnal.The estuary is under the influence of seasonal rainfall: A rainy season occurs from May to November, followed by a dry season from December to April.In the upper estuary surface water salinity varies accordingly from about 24 psu (practical salinity units) during the rainy season, to near 34 psu during the dry season (Voorhis, Epifanio, Maurer, Dittel, & Vargas, 1983).Primary productivity data from the Gulf of Nicoya indicates that this estuary may have reached the hypertrophic condition (Cloern, Foster, & Kleckner, 2014).Runoff reaches the Gulf by several rivers of which the Tempisque River at the head and the Tárcoles River near the mouth (Fig. 1A) are the more important.Surface water temperatures are near 30 o C year around, while at low tide the exposed sediments may reach temperatures near 40 o C, particularly during the sunny dry season.A strong El Niño Southern Oscillation (ENSO) event influenced the Gulf of Nicoya region during 1982-1983, as well as during 2015.In 2015 the onset of the rainy season was delayed several months and rains ended in early November due to ENSO.

Sediment and specimen collections at the
Punta Morales sand-mud flat (1984)(1985)(1986)(1987)2013): The intertidal sampling protocols used at the Punta Morales flat are outlined in Vargas (1987Vargas ( , 1988)).In summary, two sets of 14 cores (17.7 cm 2 -15 cm deep) were collected at semimonthly intervals (February, 1984to February, 1985) and one set of 14 cores at near monthly intervals (March, 1985-April, 1987;and January -December, 2013) from muddy sands during low tide and within a 400 m 2 plot 20 m apart from a sandy beach (Fig. 1A).Data from Vargas (1989) was expanded with unpublished information from March and April of 1987, for a total of 49 sampling dates.Core samples were fixed in Rose Bengal stained formalin in sea water.Preserved cores were sieved thru a 500 micron mesh.Sorted specimens of brachiopods (Glottidia spp.), sipunculans (S. nudus, A. trichocephalus), and acorn worms (Fig. 2A) were kept in vials filled with 70 % ethanol.Sediment cores were also collected for determinations of grain size (sieve analysis according to Gray & Elliot, 2010) and organic matter content (by combustion at 450 o C), dried at 65 o C, and kept in sealed polyester bags until analysis.
Specimen collection at Cocoracas sand flat (1996,2000,2015): On June 10, 1996 individuals of S. nudus (Fig. 2B, Fig. 2C, Fig. 2D) and of the brachiopod G. audebarti (Fig. 2F) were collected at the Cocorocas flat along a 100 m long transect perpendicular to the shore and running between two islets (Fig. 1A, Fig. 1C).Sediment clumps were removed with a shovel to a depth of 20 cm, broken up by hand, and brachiopods and sipunculans placed in acid-washed polyester bags.At the laboratory the specimens were vacuum dryed for international shipment.A total of 40 individuals of S. nudus, and 20 of G. audebarti were selected to perform trace metal analyses by means of Atomic Absorption Spectrometry (AAS).Analyses were accomplished by M. Abdullah at the Department of Biology of the University of Oslo (Norway) on non-depurated whole S. nudus, and in tissues, pedicles, and shells of G. audearti.
On March 8th, 2000, sediment clumps were again removed with a shovel to a depth of 20 cm.A total of 100 specimens of S. nudus were collected at Cocorocas and transported to the laboratory in a cooler filled with sea water from the site.At the laboratory the worms were kept in cooled (20 o C) sea water for two hours, rinsed with filtered sea water, blotted dry on a paper towel, measured, weighed, and separated into two groups of 50 individuals each.One group was placed in a heat-sealable, acid-washed polyester bag and stored frozen (-18 ºC) until analysis.The second group was kept in a plastic container filled with filtered sea water from the site and aerated with portable plastic aquarium pumps to allow the worms to remove sediments and other matter from their digestive tracts (Fig. 2 D).Sea water, feces, and dead organisms were discarded every 12 hours and the container filled with fresh filtered sea water.The procedure was repeated six times for a total of 72 hours.Forty-two individuals were alive (able to move) at the end of the depuration period.The depurated worms were removed by hand, blotted dry on a paper towel, measured, weighed, and stored frozen (-18 ºC) as above until analysis.Trace metal analyses were made at the Research Centre on Marine Pollution (CICA) of the University of Costa Rica.Metal determinations (Fe, Mn, Ni, Zn, Pb, Cd) were performed following methods described by Vargas, Acuña-González, Gómez, and Molina (2015).In summary: Certified reference materials and blank tests were run.Fe, Mn and Zn were analyzed by Flame AAS (Perkin Elmer ® 3300).Cd, Ni and Pb were analyzed by Graphite Furnace AAS (Perkin Elmer ® HGA 600).The results for iron are expressed as mg/g dry weight (mg/g dw), and all other concentrations are expressed in parts per million / dry weight (μg/g dw).Sediment samples for trace metal analysis were collected using an acid-washed plastic corer, placed in heat-sealable, acid-washed polyester bags, and stored frozen (-18 ºC) until analysis.
On July 6, October 2, and December 15, 2015, peanut worms (S. nudus) and brachiopods (G.audebarti and G. albida) were collected at the Cocorocas sand flat to verify its presence and estimate their relative abundances.A wooden frame (50 x 40 cm) was pressed on the sediment surface to mark the sampling area (0.20 m 2 ).A minimum of 15 (3.0 m 2 ) and a maximum of 20 (4.0 m 2 ) samples were collected during low tide, along a 200 m transect running between two islets (Fig. 1C) and perpendicular to the shoreline.Each sample was collected 10 m apart from each other.Sampling was scheduled to coincide with very low tides near noon time.Sediment clumps were removed with a shovel to a depth of 20 cm, broken up by hand, and brachiopods and sipunculans placed in plastic bags with sea water from nearby tide pools.Brachiopod specimens from each quadrat were split in two groups: those with cream-white shells (G.albida, Fig. 2E) and those with dark-green coloration on the anterior region (G.audebarti, Fig. 2F).Brachiopods were blotted dry on paper towel and shell length (center of shell tip to pedicle attachment) measured to the nearest 0.5 mm.Final preservation of organisms was in 70 % ethanol.Specimens of S. nudus were also preserved in 70 % ethanol for 24 hours and quickly blotted dry on a paper towel and weighed.Samples for grain size analysis and organic matter content were collected by hand from sediment removed by the shovel, dried at 65 o C, and kept in sealed polyester bags as described before.

Sediment compositions:
The sediment compositions of the Punta Morales (1984,2013) sand-mud flat and of the Cocorocas sand flat (2015) are included in Table 1.The sand-mud flat was characterized in 1984 by 65.5 % sand (46.4 % fine sand, 15.1 % very fine sand), the silt+clay content was 31.5 % and organic matter was 2 %.A change towards finer sediments and organic enrichment was found during the 2013 survey: 41.3 % sand (5.6 fine sand, 33.0 % very fine sand), the silt+ clay content was 57.6 % and organic matter was 8 %.Coarse + very course sands were around 1 % in both surveys (Table 1 A, B).The Cocorocas sand flat samples contained higher percentages of sand (90.7 %) of which 29.2 % was fine sand.Coarse + very coarse sands were 10.7 %.The silt+clay content was only 5.6 % and organic matter was around 3 % (Table 1 C).
The 1984-1987 survey: Examples of the acorn worms and brachiopods collected are included in Fig. 2A, Fig. 2F.Organisms of the three groups had many on a size range that may represent juvenile forms.Acorn worms in particular included many organisms around 10 mm in length (Fig. 2A).The total numbers were: 13 sipunculans (12 ind. of S. nudus, 1 ind. of A. trichocephalus) 129 brachiopods and 185 hemichordates (Fig. 3).The total area sampled was 1.83 m 2 .
The abundances of these three groups at the Punta Morales sand-mud flat for the 49 dates of the 1984-1987 survey are included in Fig. 3. Sampling started in February 1984, shortly after the impact of the severe 1982-1983 ENSO high temperatures and dry conditions had returned to near normal values in the Central American region.Red tides were frequent in the upper Gulf of Nicoya in 1985, particularly after June (Fig. 3).The periods from the start of sampling to about June 1984, and from around May 1985 to October 1985 were characterized by low to zero abundances.In spite of the fact that after March 1985 the sampling effort was smaller in space (14 cores) and in time (monthly), abundances were higher during the rainy seasons.Abundance patterns of brachiopods and acorn worms in particular were similar since both groups were absent early on 1984, begun to increase and reached high numbers around the mid rainy season of 1984, became scarcer during 1985 and early 1986, and increased slightly during the dry season of 1986-1987.The presence of a few sipunculans was restricted to the rainy seasons of 1984 and 1985 (Fig. 3).

Metal concentrations (1996, 2000):
Examples of the peanut worm Sipunculus nudus are included in Fig. 2B, Fig. 2C, Fig. 2D.The concentrations of trace metals in non-depurated S. nudus, soft parts, pedicles, and shells of Glottidia audebarti, and sediments collected in 1996 are included in Table 2. S. nudus had higher mean concentrations of Fe and Cr than G. audebarti, while the brachiopod presented higher mean concentrations of Zn and Cd (soft parts), Pb (shells) and Cu (pedicles).Pb concentrations were similar in organisms and in sediments, while those of Zn and Cd were higher in organisms than in sediments.Mar.

May May
Jun. Jun. Jul. Jul. Aug.

May
Jun. Aug.

May
Jun. Jul. Aug.

Mar.
Apr   (1984-1987). Feb. 22, 1984to Feb. 21, 1985: Semi-monthly sampling of two sets of 14 cores per date.Numbers of acorn worms and brachiopods found on each 14 core set are indicated on the top of the bars.Sipunculans were found on one of the sets.Mar. 07, 1985to Apr. 29, 1987: Monthly sampling of 14 cores per date.Core area 17.7 cm 2 , core depth 15 cm.Mesh screen: 500 microns.
Mean concentrations of Fe, Cu, and Cr were higher in sediments.
The concentrations of metals in tissues of non-depurated and depurated whole S. nudus collected in 2000 are included in Table 3A.In non-depurated worms mean concentrations of Fe and Zn were higher in 2000 than in 1996, while Pb and Cd were lower in 2000; Cu was detected in similar concentrations for both surveys (Tables 2, 3A).When median concentrations of non-depurated and depurated S. nudus are compared, these decreased for all metals except Cd (Table 3A).
The weight on non-depurated S. nudus ranged from 0.65 g to 1.67 g, with an average of 0.95 g, while the weight of depurated worms ranged from 0.61 g to 1.35 g, with a mean of 0.84 g.The weight of depurated S. nudus was significantly lower after the depuration period (Table 3B).There was no statistically significant difference in length of non-depurated vs depurated organisms The mean length of S. nudus was 45 mm, with a maximum of 65 mm (Table 3B, Fig. 2

B.C).
The 2013 survey of the sand-mud flat: During the 2013 monthly core sampling at the Punta Morales sand-mud flat no specimens of brachiopods and of S. nudus were found.However, the Sipuncula included seven individuals of A. trichocephalus (Fig. 4).The acorn worms (Enteropneusta) were represented by 16 specimens found mainly during the rainy season (Fig. 4).Jun.24 Jul. 29 Aug. 9 Sep. 9 Oct. 10 Nov. 8 Dec. 5 The 2015 survey of the sand flat: Examples of the brachiopod G. albida collected are included in Fig. 2E.Samples were collected at the Cocorocas flat in July, October, and December of 2015, a year when high temperatures and dry conditions in the Eastern Pacific region were enhanced by the ENSO event (NOAA, 2016).The presence at the sand flat of brachiopods (G.albida and G. audebarti) and S. nudus worms was confirmed during sampling conducted in 2015.These three species were found frequently living together on a single quadrat area (0.2 m 2 ).Acorn worms, however, were not found.
A total of 76 individuals of G. audebarti, 112 of G. albida and 366 of S. nudus were  found in a pooled sample area of 10.6 m 2 (Fig. 5).The size range of G. audebarti was from 9.0 mm to 38.0 mm, while for G. albida the range was from 6.0 mm to 29.0 mm (Fig. 5).The weight range of S. nudus was from 0.10 g to 2.55 g (Fig. 6).Abundances (22 -23 -31 ind.) of G. audebarti were relatively stable.The smaller and the bigger specimens were both collected in December, and mean shell lengths (22.5 -21.0 -21.8 mm) were similar during the sampling period (Fig. 5 A, Fig. 5B, Fig. 5C).On the other hand, G. albida abundances (73 -35 -4 ind.) declined sharply during the period.The smaller specimens of G. albida were found during the three sampling dates, while the bigger individuals were found in July and October, and mean shell lengths increased and then decreased (12.5 -16.0 -9.7 mm), Fig. 5 D, Fig. 5E, Fig. 5F.The abundance of S. nudus was high towards the end of the rainy season.Individuals in the ranges below 1.0 g were more common along the sampling period, but the increase in numbers of those above 1.0 g was modest (Fig. 5 A, Fig. 5B, Fig. 5C).Two individuals of the larger (approx.200 mm) sipunculan, Xenosiphon branchiatus were also found in the Cocorocas samples.

DISCUSSION
The Brachiopoda and Hemichordata are considered two separate groups each with the category of Phylum, while the Sipuncula is presently grouped by some authors within the Phylum Annelida based on molecular evidence (Brusca, Moore, & Shuster, 2016).
The lamp shells, acorn worms and peanut worms include organisms that burrow in sediments of marine and estuarine habitats worldwide.The early information on these groups was summarized by Hyman (1959) who mentioned several examples of the scarce literature available at the time for tropical regions.More recent reviews were published by Cutler (1994), Emig (1997, 2009) Cameron (2005) and Vargas and Dean (2009).These works also indicate that for tropical intertidal flats the gap of information on diversity, distribution and abundance of these benthic organisms continues to be significant.Dittmann and Vargas (2001) and Dittmann (2002) have pointed out that among key organisms structuring the fauna of tropical tidal flats are burrow building sipunculans, brachiopods and enteropneusts.Their burrowing activities and the ingestion-egestion of sediments modify the substrate in such a way, that it becomes suitable or unsuitable for the settlement of larvae or the survival of juveniles of other benthic species (Gray & Elliott, 2009).In this study we followed Cutler (1994) and Emig (1983Emig ( , 1997) ) morphological criteria for the identification of the peanut worms as Sipunculus nudus and Apionsoma trichocephalus, and the lamp shells as Glottidia audebarti and G. albida, respectively.The acorn worms need further taxonomic work and we were reluctant to assign the specimens to categories below Class level.
As its name indicates, G. albida has a white shell while G. audebarti is characterized by the green color on the anterior region of the valves.According to Dall (1921) G. albida shells are streaked with brown, but this color pattern was present in only two of the larger specimens found.The type locality of G. albida is Magdelena Bay (California, U.S.A.) in 12 m depth sandy-mud, and for G. audebarti is Guayaquil Bay (Ecuador), in intertidal sand.The presently known distribution of G. albida is from California to Costa Rica, while G. audebarti is found from Mexico to Ecuador (Emig, 2009).
Lingulides live in burrows produced by themselves with the aid of shell movements.Lingulides are euryhaline invertebrates.The organisms retract into the burrows during low tide and emerge with the incoming tide to near sediment surface and filter the water for food.The pedicle end serves as anchor to facilitate vertical movements (Emig, 1997).
The peanut worm S. nudus is a sipunculan that has been known to scientists since Linnaeus described this species in 1766 (Hyman, 1959).Specimens collected around the world have morphological features that fit the description of S. nudus.Therefore, this species has been considered as cosmopolitan (Cutler, 1994).As such, its potential for use as a worldwide indicator of pollution was advanced by Vargas and Adbullah (1997).The distribution of S. nudus as a cosmopolitan species has been questioned based on recent molecular evidence provided by Kawauchi and Giribet (2013) that points to the existence of at least four geographical lineages.However, their study did not include specimens from Costa Rica.
This sipunculan is a prey item for birds foraging at low tide in the Gulf of Nicoya (Pereira, 1996).In China and other Asian countries there is a fishery that exports live S. nudus and sells it (market name: BB worm) as fish bait, mainly in Japan (Saito, Kawai, Umino, & Imabayashi, 2014).The worms are also a food item for human consumption in Asia (Saiz-Salinas, 1993).  1 0 0 0 0 2 1 0 3 0 0 5 2 0  S. nudus burrows actively in sandy sediments feeding on them as they burrow and has integumentary specializations for gas exchange thru the body surface (Ruppert & Rice, 1995) facilitating their infaunal mode of life.Sipunculans are able to survive drastic changes in salinity and their lower limit may be around 18 psu (Murina, 1984).Salinities around 20 psu are sometimes found in the upper Gulf of Nicoya estuary during the rainy season (Voorhis, Epifanio, Maurer, Dittel, & Vargas, 1983).
The acorn worms (Enteropneusta) include about 80 species grouped into several families separated by the structures of the gills, gonads, liver sacs, and coelomic diverticula (Cameron, 2005).The need to observe internal structures in complete specimens makes their proper identification difficult for the non-specialist.In addition, benthic sampling gear (grab, corer) normally yields fragmented specimens of these fragile and usually long organisms.The study by Deland, Cameron, Rao, Ritter, & Bullock (2010) on the family Harrimaniidae of acorn worms points out that the Eastern Pacific harbors many yet unknown species.In spite of these identification difficulties, information on the presence and abundance of acorn worms is useful for future comparisons of the structure (species diversity and abundance) and functioning (energy flow) of tidal flat systems.In spite of the fact that data from the Punta Morales flat was collected more than three decades ago, it remains as the only study on the ecology of these three groups of invertebrates for this region of the Eastern Tropical Pacific.Emig & Vargas (1990) pointed out that the 102 small specimens reported by Vargas (1987) as G. audebarti may include individuals of both species, while the five larger specimens collected the following year were of G. audebarti.Thus, abundances of lamp shells for 1984-1987 are reported as Brachiopoda in Fig. 3A.Data collected during 1984Data collected during -1987 indicates that brachiopods and acorn worms were scarce in cores collected early on 1984 (dry seasonearly rainy season) and increased to a maximum towards the end of 1984 (late rainy -early dry seasons).Sampling started shortly after the end of the period of high temperatures produced by the severe 1982-1983 ENSO.The impact of this warming event on the Peruvian benthos was described by Arntz, Brey, Tarazona, and Robles (1987).In spite of the fact that sampling effort was less after February 1985, numbers of brachiopods and acorn worms were low and patchy from early 1985 to early 1986, a period when red tides were frequent in the upper Gulf of Nicoya.The impact of these red tides on other infaunal organisms such as the polychaete worms has been discussed by Vargas-Zamora, Sibaja-Cordero, Dean, and Solano-Ulate (2015).Thus, we cannot rule out the possibility that the abundances of brachiopods and acorn worms found early on 1984, and those found during 1985-1986 were in some way also influenced by the end of ENSO and red tides, respectively.In summary, our data provides evidence of the presence of acorn worms in the upper estuary in sand-mud sediments.
The relatively low abundances of S. nudus and Glottidia spp.and enteropneusts at the Punta Morales sand mud flat particularly after the red tides of 1985 may indicate that other environmental factors were acting at this time.In this context the 2013 survey data indicates that sediments of the sand-mud flat changed to finer, softer, and higher organic matter content than those found in 1984, perhaps leading to the lack of brachiopods and S. nudus in the samples.Brachiopods and S. nudus seem to prefer coarser sediments like those found at Cocorocas in 2015.The enteropneusts appear more tolerant to these changes since numbers near 10 ind./14 cores were found both in 1984 and 2013.Our data also provides information on the relative importance of these three groups in comparison with other main macro-invertebrates from the Punta Morales flat collected during the 1984-1987 survey, like cephalochordates (Vargas & Dean, 2010), mollusks (Vargas-Zamora & Sibaja-Cordero, 2011), echinoderms (Vargas & Solano, 2011), crustaceans (Vargas-Zamora, Sibaja-Cordero, & Vargas-Castillo, 2012) and polychaetes (Vargas-Zamora, Sibaja-Cordero, Dean, & Solano-Ulate, 2015).
The presence of trace elements in brachiopod shells was reported by Jope (1965) and included Fe, Mg, Mn and Cu.However, data on the composition of shells and tissues of lingulide brachiopods are scarce.Therefore, the reported concentrations of metals in tissues and shells of G. audebarti are the first informed for this species.An important observation is that concentrations were different in shells, pedicles, and soft parts indicating different rates of accumulation.Published reports on metal concentrations in other living brachiopods are rare.
The study by De Moreno, Gerpe, Moreno, and Vodopivez (1997) from a pristine site in Antarctica is a notable exception.They found mean concentrations (μg/g dw) of Cd (4.33), Cu (4.87) and Zn (26.70) in tissues of unidentified brachiopods.As a comparison, mean concentrations (μg/g dw) of these elements found in soft parts of G. audebarti from the Cocorocas sand flat were: Cd (5.20), Cu (12.88) and Zn (55.20).Thus, concentrations of Cu and Zn in G. audebarti from the tropical estuarine flat were slightly higher than those from the Antarctic pristine site.
Reports on the evaluation of trace metal concentrations in tissues of sipunculans are also scarce.De Moreno, Gerpe, Moreno and Vodopivez (1997) found in unidentified whole sipunculans from pristine sites in Antarctic mean concentrations (μg/g dw) of Cd (0.47), Cu (14.64) and Zn (69.9).Mean concentrations (μg/g dw) of these metals found in whole non-depurated S. nudus from the Cocorocas sand flat during the 1996 survey were: Cd (1.2), Cu (26) and Zn (56), and during the 2000 (Table 3) survey the concentrations were: Cd (0.6), Cu (20.7) and Zn (81).Thus, mean concentrations of these three metals in nondepurated Antarctic and tropical sipunculans were of the same orders of magnitude.Yan and Wang (2002) found that sediments are a direct source of metal accumulation in S. nudus from a bay in China.They considered that sediment concentrations of less than 0.9 μg/g of Cd, less than 49 μg/g of Cr, and less than 140 μg/g of Zn represented a clean site, and maximum concentrations (μg/g) were 9.0 for Cd, 61.6 for Cr and 428 for Zn at the contaminated sites.For comparison, maximum mean concentrations (μg/g dw) of Cd, Cr, and Zn in the Cocorocas flat sediments (1996 and 2000 surveys) were 1.1, 31.5, and 63, respectively.Thus, the Cocorocas sediments may be considered as relatively clean based on Yan andWang (2002) criteria. Moreover, Vargas, Acuña González, Gómez, andMolina (2015) concluded that metal concentrations in sediments from the Cocorocas sand flat were those expected for a region of the estuary not as yet impacted by industrial activities.Wang, Yan, and Fan (2002) found that S. nudus from the same bay in China retained ingested sediment for about a day and were able to egest unassimilated metals (Cd, Cr, Zn) for up to 72 hours.Thus, the results reported here are in agreement with those found in China in specimens of S. nudus as concentrations of metals (Fe, Mn, Ni, Zn, Pb) in 72 hour depurated S. nudus from Cocorocas were lower than those in non-depurated organisms.The fraction of non-assimilated sediments is important enough to influence weight differences between depurated and nondepurated S. nudus.Therefore, data on metal concentrations in non-depurated S. nudus must be viewed with caution.Ha, Nhuan, Ngoc and Dung (2007) found in non-depurated S. nudus from Vietnam mean ppm concentrations (μg/g dw) of: Fe (441), Mn (3.7), Ni (4.7), Zn (21.9),Pb (5.5), Cu (1.8), and Cd (0.3).With the exceptions of Fe, the Vietnamese values were of the same order of magnitude than mean values reported here (μg/g dw) for Ni (10.4), Zn (81), Pb (2.8) and Cd (0.63) found in non-depurated S. nudus from Cocorocas.Higher concentrations of Fe, Mn and Cu in non-depurated S. nudus from Cocorocas, was probably related to their abundance in ingested and unassimilated sediments.Other species of intertidal invertebrates, like certain filter feeding bivalves, are able to lower their metal concentrations after depuration in clean sea water as found by Vargas, Acuña González, Gómez, and Molina (2015) for Fe, Mn, Ni, and Zn in the razor clam Tagelus affinis from the Cocorocas sand flat.Concentrations of Pb around 20 μg/g dw in sediments, 9.3 μg/g dw in nondepurated S. nudus, and 21 μg/g dw in shells of G. audebarti collected during 1996 may be related to the use of leaded fuel in outboard engines operated by the artisanal fishing fleet at the time; a small fishing village is located at the mouth of the Lagarto River.However, Pb was found in lower concentrations during the 2000 survey than during 1996.Other metals may have been carried to the sand flat by this river, which drains agricultural lands and brings heavy loads of fine and coarse sediments during the rainy season.Other compounds are also transported by the Lagarto river.For instance, a concentration of 41.4 ng/g dw of PCB was found by Spongberg (2006) in S. nudus from Cocorocas.However, this concentration was relatively low when compared with PCB values found in industrialized estuaries.
Data collected during the 2015 survey confirmed the presence of the brachiopods G. audebarti and G. albida at the Cocorocas sand flat after both species were first collected there by Emig and Vargas (1990).The survey yielded a maximum shell length of 26 mm for G. albida, and 38 mm for G. audebarti, which are in agreement with measurements made 25 years ago.Jones and Barnard (1963) conducted a spatial survey of the abundance of G. albida in subtidal coastal sediments of the California (33 o N) coast.This species was more common at depths from 22 m to 47 m on a wide variety of sediments.The size distribution of 8 976 specimens of G. albida collected by Jones and Barnard (1963) was dominated (86 %) by individuals of less than 5 mm, and only 1 % included individuals of 21-50 mm long.In tidal flats of the Colorado River mouth (31 o N) in the Gulf of California, small individuals of G. palmeri were rare and this species was found in patches of individuals of mean shell length of about 38 mm (Kowalewski, 1996).
On the other hand, at Cocorocas individuals smaller than 15 mm long characterized G. albida.This lingulide almost vanished from the sand flat in December, 2015.In addition, as found by Jones and Barnard (1963), G. albida survives well at its Northern range in non-estuarine coastal waters cooler than 20 o C and in depths shallower than 60 m.In 1980 a grab survey of 41 subtidal (6 to 65 m) stations in the Gulf of Nicoya included 15 stations in the upper estuary.However, as reported by Vargas, Dean, Maurer and Orellana (1985) only one Glottidia sp.specimen was found during the survey (Ballena Bay, 20 m, lower estuary).Estuarine low salinities and higher intertidal temperatures at the Cocorocas sand flat enhanced by the ENSO 2015 (NOAA, 2016) conditions may be among the possible factors preventing this subtidal species to reach a relatively stable presence at the intertidal flat.Moreover, temperatures of 40 o C and 41 o C were recorded at noon in tide pools from the Cocorocas sand flat on July 6 and December 15, 2015, respectively.
The presence of G. audebarti in the Gulf of Nicoya is at the middle of the latitudinal range of this mainly intertidal species (Emig, 1997) and these two factors may influence its relatively stable presence at Cocorocas.Paine (1963) described the life cycle of G. pyramidata living in intertidal sand bars on the West coast of Florida (29 o N) in salinities ranging from 19 psu to 35 psu.At these sand bars most of G. pyramidata individuals live for about a year with some individuals living up to 20 months.However, comparative information is lacking on the life span of both Glottidia species from the Gulf of Nicoya.
Natural mortality must be important at Cocorocas, although empty shells were never observed.Paine (1963) also points out that the phosphatic shell of Glottidia and its thin organic periostracum decay rapidly in the sediments.Predation on Glottidia spp.by sandpiper shore birds has been observed by Pereira (1996) at nearby sand flats in the Gulf of Nicoya.Naticid snail boreholes are frequently found on molluscan shells at the Punta Morales region (Vargas-Zamora & Sibaja-Cordero, 2011), but no such boreholes or traces of them have been observed on brachiopod shells from the sites.
At the Cocorocas sand flat maximum mean densities of G. audebarti and G. albida were 8.6 and 18.2 ind./m 2 , respectively.Jones and Barnard (1963) reported densities of G. albida of 20 ind./m 2 in coarser sediments where the onuphid polychate Diopatra sp. was common, while higher densities (up to 500 ind./m 2 ) coincided with the presence of the ophiuroid Amphioplus hexacanthus.At both the Cocorocas and Punta Morales flats the infaunal ophiuroid, Amphipholis geminata and the onuphid Diopatra ornata are conspicuos elements of the macrofauna (Ditmann & Vargas, 2001;Vargas & Solano, 2011).In sand flats from Australia (19 o S) L. anatina reached densities of 864 ind./m 2 (Kenchingtom & Hammond, 1978), while in intertidal flats in the Gulf of California (31 o N) G. palmeri had densities of up to 300 ind./m 2 (Kowalewski, 1996).A drastic decline in abundance, from a maximum of 500 ind./m 2 (mean 25-50 ind./m 2 ) in 1969 to 0.94 ind./m 2 in 2007, has been reported for L. reevii in reef flats from Hawaii (21 o N).The decline was related to a decrease in organic matter input and algal mats covering the sediments (Hunter, Krane, & Fitzpatrick, 2008).At the Bay of Bengal (India, 14 o N) L. translucida reached densities between 4 ind./m 2 and 92 ind./m 2 (mean 46 ind./m 2 ) estimated from 0.25 m 2 quadrats collected along three transects on an intertidal beach under the influence of a brackish water river (Raughnathan & Jothinayagam, 2007).Thus, densities of both G. audebarti and G. albida at Cocorocas appear relatively modest when compared with data from higher latitudes, but similar to those of G. translucida found in a tropical habitat in India under the influence of a river.
The 2015 sampling at the Cocorocas sand flat also confirmed the presence of the sipunculan S. nudus.Individuals representing nine weight groups of this peanut worm were found, but with most of the individuals were in the 0.1 to 1.0 g range.Representatives of this range increased in abundance from July to December indicating possible growth of the individuals.Both new (lighter) and heavier (probably older) specimens were present during the sampling period, which may represent a relatively stable population.A density of 5 ind./m 2 of S. nudus is mentioned by Saiz-Salinas (1993).Therefore, densities of about 34.5 ind./m 2 found at the Cococorcas sand flat are relatively high.
However, their mean size of 45 mm is relatively small if the maximum size (340 mm) also cited by Saiz-Salinas (1993) and a photograph in Ha, Nhuan, Ngoc, and Dung (2007) are considered for comparison.These latter authors indicate that S. nudus is more abundant in sand flats from Vietnam where sediments contain more than 80 % sand, and is rarely found if the sand content is less than 60 %.These values may explain the relatively high abundance of S. nudus in Cocorocas sediments (90 % sand) when compared to their low abundance in the 65 % sand content of the Punta Morales flat.
The distribution of individuals in space is a fundamental characteristic of soft-bottom species (Thrush, 1991).G. palmeri and other lingulides occur in patches in the intertidal zone (Kowalewski, 1996).Our data indicates that spatial patchiness characterizes the abundance of the three species at the Cocorocas sand flat.The two brachiopods and S. nudus were frequently found together within a single quadrat area, as illustrated for G. audebarti and S. nudus for the 18 quadrats collected on Dec. 15, 2015.The number of individuals collected in a single quadrat ranged from 0 to 11 (G.audebarti), 0 to 18 (G.albida) and 0 to 20 (S. nudus), which may indicate subtle differences in sediment grain size and other factors (like presence of adults) attractive or dissuasive for site selection by settling larvae and later survival of juveniles (Gray & Elliot, 2009).
Spatial and temporal patchiness is also important in the context of fishing potential of certain species since relatively high mean densities at a given sediment spot or time may be misleading.The sipunculan S. nudus is not as yet exploited commercially in Costa Rica, but it has been used as fish bait for generations.The peanut worm coexists at Cocorocas with the lamp shells and the small razor clam T. affinis which is already under moderate fishing pressure.
Brachiopods of the genus Lingula are harvested commercially in Southeast Asia (Brusca, Moore, & Shuster, 2016).Future exploitation of populations of lamp shells in the Gulf of Nicoya may be expected, and management policies are urgently needed to prevent future overfishing of these species.Other tidal flat species from the Gulf of Nicoya, like Anadara spp.ark clams (Stern-Pirlot & Wolff, 2006) and the giant polychaete Americonuphis reesei (Rojas & Vargas, 2008) are examples of living resources where management strategies have been implemented with relatively good results in this tropical estuary.

Fig. 1 .
Fig. 1. A. Location of the Punta Morales and Cocorocas intertidal flats.Punta Morales Peninsula, mid upper Gulf of Nicoya estuary.Pacific, Costa Rica.(10ºN-85ºW).1. Tempisque River, 2. Tárcoles River, 3. Lagarto River.B. Sand mud flat.C. Sand flat.Two islets (outlined in A) are visible in the background.Note horizontal dark band marking high tide level on the tree islet.
is found and harvested occasionally for local consumption.

Fig. 2 .
Fig. 2. A. Anterior regions (proboscis, collar, and part of trunk) of five Rose Bengal stained acorn worms.B. Live Sipunculus nudus collected at the sand flat.C. S. nudus drawn from live specimen.D. S. nudus during the 72 h depuration process.E. G. albida.Note setae bundles.F. Rose Bengal stained Glottidia audebarti.Note tips of pedicles with agglomerated sand grains.

Fig. 3 .
Fig. 3. Number of individuals per date (49 dates) of: A. Lingulide brachiopods (Brachiopoda, light grey bars) and peanut worms (Sipuncula, black bars).B. Acorn worms (Enteropneusta, dark gray bars) found in a 400 m 2 sampling plot at the Punta Morales intertidal sand-mud flat, Gulf of Nicoya estuary, Pacific coast of Costa Rica(1984-1987).Feb. 22, 1984 to  Feb. 21, 1985:  Semi-monthly sampling of two sets of 14 cores per date.Numbers of acorn worms and brachiopods found on each 14 core set are indicated on the top of the bars.Sipunculans were found on one of the sets.Mar.07, 1985 to Apr. 29,  1987:  Monthly sampling of 14 cores per date.Core area 17.7 cm 2 , core depth 15 cm.Mesh screen: 500 microns.