Bio-Optical Characteristics of the Northern Gulf of California during June 2008

Bio-optical variables in theNorthernGulf of California were analyzed using in situ and satellite data obtained during a cruise in June 2008. The study area was divided into three bio-optical regions: Upper Gulf (UG), Northern Gulf (NG), and Great Isles (GI). Each region was characterized according to phytoplankton pigment concentration, phytoplankton and nonpigmented material spectral absorption coefficients, and spectral reflectance. Observed patterns were an indication of the shift in bio-optical conditions from north to south going from turbid and eutrophic waters tomesotrophic ones. Although there was a good agreement between satellite and in situChla (RMSE±33%), an overestimation of in situChlawas observed.Thiswas partly explained by the presence of nonalgal particles, as well as the influence of desert and continental aerosols, which is generally overcorrected in the standard processing. The UG and NG could be considered as Case 2 waters, but they did exhibit different bio-optical characteristics. This implies that both biological and optical properties should be invoked to better understand water reflectance variability in the study region and its implications for the remote sensing of Chla and biogeochemical processes.


Introduction
Observing the optical properties of the marine environment is crucial for the evaluation of satellite ocean color products and the development of spectral deconvolution techniques. Coastal zones are of particular interest, since they are influenced by mineral sediments, brought by rivers or stirred up from the bottom, as well as by colored dissolved organic matter of terrestrial origin (CDOM). In these zones, however, the algorithms (of the blue-to-green ratio type), developed for open ocean waters, generally fail, with the result of strongly overestimated chlorophyll-a concentration (Chla) [1].
The Gulf of California is a subtropical semienclosed sea characterized by a cyclonic circulation in summer and an anticyclonic circulation in winter and by a seasonally reversible upwelling system with moderate upwelling along the west coast during summer and intense upwelling along the east coast during winter [2,3]. In the northern part of the Gulf, that is, from the Colorado River Delta to 28 ∘ N (Figure 1(a)), these characteristics and also the bathymetry result in an eddy in the center of the basin and a coastal current on the mainland shelf. Direct observations show that the central eddy is ∼150 m deep, cyclonic from June to September, and anticyclonic from November to April [4,5]. Satellite observations have confirmed the cyclonic pattern during summer, but at the same time they have revealed a more complex circulation pattern with numerous plumes and eddies (some anticyclonic), traveling from coast to coast, enhancing the exchange of suspended material [6].
The Northern Gulf is mostly less than 200 m deep and tidal mixing is very important [7], especially around the sills of the midriff archipelago region (Great Isles, GI) and in 2 International Journal of Oceanography the shallow Upper Gulf (UG, north of 34 ∘ 45 ). The GI region and the sills are an area of continuous upward pumping of nutrients by tidal mixing [8,9], while in the UG region tidal mixing has the same effect. These complex circulation patterns support very high primary production rates [8] and a high-diversity, high abundance fish fauna of commercial and ecological importance [10][11][12].
In the UG region tidal mixing is also responsible for the resuspension of inorganic material generating a zone of high turbidity [11], which can eventually be observed in truecolor imagery of the area. For this reason these waters were classified as Case 2 while those to the south were classified as Case 1 waters [13]. However, the complex circulation patterns mentioned above may lead to a high spatial variability of such properties. In addition, the area studied by Pegau et al. [13] was located mostly to the south of the GI region and few stations covered the northern area.
In this context, the hypothesis is that there should be at least two distinct bio-optical regimes in the Northern Gulf of California that respond to the strong oceanographic heterogeneity of this zone. Based on the measurement of some surface bio-optical variables carried out in June 2008, the present study aims to (i) identify different bio-optical regions; (ii) evaluate the differences in spectral reflectance and spectral particulate absorption among identified biooptical regions; (iii) evaluate the role of different phytoplankton communities among regions; and (iv) quantify the influence of these differences on the remote sensing of Chla.

