Introducción

El mapeo del fondo marino genera información indispensable para la conservación y el manejo de los recursos marinos. Es fundamental para la zonificación y planeación de ordenamientos ecológicos marinos y áreas naturales protegidas. Estos instrumentos de conservación cuen tan con amplias referencias metodológicas en el medio terrestre, pero son muy escasas en el medio marino (^{Colenutt et al., 2013}; ^{Espinoza-Tenorio et al., 2014}).

El Golfo de California es reconocido por su elevada biodiversidad y productividad (^{Brusca, 2010}). Particularmente, el sureste de la bahía de La Paz (BLP) (Figura 1a) ha sido denominado como una Región Marina Prioritaria e incorpora un Área de Importancia para la Conservación de Aves (^{Arriaga-Cabrera et al., 2009}), dos sitios Ramsar, un parque nacional y 23 zonas con manglar. También incluye el Área de Protección de Flora y Fauna “Islas del Golfo” y el área de Balandra, ambas consideradas por la UNESCO como Patrimonio Natural de la Humanidad.

Figuras 1a-b a). La Bahía de La Paz y localidades. El área de estudio se muestra en el recuadro rojo, y las cuencas hidrológicas, con líneas azules. Batimetría de ^{Nava-Sánchez et al. (2001}). b) Área de estudio. Sitios de muestreo del fondo marino (puntos) y muestra independiente (triángulos). Las líneas azules representan el drenaje de las pequeñas cuencas hidrológicas. 

La BLP registra un ambiente bentónico con sedimentos mixtos (siliciclásticos-carbonatados) y carbonatados (^{Halfar et al., 2004}). Los ambientes carbonatados (>90% de CaCO₃) son poco comunes en las costas orientales del océano Pacífico (^{Halfar et al., 2004}). Estos am bientes se asocian a organismos como corales y rodolitos, que pueden formar colonias en el fondo y concentrar una biodiversidad relativa mente elevada (^{Reyes-Bonilla, 1992}; ^{Steller et al., 2003}). A pesar de las declaratorias y la presencia de estos ambientes sedimentarios, no se han propuesto mapas que clasifiquen y caractericen unidades en el fondo marino de la BLP.

Para realizar un mapeo en el fondo marino se han propuesto me todologías con base en componentes bióticos y abióticos, como la cla sificación de biotopos (^{Sergej & Ducrotoy, 2006}) y hábitats bentónicos (^{Sheppard et al., 1995}; ^{Cerdeira-Estrada et al., 2008}). Estas clasifica ciones generalmente se desarrollan con ayuda de sensores remotos sobre costas carbonatadas someras con estructuras de arrecifes cora linos (^{Giardino et al., 2016}). En estas zonas es posible capturar infor mación escénica; sin embargo, en las plataformas continentales como la BLP o cualquier zona alejada de la costa (más allá de los 20-30 m), el uso de sensores ópticos pasivos está limitado por la profundidad y la transparencia del mar (^{Hamylton, 2011}). Para zonas más profundas se han hecho clasificaciones basadas en componentes abióticos como la determinación de geoformas y estructuras del fondo (^{Walker et al., 2008}; ^{FGDC, 2012}; ^{Haris et al., 2012}; ^{Micallef et al., 2012}). Sin embar go, cuando el fondo marino carece de geoformas o estructuras bien definidas, y está compuesto principalmente por sedimentos no conso lidados, la identificación y delimitación de unidades no es sencilla. En fondos con sedimentos no consolidados se han utilizado clasificaciones basadas en un análisis textural en el que sólo se incluye el porcentaje de arena, limo y arcilla (^{FGDC, 2012}; ^{Jerosch, 2013}). No obstante, otras variables, como el contenido de CaCO₃, el tamaño medio de grano (^{Ver faillie et al., 2006}), la clasificación del sedimento (^{Ma et al., 2010}) y la profundidad de la columna de agua (^{FGDC, 2012}; ^{Jerosch, 2013}; ^{Harris et al., 2014}; ^{Snickars et al., 2014}), pueden aportar información para delimitar zonas o unidades en el fondo marino.

