Introducción

El Ordenamiento Ecológico Marino (OEM) en los mares de México se puede definir de acuerdo a la Ley General del Equilibrio Ecológico y Protección al Ambiente (LGEEPA), como "un instrumento de la política ambiental, con el propósito de inducir el aprovechamiento de los recursos naturales del medio marino y de las actividades productivas, bajo un marco de sustentabilidad que integre la protección del medio ambiente marino y costero" ^{(Diario Oficial de la Federación, DOF 1996)}. La regionalización ecológica por su lado, parte de condiciones homogéneas desde el punto de vista físico-biótico en Unidades de Gestión Ambiental (UGA), y puede implementarse en múltiples niveles jerárquicos ^{(INE, 2006)}. Las Fases para la elaboración del programa del OEM son: Caracterización, Diagnóstico, Pronóstico y Propuesta. La primera fase es de vital importancia porque en ella se regionaliza y se delimita el área a ordenar y propone la descripción de los componentes físicos y biológicos.

Los primeros intentos por ordenar el mar surgen en la década de los ´80s del siglo pasado. ^{Merino (1987} incorporó componentes oceanográficos en un primer intento de ordenamiento marino. )La Secretaría de Marina (SEMAR) realizó zonificaciones de regiones y subprovincias dentro de la Zona Económica Exclusiva (ZEE) ^{(INEGI, 2009)} definidas por sus características oceanográficas, geológicas, biológicas, físicas y químicas. ^{Arriaga et al. (1998)} y ^{Botello et al. (2000)} reconocen siete y cinco provincias costeras y oceánicas, respectivamente, en los mares mexicanos. El Instituto Nacional de la Pesca ^{(INP, 1994)} elaboró su propia regionalización en función de la explotación de los recursos bióticos-pesqueros en los mares mexicanos en el Atlas Pesquero. El World Wildlife Fund, dividió al país en cinco regiones para la conservación de zonas costeras y marinas, mientras que la Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO) diseñó una regionalización que considera áreas marinas prioritarias de conservación (^{Arriaga et al. 1998}). Se realizó un estudio por ^{Lluch-Cota y Arias-Arechiga (2000)}, sobre regionalización, basado en la distribución anual de la concentración pigmentaria a partir de imágenes de satélites en el Golfo de California. Se formalizó una propuesta de regionalización marina por ^{Bezaury et al. (1996)} y ^{Gutierrez & Bezaury (2001)} como parte del "Proyecto Mapeo de Ecosistemas Marinos y Estuarinos de Norte América (B2B)" que también incluyó el trabajo de "Eco-regiones Marinas de Norte América" para la Comisión de Cooperación Ambiental (^{Wilkinson et al., 2003}).

En el plano internacional, el área del Pacífico Norte mexicano (PNM) se regionalizó con múltiples criterios de clasificación, como por ejemplo: "California Upwelling" propuesto por ^{Caddy & Bakun (1994)} basados en procesos dominantes de suplemento de nutrientes; "Mexican Temperate Pacific y Magdalena Transition" como regiones biogeográficas costeras propuesto por ^{Sealey y Bustamante (1999)}; "Provincia CC" (CALC) propuesta por ^{Longhurst (2007)} y "LME, California Current", dentro del esquema de los "Large Marine Ecosystems" de gran aceptación según ^{Sherman (1994)}, ⁽¹⁹⁹⁶⁾ y ^{Sherman & Tang (1999)}.

Entre el 2004 y el 2007 el entonces Instituto Nacional de Ecología (INE), ahora Instituto Nacional de Ecología y Cambio Climático (INECC) y la Secretaría de Medio Ambiente y Recursos Naturales (SEMARNAT), realizaron tres Talleres de Expertos en materia de OEM durante el 2004, 2006 y 2007 (^{GEA, 2007}; ^{Córdova et al., 2009}), para discutir e intentar consensuar elementos técnicos y metodológicos para la regionalización marina de los mares mexicanos. Durante el I Taller de OEM en 2004, la Mesa Oceánica aprobó presentar una regionalización espacial dinámica, con temperaturas superficiales correlacionadas con la clorofila y la productividad primaria ^{(Córdova et al., 2006)}. En los procesos de caracterización oceánica, el uso de la producción del carbono orgánico nuevo se explica sobre la base de su potencial de aplicaciones: el suplemento de nutrientes, la función básica de la física de los ecosistemas y cómo éstos se relacionan con su ciclo biogeoquímico y la física del océano (^{Hernández-de la Torre, et al., 2011}).