Methodology
2.1. Sampling Procedure. During June 2 to 17, 2008, 48 stations were made in the Northern Gulf of California, from 28 ∘ to 32 ∘ N and 112 to 115 ∘ W (Figure 1(a)). Water transparency was measured using Secchi disk from which we calculated the attenuation coefficient (Kd) and the first optical depth (1OD) [14]. Surface water samples (∼0.40 m depth) were taken using Niskin bottles to measure phytoplankton pigments concentration and light absorption coefficient by particulate material.
In some stations, sampled between 10 a.m. and 3 p.m. and with clear sky conditions, marine reflectance ( ( )) was measured using a SIMBADA radiometer, an improved version of the SIMBAD radiometer [15] that measures in spectral bands centered at wavelengths of 350, 380, 412, 443, 490, 510, 560, 620, 670, 750, and 870 nm [16]. The instrument was calibrated in the laboratory using an integrating sphere and in the field using a Spectralon plaque. Procedures described in Bécu [16] were applied to process the data into marine reflectance.
For pigment analysis, one to two liters (depending on station location) of seawater was filtered through 25 mm diameter Whatman GF/F filters with positive pressure. Filters were wrapped in aluminum foil and stored in liquid nitrogen until laboratory analysis. One or two liters of seawater was also filtered through Whatman GF/F filters for measurement of particulate light absorption. Filters were placed in Histoprep tissue caps and stored in liquid nitrogen until laboratory analysis.

Pigment Concentration.
Pigments were extracted in 100% acetone for 24 h with trans-beta-apo-carotene as internal standard in a freezer at −20 ∘ C, with a previous threeminute sonication. Pigment extract was filtered with 0.2 m acrodisc filters and injected into a Varian HPLC system with a degasser and Adsorbosphere C8 column (150 × 4.6 mm, 3.5 m particle size). A two-solvent gradient, methanol and ammonium acetate (70 : 30 v : v), was used following the methodology described by Barlow et al. [17]. The HPLC pigment standards were purchased from Sigma-Aldrich and DHI. They were quantified with a spectrophotometer using published extinction coefficients [18] and used to identify pigment peaks and calibrate pigment concentrations based on the peak areas.
Pigment information was used to calculate phytoplankton size fractions (micro, nano, and picoplankton) according to Uitz et al. [19].

Particulate Absorption Coefficient.
In the laboratory, filters were saturated with filtered seawater, which had been irradiated with UV lamps (25 W) and the optical density (OD) was measured on a Perkin-Elmer Lambda 10 spectrophotometer with integrating sphere, following the procedure in Mitchell et al. [20]. OD was measured from 400 to 750 nm with a resolution of 1 nm before and after rinsing the filters with hot methanol [20] for 15 min twice. Absorption coefficients for a particles and not pigmented particles, that is, ( ) and ( ), were determined using the equation: where is the filter clearance area, is the filtered volume, and 0.4068 and 0.368 are the coefficients to correct the increase in path length caused by multiple scattering in the glass-fiber filter, which were previously determined for this spectrophotometer. Phytoplankton absorption ℎ ( ) was determined by the difference between absorption by total particulate matter, ( ), and absorption by nonpigmented material, ( ). The specific absorption coefficient by phytoplankton ( * ℎ ( )), with units m 2 (mg Chla) −1 , was obtained by dividing ℎ ( ) (m −1 ) by Chla (mg m −3 ) measured by HPLC. To compare curve shapes, normalized absorption (A, dimensionless) was plotted, which was calculated by normalizing the spectrum ℎ ( ) to the absorption maxima in the blue (around 440 nm).

Satellite Images.
Sea-viewing Wide Field-of-View Sensor (SeaWiFS) and Moderate resolution Imaging Spectroradiometer (MODIS/Aqua) images were processed at 1 km spatial resolution for the entire study period using SeaDAS V5.5 to obtain Chl SAT . For comparison between the in situ data and satellite products and to minimize georeference errors, median values were computed for all products in a 3-by-3 pixel window centered on the locations of the oceanographic stations. Furthermore, the coefficient of variation (CV) was computed for the Chla within each 3-by-3 window, and the retrieved value was excluded if CV > 0.2. This process was carried out to avoid strong variation from nonhomogeneous regions within each window. For the temporal threshold coincidence between the satellite and in situ measurements, a ±3 hour window around the satellite overpass was considered, which follows the NASA criterion [21]. The calculation of Chl SAT from SeaWiFS data is based on OC4V4 algorithm while from MODIS data it is based on OC3 algorithm [22], and Chl SAT was evaluated against in situ Chla. In the comparisons, SeaWiFS and in situ Chla ( = 26), MODIS and in situ Chla ( = 22), and the average Chla from SeaWiFS and MODIS and in situ Chla ( = 28) were used.
The OC4V4 algorithm was applied to SIMBADA data considering that SIMBADA has a band centered at 560 nm while SeaWiFS is 555 nm; that is, the effects of different SeaWiFS and SIMBADA spectral bands were ignored. These effects are expected to be much smaller than algorithm performance.