El contenido porcentual de CaCO₃ y la composición de bioclas tos indican el origen y la asociación de organismos marinos (^{Morse & Mackenzie, 1990}), así como la relación entre ambientes sedimentarios mixtos. El tamaño medio de grano se asocia a los niveles de energía de las condiciones hidrodinámicas del ambiente marino (^{Middleton, 1976}; ^{Syvitski, 1991}). La clasificación del sedimento ayuda a interpretar pro cesos de ambientes sedimentarios (^{Boggs, 2012}), y la profundidad de la columna de agua es una variable frecuentemente utilizada en estu dios ecológicos marinos (^{Sheppard et al.,1995}; ^{McLachlan & Brown, 2006}; ^{Snickars et al., 2014}). En este sentido, este trabajo propone un método multivariante para clasificar unidades en el fondo marino de zonas costeras y plataformas continentales donde se incorpora la va riabilidad espacial del contenido de CaCO₃, el tamaño medio de grano, la clasificación del sedimento y la profundidad de la columna de agua. Particularmente, en este estudio: 1) se modela la distribución espacial de las cuatro variables citadas con base en un algoritmo de interpola ción y 2) se clasifica y caracteriza las unidades del fondo marino a partir de la distribución espacial de dichas variables utilizando un análisis de agrupamiento jerárquico.

Marco geológico y sedimentos marinos. El fondo de la BLP (Figura 1a) está conformado por plataformas estrechas con depresiones y bloques elevados (^{Nava-Sánchez et al., 2001}). La porción noroeste de la BLP es más profunda e incluye a la cuenca Alfonso, con 420 m de profundidad (^{Cruz-Orozco et al., 1996}). La depresión recibe sedimentos terrígenos, eólicos y hemipelágicos (^{Silverberg et al., 2014}). Hacia el sureste de la BLP la profundidad disminuye desde los 220 m hasta la línea de cos ta, configurando una plataforma con pendiente suave donde el fondo marino contiene arenas finas de origen terrígeno (^{Cruz-Orozco et al., 1996}). Sobre la línea de costa, el transporte litoral se lleva a cabo prin cipalmente en dirección sur-sureste hasta configurar una barra arenosa conocida como El Mogote (^{Nava-Sánchez & Cruz-Orozco, 1989}). Esta barra disminuye la energía del mar enmarcando a la laguna de La Paz, la cual está conectada con la BLP por medio de un canal de mareas donde se registra la profundidad máxima de 10 m. La laguna recibe sedimentos por aporte fluvial (arroyos) y por transporte litoral y eóli co (^{Nava-Sánchez & Cruz-Orozco, 1989}); ^{Choumiline et al., 2009}). Los sedimentos de la laguna de La Paz se clasifican como arenas gruesas hasta limos muy finos (^{Álvarez-Arellano et al., 1997}). En el canal de San Lorenzo (Figura 1b) los sedimentos tienen un tamaño de arenas gruesas a muy gruesas, con una composición mineralógica mayor que 90% de CaCO₃. Los bioclastos corresponden a fragmentos de moluscos, cora les, rodolitos y foraminíferos (^{Halfar et al., 2000}; ²⁰⁰⁴).

Materiales y métodos

Trabajo de campo y laboratorio. Dentro del área de estudio (Figura 1b), se programaron 85 sitios de muestreo (70 aleatorios y 15 dirigidos) a través del sistema de información geográfica (SIG) Esri ArcGIS 10.3. También se generó una muestra independiente de 18 sitios para evaluar los modelos de interpolación. El trabajo de campo se realizó mediante embarcaciones equipadas con sistema de navegación satelital y con dragas de tipo Van Veen, desde febrero hasta mayo de 2016. Este perio do se eligió previo a la temporada de ciclones tropicales para evitar que las muestras estuvieran bajo la influencia del efecto hidrodinámico que éstos ocasionan. Las muestras se tomaron del sedimento superficial del fondo marino y fueron trasladadas al Laboratorio de Geología del IPN-CICIMAR. Adicionalmente, en cada sitio se tomó video submarino SCUBA y videos asistidos con una carcasa desde la borda para regis trar la presencia de organismos conspicuos y estructuras bentónicas. Tanto las muestras de sedimento como los videos se obtuvieron en un intervalo entre 1 a 70 m de profundidad. Para determinar el origen biogénico de los sedimentos gruesos, se realizaron observaciones con microscopios estereoscópicos directamente sobre las muestras.