Justificación. ¿Por qué utilizar la Pnueva para establecer una regionalización marina? La Pnueva es la fracción de la producción primaria total que se genera a partir de una fuente de nitratos alóctonos (^{Dugdale & Goering, 1967}) presentes por debajo de la termoclina estacional y que son transportados hacia las capas superficiales del océano. La Pnueva es considerada como la porción de la producción fitoplanctónica que es exportada por debajo de la zona eufótica (ZE) como materia orgánica particulada sin que el ecosistema pelágico se colapse (^{Eppley, 1989}). El transporte vertical de nitratos hacia la ZE es el factor principal en la regulación de la biomasa y producción del fitoplancton en las aguas costeras de la Corriente de California (CC), donde la distribución vertical y la profundidad del inicio de la nutriclina están correlacionadas con la razón de producción (Eppley & Peterson, 1979). Las pesquerías están controladas por la producción de carbono orgánico nuevo a través de su consumo por organismos de niveles tróficos mayores en la ZE y de su transporte vertical hacia el fondo marino (^{Chen-Tung, 2003}). Este transporte de material orgánico desde aguas superficiales sustenta a los organismos meso y batipelágicos e influye en el intercambio de gases con la atmósfera, tales como el dióxido de carbono y oxígeno. La Pnueva no se puede medir de forma directa a través de la percepción remota por no existir una expresión electromagnética que pueda ser correlacionada directamente con ella. Diferentes grupos de investigadores han calculado la Pnueva por imágenes de sensores remotos a través de la relación temperatura-nitratos derivada de mediciones in situ y temperatura de satélite (^{Goes et al., 2000}, ²⁰⁰¹, ²⁰⁰⁴); por imágenes de la temperatura superficial del mar (TSM) del AVHRR o del MODIS con modelos que determinan el máximo profundo de clorofila-a (Chl-a), (e. g. ^{Sathyendranat et al. 1991}); o por la relación de Chl-a con la Pnueva (^{Kudela & Chavez, 2002}). El cálculo de la Pnueva por imágenes de sensores remotos a bordo de satélites ha ganado importancia y aceptación, ya que es una herramienta útil para determinar grandes áreas de producción fitoplanctónica (^{Kudela & Dugdale, 1996}; ^{Dugdale et al., 1997}; ^{Watts et al., 1999}; ^{Carr, 2001}; ^{Henson et al., 2003}; y ^{Coles et al., 2004}). El estudio de la variabilidad temporal de la producción primaria nueva (Pnueva), es uno de los aportes más importantes de la oceanografía biológica puesto que permite semicuantificar la producción exportada y por tanto, estimar el potencial de los océanos en remover el carbono atmosférico y capturarlo a las grandes profundidades.

Para implementar la metodología que se propone dentro del OEM en el PNM, es indispensable contar con la delimitación y regionalización de cada área específica en el océano, -a escalas adecuadas- a través de una caracterización marina que describa los componentes físicos y biogeoquímicos de mayor importancia. Igualmente, identificar y priorizar aquellos atributos ambientales que permitan efectuar acciones para el desarrollo de actividades productivas y la conservación de los ecosistemas mediante las UGAs. Una UGA (Unidad de Gestión Ambiental) es la unidad mínima territorial donde se aplican lineamientos y estrategias ambientales, con esquemas de manejo de recursos naturales y orientados a un desarrollo que transite a la sustentabilidad (^{Rosete et al., 2006}).

La introducción de novedosas técnicas aplicadas al monitoreo e integración de datos ambientales en el océano, pudiera ser una alternativa para llegar a entender la respuesta del ecosistema marino ante perturbaciones de corto y mediado período en el establecimiento de un orden marino. Por consiguiente el objetivo del presente trabajo es proponer una metodología que estime temperatura y tasas de Pnueva en profundidad partir de imágenes satelitales de temperatura superficial y relaciones empíricas temperatura-nitratos, con vistas a incorporar los principales procesos físicos y biogeoquímicos en la costa y el océano frente a la Península de Baja California. Ello contribuirá a la elaboración de la propuesta de regionalización de la costa y océano en el PNM -ya que en la práctica no existe una regionalización basada en lo estrictamente marino- sino que están basados en métodos propuestos de Ordenamientos Ecológicos Terrestres, que contienen y circunscriben numerosos aspectos sociales y económicos.

Material y métodos

