Coastal development in Costa Rica increased drastically towards the end of the XXth Century, and information on water conditions, water dynamics, and concentrations of nutrients and pollutats was urgently needed. Thus, an ecological assesssment was started in 2000 at four coastal marine embayments of Costa Rica, and the sudy of water characteristics, concentration of pollutants and primary productivity estimates, were part of the study. Three embayments were located along the Pacific coast and the fourth on the Caribbean coast. The results of this research have been published mostly by the Revista de Biología Tropical, as part of the series Ecosystems of Costa Rica, as in this case.

The tidal currents in the Gulf of Nicoya were studied by Lizano & Alfaro (2004). Data on concentrations of Poly-Chlorinated Biphenyls (PCBs) in sediments and tissues of sipunculan worms were discussed by Spongberg (2004 a, b, c; 2006); a review of the literature on PCB concentrations in tropical media was published by Spongberg & Witter (2008). Concentrations of petroleum hydrocarbons in waters were determined by Acuña-González, Vargas-Zamora, Gómez-Ramírez, & García-Céspedes (2004). The presence of coliform bacteria in waters and beach litter at the four sites were reported by García, Acuña-González, Vargas-Zamora, & García-Céspedes (2006). The impact of certain pollutants as potential endocrine disruptors of coastal organisms was reported by Cheek (2006), and Valverde, Selcer, Lara, & Sibaja-Cordero (2008). Gravel, Johanning, McLachlan, Vargas, & Oberdörster (2006) published the first report for Central America on imposex in the intertidal snail Acanthais brevidentata. As a follow-up of the evaluation of trace compounds, Spongberg et al. (2011) reported the presence of Pharmaceutical and Personal Care Products (PPCPs) in rivers flowing into the four embayments and other coastal sites.

Data on the ecology of selected coastal organisms from the embayments were published by Sibaja-Cordero & Vargas-Zamora (2006) and by Rojas-Figueroa & Vargas-Zamora (2008). As part of the study, the lists of polychaete and sipunculan worms found in Costa Rica were updated with the goal of identifying potential pollution-indicator species (Dean, 2004; 2009; Vargas & Dean, 2009). The use of polychaetes as indicators of pollution in tropical sites was further discussed by Dean (2008). Data on trace metals in sediments from the four embayments were summarized by García-Céspedes, Acuña-González, & Vargas-Zamora (2004), while Vargas, Acuña-González, Gómez, & Molina (2015) published data on trace metals from snails (including A. brevidentata) and clams collected at the four sites. More recently Vargas, Acuña-González, Vásquez, & Sibaja-Cordero (2016) reported on the abundances and concentrations of trace metals in tissues of the sipunculan, Sipunculus nudus and the lingulide brachiopods, Glottidia audebarti and Glottidia albida from a sand flat in the upper Gulf of Nicoya.

Thus, the objective of this report is to make accessible data collected on water parameters, their temporal and spatial distribution, and on primary productivity, gathered during the 2000-2002 environmental assessment of four coastal sites of Costa Rica.

MATERIALS AND METHODS

Study sites: The four coastal embayments sampled in this study differ in the extent of their geographical area. The Gulf of Nicoya has the greatest area, followed by Culebra Bay, Golfito Bay, and Moín Bay. Three of the sampling sites (Culebra Bay, Gulf of Nicoya, and Golfito Bay) are located along the Pacific coast of Costa Rica (Fig. 1A), where mean tidal range is near 3 m. The fourth site, Moín Bay, is on the mid Caribbean coast (Fig. 1A), where mean tidal range is usually less than 0.5 m.

Culebra Bay is an embayment on the Northern Pacific coast of Costa Rica and its mouth opens into the Gulf of Papagayo (Fig. 1B), an upwelling region. No major river or creek flows into Culebra Bay. The shores are the sites for many tourism facilities that benefit from numerous white sandy beaches and coral patches (Jiménez, 2001). The Bay has a maximum depth of 30 m at its mouth. Three sampling stations were located inside the Bay in front of Esmeralda, Panamá, and Nacascolo beaches (Fig. 1B).

The Gulf of Nicoya is an estuary on the central Pacific coast of Costa Rica (Fig. 1C). Since the Mid-20th Century the Gulf has been the main fishing ground of the country for finfish and shellfish. Puntarenas and Caldera, located on the middle Gulf, are the most important ports on the Pacific coast. The estuary is under the influence of seasonal rainfall (May to November) followed by a dry season (December to April). Wind is also important as a water mixing force mainly during the early dry season. Depths at low tide in the upper Gulf range from about 40 m between Chira Island and the Port of Puntarenas, to less than 2 m near the mouth of the Tempisque River (Fig. 1C). As for this study the Gulf was divided into four regions: upper (near the mouth of the Tempisque River, including Chira Island and Jicaral estuary), middle right (mear Cortezas islands and the Punta Morales port), middle-left (near Caballo Island) and lower region, near the port city of Puntarenas and San Lucas Island (Fig. 1C).

Golfito Bay is a small, shallow (less than 5 m at low tide) and semi-enclosed embayment on the Northern shore of Golfo Dulce, the main gulf in the Southern Pacific coast of Costa Rica (Fig. 1D). The Golfo Dulce region has heavy rains year around. Golfito Bay is the site of the main port of Golfo Dulce and the bay is bordered mostly by red mangroves along most of its shores. Golfito Bay receives freshwater input from many small creeks and the sewerage system of the port city. Two sampling stations were deployed at the mouth of Golfito Bay and one station near the Pelícano islet on the upper bay region (Fig. lD).

Moín Bay is on the central Caribbean coast of Costa Rica and is directly exposed to the Caribbean Sea (Fig. 1E). The Moín River runs parallel to the shore and is bordered by mangroves. Brackish waters are found at the river mouth (Mata, Acuña, Murillo, & Cortés, 1987). The Moín port houses the main Caribbean facilities for docking of small fuel tankers and container ships, and is under the process of further development. Depths range from about 10 m inside the port to near 25 m at the port entrance. Sampling stations included one near the main dock, and two along the navigation channel (green and red buoys in Fig. 1E).

Water parameters: A total of 128 water stations were taken at the four embayments: 68 on the Gulf of Nicoya estuary, 30 on Moín Bay, 19 on Golfito Bay, and 11 on Culebra Bay. The term station is used here to mean a geographical sampling point (located by GPS coordinates) where water samples were taken at a given date and time of the day. For each station, 13 parameters were recorded: 1. Light penetration, 2. Water temperature, 3. Salinity, 4. Dissolved oxygen, 5. Oxygen saturation, 6. Chlorophyll-a, 7. Phaeopigments, 8. Suspended matter, 9. Ammonium, 10. Nitrite, 11. Nitrate, 12. Phosphate and, 13. Silicate. The stations were reached by an inflatable boat, and a hand operated winch was used to handle the sampling equipment.

The penetration of light in the waters was estimated by the use of a 0.3 m diameter Secchi disk and following Holmes (1970). Water samples were collected with a 5 L oceanographic General Oceanics® Niskin bottle at 1 m depth. The water-filled closed Niskin bottle was brought aboard and temperature and dissolved oxygen were measured in the bottle with a probe of a calibrated Yellow Springs® Instrument Model 57 Oxygen meter, after compensating for salinity. Percent oxygen saturation was also computed. Water salinity was determined from drops of seawater placed on a calibrated (Standard Sea Water) American Optical® refractometer. Water samples from the Niskin bottle were transferred to dark plastic 2.5 L bottles and placed in a cooler filled with ice and transported to the nearby ground facility. Water samples were filtered at the field laboratory and the glass microfibre filters were folded, placed in plastic bags and preserved on ice for chlorophyll-a and suspended matter determinations following Strickland & Parsons (1972). A Shimadzu® UV-1700 spectrophotometer was used for the determination of chlorophyll-a and phaeopigments. Filtered water samples for nutrient analyses were kept in dark plastic 0.5 L bottles and stored in a freezer (-18 oC) until analysis. Spectrophotometric (Shimadzu UV-1700) determination of ammonium, nitrite, nitrate, phosphate, and silicate concentrations were carried out in triplicate (10 mL samples) following methods based on Strickland & Parsons (1972) and validated protocols in use at the UCR-CIMAR Chemical Oceanography Laboratory. Before each batch of samples was analyzed, calibration curves were run for each analyte and the detection limits were obtained according to Meier & Zünd (1993).