En el laboratorio, con una balanza analítica se pesaron entre 3 y 4 g de cada muestra, las cuales fueron tratadas con ácido clorhídrico diluido al 10%. Luego, por diferencias de pesos se calculó el contenido porcentual del CaCO₃ (^{Halfar et al., 2000}). Por otro lado, una muestra de aproximadamente 27 g fue fraccionada en tres grupos: sedimento muy grueso (-2.25 a -0.75 Φ), grueso (-0.75 a 4 Φ) y fino (>4 Φ), con tami ces estandarizados. La unidad Φ es una transformación del diámetro de las partículas del sedimento (Φ = -log₂d, donde d es el diámetro de la partícula en mm) (^{Blott & Pye, 2001}). La fracción muy gruesa fue tratada por método de tamizado ro-tap, con intervalos de 1/4 Φ. Las fracciones gruesa y fina se trataron mediante un módulo seco y otro húmedo utilizando un analizador de partículas de difracción láser, mo delo Beckman Coulter LS13320, en intervalos de 1/8 Φ. Con los datos se calculó el tamaño medio de grano y la clasificación según el método logarítmico de los momentos con el programa Gradistat (^{Blott & Pye, 2001}). A partir de los datos obtenidos para cada sitio y sus coordena das, se procedió a modelar la distribución espacial de cada variable.

Modelación de la distribución espacial (interpolaciones). Los mo delos se construyeron con el algoritmo de funciones de base radial del SIG (^{Buhmann, 2003}). Los modelos fueron evaluados mediante la com paración de los valores pronosticados contra los medidos en la muestra independiente. En la comparación se calculó la media del error (ME) y la raíz cuadrada del error cuadrático medio (RCEC), (^{Luo et al., 2008}). Adicionalmente, se estimó un intervalo de confianza de la diferencia de medianas entre cada modelo y las observaciones de la muestra independiente, basado en la prueba de rango con signo de Wilcoxon (^{Hollander et al., 2014}). El modelo de profundidad se elaboró con la información disponible en cartas batimétricas escala 1:20,000 (^{DGID, 1996}; ²⁰⁰⁸). Todas las variables se manejaron en formato ráster y tie nen la misma alineación y tamaño de pixel, de 100 x 100 m. El formato de los modelos permite caracterizar el comportamiento espacial de cada variable en el fondo de la BLP.