Se utilizó la red de estaciones del programa CalCOFI, http://www.calcofi.org/ (Figura 1a) , para la obtención de los datos de temperatura (ºC) y nitratos (NO₃) de enero de 1952 a diciembre de 2002 y los datos de la red IMECOCAL (Investigaciones Mexicanas de la Corriente de California) http://imecocal.cicese.mx de 1997 hasta 2008. La red va desde la Línea 90.28 (33º55´N y 117º63´O) hasta la 90.120 (30º41´N y 124º00´O) y por el sur, desde la Línea 157.15 (22º53´N y 109º43´O), hasta la 157.90 (20º19´N y 113º74´O), (Figura 1), con más de 400 km de la costa hacia el océano y de forma perpendicular al eje de la península. Se consideró agrupar los datos de la temperatura superficial del mar en clases de cada 0.5 oC, mientras que los perfiles de temperatura se agruparon en profundidades estándar (0, 10, 30, 50, 75 y 100 m). Se eliminó todo perfil menor a 200 metros y los mayores se recortaron a esta última profundidad. La información oceanográfica se agrupó y se clasificó de acuerdo al Índice Multivariado del ENSO (El Niño/Southern Oscillation) (MEI, ^{Wolter & Timlin 1993}, ¹⁹⁹⁸) en años "El Niño", "La Niña" y "Año Normal". Se sacó la media, desviación estándar, máximo y mínimo para cada profundidad por rango de temperatura, a los cuales se les denominó perfiles verticales climatológicos (PVC) mismos que se adjudicaron a las imágenes compuestas mensuales de temperatura del sensor AVHRR del 2003, para crear nuevos mapas en profundidad. Se tomó como base el diapasón térmico de la serie 1952-2002 para la región y se decidió catalogar las temperaturas entre 10º y 14.9 ºC como aguas frías; entre 15º y 19.9 ºC medianamente frías; entre 20º y 24.9 ºC medianamente cálidas y ≥ 25 ºC aguas cálidas. Se ajustaron modelos estadísticos lineales de la relación temperatura-nitratos separados en regiones costeras y oceánicas según ^{Hernández-de la Torre et al. (2003)} . Estos últimos fueron trasferidos a los PVC de temperatura y se interpolaron finalmente los parámetros de la regresión en forma latitudinal y longitudinal, y se confeccionaron mapas de nitratos para la zona de estudio.

Figuras 1a-b: a) Península de Baja California y área del Pacífico Norte mexicano. b) Líneas y estaciones de la red CalCOFI e IMECOCAL. Abarca más de 400 km fuera de la costa y de forma longitudinal (Tomado de CalCOFI Basic Station Plan inside back cover). La líneas punteada representa el límite de la Zona Económica Exclusiva (ZEE). 

La razón-ƒ_z es un cociente entre la Pnueva y la producción primaria total (Ptotal), (razón-ƒ_z = P_nueva/P_total) (1) y varía entre 0 y 1 (^{Dugdale & Goering, 1967}). Con base en los datos reportados por ^{Eppley y Peterson (1979)} para el Southern California Bight (33ºN; 119ºO) y el modelo empírico propuesto por ^{Harrison et al. (1987}) para la misma zona, se calculó la razón-ƒ_z mediante la ecuación:

donde:

ƒ_max es el valor máximo esperado y

m es la pendiente inicial de la curva (NO₃ vs ƒ).

Debido a que los valores altos de NO₃ se presentan con mayor frecuencia en las estaciones costeras, se tomó ƒ_max = 0.67 y m = 13.7 (Tabla II del modelo empírico de ^{Harrison et al., 1987}), mientras que con base en el conocimiento de que en las estaciones oceánicas se presenta una disminución progresiva de los NO₃ a medida que se aleja de la costa, se tomó ƒ_max = 0.61 y m = 10.5 (Tabla II del modelo empírico de ^{Harrison et al., 1987}).

La producción primaria total integrada (PPi, mgC m^-2 d^-1) se calculó mediante imágenes compuestas mensuales para el 2003 de TSM del sensor AVHRR de 4 x 4 km en formato HDF del sensor Sea-viewing Wide Field-of-View Sensor (SeaWiFS) (^{Behrenfeld & Falkowski, 1997}), http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/