Bio-Optical
Regions. The 1OD was used as an index to define patterns in the optical properties of waters from the study area. Three regions were identified (Figure 1(b)) and named according to their geographical location. The UG region (1OD < 3 m) included stations located in the shallow Upper Gulf. The Northern Gulf region or NG (1OD > 5 m) included stations mostly located in the central part of the study area, although there were three stations in the southernmost zone. Finally, the GI region (3 m < 1OD < 5 m) had stations located in the Great Isles zone. Five stations to the north and between the UG and NG regions could be classified as belonging to the GI region considering its 1OD values, but we decided to consider this area as a transitional zone between the UG and NG. The selection of the optical depths (3 and 5 m) was based on a cluster and a stepwise discriminant analysis. The cluster analysis divided stations into groups and the second was used to confirm which station belongs to each group.

In Situ Chla and
Pigments. In situ Chla shows a strong spatial variability with values from 0.53 to 3.28 mg/m 3 (Figure 2(a)), with the lowest values on station 36 in the NG and the highest ones on station 22 in the UG region. Most stations were dominated by microplankton size fraction ( Figure 2(b)), with exception of station 27 where picoplankton contributed by more than 80% to Chla.
Some phytoplankton pigments can be used as markers of certain microalgae groups [17]. We used pigment composition to evaluate phytoplankton community composition according to Table 1, modified from Jeffrey et al. [18].
The UG region presented the highest concentration of the pigment Fucoxanthin, indicating the presence of diatoms, followed by Chlorophyll b, Hexanoyloxyfucoxanthin, Zeaxanthin, and Alloxanthin. The high ratio of Fucoxanthin/Chla coincides with Chla higher than 1.5 mg/m 3 and the highest Chla of the entire area was observed on station 22 (Figure 2(a)), where this ratio was almost one. Santamariadel-Angel et al. [23] mentioned that diatoms are the dominant phytoplankton group in this region, and also that Chla/cell is higher due to their acclimation to the high turbidity characteristic of this zone.
The transition zone corresponds to stations 16, 17, 18, 20, and 21 ( Figure 2). From these we only obtained pigment data  Figure 2: (a) The ratio pigment : Chla, expressed in percent, is indicated in vertical bars (left axis). The bar color corresponds to some specific pigment as indicated in the legend. Pigment abbreviation is indicated in The NG regions were characterized by Chla lower than 1 mg/m 3 , except station 32 located close to the east coast of the Gulf. In addition, there is a general increase of Zeaxanthin, Hexanoyloxyfucoxanthin, and Peridinin. At station 27 Zeaxanthin/Chla ratio was almost 0.9, while at station 36 Neoxanthin/Chla was high (∼0.5).
In region GI Chla is higher although values are lower than in region UG, and Zeaxanthin is barely observed. Hexanoyloxyfucoxanthin is found in increased proportion, while the Peridinin ratio is up to 0.7 at station 48 ( Figure 2).   (Figures 3(a), 3(c), and 3(e)). The individual contribution of nonpigmented material and phytoplankton for light absorption was also analyzed ( Figure 4) and a gradual decrease of (440) contribution from the UG to the GI with average values of 56% in the UG region, 41% in the NG region, and 18% in the GI region was observed.

Light Absorption
The shape of the absorption spectra ( Figure 5) has the two Chla absorption maxima around 440 and 675 nm and an additional broad maximum at wavelengths between 450 and 500 nm, which is the highest in the NG region and slightly lower in the GI region.
The absorption coefficients ℎ (440) and ℎ (675) were compared to Chla and a linear fitting between log-transformed data provided a low 2 of 0.24 for ℎ (440) ( = 30) and 0.33 for ℎ (675). When compared with the fit proposed by Bricaud et al. [24] (data not shown), a lot of dispersion is observed, especially at 440 nm. The chlorophyll-specific phytoplankton absorption at 440 nm [ * ℎ (440)] decreased with Chla and ranged from 0.017 to 0.155 m 2 mg Chla −1 , while the chlorophyll-specific phytoplankton absorption at 675 nm [ * ℎ (675)] also decreased varying from 0.006 to 0.065 m 2 mg Chla −1 . The general tendency is a decrease with the Chla increase, as expected, although with a high dispersion as observed above for ℎ (440) and ℎ (675) (Figure 6).
Light absorption by nonpigmented material [ (440)] did not covary with Chla (Figure 7(a)) with the exception of the UG where an inverse relationship was observed. This indicates that in this region stations with high Chla had generally less inorganic particles. A regression line was adjusted to these data and the resulting equation ( 2 = 86%, 95% confidence) was The attenuation coefficient was also compared to Chla (Figure 7(b)) and an inverse relationship was also observed for UG data, while a positive but weak one was observed for the other sites. A regression line was adjusted for those UG data and the resulting equation ( 2 = 65%, 95% confidence) was = 1.93 − 0.573Chl . (3)