Clasificación y caracterización de unidades del fondo marino. Los procedimientos para la clasificación de unidades del fondo marino se resumen en la Figura 2. Para cada modelo, los valores de los pixeles fueron agrupados según los intervalos reportados en la literatura (Tabla 1). Para el tamaño medio de grano y la clasificación se utilizaron los límites de ^{Folk y Ward (1957}) y ^{Friedman y Sanders (1978}). Para la profundidad se utilizaron los contornos usados por ^{Harris et al. (2014}), donde el límite de 100 m representa la profundidad máxima de la pla taforma continental (^{Nava-Sánchez et al., 2001}). El modelo de CaCO₃ se dividió en cinco clases equidistantes con base en los siguientes cri terios: 1) el número de clases debe ser similar al resto de las variables para evitar sesgo, 2) todos los valores posibles fueron incorporados, esto es de 0 a 100% de CaCO₃, y 3) los valores de los límites de las cinco clases y su distribución geográfica coinciden razonablemente con los valores del método de rupturas naturales Jenks (^{De-Smith et al., 2015}). Estas divisiones permitieron establecer los límites espaciales en el área de estudio. Los límites se editaron en el SIG para generar polígonos con los valores clasificados. Los polígonos de cada modelo se sobrepusieron y se cortaron unos con otros para obtener una matriz de n polígonos con p variables. Posteriormente, se realizó un análisis de agrupamiento jerárquico (multivariado) con el software Minitab v.17. El análisis requiere que las variables tengan el mismo rango para evitar sesgos; por lo tanto, se utilizaron los valores transformados entre 1 y 5 (Tabla 1). En el análisis se utilizó la distancia euclidiana al cuadrado y el método de aglomeración de Ward (^{Kaufman & Rousseeu, 2005}). A continuación, se agruparon los polígonos en clases (conglomerados), en función de la distancia multivariada. Las unidades del fondo marino se determinaron por medio de la inspección de los conglomerados en el dendrograma y, según la configuración espacial de los mismos, en diferentes niveles de agrupación. Finalmente, el perímetro de cada uni dad agrupó un conjunto de pixeles de cada modelo. Los valores de los pixeles permitieron caracterizar numéricamente las unidades.

Figura 2 Procedimientos para la clasificación de unidades del fondo marino a partir de la distribución espacial de cuatro variables (contenido de CaCO₃, tamaño medio de grano, clasificación del sedimento y profundidad). 

Tabla 1 Límites de los intervalos de clases de las variables utilizadas para la clasificación de unidades del fondo marino. 

T=variables transformadas jerárquicamente de 1 a 5. LMP= límite máximo de plataforma.

Resultados

Modelos de distribución espacial y cálculos del error. La distribu ción del contenido de CaCO₃, el tamaño medio de grano, la clasificación del sedimento y la profundidad (batimetría) se presentan en la Figura 3. El modelo de CaCO₃ presenta valores entre 2.7 y 95% (Figura 3a). Los valores menores que 20% se distribuyen ampliamente sobre la pla taforma interna de la BLP, al norte de la barra arenosa El Mogote. Los valores mayores que 60% se encuentran bordeando estrechamente la zona de Punta Coyote. En el canal de San Lorenzo se reportan valores mayores al 80%. Los valores mínimos y máximos del modelo de tama ño medio de grano (Figura 3b) fueron de -0.9 a 7.8 Φ, esto es, entre arena muy gruesa hasta limo fino. Los valores bajos (tamaños gruesos) se presentan en el canal de mareas de la Laguna de La Paz (-0.5 a 1.6 Φ) y en algunas áreas bordeando la zona de Punta Coyote, sobre todo en el canal de San Lorenzo, donde alcanza el tamaño de arena muy gruesa. Los valores más altos (tamaños finos) se encuentran en la parte interna de la laguna de La Paz y al noroeste del área de estudio. La clasificación del sedimento varía desde 0.1 hasta 4 Φ (Figura 3c). Los valores bajos se distribuyen a lo largo de la línea de costa, al norte de la barra arenosa El Mogote (0.5 a 0.7 Φ), en la plataforma interna y también en el canal de San Lorenzo (0.7 a 1 Φ). Los valores altos se ubican en la parte central de la laguna de La Paz y al noroeste del área de estudio.

Figuras 3a-d Modelos de distribución espacial. a) Contenido de CaCO3. b) Tamaño medio de grano. Los valores mínimos representan los tamaños gruesos, y los máxi mos, los tamaños finos. c) Clasificación del sedimento. Los valores mínimos representan una buena a moderada selección (baja variabilidad), y los valores máximos representan una mala a muy mala selección (alta variabilidad). d) Profundidad (batimetría). Las líneas representan los límites de la Tabla 1. 

La evaluación de los modelos a través de la media del error (ME) y la raíz cuadrada del error cuadrático medio (RCEC) se muestra en la Figura 4. La estimación de la diferencia de medianas entre los valores pronosticados y los medidos con la muestra independiente fue de 0.5 a 6.4% (contenido de CaCO₃), -0.5 a 0.2 Φ (tamaño medio de grano) y -0.4 a 0.0 Φ (clasificación del sedimento) para intervalos de confianza del 95%.