Primary productivity: Water column productivity was estimated at selected stations at the four sites using the light and dark bottles method, determining changes in dissolved oxygen concentrations, as in Riley & Chester (1971), and Strickland & Parsons (1972). Oxygen chemical determinations were performed according to the modified Winkler method (Gocke, Cortés, & Villalobos, 1990). Plastic rims were designed to hold the set of bottles for a selected depth, and the set of rims was kept in place by an anchor weight at the bottom and a buoy at the surface. Depths for the incubation of two dark and two clear volume-calibrated 300 mL BOD bottles, plus two clear bottles as reference, were selected based on the depth of the Secchi disk at the site, and on the light extinction exponential equation using an extinction coefficient k = 1.7, as indicated by Holmes (1970). Also, depending on the bottom depth, two to five incubation rims of light and dark BOD bottles representing from 100 % to 1 % of the incident light were selected. The incubation period lasted three hours on the average, and usually one period in the morning and another in the afternoon were accomplished. To convert data from mg of oxygen to mg of carbon, a photosynthetic coefficient of 1.2 and a respiration coefficient of 0.83 were used based on Gocke et al. (1990).

Statistical analyses: All statistical analyses were carried out with software PAST (Hammer, Harper, & Ryan, 2001).

A. The four sites: A between groups Principal Component Analysis (PCA) was computed with the standardized (z scores) values for the 13 parameters of stations sampled at the four embayments, with a star representing the centroid for the group of stations from an embayment (Culebra Bay, Gulf of Nicoya, Golfito Bay, Moín Bay). Ellipses were computed to include stations from an embayment within 95 % confidence limits. A PCA Gabriel biplot was also computed to illustrate the relative importance of each parameter on the grouping of stations representing the embayments (Dolédec, & Chessel, 1991; Hammer et al., 2001; Legendre, & Gallagher, 2001). To evaluate for spatial differences in the environmental parameters found by PCA, two PERMANOVA (Anderson, 2001) tests were computed with a posteriori comparisons and 9999 permutations. Two matrices were employed: An Euclidean distance matrix based on the standardized environmental values for Gulf of Nicoya only, and another matrix of standardized values of the parameters of the four embayments.

B. The Gulf of Nicoya: Two groups of similar number of stations taken during the dry and rainy seasons in regions 1, 2 and 3 (Fig. 1C) of the Gulf of Nicoya estuary were selected to test for seasonality of the 13 parameters measured. Seasonality was evaluated by means of a Mann-Whitney U test comparing medians (Conover, 1971). A between groups PCA was computed as a graphic display of stations sampled in the Gulf. In order for the software to run properly, all data on concentrations with values below the detection limit were included in the data matrix as 0.00. Data was standardized to z scores before computing PCA. The centroid for each group of stations from a region (Upper, Middle right, Middle left, Lower) was represented by a star, and a hand-drawn polygon enclosed all stations for a region. A PCA Gabriel biplot was computed to illustrate the relative importance of each parameter on the group of stations representing a region (Dolédec & Chessel, 1991; Hammer et al., 2001; Legendre & Gallagher, 2001).

C. Primary productivity: A Spearman rho rank correlation (Conover, 1971) was computed between depth-integrated net primary production per day (Net PP/day) and each of the 13 environmental parameters.

RESULTS

A. The four sites: Code numbers for the 128 stations taken at the four sites and date of collection are included in Table 1. Results of the PCA performed on the total of 128 stations sampled at the four embayments are included in Fig. 2A. Three stations from the Gulf of Nicoya (14, 92, 102) are outside the 95 % confidence limits. The PCA separation of the four embayments included three (Gulf of Nicoya, Culebra Bay and Golfito Bay) with stations that stand alone in one of the sides of their ellipses (and maximizing the differences), while the centroid for the Moín Bay stations was near the other three centroids (Fig. 2A). Moreover, Moín was separated by a third axis (PC-3, not plotted) that explains 6 % of the variation. The PERMANOVA test confirmed that the four embayments could be separated based on their environmental parameters (F = 9.0, p < 0.001). However, the four sites also presented several stations that overlapped in their values (Fig. 2A).

Results of the PCA Gabriel biplot for the 128 stations are included in Fig 2B. Some Culebra Bay stations appeared to cluster more along the line representing higher Secchi depth and salinity, while some stations from the Gulf of Nicoya clustered along the lines representing higher nitrate, phosphate, phaeopigments and chlorophyll-a. Several stations from Golfito Bay seemed to follow the lines for increasing temperature and concentrations of oxygen, silicate, ammonium. Stations from Moín bay, however, were foundto group with higher values of salinity and Secchi depth and lower chlorophyll-a and phaeopigments than the Gulf of Nicoya.

Examples of stations (26, 27, 28, 32, 44, 47, 45, 46, 48, 49, 59, 62, 63, 66, 67, 64, 119) from the other three bays (Golfito, Culebra, Moín) are listed in Table 2. Parameters for stations 92 102 and others from the Gulf of Nicoya are included in Table 3.

B. Gulf of Nicoya: The list of the 68 stations taken on the Gulf of Nicoya is included in Table 1. The maximum and minimum values for the 13 parameters, their means, medians, standard errors, standard deviations, and the 25 and 75 percentiles, as well as selected examples of representative stations, are included in Table 3.

Results of the PCA performed on the 68 stations from the Gulf of Nicoya indicate that PC-1 and PC-2 explained 90 % of the variance (Fig. 3A). A gradient of change was found for the environmental parameters from Upper to Lower region following the PC-1. PC-2 separated the Middle region from the Upper and Lower regions of the estuary.

Results of the PCA Gabriel biplot for the Gulf are included in Fig 3B, and indicate three main parameter trends for station separation: 1- Secchi depth, salinity and nitrate (with higher values in the Lower region), 2- Percentage of dissolved oxygen saturation, chlorophyll-a and phaeopigments (with higher values at the Middle of the Gulf) and, 3- Temperature, suspended matter, phosphate, silicate, and nitrite (higher values in the Upper region).

Results of the Mann-Whitney U test, performed on a set of 49 stations (25 during the dry season and 24 during the rainy season) from the Gulf of Nicoya are included in Table 4. The U test detected seasonal differences (dry vs rainy, p < 0.05) for: salinity, dissolved oxygen, and the percentage of dissolved oxygen saturation, with higher values during the dry seasons. Phosphate, silicate, and nitrite, had higher concentrations during the rainy seasons. No seasonal differences were found for the other seven parameters (Table 4).

There was a gradual change in the parameter values within the estuary (PERMANOVA, F = 3.39, p = 0.0001). The Lower region was similar to the Middle right region (p = 0.105), but it differed from the rest (p < 0.005). Both Middle regions were similar (p > 0.281), but only the left region was similar to the Upper region (p = 0.080).

C. Primary productivity: Results of the depth-integrated Net Primary Productivity estimates at the four embayments are included in Table 5. Net productivity varied within and between days. Oxygen consumption was greater than production (negative in Table 5) ) at several dates and times in Culebra Bay, Golfito Bay and Moin Bay (Table 5, A, C, D). Oxygen production was positive for all dates in the Gulf of Nicoya (Table 5B). However, the maximum net primary productivity estimate for the four embayments was 6.044 g C/m2/day (2 206 g C/m2/year) from Culebra Bay on June 29, 2000 (Table 5A). This highest productivity estimate was preceeded by 1 218 g C/m2/year at the site during the afternoon of the previous day, and a negative production in the morning of the same day. The lowest net primary production was recorded at Moín Bay on October 5, 2000 (30 g C/m2/year).

Results of the Spearman rho rank correlation test between net primary productivity (g C/m2/day) and the environmental parameters are included in Table 5 E. Only temperature and silicate had p values below 0.05.

DISCUSSION

Our initial hypothesis was that the four coastal embayments were different enough in their area, depth, tidal amplitude, exposure, rainfall, and morphology, so they could be clearly separated from each other based on the data set. Although the results of the PCA test confirmed that the four sites could be separated based on the 13 parameters, they had several overlapping stations. This latter result points out to the need of deploying as many stations as possible to account for spatial and temporal variability. This is important if future studies dealing with within-site and between-sites comparisons are conducted aimed at the identification of trends at different scales, like eutrophication and climate change.

The four embayments sampled during the 2000-2002 are under the influence of seasonal rainfall. This seasonality is a key factor on the estuarine dynamics of the Gulf of Nicoya. Seasonal rainfall is also important in Golfito and Moín bays, and less so in Culebra Bay, where upwelling events may be more important. The key role played by rainfall as a source of nutrients via runoff to temperate estuaries has been emphasized by Mallin, Pearl, Rudek, & Bates (1993). They also point out that changes in global or local weather patterns influencing rainfall are likely to have a significant impact in estuarine nutrient dynamics and productivity.