Para el procesamiento de las imágenes de clorofila, temperatura y radiación fotosintéticamente activa (luz PAR) del sensor SeaWiFS, se utilizó el programa Windows Image Manager (WIM-WAM, versión 6.56, Junio 2009; ^{Kahru, 2007}). Para pasar la PPi (mgC m^-2 d^-1) a Producción Primaria Total en profundidad (PPz, mgC m^-3 d^-1), se utilizaron los perfiles promedio de producción primaria de mediciones in situ durante los cruceros IMECOCAL en 2003. Se realizó una integración con base en la irradiancia del área bajo la curva para obtener la proporción o el porcentaje de área de PPz correspondiente a cada profundidad, de acuerdo al método descrito por ^{Hernández-de la Torre et al. (2003)}. Finalmente, la matriz de datos de la razón-ƒ_z se multiplicó por la matriz de datos de la imagen de PPz (Pnueva = PPz* razón-ƒ_z) y se obtuvo una imagen final que estima tasas de Pnueva a 0, 10, 30, 50, 75 y 100 m de profundidad. Se tomó como base el diapasón trófico de la serie 1952-2002 para la región de estudio, por lo cual se decidió catalogar la Pnueva <5.99 (mgC m^-3 d^-1) aguas oligotróficas; entre 6 y 14.99 (mgC m^-3 d^-1) aguas mesotróficas y >15 (mgC m^-3 d^-1), aguas eutróficas.

Mapas en profundidad. Para confeccionar los mapas de temperatura en profundidad se realizó un programa elaborado en MatLab el cual lee el archivo de datos contenido en la imagen de TSM del sensor AVHRR con resolución espacial de 4 x 4 km en formato HDF. Después de leer la información por pixeles de la imagen seleccionada, se adjudicaron los PVC correspondientes al año que corresponda ("El Niño", "La Niña" o "Año Normal") y se confeccionó una nueva imagen para cada profundidad estándar de temperatura y de nitratos. A las matrices de datos de la imagen de nitratos, se les aplicó los coeficientes del modelo de ^{Harrison et al. (1987)} para estimar la razón-ƒ_z por profundidad. Por último, la matriz de datos de la razón-ƒ_z se multiplicó por la matriz de datos de la imagen de producción primaria total, previamente elaborada por el programa WIM-WAM (obtenida de un compuesto mensual de TSM, clorofila a y luz PAR del sensor SeaWiFS), y se logró una imagen final que estima los valores de Pnueva a 0, 10, 30, 50, 75 y 100 m de profundidad. Los promedios climatológicos mensuales (1970-2008) de las surgencias costeras frente a la Península de Baja California, provienen del sitio Web: http://www.pfeg.noaa.gov/products/PFEL/modeled/indices/upwelling/upwelling.html

Para caracterizar en el plano horizontal el área de estudio, se utilizó la propuesta de "esquemas de franjas" de Grandes Ecosistemas Marinos (GEMs 3, Corriente de California) de ^{Escofet (2004)}, mientras que la primera propuesta de caracterización en profundidad, es la que se deriva del presente trabajo.

La escala de los mapas de 0-100 m es de 1:2 500 000. Con esta información se regionalizó el PNM con criterio de las UGAs. Se realizó un análisis de conglomerados basado en las escalas de enumeración continuas entre Pnueva y la temperatura. De tal modo, el espacio de la regionalización de las UGAs para la región del PNM queda delimitado hacia el océano por la línea exterior de la ZEE y hacia el continente por la línea de costa o franja intermareal.

Resultados