Figuras 4a-c Evaluación de los modelos. Comparación entre los valores pro nosticados y los medidos con la muestra independiente. a) Contenido de CaCO₃. b) Tamaño medio de grano (en Φ). c) Clasificación del sedimento (en Φ). R²= coeficiente de determinación. ME= media del error, RCEC= raíz cuadrada del error cuadrático medio. La línea discontinua representa una estimación no ses gada 1:1. 

Clasificación y caracterización de unidades del fondo marino. Al sobreponer los límites de los cuatro modelos de distribución (Figura 3, Tabla 1) se estableció un total de 413 polígonos. Los polígonos se agrupa ron como se muestra en la Figura 5a. Se eligió la agrupación para cuatro clases (conglomerados) con base en el nivel de distancia euclidiana, con valores de 77.8, 108.7, 175.8 y 197.5. Las cuatro clases fueron sub-divididas, a su vez, en 14 unidades, según su separación geográ fica (Figura 5b). La agrupación de polígonos estableció cuatro principales unidades (1a, 2a, 3a y 4a), que representan el 91% del área de estudio en el fondo marino de la BLP. Los valores de los pixeles agrupados para cada variable dentro de cada unidad se resumen en la Tabla 2.

Figuras 5a-b Clasificación de unidades del fondo marino. a) Agrupación de polígonos. b) Configuración espacial de las 14 unidades del fondo marino. Las cuatro clases fueron divididas en unidades según el análisis de agrupamiento y su configuración geográfica. 

Tabla 2 Caracterización de unidades. Valores del conjunto de pixeles obtenidos en los modelos dentro de cada unidad clasificada. 

ID= identificador de la unidad del fondo marino (ver Figura 5); se presenta el promedio (μ), la desviación estándar (σ) y el rango (R) del conjunto de pixeles agrupados según los límites espaciales de cada clase.

La unidad 1a se caracteriza por tener una profundidad promedio de 61 ± 25 m y un tamaño medio de grano de 4.5 ± 0.6 Φ, corres pondiente a limo grueso y limo muy grueso. La unidad 2a presenta una profundidad de 19 ± 10 m, una pobre a moderada clasificación del se dimento (1.1 ± 0.3 Φ) y muy bajo contenido de CaCO₃ (13 ± 4.2%). En las unidades 1a y 2a es común observar orificios en el sedimento realizados por organismos de la infauna bentónica. La unidad 3a se ca racteriza por una muy mala clasificación del sedimento (2.2 ± 0.7 Φ) y baja profundidad (2.2 ± 1.9 m). En el fondo marino se observaron sitios con algas (Sargassum spp.) adheridas a valvas de moluscos bivalvos o rocas (guijarros) y, en ocasiones, sobre corales muertos. Algunos bival vos observados en el fondo fueron Chione californiensis, Chione sp., Abra tepocana, Laevicardium sp., Nuculana sp. y los gasterópodos Nassarius luteostoma y Cerithium stercusmuscarum. La unidad 4a se encuentra a 13 ± 10 m de profundidad y presenta el mayor contenido de CaCO₃ (72 ± 20%) y el mayor tamaño medio de grano, con arenas finas a gruesas (1.4 ± 0.8 Φ). La unidad 4a presenta una pobre a moderada clasificación del sedimento (1 ± 0.3 Φ). Dentro de la unidad, en los sitios 62, 18 y 0-independiente (Figura 1b) se observó la presencia de mantos de rodolitos, mientras que en los sitios 30, 43 y 49 se observaron colonias de corales.

Discusión