Culebra Bay was the site with the clearest water of the four embayments (Secchi depth: 12.6 m) and the temperature had the widest range (24.0 to 32.5 oC). It was also a bay with low values for suspended matter (4.12 mg/L), and non detected phaeopigments and nitrite. Maximum chlorophyll-a concentration was the lowest (2.23 mg/m3) of the four embayments. In a year study (2004-2005), Fernández-García, Cortés, Alvarado, & Nivia-Ruiz (2012) reported maximum / minimum values (Units as in Table 2) of: temperature (29.6 / 22.9), salinity (34.3 / 30.0), chlorophyll-a (3.8 / 0.4), suspended sediments (10.2 / 3.3), phosphate (1.45 / 0.03), silicate (5.88 / non detected), nitrite (1.18 / non detected), nitrate (15.60 / non detected), and ammonium (0.52 / non detected). These parameters make Culebra Bay waters characteristic of open ocean conditions. However, cooling conditions occcur from time to time (Alfaro & Cortés, 2012; Rixen, Jiménez, & Cortés, 2012; Rodríguez-Sáenz & Morales-Ramírez, 2012) Water warming related to El Niño Southern Oscillation (ENSO) is also important in the region and has been the cause of coral bleaching and mortality at Culebra Bay (Jiménez, Cortés, León, & Ruiz, 2001).