Caracterización marina con base en la temperatura. La Figura 2a muestra la caracterización física en profundidad (0-100 m) de la temperatura en la primavera de 2003. Se destaca la presencia de una zona con temperatura medianamente frías a frías (16º y 13 ºC) desde superficie hasta los 75 m en la zona costera, desde la frontera con el sur de California, Estados Unidos, hasta el norte del Golfo de Ulloa, que funge como línea oceanográfica divisoria entre el norte y el sur de la Península de Baja California. A 75 y 100 m se observa un área de surgencias costeras desde Punta Santo Tomás hasta Punta Baja y tres regiones oceánicas perpendiculares a la costa principalmente alrededor de los 75 m de profundidad. Al suroeste de Los Cabos, en los primeros 30 m, la zona es considerada como medianamente cálida (21º-23 ºC). La Figura 2b muestra la caracterización para el verano: de 0 a 50 m la región de la Isla Navidad se destaca como la línea oceanográfica divisoria entre el agua medianamente fría al norte y medianamente cálidas del sur. Un área de aguas medianamente cálidas proveniente de la región ecuatorial llega hasta la zona costera de Punta Abreojos hasta el Golfo de Ulloa desde superficie hasta 50 m. La temperatura a 75 m fue medianamente fría para toda la región y a 100 m de profundidad presentó valores entre 13º y 14 ºC (fría) al norte hasta Bahía Magdalena, y al sur de la misma, se comporta medianamente fría (15º-17 ºC). La Figura 2c señala la caracterización para otoño, con temperaturas medianamente frías (entre 19º y 16 ºC) desde superficie hasta 75 m al norte de Isla Navidad. Aguas cálidas (entre 25º y 30 ºC) se presentan al sur de Isla Navidad por toda la ZEE hasta 30 m de profundidad. A partir de 50 m desde Punta Abreojos hacia el sur, las temperaturas son medianamente cálidas. A 75 m las temperaturas son medianamente frías en toda la región, mientras que a 100 m oscila entre frías al norte y medianamente frías al sur de Isla Navidad. La Figura 2d muestra la caracterización física en profundidad con base en la temperatura en invierno del 2003. Al norte de 25º N desde la superficie hasta 50 m, la temperatura del océano presentó un comportamiento medianamente frío (15º-19 ºC), excepto al sur de esta misma latitud que se muestra con aguas medianamente cálidas (21º-24 ºC) desde Bahía Magdalena hacia el suroeste. Entre 75 y 100 m se observó penetración de aguas frías (14º-16 ºC) por el norte que llegan hasta la zona de Bahía Magdalena, mientras que al sur de la misma se reportó aguas medianamente frías (15º-17 ºC).

Figuras 2a-d: Caracterización física en profundidad de la temperatura del mar en la ZEE del Pacífico Norte mexicano por época del año para 2003. Regiones Frías (F) aquellas que presentan una dominancia de temperatura entre 10º y 14.9 ºC; Regiones medianamente frías (MF) aquellas que presentan una dominancia de temperatura entre 15º y 19.9 ºC; Regiones medianamente cálidas (MC) las que presentan una dominancia de temperatura entre 20º y 24.9 ºC; y Regiones Cálidas (C) son las que presentan una dominancia de temperatura entre 25º y 30.0 ºC. 

Caracterización marina con base en la Pnueva. La Figura 3a muestra la caracterización biogeoquímica en profundidad con base en la Pnueva en la primavera del 2003. Se destacó una región mesotrófica costera en los primeros 10 m desde Tijuana hasta Isla Navidad; y al sur, desde el Trópico de Cáncer hacia la ZEE. En la parte oceánica se presentó una región oligotrófica (0-2, mgC m^-3 d^-1) al norte de 25º N, posiblemente relacionada con el aporte de aguas pobres en nutrientes procedente del Pacífico Central hasta los 100 m. Entre 30 y 50 m se presentó región mesotrófica en el área cerca de Punta Santo Tomás y llegó hasta el sur de la ZEE del PNM. Esta región presentó giros anticiclónicos al oeste de Isla Navidad con una región eutrófica a 50 m de profundidad que circundó el Golfo de Ulloa. A 75 m se observó una franja al oeste del Golfo de Ulloa hasta la ZEE en forma de giros con carácter mesotrófico (6-14, mgC m^-3 d^-1). Para el verano, el panorama biogeoquímico cambió drásticamente a condiciones oligotróficas (0-5, mgC m^-3 d^-1) en los primeros 30 m de toda el área de estudio (Fig. 3b). Al sur de la Isla Navidad la Pnueva se presentó con condiciones mesotróficas en 50 m, entre (6-9, mgC m^-3 d^-1), mientras que al norte permanecieron oligotróficas. La Figura 3c muestra escenarios oligotróficos en toda el área de estudio durante el otoño. En invierno (Fig. 3d), se presentó región oligotrófica de 0-10 m en toda el área, y permaneció en el norte de la Península de Baja California hasta los 50 m. Un área mesotrófica se presentó al sur de los 25º N entre 30-50 m, mientras que entre 75 y 100 m se presentó condiciones oligotróficas.

Figuras 3a-d: Caracterización biogeoquímica con base en la Pnueva del mar para el Pacífico Norte mexicano por época del año. Regiones Oligotróficas (O) son aquellas que presentan una dominancia de Pnueva </= 0-5 mgC m^-3 d^-1; Regiones Mesotróficas (M) aquellas que presentan una dominancia de Pnueva entre 5 y 15 mgC m^-3 d^-1; Regiones Eutróficas (E) aquellas que presentan una dominancia de Pnueva >15 mgC m^-3 d^-1. 