Evaluación de los modelos de distribución espacial. En la BLP es escasa la información publicada relacionada con las variables modela das en este estudio; sin embargo, se pueden hacer comparaciones con estudios puntuales. En el canal de San Lorenzo (Figura 1b), ^{Halfar et al. (2000}; ²⁰⁰⁴) reportan que los valores de CaCO₃ en el sedimento marino fueron mayores a 90% a lo largo de dos transectos. Estos valores son parcialmente consistentes con los resultados de este estudio, el cual reporta entre 75 y 95%. Las diferencias pueden ser atribuidas a los diferentes enfoques metodológicos utilizados, aunque los autores men cionados utilizaron transectos y no propiamente un mapeo (variabilidad espacial). La ventaja de los modelos basados en interpolaciones es que predicen valores entre puntos conocidos en el espacio, lo que permite determinar una distribución. Por otro lado, el tamaño medio de grano y la clasificación en cuatro sitios muestreados por ^{Cruz-Orozco et al. (1996}) dentro del área de estudio son consistentes con los modelos desarrollados. Los autores presentan un esquema conceptual basado en una densidad de muestreo de 0.01 muestras km^-2, mientras que el presente estudio, basado en 0.2 muestras km^-2, aporta mayor detalle de la variabilidad espacial y utiliza un algoritmo de interpolación.

La información generada en los modelos y el uso de un SIG permi ten establecer los límites o intervalos numéricos (Tabla 1) en el espacio geográfico. La calidad de estos límites depende del tamaño de muestra y del análisis del error de los modelos. En este aspecto, el error de los modelos propuestos ha sido calculado mediante una muestra in dependiente. Esto permite la comparación de futuros modelos, incluso con diferentes algoritmos de interpolación. El modelo del contenido de CaCO₃ sobreestima a los valores observados con un promedio de 3.5% (ME). En valores absolutos, cada punto tiene una diferencia promedio de 6.4% (RCEC, Figura 4). Los modelos del tamaño medio de grano y clasificación subestimaron a los valores observados con ME = -0.12 Φ y -0.19 Φ, respectivamente. En valores absolutos de cada punto con trastado, se estima una diferencia promedio (RCEC) de 0.68 Φ (tamaño medio de grano) y 0.39 Φ (clasificación del sedimento). Los modelos predicen valores cercanos en promedio a las observaciones realizadas con la muestra independiente. Este trabajo aporta los primeros modelos de la distribución espacial y el cálculo de su error (ME, RCEC) del conte nido de CaCO₃, tamaño medio de grano y clasificación del sedimento en el fondo de la BLP, con base en algoritmos de interpolación.

Unidades del fondo marino. El análisis de agrupamiento jerárquico permitió agrupar los 413 polígonos en 14 unidades; sin embargo, cua tro unidades representan el 91% del área estudiada. La agrupación geográfica preferencial de los polígonos en cuatro unidades se debe a que la distancia multivariante (o similitud) entre éstos es próxima. Una posible explicación es que en el fondo de la BLP prevalecen prin cipalmente cuatro conjuntos de condiciones ambientales, las cuales se infieren a partir de la caracterización de los valores de las cuatro variables modeladas y su distribución espacial. La unidad 4a agrupó polígonos con pixeles de alto contenido en CaCO₃ (72 ± 20%). Los va lores elevados pueden ser atribuidos a la tasa baja de aporte de sedi mentos siliciclásticos en la zona de Punta Coyote y, por consiguiente, a la alta producción de exoesqueletos carbonatados de los organismos que habitan dicha unidad (^{Halfar et al., 2004}). El aporte de terrígenos proviene de la erosión de rocas en la línea de costa y, principalmente, por aporte de sedimentos desde las cuencas hidrológicas adyacentes (^{Halfar et al., 2000}, ²⁰⁰⁴; ^{Velasco-García, 2009}) (Figura 2b). Estas cuen cas son pequeñas y sólo en eventos de precipitación alta (asociado con ciclones tropicales) se consigue un aporte de sedimentos siliciclásticos. El aporte de carbonatos en el fondo marino de la unidad 4a se deriva de los bioclastos de fragmentos de conchas (moluscos), corales, algas rojas (rodolitos), foraminíferos y balanos. Estos componentes han sido reportados previamente por ^{Halfar et al. (2000}). Específicamente en el canal de San Lorenzo se ha reportado que el sedimento marino está dominado por mantos de rodolitos (^{Halfar et al., 2004}); sin embargo, en la inspección del sedimento en el fondo, se observaron mantos de rodolitos sólo en los sitios 62 y 0-independiente. Adicionalmente, se en contró un manto de rodolitos en el sitio 18, no reportado en la literatura, en una zona con valores entre 65 y 70% de CaCO₃. El tamaño de grano de arenas finas a gruesas (1.4 ± 0.8 Φ) en el fondo indica niveles de energía en el ambiente (por corrientes, oleaje y marea) más altos que el resto de las unidades estudiadas. La moderada a pobre clasifica ción del sedimento (1 ± 0.3 Φ) se asocia al poco aporte de terrígenos y dominancia de bioclastos carbonatados, los cuales mantienen una dispersión moderada del tamaño medio de grano. La unidad 4a por su contenido alto en CaCO₃ asociado a moluscos, corales y mantos de ro dolitos podría destacar en biodiversidad sobre el resto de las unidades estudiadas.