The information on water characteristics of the Gulf of Nicoya comes mainly from cruises aboard the RV Skimmer in the early 1980’s (Vargas, 2016). The reports by Voorhis, Epifanio, Maurer, Dittel, & Vargas (1983) and Epifanio, Maurer, & Dittel (1983) clearly identified the estuarine nature of the embayment. The influence of seasonal (May to November) freshwater input due to rains falling directly on the estuary, or carried there mainly by the Tempisque river as runoff from its extensive watershed, has been emphasized by Vargas & Mata (2004) and Vargas (2016). The results of this study confirm the seasonality of certain water parameters in the estuary. Higher salinities are expected during the dry season and lower values during the rainy season. Moreover, the influx of certain nutrients, like phosphate and silicate, during the rainy season were also expected. On the other hand, during the dry season (December to April) the estuary is under the influence of Easterly winds that increase water mixing. Thus, higher values for dissolved oxygen and percent saturation were expected during the dry season. Epifanio et al. (1983) found that, with the exception of phosphate, concentrations of ammonium, nitrite and nitrate varied seasonally. During 2000-2001, Palter, León-Coto, & Ballestero (2007) collected data on dissolved oxygen, nutrients and chlorophyll-a in the Gulf of Nicoya. They found that nutrients were more concentrated in the upper Gulf during the rainy season and concluded that the Tempisque River controls spatial and temporal variability of nutrients in the upper Gulf. The maximum values found by Palter et al. (2007) in the Gulf (Stations taken along the estuary at a depth range of 0 to 35 m (Units as in this study), were: Dry / Rainy seasons (Secchi depth: 0.4 / ND; Temperature: 28.4 / 29.9; Salinity 29.2 / 26.9; Percentage of oxygen saturation: 130 / 131; Chlorophyll-a 3.3 / 12.63; Phosphate 1.75 / 3.35; Silicate 127.7 / 165.3: Nitrite 0.77 / 1.71; Nitrate 9.10 / 10.27. With the exception of higher phosphate (9.50) and lower nitrate (4.59) found in our study, the values found in their survey for other parameters are similar to those found in our study.

The estuarine nature of Golfo Dulce, where Golfito Bay is located, has been studied by Svendsen et al. (2006). However, Golfo Dulce differs from the Gulf of Nicoya on its deep bathymetry (200 m) and the presence of a sill (60 m) at the entrance, which makes Golfo Dulce similar to an European fjord. This similarity has encouraged comparative research leading to significant discoveries, such as the anammox reaction (Dalsgaard, Canfield, Petersen, Thamdrup, & Acuña-González, 2003). In contrast with Golfo Dulce, Golfito Bay is a small shallow semi-enclosed embayment connected to Golfo Dulce by a tidal channel. Comparative data on dissolved oxygen and nutrient profiles for Golfito bay is available from a 1998 survey near the mouth of a creek that flows near the Northeastern end of the Bay (Silva & Acuña-González, 2006). They found maximum values (Units as in this study) of: Secchi depth (2.0), salinity (33), temperature (33.0), suspended matter (96.0). Maximum nutrient concentrations, were: nitrate (13.52), nitrite (0.88), phosphate (0.86), and silicate (110.4). Maximum chlorophyll-a concentration was 7.9 mg/L. The data reported by Silva & Acuña-González (2006) is similar to data found during our 2000-2002 study, except for a lower maximum concentration (1.33) of nitrate. Moreover, suspended matter in Golfito Bay had the highest concentration (156.04 mg/L) of the four embayments. The Golfito Bay region is under the influence of heavy rains during most of the year and runoff flows directly into the bay. This runoff carries nutrients and suspended matter from the nearby rain forest, the mangrove forest, and also brings into the bay loads of untreated sewage and other pollutants from the city and port activities (Spongberg & Davis, 1998; García et al., 2006).

Moín port is within a bay exposed to the Caribbean Sea. The Moín River carries sediments and nutrients from mangrove forests and agricultural lands. The estuary also houses a fluvial port and an oil refinery further upstream (Acuña-González et al., 2004). Data on nutrients from the Caribbean coast comes mainly from a 1997 survey by Muller-Parker & Cortés (2001). They found minimum-maximum concentrations (µmol/L) of: phosphate (0.06 to 0.54), ammonium (0.3 to 15.4), nitrite (0.01 to 0.13) and nitrate (0.20 to 3.56). Minimum-maximum concentrations (µmol/l) found at Moín Bay in our 2000-2002 survey, were: phosphate (0.02 to 2.91), ammonium (0.07 to 2.24 µmol/L), nitrite (0.00 to 0.40 µmol/L), nitrate (0.01 to 1.63). Their results are within the range of those reported in this study. These relatively low values are similar to those found in Culebra Bay. Both Culebra and Moin bays are under the direct influence of Pacific and Caribbean conditions, respectively.

In this study we followed the definition of Cloern, Foster, & Klechner (2014) of productivity as the rate of production per unit time. Based on early reports like Riley & Chester (1971), our hypothesis was that higher productivity was to be expected in estuarine embayments like the Gulf of Nicoya and Golfito Bay, and lower values in near open-ocean conditions, such as Culebra and Moin bays. However, the highest net primary productivity (2 206 g C/m2/year) found in this study was mesusred in Culebra Bay. This value followed lower and negative productions during the afternoon and morning of the previous day, respectively. This high variability may be indicative of patchy environmental conditions acting within hourly and daily time frames. Red tide outbreaks have been observed occasionally in Culebra Bay (Morales-Ramírez, Víquez, Rodríguez, & Vargas, 2001), but no patches were observed during our survey.

Primary productivity estimates using the same methods (Oxygen change after incubation) as in this study have been conducted in the Gulf of Nicoya by Córdoba-Muñoz (1998) and Gocke, Cortés, & Murillo (2001a). Both studies were performed near the Cocorocas islands. Córdoba-Muñoz (1998) found a maximum value of 1.81 g C/m2/day (662 g C/m2/year), while Gocke et al. (2001a) reported maximum net primary production of 2.76 g C/m2/day (1 007 g C/m2/year) at a mangrove tidal channel, and higher net PP values were influenced by red tides. Red tides have been reported frequently on the Gulf of Nicoya (Viquez & Hargraves, 1995) but no patches were observed during our 2000-2002 survey. Net primary productivity estimates in the Gulf of Nicoya estuary found during this study had all values as positive, with a maximum of 1 601 g C/m2/year and a minimum of 191 g C/m2/year. These estimates of net PP are similar to those previously reported for the Gulf of Nicoya.

Net primary productivity in Golfito Bay and Moín Bay had maximum values of 1 371 and 1 100 g C/m2/year, respectively. These values are similar to those found in the other two embayments.

If the above mentioned maximum values of Net Primary Production (Net PP) found at the four embayments are considered as representing unusual conditions and therefore are excluded when computing the average values, the following trend of decreasing estimates of Net PP (g C/m2/year) emerges: Gulf of Nicoya (619) > Culebra Bay (456) > Golfito Bay (180) > Moín Bay (92), which fits into the general acceptance of estuaries and upwelling regions as highly productive ecosystems. Moín Bay is close to low productivity tropical ocean conditions. Golfito Bay is a special case due to its small area, restricted circulation, and the impact of urban and forest runoffs.

It has been argued that primary production in coastal areas may be nutrient-limited (particularly nitrogen and phosphorus) mainly in temperate estuaries (Epifanio, Maurer, & Dittel (1983). However, information is scarce for tropical ecosystems. This problem is becoming more relevant in the context of eutrophication of coastal ecosystems and climatic change (Chávez, Messie, & Pennington, 2011). The recent review by Cloern, Foster, & Kleckner (2014) found a range from negative to 1 890 g C/m2/year in 1 148 estimates of primary production in estuarine systems (including the Gulf of Nicoya). They pointed out that PP may vary up to ten fold within ecosystems. Based on their figure of 500 g C/m2/year separating eutrophic from hypertrophic conditions, estuaries like the Gulf of Nicoya may reach eutrophication frequently (mean 619 g C/m2/year) or become hypertrophic at times. Future studies may provide information to evaluate this trend towards hypertrophic conditions.

Primary productivity estimates are based on the relationship between oxygen produced and carbon fixed thru the photosynthetic quotient (PQ). In this study we used a value of PQ of 1.2. However, PQ may vary from 1.0 to 1.4 (Cloern et al., 2014). Estimates of primary productivity, based on the change of dissolved oxygen concentration (determined by the modified micro Winkler method) after a period of incubation of a water sample, have been reported for decades worldwide. However, the results must be viewed with caution when the method is employed either in oligotrophic waters, or in waters with high bacterial counts (Riley & Chester, 1971). In our study, oligotrophic conditions were expected in Moín Bay year around, and in Culebra Bay during non-upwelling periods. During the 2000 survey Culebra Bay was found the least polluted (coliform bacterial counts) of the four embayments, while Golfito Bay had the highest bacterial counts (García et al., 2006). Moreover, the role of phytoplankton in net primary productivity is crucial. However, phytoplankton in general has a patchy distribution in space and time. Thus, which type and how many microalgal cells are finally incubated in net PP estimates is a guess, and no correlation of net PP with chlorophyll-a is frequently found, as in this study. To what extent the results of the net PP estimates for the four embayments were influenced by the above mentioned factors is unknown.

Of the 1 148 estimates of Annual Phytoplankton Primary Productivity (APPP) compiled by Cloern et al. (2014) only 15 were from the South and Central Americas, and included only one report (Gocke, Cortés, & Murillo, 2001b) from Costa Rica (Gulf of Nicoya). Thus, the need for future estimates of net primary productivity as APPP is growing to facilitate comparisons with other environments worldwide, and to assess trends of productivity influenced by coastal development, pollution, and global climate change. Data included in this study are a key reference for the Central American region.

ACKNOWLEDGMENTS

The collection and analyzes of samples was made possible by grants from the Costa Rica-United States of America Foundation for Cooperation (CR-USA), and the University of Costa Rica (UCR) Project 808-AO-506. We thank Davis Morera and Eleazar Ruiz for their help in the field and the laboratory.

RESUMEN

Parámetros del agua y productividad primaria en cuatro bahías marinas de Costa Rica (2000-2002). La disponibilidad de datos sobre parámetros acuáticos y productividad primaria es esencial en estudios comparativos orientados a la identificación de tendencias ambientales, como la eutroficación, en ambientes costeros. En este contexto, el objetivo de este trabajo es hacer accesible datos, colectados entre el 2000 y 2002, de 13 parámetros acuáticos (Profundidad de Secchi, temperatura, salinidad, oxígeno disuelto, saturación de oxigeno, material suspendido, clorofila-a, faeopigmentos, fosfato, silicato, amonio, nitrito y nitrato) su variación temporal y distribución espacial, y productividad primaria, en cuatro sitios costeros de Costa Rica. El estudio multi-paramétrico fue realizado con base en muestras superficiales (1 m profundidad) de agua recolectadas en 128 estaciones en tres bahías costeras en la costa del Pacífico y una en la del Caribe de Costa Rica: una bahía (Bahía Culebra) bajo la influencia de afloramiento costero, un estuario de marea (Golfo de Nicoya), una pequeña bahía semi-cerrada rodeada de manglares (Bahía de Golfito), y una bahía y puerto del Caribe (Bahía Moín). Mediante Análisis de Componentes Principales (PCA) se obtuvo una distribución grafica de las estaciones de los cuatro sitios en general y para el Golfo de Nicoya en particular. PCA Gabriel biplots fueron utilizados para ilustrar la importancia relativa de cada parámetro en la separación de las estuaciones. Los cuatro ambientes fueron separados por el PCA con base en los 13 parámetros. Los valores máximos encontrados fueron: Secchi: 12.6 m; salinidad: 35 psu, temperatura: 32.5 oC, oxígeno disuelto: 11.4 mg/L; porcentaje de saturación de oxígeno disuelto: 188 %; material suspendido: 156.04 mg/L; clorofila-a: 16.21 mg/m3; faeopigmentos: 6.37 mg/m3; fosfato: 9.53 µmol/L; silicato: 156.11 µmol/L; amonio: 2.66 µmol/L; nitrito: 2.47 µmol/L y nitrato: 4.59 µmol/L. La salinidad y algunos nutrimentos variaron estacionalmente (seca vs lluviosa) en el estuario del Golfo de Nicoya. Estimaciones de la productividad primaria neta (P.P. neta, método de oxígeno de Winkler) produjo un máximo de 2 206 g C/m2/año en Bahía Culebra. Los máximos de PP neta (g C/m2/año) en los otros tres sitios, fueron: 1 601 (Golfo de Nicoya), 1 371 (Bahía de Golfito), y1 100 (Bahía Moín). Los valores de PP neta mostraron una alta variabilidad espacial y temporal, y oscilaron desde negativos hasta aquellos asociados a condiciones oligotróficas, eutróficas e hipertróficas. En el Golfo de Nicoya 60 % de los valores de PP neta excedieron el valor de referencia (500 g C/m2/año) para separar condiciones eutróficas de las hipertróficas.

Palabras clave: nutrientes, clorofila-a, eutrofización, estuario, afloramiento, Pacífico Este, Caribe.

REFERENCES

Anderson, M. J. (2001). A new method for non-parametric multivariate analysis of variance. Austral Ecology, 26, 32-46.

Acuña-González, J., Vargas-Zamora, J. A., Gómez-Ramírez, E., & García-Céspedes, J. (2004). Hidrocarburos de petróleo, disueltos y dispersos, en cuatro ambientes costeros de Costa Rica. Revista de Biología Tropical, 52(Supl. 2), S43-S50.

Alfaro, E. J., & Cortés, J. (2012). Atmospheric forcing and cool subsurface water events in Bahia Culebra, Gulf of Papagayo, Costa Rica, Revista de Biologia Tropical, 60(Suppl. 2), S173-S186.

Chávez, F., Messie, M., & Pennington, J.T. (2011). Marine Primary Production in relation to climate variability and change. Annual Review of Marine Science, 3, 227-260.

Cheek, A.O. (2006). Subtle sabotage, endocrine disruption in wild populations. Revista de Biología Tropical, 54(Suppl. 1), S1-S19.

Cloern, J.E., Foster, S.O., & Klechner, A. E. (2014). Phytoplankton primary production in the worlds’s estuarine-coastal ecosystems. Biogeosciences, 11, 2477-2501.

Conover, W.J. (1971). Practical non-parametric statistics. New York, USA: John Wiley & Sons, Inc.

Córdoba-Muñoz, R. (1998). Primary productivity in the water column of Estero Morales, a mangrove system in the Gulf of Nicoya, Costa Rica. Revista de Biología Tropical, 46(Suppl. 6), S257-S262.

Dalsgaard, T., Canfield, D.E. Petersen, J. Thamdrup, B., & Acuña-González, J. (2003). N2 production by the anammox reaction in the anoxic water column of Golfo Dulce, Costa Rica. Nature, 422, 606-608.

Dean, H. K. (2004). Marine biodiversity of Costa Rica; Class Polychaeta (Annelida). Revista de Biología Tropical, 52(Suppl. 2), S131-S204.

Dean, H. K. (2008). The use of polychaetes (Annelida) as indicator species of marine pollution, a review. Revista de Biología Tropical, 56(Suppl. 4), S11-S38.

Dean, H. K. (2009). Polychaetes and echiurans. In I. S. Wehrtmann, & J. Cortés (Eds.). Marine Biodiversity of Costa Rica, Central America (pp. Text: 181-191, Species List: Compact Disc pp. 122-159). Berlin: Springer + Business Media B.V.

Dolédec, S., & Chessel, D. (1991). Recent developments in linear ordination methods for environmental sciences. Advances in Ecology, India, 1, 133-155.

Epifanio, C. E., Maurer, D., & Dittel, A. I. (1983). Seasonal changes in nutrients and dissolved oxygen in the Gulf of Nicoya, a tropical estuary on the Pacific coast of Central America. Hydrobiologia, 101, 231-238.

Fernández-García, C., Cortés, J., Alvarado, J. J., & Nivia-Ruiz, J. (2012). Physical factors contributing to the benthic dominance of the algae Caulerpa sertularioides (Caulerpaceae, Chlorophyta) in the upwelling Bahía Culebra, north Pacific of Costa Rica. Revista de Biología Tropical, 60(Suppl. 2), S93-S107.

García, V., Acuña-González, J., Vargas-Zamora, J. A., & García-Céspedes, J. (2006). Calidad bacteriológica y desechos sólidos en cinco ambientes costeros de Costa Rica. Revista de Biología Tropical, 54(Supl. 1), S35-S48.

García-Céspedes, J., Acuña-González, J. A., & Vargas-Zamora, J. A. (2004). Metales traza en sedimentos de cuatro ambientes costeros de Costa Rica. Revista de Biología Tropical, 52(Supl. 2), S51-S60.

Gocke, K., Cortés, J., & Murillo, M. M. (2001a). Planktonic primary production in a tidally influenced mangrove forest on the Pacific coast of Costa Rica. Revista de Biología Tropical, 49(Suppl. 2), S279-S288.

Gocke, K., Cortés, J., & Murillo, M.M. (2001b). The annual cycle of primary productivity in a tropical estuary, The inner regions of the Golfo de Nicoya, Costa Rica. Revista de Biología Tropical, 49(Suppl. 2), S289-S306.

Gocke, K., Cortés, J., & Villalobos, C. (1990). Effects of red tides on oxygen concentration and distribution in the Gulf of Nicoya, Costa Rica. Revista de Biología Tropical, 38, 401-407.

Gravel, P., Johanning, K., McLachlan, J., Vargas, J.A., & Oberdörster, E. (2006). Imposex in the intertidal snail Thais brevidentata (Gastropoda, Muricidae) from the Pacific coast of Costa Rica. Revista de Biología Tropical, 54(Suppl. 1), S21-S26.

Hammer, Ø., Harper, D. A. T., & Ryan, P. D. (2001). PAST, Paleontological Statistics Software Package for Education and Data Analysis. Palaeontologia Electronica, 4, 1-9.

Holmes, R. W. (1970). The Secchi disk in turbid waters. Limnology, & Oceanography, 15, 688-694.

Jiménez, C. (2001). Arrecifes y ambientes coralinos de Bahía Culebra, Pacífico de Costa Rica, aspectos biológicos, económico-recreativos y de manejo. Revista de Biología Tropical, 49(Supl. 2), S215-S231.

Jiménez, C., Cortés, J., León, A., & Ruiz, E. (2001). Coral bleaching and mortality associated with the 1997-98 El Niño in an upwelling environment in the eastern Pacific (Gulf of Papagayo, Costa Rica). Bulletin of Marine Science, 69, 151-169.

Legendre, P., & Gallagher, E. D. (2001). Ecologically meaningful transformations for ordinations of species data. Oecologia, 129, 271-280.

Lizano, O. G., & Alfaro, E. (2004). Algunas características de las corrientes marinas en el Golfo de Nicoya, Costa Rica. Revista de Biología Tropical, 52(Supl. 2), S77-S94.

Mallin, M. A., Pearl, H. W., Rudek, J., & Bates, P. W. (1993). Regulation of estuarine primary production by watershed rainfall and river flow. Marine Ecology Progress Series, 93, 193-203.

Mata, A., Acuña, J. A., Murillo, M. M., & Cortés, J. (1987). La contaminación por petróleo en el Caribe de Costa Rica, 1981-1985. Caribbean Journal of Science, 23, 41-49.

Meier, P. C., & Zünd, R. E. (1993). Statistical Methods in Analytical Chemistry. New York, USA: John Wiley & Sons, Inc.

Morales-Ramírez, A., Víquez, R., Rodríguez, K., & Vargas, M. (2001). Marea roja producida por Lingulodinium polyedrum (Peridiniales, Dinophyceae) en Bahía Culebra, Golfo de Papagayo, Costa Rica. Revista de Biología Tropical, 49(Supl. 2), S19-S23.

Muller-Parker, G., & Cortés, J. (2001). Spatial distribution of light and nutrients in some coral reefs of Costa Rica during January 1997. Revista de Biología Tropical, 49(Suppl. 2), S251-S263.

Palter, J., León-Coto, S., & Ballestero, D. (2007). The distribution of nuttrients, dissolved oxygen and chlorophyll a in the upper Gulf of Nicoya, Costa Rica, a tropical estuary. Revista de Biología Tropical, 55, 427-36.

Riley, J. P., & Chester, R. (1971). Introduction to Marine Chemistry. New York, USA: Academic Press.

Rixen, T., Jiménez, C., & Cortés, J. (2012). Impact of upwelling events on the sea water chemistry in the Gulf of Papagayo (Culebra Bay), Costa Rica. Revista de Biología Tropical, 60(Suppl. 2), S187-S195.

Rojas-Figueroa, R. E., & Vargas-Zamora, J. A. (2008). Abundancia, biomasa y relaciones sedimentarias de Americonuphis reesei (Polychaeta, Onuphidae) en el Golfo de Nicoya, Costa Rica. Revista de Biología Tropical, 56(Supl. 4), S59-S82.

Rodríguez-Sáenz, K., & Morales-Ramírez, A. (2012). Composición y distribución del mesozooplancton en una zona de afloramiento costero (Bahía Culebra, Costa Rica) durante La Niña 1999 y el 2000. Revista de Biología Tropical, 60(Supl. 2), S143-S157.

Sibaja-Cordero, J. A., & Vargas-Zamora, J. A. (2006). Zonación vertical de epifauna y algas en litorales rocosos del Golfo de Nicoya, Costa Rica. Revista de Biología Tropical, 54(Supl. 1), S49-S67.

Silva, A. M., & Acuña-González, J. (2006). Caracterización físico-química de dos estuarios en la bahía de Golfito, Golfo Dulce, Pacífico de Costa Rica. Revista de Biología Tropical, 54(Supl. 1), S241-S256.

Spongberg, A. L. (2004a). PCB contamination in surface sediments in the coastal waters of Costa Rica. Revista de Biología Tropical, 52(Suppl. 2), S1-S10.

Spongberg, A. L. (2004b). PCB concentrations in sediments from the Gulf of Nicoya estuary, Pacific coast of Costa Rica. Revista de Biología Tropical, 52(Suppl. 2), S11-S22.

Spongberg, A. L. (2004c). PCB contamination in marine sediments from Golfo Dulce, Pacific of Costa Rica. Revista de Biología Tropical, 52(Suppl. 2), S23-S32.

Spongberg, A. L. (2006). PCB concentrations in intertidal sipunculans (Phylum Sipuncula) from the Pacific of Costa Rica. Revista de Biología Tropical, 51(Suppl. 1), S27-S33.

Spongberg, A., & Davis, P. (1998). Organochlorinated pesticide contaminants in Golfo Dulce, Costa Rica. Revista de Biología Tropical, 46(Suppl. 6), S111-S124.

Spongberg, A. L., & Witter, J. D. (2008). A review of PCB concentrations in tropical media, 1996-2007. Revista de Biología Tropical, 56(Suppl. 4), S1-S9.

Spongberg, A. L., Witter, J. D., Acuña, J., Vargas, J., Murillo, M., Umaña, G., … Perez, G. (2011). Reconnaissance of selected PPCP compounds in Costa Rican surface waters. Water Research, 45, 6709-6717.

Strickland, J. D. H., & Parsons, T. R. (1972). A Pratical Handbook of Sewater Analysis. Fisheries Research Board of Canada Bulletin, 167, 1-310.

Svendsen, H., Rosland, R., Myking, S., Vargas, J. A., Lizano, O. G., & Alfaro, E. J. (2006). A physical-oceanographic study of Golfo Dulce, Costa Rica. Revista de Biología Tropical, 54(Suppl. 1), S147-S170.

Valverde, R. A., Selcer, K. W., Lara, L. R., & Sibaja-Cordero, J. A. (2008). Lack of xenoestrogen-induced vitellogenin in male olive ridley sea turtles (Lepidochelys olivacea) from the Pacific coast of Costa Rica. Revista de Biología Tropical, 56(Suppl. 4), S49-S57.

Vargas, J. A. (2016). The Gulf of Nicoya estuarine ecosystem. In M. Kapelle (Ed.). Ecosystems of Costa Rica (pp. 106-124). Chicago, USA: University of Chicago Press.

Vargas, J. A., Acuña-González, J., Gómez, E., & Molina, J. (2015). Metals in coastal mollusks of Costa Rica. Revista de Biología Tropical, 63(4), 1007-1019.

Vargas, J. A., Acuña-González, J., Vásquez, F., & Sibaja-Cordero, J. A. (2016). Brachiopods, sipunculans, enteropneusts and metals from two estuarine tidal flats, Pacific, Costa Rica. Revista de Biología Tropical, 64(3), 1311-1331.

Vargas, J. A., & Dean, H. K. (2009). Sipunculans. In I. S. Wehrtmann, & J. Cortés (Eds.). Marine Biodiversity of Costa Rica, Central America (pp. Text: 175-180, Species List: Compact Disc pp. 119-121). Berlin: Springer + Business Media B.V.

Vargas, J. A., & Mata, A. (2004). Where the dry forest feeds the sea, the Gulf of Nicoya Estuary, In G. W. Frankie, A. Mata, & S. B. Vinson (Eds.). Biodiversity Conservation in Costa Rica, Learning the Lessons in a Seasonal Dry Forest (pp. 126-135). California, USA: University of California Press.

Víquez, R., & Hargraves, P. E. (1995). Annual cycle of potentially harmful dinoflagellates in the Gulf of Nicoya, Costa Rica. Bulletin of Marine Science, 57, 467-475.

Voorhis, A., Epifanio, C. E., Maurer, D., Dittel, A. I., & Vargas, J. A. (1983). The estuarine character of the Gulf of Nicoya, an embayment on the Pacific coast of Central America. Hydrobiologia, 99, 225-237.

Fig. 1. A: Location of the four coastal embayments on the Pacific and Caribbean coasts of Costa Rica, 2000-2001-2002. B. Gulf of Papagayo, Culebra Bay (10o 62’ N - 85o 66’ W) Panamá beach, Nacascolo beach, Esmeralda beach. C. Gulf of Nicoya estuary (10o 00’N - 85o 00’W): 1- Upper region, 2- Cortezas islands, 3- Caballo island, 4- Lower region. D. Northern shore of Golfo Dulce, Golfito Bay (83o 10’ N - 8o 35’ W) Bay channel, Cacao beach, Inner islet. E. Container ship port of Moín (10º 00’ N - 83º 07’ W) Green navigation buoy, Red navigation buoy, Inner Moín Port B-E. Scale bars in km.

Table 2

Golfito Bay, Culebra Bay, Moin Bay: Number of stations per parameter (N), Maximum and Minimum, Mean, Standard error, Standard deviation, Median, 25 percentile, 75 percentile, of 13 water parameters at surface (1 m depth) of a total of 60 stations. Six selected examples of stations per embayment are included, and maximun values are in bold

TABLE 2 (Continued)
	Secchi depth	S ‰	T ºC	Diss. Oxygen	% O2 saturation	Susp. matter	Chlorophyll-a	Phaeopigments	PO43-	SiO(OH)3-	NH4+	NO2-	NO3-
	Secchi depth	S ‰	T ºC	Diss. Oxygen	% O2 saturation	Susp. matter	Chlorophyll-a	Phaeopigments	PO43-	SiO(OH)3-	NH4+	NO2-	NO3-
Culebra Bay
N	11	11	11	11	11	11	11	11	11	11	11	11	11
Min	4.2	23	24.0	5.16	85.2	4.12	0.36	0.00	0.11	0.36	0.36	0.00	0.02
Max	12.6	35	32.5	8.20	141.6	17.43	2.23	0.71	0.36	4.37	0.57	0.08	0.69
Mean	7.9	30	28.7	6.89	116.5	7.52	0.95	0.21	0.19	1.02	0.51	0.01	0.10
Standard	0.8	1	0.6	0.31	5.6	1.08	0.22	0.08	0.02	0.35	0.03	0.01	0.06
Standard dev	2.7	3	2.1	1.02	18.5	3.57	0.67	0.24	0.08	1.15	0.09	0.02	0.20
Median	8.0	29	29.0	7.10	120.2	6.67	0.69	0.14	0.16	0.78	0.55	0.01	0.04
25 percentile	5.4	29	28.0	5.66	98.4	5.82	0.42	0.00	0.13	0.36	0.36	0.00	0.04
75 percentile	10.0	31	30.0	7.61	132.8	7.63	1.56	0.35	0.27	0.99	0.57	0.01	0.05
Sta 59 Panamá	12.6	30	29.0	5.52	94.1	7.63	0.36	0.00	0.16	0.36	0.57	0.01	0.04
Sta 62 Panamá	4.5	29	28.5	7.30	122.5	9.54	2.23	0.23	0.21	4.37	0.55	0.00	0.05
Sta 63 Nacascolo	5.4	31	24.0	6.98	108.1	6.90	1.61	0.14	0.14	1.06	0.55	0.00	0.05
Sta 67 Nacascolo	10.0	29	29.3	7.07	120.2	6.10	0.55	0.38	0.36	0.36	0.57	0.01	0.04
Sta 66 Esmeralda	10.0	29	29.0	7.10	120.3	4.93	0.36	0.00	0.28	0.99	0.36	0.08	0.04
Sta 64 Esmeralda	7.0	32	27.0	7.10	117.6	6.67	1.51	0.71	0.13	0.58	0.55	0.00	0.02
Golfito Bay
N	19	19	19	19	19	18	18	18	19	19	19	19	19
Min	0.9	10	28.0	6.58	111.3	6.10	0.00	0.00	0.05	1.60	0.12	0.00	0.03
Max	4.5	30	32.0	9.60	158.4	156.04	3.52	0.68	0.93	144.94	2.24	0.82	1.33
Mean	2.6	24	29.8	7.68	128.3	26.25	1.32	0.34	0.26	35.53	0.93	0.16	0.43
Standard error	0.2	1	0.3	0.24	3.4	11.09	0.23	0.06	0.05	9.15	0.18	0.07	0.10
Standard dev	1.1	6	1.2	1.06	14.7	47.03	0.96	0.26	0.20	39.90	0.78	0.29	0.46
Median	2.3	26	29.6	7.18	122.6	9.95	1.44	0.38	0.26	21.57	0.97	0.05	0.27
25 percentile	2.0	25	29.0	6.79	116.1	8.61	0.59	0.03	0.10	8.19	0.12	0.02	0.09
75 percentile	3.5	28	31.0	8.83	143.5	12.16	1.98	0.58	0.33	51.64	1.41	0.06	0.86
Sta 26 Channel	0.9	10	28.0	9.20	137.5	14.30	2.10	0.08	0.55	81.51	0.97	0.82	0.03
Sta 119 Channel	4.0	26	28.5	7.36	121.4	156.04	0.86	0.58	0.26	11.23	1.41	0.05	0.14
Sta 32 Cacao	4.5	25	31.5	9.10	158.4	8.97	2.63	0.32	0.14	8.82	2.00	0.02	0.09
Sta 35 Cacao	0.9	9	28.0	9.60	143.5	18.14	3.52	0.19	0.93	112.61	0.97	0.80	1.33
Sta 27 Islet	1.0	12	28.5	9.30	142.0	11.45	1.97	0.45	0.20	144.94	1.09	0.78	0.03
Sta 28 islet	2.0	28	31.0	6.90	120.9	9.47	0.00	0.00	0.09	69.00	0.16	0.06	1.16
Moín Bay port
N	30	30	30	27	27	30	28	28	30	30	30	30	30
Min	0.0	10	26.5	5.55	91.8	5.84	0.15	0.00	0.02	0.19	0.07	0.00	0.01
Max	12.0	35	31.5	8.29	142.5	34.26	2.44	0.88	2.91	143.32	2.24	0.40	1.63
Mean	3.3	27	28.7	6.70	111.9	11.06	0.94	0.40	0.36	25.78	0.63	0.08	0.30
Standard error	0.5	1	0.2	0.13	2.5	1.09	0.13	0.05	0.11	6.28	0.12	0.02	0.08
Standard dev	2.5	6	1.1	0.66	13.0	5.95	0.69	0.25	0.58	34.42	0.65	0.12	0.43
Median	2.8	28	28.5	6.68	113.6	9.06	0.67	0.38	0.15	17.02	0.38	0.03	0.08
25 percentile	1.7	26	28.0	6.28	101.2	7.16	0.43	0.20	0.07	0.79	0.22	0.01	0.03
75 percentile	4.0	30	29.5	7.14	121.6	13.59	1.49	0.58	0.44	34.34	0.75	0.08	0.58
Sta 44 Port	7.0	32	27.5	5.55	92.8	8.83	0.27	0.32	0.06	4.82	0.25	0.03	0.29
Sta 120 Port	5.0	35	29.5	6.50	115.0	6.99	0.29	0.48	0.05	0.79	0.38	0.14	0.64
Sta 39 Red	2.0	26	28.0	6.35	103.7	11.84	0.86	0.00	0.15	21.75	2.24	0.02	0.06
Sta 48 Red Buoy	2.0	27	27.0	5.70	91.8	13.46	0.53	0.18	0.31	23.06	0.42	0.01	0.01
Sta 46 Green Buoy	12.0	18	27.5	6.30	97.4	9.40	2.27	0.77	0.68	89.98	0.50	0.06	0.15
Sta 49 Green Buoy	2.0	26	26.5	6.00	95.1	7.95	0.84	0.57	0.22	39.39	0.42	0.01	0.01

Secchi depth (m), Salinity (psu), Temperature (ºC), Dissolved oxygen (mg/L), Oxygen saturation (%), Suspended matter (mg/L), Chlorophyll-a (mg/m3), Phaeopigments (mg/m3), Phosphate (µmol/L), Silicate (µmol/L), Ammonium (µmol/L), Nitrite (µmol/L) and Nitrate (µmol/L).

Table 3

A. Gulf of Nicoya estuary: Number of stations per parameter (N), Maximum and Minimum, Mean, Standard error, Standard deviation, Median, 25 percentil, 75 percentil,

of 13 water parameters at surface (1 m depth) water stations in four regions (Upper Gulf, Caballo island, Cortezas islands, Lower Gulf), for a total of 68 stations.

Pacific coast of Costa Rica. 2000-2002. B. Examples of stations with a parameter maximum, except salinity minimum (Sta. 102).

C. Examples of stations with a parameter minimum, except phosphate maximum (Sta. 112)

TABLE 3 (Continued)
	Secchi depth	S ‰	T ºC	Diss. Oxygen	% O2 saturation	Susp. matter	Chlorophyll-a	Phaeopigments	PO43-	SiO(OH)3-	NH4+	NO2-	NO3-
A	Secchi depth	S ‰	T ºC	Diss. Oxygen	% O2 saturation	Susp. matter	Chlorophyll-a	Phaeopigments	PO43-	SiO(OH)3-	NH4+	NO2-	NO3-
Upper Gulf
N	24	24	24	24	24	24	24	24	24	24	24	24	24
Min	0.1	5	24.8	3.60	55.8	12.26	0.78	0.00	0.34	0.99	0.03	0.00	0.05
Max	1.1	28	31.1	10.50	180.4	107.90	15.60	6.37	9.53	156.11	1.45	2.47	4.06
Mean	0.6	23	29.5	6.69	110.7	35.63	5.78	1.65	2.10	42.55	0.52	0.32	0.69
Standard error	0.0	1	0.2	0.34	5.9	4.62	0.80	0.29	0.48	7.29	0.09	0.11	0.21
Standard dev	0.2	5	1.2	1.68	28.7	22.63	3.93	1.43	2.36	35.72	0.43	0.53	1.02
Median	0.6	24	29.5	6.90	111.3	30.36	4.96	1.50	1.45	33.64	0.36	0.17	0.28
25 percentile	0.5	22	29.2	5.35	92.2	18.06	2.79	0.53	1.05	20.50	0.26	0.02	0.11
75 percentile	0.8	25	30.0	7.55	124.6	48.77	7.77	2.23	1.70	61.42	0.77	0.32	0.87
Caballo island
N	6	6	6	6	6	6	6	6	6	6	6	6	6
Min	0.8	23	25.5	5.15	82.7	15.80	2.29	0.43	0.54	0.18	0.03	0.00	0.06
Max	1.1	32	28.9	11.10	185.2	32.95	9.86	5.77	1.21	36.52	0.83	0.33	3.47
Mean	0.9	27	27.9	7.36	120.8	24.94	6.60	2.08	0.88	9.30	0.50	0.13	0.89
Standard error	0.1	1	0.5	0.81	14.0	3.23	1.25	0.79	0.12	5.59	0.12	0.06	0.54
Standard dev	0.1	3	1.2	1.99	34.3	7.23	3.05	1.92	0.29	13.68	0.28	0.14	1.33
Median	0.9	27	28.4	7.08	113.3	28.15	7.48	1.68	0.92	4.71	0.54	0.10	0.27
25 percentile	0.8	25	27.4	6.20	102.4	17.39	3.18	0.66	0.55	0.79	0.28	0.00	0.08
75 percentile	1.0	29	28.6	8.16	137.4	30.90	9.23	3.09	1.15	15.64	0.74	0.27	1.75
Cortezas islands
N	27	27	27	27	27	26	27	27	27	27	27	27	27
Min	0.8	24	26.6	4.01	66.2	9.52	0.86	0.03	0.28	0.18	0.03	0.00	0.03
Max	2.0	31	30.0	11.40	188.3	51.55	16.21	4.20	4.10	38.54	2.66	0.88	4.59
Mean	1.2	28	28.5	6.19	103.2	23.08	5.64	1.49	1.17	13.59	0.58	0.32	1.68
Standard error	0.1	0	0.1	0.30	5.1	2.48	0.65	0.17	0.18	2.76	0.10	0.06	0.27
Standard dev	0.4	2	0.7	1.57	26.3	12.65	3.36	0.89	0.93	14.36	0.53	0.31	1.41
Median	1.1	28	28.5	6.25	105.2	18.23	5.23	1.29	0.93	7.46	0.50	0.18	1.34
25 percentile	0.9	26	28.0	5.20	87.4	13.17	4.51	0.90	0.58	0.26	0.26	0.02	0.54
75 percentile	1.4	30	29.0	7.00	116.2	31.56	5.83	2.10	1.57	27.15	0.71	0.59	3.00
Lower region
N	11	11	11	11	11	11	11	11	11	11	11	11	11
Min	0.6	23	27.8	5.05	81.8	7.31	2.47	0.00	0.26	0.18	0.17	0.00	0.03
Max	2.0	35	32.0	7.50	128.8	60.97	5.00	4.81	3.40	50.07	0.92	0.26	4.18
Mean	1.4	30	29.2	6.21	106.5	32.49	3.71	1.23	1.05	9.88	0.51	0.06	1.43
Standard error	0.2	1	0.4	0.23	4.2	6.20	0.28	0.40	0.27	4.30	0.09	0.03	0.41
Standard dev	0.5	4	1.4	0.77	14.0	20.55	0.93	1.34	0.91	14.25	0.30	0.08	1.36
Median	1.5	30	29.0	6.25	109.5	37.77	3.89	0.72	0.69	4.89	0.36	0.02	1.38
25 percentile	1.0	28	28.0	5.50	100.9	9.55	2.62	0.43	0.49	0.99	0.30	0.00	0.31
75 percentile	2.0	34	29.5	6.79	113.5	51.25	4.79	1.74	1.41	11.60	0.92	0.13	2.32
B
Station 05 Cortezas	1.1	26	29.5	6.25	105.2	17.03	1.39	1.52	0.82	7.46	2.66	0.28	0.56
Station 06 Cortezas	2.0	28	28.0	11.40	188.3	15.49	5.23	1.30	0.28	6.46	0.57	0.18	3.59
Station 10 Caballo	1.0	28	28.5	11.10	185.2	18.97	7.83	5.77	0.54	4.99	0.57	0.18	0.40
Station 15 Cortezas	1.5	25	28.2	7.20	117.5	18.71	16.21	3.45	0.66	0.26	0.09	0.06	0.73
Station 73 Cortezas	0.8	27	28.0	4.30	70.6	46.47	4.45	1.26	1.36	27.15	0.71	0.70	4.59
Station 92 Jicaral	0.4	27	31.1	5.20	90.8	107.9	15.6	6.37	2.15	31.38	0.33	0.02	4.06
Station 102 T. R. mouth	0.5	5	24.8	7.40	100.5	46.38	1.74	2.29	0.36	156.11	0.03	0.02	0.22
Station 111 T. R. mouth	0.5	20	29.1	4.00	64.6	34.86	4.23	1.45	7.65	3.93	1.45	2.47	1.21
C
Station 01 Cortezas	1.5	28	30.0	6.60	113.4	9.52	5.83	0.95	0.44	9.47	0.17	0.00	0.03
Station 11 Caballo	0.8	26	28.9	6.55	108.9	15.80	9.86	1.59	0.55	0.18	0.50	0.00	1.17
Station 81 Puntarenas	1.0	33	32.0	5.50	100.9	56.94	2.62	0.00	1.06	2.68	0.30	0.03	1.72
Station 84 San Lucas	2.0	28	29.5	6.25	106.4	9.55	3.89	0.42	0.54	11.60	0.17	0.02	0.03
Station 91 Puntarenas	2.0	35	27.8	6.41	109.5	37.77	3.35	0.73	0.49	0.99	0.36	0.00	1.38
Station 99 Cortezas	0.8	27	28.5	7.00	116.2	28.15	2.29	2.20	0.78	4.42	0.03	0.02	0.13
Station 101 W. Chira	0.5	24	28.2	7.60	123.3	12.60	0.78	0.60	1.52	38.39	0.03	0.02	0.36
Station 112 T. R. mouth	0.1	13	29.0	3.60	55.8	49.57	3.25	1.04	9.53	19.60	1.45	0.17	0.44

Secchi depth (m), Salinity (psu), Temperature (ºC), Dissolved oxygen (mg/L), Oxygen saturation (%), Suspended matter (mg/L), Chlorophyll-a (mg/m3), Phaeopigments (mg/m3), Phosphate (µmol/L), Silicate (µmol/L), Ammonium (µmol/L), Nitrite (µmol/L) and Nitrate (µmol/L).

Table 4

Gulf of Nicoya estuary: Results of the Mann-Whitney U test comparison of medians of stations sampled during the dry and rainy seasons in the Upper and Middle Gulf regions. Pacific coast of Costa Rica 2000-2002

Parameter	Season	N	Min	Max	Median	U test	p
Secchi	Dry	25	0.4	2	1	280.50	0.7013
	Rainy	24	0.1	1.5	1
Salinity	Dry	25	22	32	27	157.50	0.004
	Rainy	24	13	31	24
Temperature	Dry	25	26.6	31	28.5	215.00	0.089
	Rainy	24	28	31.1	29.05
Oxygen	Dry	25	5.2	11.4	7.18	93.50	<0.001
	Rainy	24	3.6	7.8	5.27
% O2 sat.	Dry	25	87.4	188.3	118	92.50	<0.001
	Rainy	24	55.8	131	89.5
Susp. matter	Dry	22	10.33	68.22	19.32	196.00	0.138
	Rainy	24	9.52	107.89	27.89
Chlorophyll-a	Dry	25	0.78	16.21	5.59	283.00	0.741
	Rainy	24	3.03	15.55	5.49
Phaeopigm.	Dry	25	0	5.77	1.29	240.00	0.234
	Rainy	24	0.43	6.37	1.485
PO43-	Dry	25	0.28	1.82	0.78	126.50	<0.001
	Rainy	24	0.44	9.53	1.29
SiO(OH)3-	Dry	25	0.18	63.64	4.42	114.00	<0.001
	Rainy	24	3.93	107.79	34.06
NH4+	Dry	25	0.03	2.66	0.5	225.00	0.134
	Rainy	24	0.17	1.45	0.71
NO3-	Dry	25	0.06	3.59	0.4	224.50	0.133
	Rainy	24	0.03	4.59	0.72
NO2-	Dry	25	0	0.51	0.08	132.50	<0.001
	Rainy	24	0	2.47	0.56

T. River mouth: dry (88, 101, 109), rainy (94, 97, 112) W. Chira Island: dry (87, 103, 105, 108), rainy (82, 93, 96, 111) Cortezas islands: dry (04, 05, 06, 13, 14, 15, 16, 69, 71,72, 90, 99), rainy (01, 02, 03, 07, 17, 18, 73, 76, 77, 95, 113, 114) Caballo Island: dry (08, 10, 12), rainy (19) Jicaral: dry (86, 104, 107), rainy (92, 110, 95). Total dry = 25. Total rainy = 24.

Table 5

The four embayments: Estimates of Net Primary Productivity (grams of Carbon / unit area / unit time)

by the light vs dark BOD bottle oxygen change method at: A. Gulf of Nicoya, B. Golfito Bay, C. Culebra Bay,

D. Moín Bay. E. Spearman rho rank correlation between Net PP / day and 13 parameters. 2000-2002. Costa Rica

Locality	Date	Incubation	Secchi	T ºC	Chlorophyll-a	SiO(OH)3-	Net PP / day	Net PP / year
A
Cortezas islands	Mar. 22-00	14:30-17:30	1.5	28.2	16.2	0.26	4.381	1 601
Cortezas islands	Apr. 12-00	11:20-14:15	2.0	28.6	5.59	0.18	3.261	1 190
Cortezas islands	Apr. 12-00	11:55-14:15	0.8	27.7	5.66	0.18	2.658	970
Puntarenas	Apr. 11-00	11:40-14:00	2.0	28.0	5.00	0.18	1.597	583
Cortezas islands	Sep. 01-00	11:45-14:05	1.1	29.0	6.51	35.18	1.673	611
Caballo Island	Sep. 01-00	11:20-13:50	1.0	28.0	9.02	36.52	1.298	474
Puntarenas	Aug. 31-00	12:03-14:30	1.0	28.0	4.79	15.87	0.524	191
Cortezas islands	Apr. 06-01	09:40-12:40	1.5	28.0	5.23	0.99	2.047	747
Caballo Island	Apr. 05-01	08:55-12:00	1.0	28.2	7.13	0.99	0.953	348
Puntarenas	Apr. 04-01	09:00-12:00	2.0	27.8	3.35	0.99	1.253	457
B
Golfito islet	May 10-00	09:22-12:00	2.8	29.6	0.65	6.91	Negative	Negative
Golfito Cacao B.	May 10-00	11:00-13:30	2.5	29.5	1.69	8.63	0.187	68
Golfito channel	May 10-00	13:00-15:20	2.2	29.5	1.71	6.49	1.905	695
Golfito islet	Nov. 24-00	11:53-14:18	2.2	29.0	1.51	32.11	1.189	434
Golfito Cacao B.	Nov. 23-00	11:37-14:00	4.0	28.5	0.86	11.23	0.171	62
Golfito channel	Nov. 22-00	12:26-15:00	4.0	29.0	0.63	21.57	3.757	1 371
Golfito islet	Feb. 07-02	09:22-11:42	2.3	32.0	0.93	14.62	Negative	Negative
C
Culebra Panamá	Jun. 28-00	11:10-13:30	12.6	29.0	0.36	0.36	Negative	Negative
Culebra Nacascolo	Jun. 28-00	13:00-15:40	10.0	29.3	0.55	0.36	3.337	1 218
Culebra Esmeralda	Jun. 29-00	08:42-11:42	9.0	28.0	0.48	0.36	6.044	2 206
Culebra Panamá	Oct. 19-00	08:00-11:00	4.5	28.5	2.23	4.37	Negative	Negative
Culebra Nacascolo	Oct. 18-00	12:00-14:20	4.2	28.5	1.61	0.78	0.330	120
Culebra Esmeralda	Oct. 19-00	08:37-11:38	7.0	27.0	1.51	0.58	2.884	1 053
Culebra Esmeralda	May. 24-01	09:12-12:12	10.0	29.0	0.36	0.99	2.203	804
D
Moín port	May. 25-00	12:25-14:45	3.3	29.0	0.21	0.79	1.227	448
Moín Red buoy	May. 25-00	12:02-14:16	5.0	29.5	0.20	0.79	Negative	Negative
Moín Green buoy	May. 25-00	10:50-13:10	6.0	29.5	0.15	0.79	0.455	166
Moín port	Oct. 05-00	12:00-14:28	3.0	29.5	1.52	0.79	3.015	1 100
Moín Green buoy	Oct. 05-00	08:40-11:10	4.1	30.0	1.41	5.86	0.081	30
Moín Red buoy	May. 10-01	09:25-12:25	6.5	29.0	0.58	26.13	Negative	Negative
Moín port	May. 10-01	11:45-15:00	3.0	30.0	0.59	29.23	Negative	Negative

E
Net PP / day	Secchi dept	Salinity	T ºC	O2	% Sat.	Susp. Matt.	Chlorophyll	Phaeopigmets	PO43-	SiO(OH)3-	NH4+	NO2-	NO3-
rs	-0.10	-0.04	-0.44	-0.07	-0.12	0.04	0.21	0.35	0.31	-0.41	0.10	0.02	0.28
p	0.585	0.850	0.013	0.702	0.508	0.820	0.254	0.057	0.094	0.022	0.578	0.905	0.131

Locality. Date (Month-Day-Year), Incubation time interval, Secchi depth (m), Chlorophyll-a (mg/m3), Silicate (µmol/L), Net PP / day (gC/m2/day), Net PP / year (gC/m2/year). Negative = oxygen consumption greater than production.