Regionalización marina. En la Tabla 1 se muestra la propuesta de regionalización oceánica y costera para el PNM en nueve UGAs costeras y oceánicas.

Tabla 1: Propuesta de regionalización de UGAs costeras (c) y oceánicas (o) para el Ordenamiento Ecológico Marino del Pacífico Norte mexicano. 

El análisis de agrupamiento de las UGAs originó resultados especialmente coherentes con la ecología del área de estudio. En la Figura 5a se muestra el clúster para la temperatura con índice de afinidad para las UGAS I, II, III, VII y IV, (todas referidas a la zona norte de la Península de Baja California), mientras que para la zona sur las UGAs V, IX, VI y VIII mostró igualmente afinidad. Las UGAs IV y V evidenció una región intermedia costera para la temperatura. En la Figura 5b se muestra el vínculo de la Pnueva para las UGAs I, II, III, VII y V referidas al área norteña de la Península de Baja California, mientras que para la zona sur las UGAs VIII, VI y IX representó igualmente acoplamiento. Al igual que en la temperatura, las UGAS IV y V evidenciaron la zona intermedia costera para la Pnueva con la UGA VIII oceánica. Toda la franja, tanto costera como oceánica funge como "zona de transición" o "parteaguas marino" en las siguientes condiciones anuales medias: a) ante un fortalecimiento de incremento de la Corriente de California acontecidas por eventos extremos en el Pacífico nororiental (frío o cálido), las condiciones físicas y biogeoquímicas imperantes en esa "zona de transición" se mueve hacia el sur. Por el contrario, ante perturbaciones ecuatoriales como las condiciones de "El Niño o La Niña", esa "área de transición" se traslada o se mueve más hacia el norte.

Figura 4: Propuesta de Regionalización Oceánica y Costera para el Pacífico Norte mexicano con base en la caracterización física y biogeoquímica por profundidad. 

Figuras 5a-b: Análisis de agrupamiento de las UGAS para el Pacífico Norte mexicano. La figura (a) muestra el clúster para la temperatura con índice de afinidad para las UGAS I, II, III, VII y IV, referidas a la zona norte de la Península de Baja California. Para la zona sur las UGAS V, IX, VI y VIII muestran igualmente afinidad. Las UGAS IV y V evidencia la zona intermedia costera. La figura (b) señala el vínculo de la Pnueva para las UGAS I, II, III, VII y V referidas al área norteña de la Península de Baja California, mientras que para la zona sur las UGAS VIII, VI y IX muestran igualmente acoplamiento. Las UGAS IV y V indican la zona intermedia costera con la UGA VIII oceánica con respecto de la Pnueva. 

La Tabla 2 muestra un resumen de la ubicación geográfica de cada UGA. Señala la cercanía de las Líneas que conforman la red CalCOFI e IMECOCAL dentro del área de estudio; las condiciones promedios de la temperatura y Pnueva para 30 y 50 m en primavera, verano, otoño e invierno; las pendientes de la relación temperatura-nitrato obtenidas en cada UGA, así como la climatología (1970-2008) de las surgencias costeras en el área.

Tabla 2: Resumen de la ubicación geográfica de cada UGA. Se señala una aproximación de las líneas que conforma la red CalCOFI e IMECOCAL dentro del área de trabajo, las condiciones medias en primavera, verano, otoño e invierno de la temperatura y Pnueva para 30 y 50 m; las pendientes de la relación temperatura-nitrato obtenidas en cada UGA, así como la climatología (1970-2008) de las sugerencias costeras en el área de estudio. (C=costeras O=oceánicas). Temperatura 10º -15ºC Región fría (F); 15º-20ºC Región medianamente fría (MF); 20º-25º Región medianamente cálida (MC); >25ºC Región cálida (C). Pnueva: Región oligotrófica (0): </=(0-5, mgC m^-3d^-1); Región mesotrófica (M): (5-15mgC m^-3 d^-1); Región eutrófica (E): (15-25, mgC m^-3 d^-1). 

Discusión