Los polígonos agrupados en la unidad 3a coinciden, en su mayor parte, con la zona de la laguna de La Paz (Figura 5b). La unidad registra una muy mala clasificación del sedimento (2.2 ± 0.7 Φ); es decir, una amplia variabilidad del tamaño medio de grano, que se ve reflejada en una variedad de tipos de sedimento (desde limos gruesos hasta arenas medias). La variabilidad es consistente con la reportada por ^{Álvarez-Arellano et al. (1997}) y se puede explicar por las diferentes fuentes de partículas que llegan a la unidad. La laguna de La Paz recibe el aporte de materiales limosos y arenosos a través de los arroyos, desde diferentes cuencas hidrológicas (Figura 1a) y diferentes fuentes litológi cas, así como por aporte litoral y eólico (^{Nava-Sánchez & Cruz-Orozco, 1989}; ^{Choumiline et al., 2009}). Las descargas hacia la laguna traen consigo materiales finos y nutrientes que son distribuidos en la ense nada (^{Mendoza-Salgado et al., 2005}) y el fondo marino (^{Green-Ruiz & Páez-Osuna, 2003}; ^{Pérez-Triboullier, 2014}). La poca profundidad de la unidad 3a (2.2 ± 1.9 m) permite que el efecto de la energía del viento (el oleaje interno) en la laguna tenga una fuerte influencia sobre ella. Este fenómeno genera la redistribución y suspensión de los sedimen tos. Aunque la unidad contiene en promedio 24 ± 7.8 de CaCO₃, los mo luscos (bivalvos y gasterópodos) que habitan en el sedimento aportan bioclastos mayores que 4 mm. Estos bioclastos tienen un efecto en el incremento del tamaño medio y en la dispersión del tamaño de grano. La amplia variedad de tipos de sedimento en el fondo de la unidad 3a establece las condiciones para el hábitat de diversos macromoluscos, como bivalvos y gasterópodos (^{Pérez-Nevarez, 1995}; ^{Holguín-Quiñones & García-Domínguez, 1997}; ^{Félix-Pico et al., 2011}). En la laguna de La Paz (unidad 3a) se ha reportado la presencia de diversos grupos de moluscos; algunos son más comunes en esta zona lagunar, en compa ración con el resto de la bahía de La Paz.