Variability of Reflectance Spectra.
The reflectance spectra variability is indicated in Figure 8. Unfortunately, in the UG region we only have data for station 15, which has the highest ( ), especially at 560 nm, that can be related to high concentrations of total suspended matter. The NG region was International Journal of Oceanography the one with the highest number of data, 12 stations, while there were only data at two stations for region GI (stations 47 and 48). Due to variability in water composition, no distinct difference is apparent between the reflectance spectra of the NG and GI regions, which therefore can hardly be classified based on marine reflectance.

Chl
and Bio-Optical Algorithm Evaluation. Chl SAT images during the study period demonstrated a strong variability in the study area (Figure 9) with the appreciation of several structures like eddies, filaments, and meanders. In general Chl SAT is higher to the north related to the UG zone, to the south around the Great Isles region (GI), and along the west side of the NG. However, a decrease of Chl SAT is noted in this last zone, especially after June 9.
The performance of the empirical algorithms OC4V4 and OC3 was evaluated through the comparison of Chl SAT with in situ Chla (Figures 10(a) and 10(b)). For SeaWiFS Chl SAT (Figure 10(a)), the calculated RMSE was 36% while the BIAS was 28% ( = 26). For MODIS Chl SAT (Figure 10 RMSE and BIAS were similar, that is, 29% and a 22% ( = 22), respectively. Finally, we calculated the average between MODIS and SeaWiFS Chla and these data were compared to in situ Chla, which resulted in a RMSE of 33% and a BIAS of 26% ( = 28). The satellite and in situ maximum reflectance ratios (MRRs) were also compared ( Figure 10(c)). They were calculated using the highest reflectance value in the blue ( (443), (490), or (510)) divided by (555) or (560). Results show that satellite values are lower than the SIMBADA values, explaining the higher (overestimated) satellite estimates of Chla (Figures 10(a) and 10(b)). The lower satellite values are due to lower reflectance in the blue.
3.6. Statistical Analysis. The differences among the proposed bio-optical regions were firstly examined through the comparison of ( ), ℎ ( ), and the shape A, using a nonparametric Kruskall-Wallis test for independent samples. Results indicated that regions UG, NG and GI are statistically different (95% confidence level) from each other. However, when comparing the reflectance spectra using the same test the NG and GI regions were not significantly different, although both are different from the UG. This was mentioned qualitatively in Section 3.4. The differences among regions are summarized in Table 2.

Bio-Optical
Regions. In situ Chla indicates the strong spatial variability of the study area, which agrees with the proposed bio-optical classification, that is the highest values in the UG, a decrease in the NG, and a new increase in the GI (Figure 2). The difference between the UG and NG is especially remarkable, as previously noted by other studies based on some hydrographic and biological aspects [10,11]. These changes in Chla are also linked to changes in phytoplankton community structure (size and composition) which could be inferred from pigment composition [19,25].
In this study, the dominance of microplankton could be observed in the entire study area, with a slight increase in the proportion of picoplankton in the NG region (Figure 2(b)), which is characteristic of zones with more oligotrophic waters [24,26]. The index used to estimate size fraction [19] considers Fucoxanthin and Peridinin to infer the proportion of microplankton, which indicates the presence of diatoms and/or dinoflagellates, respectively (Table 1). These pigments were dominant in our data (Figure 2(a)) but there was observed a higher proportion of Fucoxanthin in the UG while in the NG and GI region Peridinin increased its proportion. It is also important to observe the increase in the pigment Zeaxanthin in the NG, which attains its maximum value at station 27, which results in the dominance of picoplankton in this station (Figure 2(b)). In a natural environment, shifts in the hydrological features of the water column are related to changes in nutrient concentration, which in turn is associated with the modification of phytoplankton composition and their mean size [19]. As previously mentioned, diatoms are the dominant group in the UG region [23] and a high biomass of this group has been also reported in the GI region [27,28] and confirmed in our data.
Larger cells suffer with a higher package effect and lower specific absorption coefficients [24]. It is assumed that higher Chla are related to larger cells [24] and as a consequence an inverse relationship between Chla and the specific absorption coefficients * ℎ (440) and * ℎ (675) is expected [26]. In our data Chla varied between 0.53 and 3.28 mg/m 3 (Figure 2(a)) but there was almost no change in size fraction (Figure 2(b)) which can explain the high dispersion in the observed relationship ( Figure 6). This dispersion was especially high for the UG data and at 440 nm. Millán-Núñez et al. [27] also observed more variability at 440 nm than at 675 nm when comparing * ℎ ( ) with Chla for some stations in our study area. They attributed this pattern to chlorophyll per cell volume which was not independent of cell size, as assumed in models that explain the inverse relationship between * ℎ ( ) and cell size [24]. In addition, the inverse relationship between these parameters may not necessarily hold in coastal waters because of events such as upwelling, advection, fronts, and freshwater inputs, as well as the complexity of light field conditions [29,30], processes that are frequent in our study area. In particular the UG is a region with very complex hydrodynamics, high concentrations of suspended matter, and consequently very peculiar light conditions [11,23] in comparison with the southernmost regions. The concentration of suspended matter inferred from (440) had the highest contribution to light absorption in the UG, while it had the lowest contribution in the GI region ( Figure 4). Taking into account the observed differences among the three bio-optical regions, we evaluated the relationship between  (440) and Chla for each one individually given that for the UG a power low fit resulted in a 2 of 9%, while for the NG it was 53% ( * ℎ (440) = 0.0472Chla −1.006 ) and for the GI it was 68% ( * ℎ (440) = 0.0795Chla −0.631 ). This pattern is an indication of the shift in bio-optical conditions from north to south going from turbid and eutrophic to mesotrophic waters. In addition, different associations were observed between variables in the UG; an inverse relationship between (440) and Chla, and between and Chla (Figure 7). An inverse relationship between (440) and Chla is an indication that Chla increases when nonpigmented particles are found in low concentration, that is, in clearer waters. Besides, the tendency observed for can be explained by the high concentration of nonalgal and mineral particles in the UG (Figures 3 and 4).
Phytoplankton absorption spectra also show important differences among regions with a tendency to lower values in the NG (Figure 3(d)) and at most stations of the UG (Figure 3(b)). In this region, only stations 9 and 15 presented values higher than 0.1 m −1 . In the GI region (Figure 3(f)), only station 47, located to the north and in the limit with the NG, presented values lower than 0.1 m −1 . At this station the proportion of microplankton was the highest (Figure 2(b)), while at stations 5 and 6 the proportion of nanoplankton increased. When we analyzed the average shape of the absorption spectra ( Figure 5), the differences among regions were more evident, and the absorption increase around 480 nm from the UG to the NG and a new decrease in the GI region are remarkable. Indeed, the significance of these differences (magnitude and shape) was statistically tested and confirmed. An absorption increase at wavelengths between 440 and 550 nm indicates a higher proportion of carotenoids: Chla or photoprotective pigments [19,23,30]. A low concentration of photoprotective pigments can be inferred in the UG region as an indication of low-light conditions associated with the high turbidity.
The reflectance spectra showed a small variability in the study area (Figure 8) with the exception of station 15, the only one in the UG, which has a strong signal at 560 nm related to high scattering by nonpigmented and mineral particles [30], which was confirmed by high nonpigmented particles absorption values (Figure 3(a) and Table 2). The statistical analysis confirmed the observed difference between the UG station and the others.

Chl
and Bio-Optical Algorithm Evaluation. The large spatial and temporal variability of the study area is manifested in Chl SAT (Figure 9) with gradients associated with the defined bio-optical regions. In general, Chl SAT is higher to the north related to UG region, decreases in the NG region, and increases in the GI region. However, it has to be noted that Chl SAT decreased in this last zone, especially after June 13 when the extension of the central oligotrophic region increased. In fact, a frontal zone is clearly visible in the images (Figure 9), as also noted by Sánchez-Velasco et al. [11], in which are located those stations classified as in Table 2: Average (minimum and maximum in parenthesis) values for each bio-optical region of first optical depth (1OD), chlorophyllconcentration (Chla), phytoplankton light absorption coefficient at 440 nm ( ℎ (440)) and 675 nm ( ℎ (675)), nonpigmented particles light absorption coefficient at 440 nm ( ℎ (440)), percent contribution of nonpigmented particles to light absorption (% ), reflectance at 560 nm (R(560)), and maximum reflectance ratio (MRR). In column "Region" number in parenthesis is the number of data used for average calculations for all data except R(560) and MRR. For these last two variables the number of data used for calculation is indicated at the right of the slash. a transitional zone (stations 16, 17, and 21) that were sampled in June 8 and 9. After these two days the frontal zone moved farther north and Chl SAT in the central zone (NG) of the study area gradually decreased until the end of the sampling period and in coincidence with the in situ sampling. It has to be noted that the west coast off Baja California peninsula presented higher Chl SAT with filaments that penetrate into the central part of the gulf. This zone was classified as NG despite the differences in Chl SAT from coast to coast and it was probably because stations that defined this region were sampled when this area had the lower Chl SAT and those high concentrations from the west coast strongly decreased (June 11 to 16). When a synoptic map is created joining the information from all stations sampled during a cruise, it has to be considered that during the cruise oceanographic conditions are changing. As previously noted, the study area is strongly affected by tidal mixing, especially around the sills of the midriff archipelago region (GI) and in the shallow UG and by a complex circulation pattern that is visible in Chl SAT .
The satellite data product accuracy generally obtained by the international missions is ±35% for Chla in open ocean waters [31], which is compatible with performance statistics reported in this work (Figures 10(a) and 10(b)), although with differences between sensors. However, a positive BIAS was observed for both sensors (SeaWiFS and MODIS), indicating an overestimation of in situ Chla. This is to be expected if we consider that the region is clearly influenced by high concentrations of nonpigmented material, especially to the north (UG and some stations of the NG) ( Figure 4 and Table 2). The relatively good performance of the algorithms can be explained by some compensation among variables responsible for ( ) variability; that is the backscattering coefficient ( ( )) can be higher in the green but it will be also high in the blue, but explaining how does this compensation occurs is difficult because we do not have measurements of ( ) or even of CDOM absorption coefficient. In addition, the satellite ( ) estimates in the blue are usually too low, resulting in too low satellite MRR ( Figure 10); therefore satellite Chla is too high. This may be due not only to absorption by nonalgal material, in particular sediments, and CDOM, but also to the influence of desert and continental aerosols (i.e., absorbing aerosols) over the Gulf region.

Conclusion
Our results indicate that the Northern Gulf of California can be classified into three bio-optical regions associated to previously identified surface circulation patterns and bottom topography. Besides, a transition zone was identified to the north, and the three southernmost stations (1, 2, and 3) are probably part of another region with different characteristics compared with the hither GI region. Pegau et al. [13] examined some inherent optical properties of these waters, focusing mainly on the south of our study area (Guaymas Basin) although sampling some stations also in the GI and NG regions. They considered the NG as belonging to Case 2 waters while those to the south belong to Case 1 (including the GI region). In our study, the differences between the NG and GI regions were considered significant when comparing the absorption spectra but they were not when comparing the reflectance spectra. Marine reflectance varies as a function of the ratio of the backscattering coefficient ( ) and the absorption coefficient ( ) [32], and we have measured only the last one. In addition, the absorption coefficient ( ) is a result of the combined contribution of water components like phytoplankton, nonpigmented particles, and colored dissolved organic matter ( CDOM ( )), of which we did not measure the last one, even though it was previously considered as important in the study area [13]. Thus our observations appear to confirm the differences in both ( ) and CDOM ( ) between the NG and GI explain the similarities between ( ) in these two regions. In addition, in situ absorption data were from surface samples while ( ) comes from the first optical depth, which ranged from 3 to 5 m at the GI region and it was higher than 5 m at the NG. Furthermore, the proportion of (440) to light absorption ( Figure 4 and Table 2) in the GI region was significantly lower than in the NG, although being very similar to the proportion at those stations to the south. It means that our observations probably confirm those of Pegau et al. [13]; that is, the NG region (and surely the UG) belongs to Case 2 while those to the south to Case 1. However, even when the UG and NG waters are considered as Case 2 they do have different biooptical characteristics. This implies that classifying this area as Case 1 or Case 2 is incomplete and instead both biological and optical properties should be considered to gain a better understanding of the ( ) variability in the Gulf region and the underlying biogeochemical processes. The bio-optical spatial variability of the study area has important implications for interpreting variations in ocean color. Although reasonable agreement were found between measured and satellite-estimated Chla, regional refinements of the empirical algorithm that take into account differences between bio-optical regions should improve the satellite estimates of Chla and also of the inherent optical properties for the study area. It should be noted that our findings are based on both (1) two weeks of observations and limited parameters and (2) one season; therefore it cannot provide a complete picture of bio-optical spatial and seasonal variability. Complementary field experiments and further research, including new cruises, are underway to understand