^{Pielou (1979)} , ^{Hayden et al. (1984)} y ^{Zacharías & Roff (2000)}, coinciden en que la regionalización basada en el entorno físico es generalmente una buena forma de predecir respuestas de las regiones bióticas. ^{Vallega (2002)} afirmó que con la creación de regiones oceánicas no se pretende ser estático, sino dinámico. ^{Bocco et al. (2001)} por su parte, señalaron que la vegetación y los suelos tienden a variar en forma predecible dentro de una unidad de relieve y señalan que éstas son clasificadoras integradas del paisaje. Si se toma el color del mar (pigmentos derivados de la abundancia del fitoplancton marino) como una analogía de la vegetación terrestre y la temperatura del mar como el suelo (sustrato físico característico) donde crece la vegetación, e interpretamos el "paisaje" como "paisaje marino" dentro del concepto de ordenamiento, entonces podemos justificar la cimentación de la propuesta de regionalización marina basada en la temperatura y la Pnueva por las dos razones siguientes: primero, el color del océano representa a la biomasa del fitoplancton derivada de la concentración de los pigmentos fotosintéticos y en última instancia, es la que caracteriza las diferencias de las imágenes obtenidas por sensores remotos. Segundo, la temperatura de la superficie del mar se reconoce como la 'huella digital térmica' que se imprime en diversos procesos mecánicos y termodinámicos que ocurren en la interface océano-atmósfera (^{Gallegos-García et al., 2003}), los que a su vez inducen en el mar efectos y fenómenos cuya expresión en términos de ubicación geográfica, forma, regularidad, duración, intensidad, tamaño y variabilidad determinan una regionalización natural, sustantiva para un OEM ^{(Gallegos-García et al., 2009)}. Es por ello que el uso potencial del carbono marino, es un instrumento operativo para hallar respuestas de regionalización en la administración de los ecosistemas marinos mexicanos.

Las diferencias espaciales y temporales de norte a sur en las anomalías de temperatura, sugieren que la región IMECOCAL está separada en dos provincias dinámicamente distintas con respuestas disímiles a los forzamientos físicos (^{Durazo, 2009}). Los vientos a lo largo de la costa de la Península de Baja California son más intensos durante la primavera (^{Pérez-Brunius et al., 2007}, ^{Castro & Martínez, 2010}) y la advección norte-sur asociada a estos vientos, produce dominancia de aguas relativamente bajas en temperatura en la región (^{Durazo et al., 2010}). El viento paralelo a la costa origina surgencias costeras y éstas son una característica oceanográfica principal de la zona costera en la región occidental de la Península de Baja California. Los afloramientos de aguas costeras se identifican por presentar en sus inicios temperaturas frías, altos nutrientes y bajas concentraciones de clorofila cerca de la costa e inclinan la nutriclina hacia la superficie en zonas cercanas a la costa (^{Gaxiola-Castro et al., 2010}). Al norte de la península las surgencias se presentan durante todo el año, pero con mayor intensidad durante la primavera; mientras que en la zona sur ocurren principalmente durante primavera y verano (^{Durazo et al., 2010}). Debido al efecto combinado de las mareas y el forzamiento por el viento, las aguas originadas durante las surgencias con altas concentraciones de nutrientes inorgánicos disueltos y material orgánico, es advectada hacia el interior de las bahías y lagunas bajacalifornianas (^{Camacho, 2009}). Los valores más altos de producción primaria nueva se obtuvieron a 50 m en la región sur frente a Punta Abreojos dentro del Golfo de Ulloa (Fig. 3a). Los aportes de nutrientes y la mayor cantidad de luz en primavera, determinan un perfil de condición mesotrófica en los primeros metros de profundidad para las áreas del norte de Isla Navidad, con características eutróficas entre 30 y 50 m para el resto de las áreas costeras del centro y en el sur de las áreas oceánicas (Figs. 3a , ⁴).