La unidad 2a forma parte de la plataforma interna de la BLP y el ca nal de marea de la laguna de La Paz. La unidad está influenciada por la energía del oleaje, sobre todo en la parte próxima hacia la barra areno sa El Mogote. En esta franja el transporte litoral distribuye el sedimento con el nivel de energía necesario para sostener el tamaño de grano de arenas muy finas (^{Nava-Sánchez & Cruz-Orozco, 1989}). La moderada clasificación del sedimento (1.1 ± 0.3 Φ) sugiere que las condiciones hidrodinámicas se mantienen regulares a lo largo de la barra. El tamaño medio de sedimento (Figura 3b) incrementa en relación con la alta energía hidrodinámica de las corrientes de marea hacia el canal que conecta la BLP con la ensenada de la laguna de La Paz (^{Gómez-Valdés et al., 2003}). Los valores bajos (<20%) de CaCO₃ en esta unidad indican la poca producción de bioclastos. Esto contribuye a que el sedimento ma rino esté dominado por material siliciclástico aportado por las cuencas hidrológicas de la margen occidental de la BLP (Figura 1a), de acuerdo con ^{Nava-Sánchez y Cruz Orozco (1989}) y ^{Cruz-Orozco et al. (1996}). La arena muy fina permite el desarrollo de organismos bentónicos infau nales, los cuales generan orificios y galerías en el fondo marino. Por lo tanto, esta unidad es potencialmente un hábitat para comunidades de macrofauna bentónica (^{Gray, 2002}).

Los polígonos agrupados en la unidad 1a corresponden a la parte profunda (61 ± 25 m) de la plataforma de la BLP dentro del área de estudio. Es la unidad con menor tamaño medio de grano (4.5 ± 0.6 Φ) y se puede atribuir a que la zona tiene la menor o nula influencia de la energía del oleaje. Por consiguiente, representa un área para el depósito de sedimentos terrígenos (^{Nava-Sánchez et al., 2001}) o por hundimiento de partículas en suspensión en la columna de agua (^{Sil verberg et al., 2014}). La unidad 1a registra un contenido de CaCO₃ entre 20 y 40% influenciado por la proximidad de la unidad 4a, el aporte de bioclastos de foraminíferos, fragmentos de conchas de moluscos bentónicos y nanoplancton calcáreo (^{Urcádiz-Cázares, 2005}). En esta unidad también se observaron orificios y galerías hechas por organis mos de la infauna bentónica, pero sólo hacia la parte más profunda. La unidad representa un hábitat para micro y macroinvertebrados marinos típicos de una plataforma continental (^{Gray, 2002}).

Este estudio representa una de las primeras propuestas metodo lógicas para la clasificación de unidades en el fondo marino dominado por sedimentos no consolidados en la BLP. A diferencia de otros méto dos, se toma en cuenta la distribución espacial de cuatro variables, las cuales son indicadores directos e indirectos de condiciones ambienta les del fondo marino. Los resultados, aunque no definitorios, permiten establecer para la BLP 14 unidades en el fondo marino, de las cuales cuatro corresponden al 91% de la superficie estudiada. Estas unidades sugieren que el fondo de la BLP está dominado por al menos cuatro principales conjuntos de condiciones ambientales que se ven reflejados en la distribución y las características propias de cada unidad.

Una de las limitaciones de este estudio es que los organismos que habitan en cada una de las unidades no fueron considerados de for ma sistemática. Por lo tanto, fue asumido que algunos grupos marinos podrían asociarse a las unidades propuestas, como lo muestran las observaciones puntuales realizadas durante la toma de muestras de sedimento. Por ello, recomendamos estudios ecológicos específicos, para validar si existen grupos de organismos marinos relacionados con las unidades y sus condiciones ambientales, considerando una escala espacial apropiada (1:20,000 a 1:50,000) y la variación estacional de la abundancia de cada especie. De igual manera, debe ser evaluado el efecto de los ciclones tropicales sobre la distribución de las variables estudiadas y en su caso, las correspondientes unidades del fondo ma rino.

Finalmente, esta propuesta debe tomarse como una alternativa robusta respecto a las metodologías utilizadas en el mapeo de hábi tats bentónicos, tipos desedimento y geomorfología del fondo marino. Se recomienda como línea base en estudios relacionados con orde namientos ecológicos marinos y planes de manejo. La propuesta es accesible económicamente y puede utilizarse en las zonas costeras y las plataformas continentales del Golfo de California y otras regiones marinas.