Una de las principales características del océano frente a la Península de Baja California es el intercambio de agua con altas concentraciones de nutrientes que se desprenden de Punta Concepción, California y los aportes de agua oligotrófica proveniente del Océano Pacífico central. Este desprendimiento tiende a formar chorros y separarse de la costa por la presencia de una cordillera submarina que no permite totalmente su ingreso al Southern California Bight, por lo que en su movimiento hacia el sur, arrastra agua oligotrófica del Pacífico central que afectan la región oceánica norte de la Península de Baja California principalmente entre los 30 y 75 m (Figs. 3c-3d). En la región oceánica central de la península se entremezclan la influencia de dos tipos de aguas: tropical-subtropical y aguas de origen subártica (^{Durazo et al., 2010}). Este patrón promedio anual se manifiesta como un "parteaguas" frente a la región de Isla Navidad y se caracteriza por la mayor influencia de agua subártica al norte, mientras que al sur, combina éstas en temporadas frías y aguas subtropicales el resto del año manifestándose en los diagramas T-S (^{Durazo et al., 2010}). Estos mismos autores confirmaron la presencia de la contracorriente en profundidades intermedias dirigida hacia el polo, localizada en la región sur de la Península de Baja California. Esta contracorriente transporta agua de origen subtropical, con un patrón de circulación estacional formado por dos giros ciclónicos (a favor de las manecillas del reloj) separados a la altura de Isla Navidad. Algunos indicadores biológicos como la anchoveta norteña, la sardina y la merluza del Pacifico, muestran una separación poblacional entre las dos regiones lo que sugiere la existencia de dos provincias (^{Hewitt, 1981}). Adicionalmente se ha reportado una diferenciación en la asociación faunal y la variabilidad oceánica del zooplancton frente a Baja California (^{Lavaniegos et al., 2002}, ^{Lavaniegos, 2009}).

En la oceanografía biológica es muy importante tener un acercamiento cuantitativo de la fertilidad debajo de la superficie marina, que auxilie a diferentes usuarios del Sector Pesquero, Turismo y aquellas secretarías del Gobierno Federal y Estatal involucrados en el OEM. En consecuencia, debido a que la producción de fitoplancton marino constituye la base de los recursos pesqueros, se debe regular para aprovecharlos de manera sustentable, en congruencia con los intereses de la sociedad. La Pnueva es una "bolsa general" de alimento nuevo que se desarrolla a medida que crece el fitoplancton, donde participa el zooplancton como consumidor de esta biomasa, para formar un equilibrio entre las poblaciones vegetales y animales que componen la base de la trama trófica marina. Las regiones mesotróficas y eutróficas en primavera frente a Baja California, propiciaron a una mayor presencia de organismos desovadores, que debido al "frenesí alimenticio", construyeron un nicho ecológico cuya duración persistió en dependencia de las abundancias de los productores secundarios. La biomasa del zooplancton en la región IMECOCAL es mayor en verano y otoño, compuesta principalmente por copépodos eufáusidos y otros grupos con menores abundancias (^{Lavaniegos et al., 2002}). La región oceánica central de la Península de Baja California se caracteriza por presentar altos valores de Pnueva, con regiones meso-eutróficas hasta cerca de 25º N (Fig. 3), principalmente en primavera y verano. Las áreas de Pnueva en primavera a 30m de profundidad son consistentes con la distribución de biomasa fitoplanctónica, concentraciones de macrozooplanton y densidad promedio de larvas de peces neríticos, pelágicos y bentónicos según ^{Lluch-Belda (2000)} para la zona adyacente a Punta Baja, Bahía de Sebastián Vizcaíno y el Golfo de Ulloa. Estas áreas también fungen como zonas marinas de "resiliencia natural" convirtiéndose en sumideros de carbono extremadamente eficiente. La Pnueva subsuperficial puede ser utilizada para determinar el hábitat favorable de desovadores y entender la conexión en profundidad de recursos de importancia comercial y ecológica como sardina, anchoveta y otros pelágicos menores.

Por primera vez se expone una metodología de la regionalización marino-costera a partir de los conceptos de regionalización territorial, definidos por modelos con las variables de mayor estabilidad de la Ecología Marina: la Oceanografía Física (temperatura) y la Biogeoquímica Marina (Pnueva). Se demostró que los procesos físicos entre la atmósfera y el océano representados por la variabilidad espacio-temporal de la temperatura, promueven el transporte vertical de nitratos hacia la zona eufótica, los que determinaron los cambios en la producción de carbono nuevo en el área del sur de la Corriente de California frente a Baja California. Con la metodología implementada en este trabajo, se propone separar el ecosistema epipelágico frente a la Península de Baja California en seis UGAs costeras y tres oceánicas. De las UGAs oceánicas la VIII y las vecinas UGAS costeras IV y V son consideradas "zonas de transición" entre las condiciones oceanográficas norteñas y sureñas del océano frente a Baja California. La presente propuesta metodológica da un paso más en la consecución de mejoras a nivel técnico del OEM, -que es el instrumento de política ambiental en México que recomienda las mejoras en las actividades marinas- y un impulso académico para coadyuvar a cambiar el esquema de que el océano cercano a la costa es muy dinámico y no se puede ordenar desde el punto de vista ecológico.