En estudios micropaleontológicos los métodos estadísticos permiten establecer relaciones y ordenar conjuntos amplios de variables (especies o muestras) que pueden ser interpretadas en función de parámetros o variables paleoambientales (se remite a Murray, 1991, 2006; Shi, 1993; Martin, 2000 para una revisión del tema). De entre los muchos tipos de análisis estadísticos que se pueden implementar, en este trabajo se ha optado por utilizar un análisis de correspondencia (Benzécri, 1992), ya que éste permite visualizar en un mismo gráfico bidimensional el conjunto de datos obtenidos. El análisis de correspondencia sin tendencia (DCA, por sus siglas en inglés Detrended Correspondence Analysis), también denominado rectificado o segmentado, es habitualmente empleado en ecología (Jongman et al., 1995) e incluso en reconstrucciones paleobatimétricas basadas en foraminíferos (e.g. Herkat, 2007; Herkat y Ladjal 2013). Este tipo de análisis se utiliza básicamente para la ordenación de las especies (categorías taxonómicas) y muestras en un mismo gráfico bidimensional. Particularmente, éste se emplea cuando parte de las asociaciones están caracterizadas por frecuencias muy bajas o muchos ceros (Hammer y Harper, 2006). Básicamente el algoritmo DCA genera ejes o factores no correlacionados que maximizan la dispersión de los valores de las especies. La principal ventaja que ofrece esta técnica estadística es que permite visualizar en una misma gráfica las variables y los casos, es decir, las categorías taxonómicas y las muestras. Así pues, el conjunto de datos puede ser analizado de forma integral, estableciendo similitudes y diferencias entre el total de muestras y categorías taxonómicas consideradas.

La pasada década ha sido testigo de un avance considerable en el conocimiento estratigráfico, sedimentológico y micropaleontológico de la Cuenca del Bajo Segura, una de las cuencas neógenas más representativas del Mediterráneo occidental en cuanto a intervalo temporal registrado (desde el Tortoniense al Cuaternario) y variedad de paleoambientes marinos y medios de transición representados (principalmente batial, circalitoral, infralitoral y laguna litoral) (Soria et al., 2008a, 2008b, Caracuel et al., 2011, García-García et al., 2011; Soria et al., 2014; Corbí et al., 2016; Corbí y Soria, 2016; Gamonal et al., 2017; Martínez-Martínez et al., 2017; Corbí et al., 2018, entre otros). Sin embargo, hasta la fecha los análisis estadísticos multivariantes basados en datos cuantitativos de asociaciones fósiles no han sido explorados en esta cuenca. Este tipo de análisis estadístico supone una excepcional herramienta para abordar el análisis paleoecológico y paleoambiental del registro micropaleontológico, tal como ponen de manifiesto las investigaciones realizadas en otros registros neógenos del Mediterráneo (e.g. Gliozzi y Grossi, 2008, Gennari et al., 2013, Walid et al., 2017). Es por ello, que en el presente trabajo se aborda por primera vez el análisis estadístico multivariante (en particular el DCA, mencionado anteriormente) del registro estratigráfico marino y medios de transición, básicamente lagunas litorales de edad Messiniense y Plioceno de la Cuenca del Bajo Segura, tomando como punto de partida las investigaciones previas anteriormente comentadas. La finalidad de este artículo es, por tanto, examinar estadísticamente el registro estratigráfico desde la perspectiva de las unidades aloestratigráficas o sintemas, sistemas de depósito y asociaciones de foraminíferos registrados en dicha cuenca.

Los materiales analizados pertenecen a la Cuenca del Bajo Segura (Figura 1), una de las cuencas post-orogénicas de la Cordillera Bética (Viseras et al., 2004) que incluye un registro sedimentario desde el Tortoniense hasta la actualidad (Corbí y Soria, 2016). Esta cuenca, ubicada geológicamente en el extremo oriental de la Cordillera Bética, se localiza sobre el contacto entre los dos dominios principales de esta cordillera, las Zonas Externas y las Zonas Internas. El límite norte de la cuenca coincide con la Falla de Crevillente mientras que el sur no está definido dado que esta cuenca se encuentra en continuidad con la de Cartagena. Hacia el oeste los relieves de la Sierra de Orihuela definen el límite con la Cuenca de Fortuna, mientras que, hacia el este, la cuenca se abre hacia el Mediterráneo.

En esta cuenca es posible distinguir los sectores, norte y sur, claramente definidos desde el punto de vista tectónico y paleogeográfico (Figuras 1 y 2). Ambos sectores están separados por el valle fluvial del río Segura, cuya génesis esta relacionada con la actividad reciente de la Falla del Bajo Segura. En el sector norte el basamento de la cuenca lo forman los materiales de las Zonas Externas, mientras que en el sector sur está constituido por materiales de las Zonas Internas. La estructura de estos sectores presenta notables diferencias, de modo que en el sector norte los materiales definen un dispositivo monoclinal con frecuentes discordancias progresivas asociadas al movimiento de la Falla de Crevillente y al levantamiento de los relieves de las Zonas Externas. En cambio, en el sector sur, los materiales se encuentran afectados por pliegues, entre los que destaca, por sus proporciones, el Anticlinal de Torremendo. Desde el punto de vista paleogeográfico estos dos sectores también muestran ciertas diferencias. En el sector norte las unidades aloestratigráficas presentan, el líneas generales, asociaciones de facies típicas de medios marinos someros y continentales; en cambio, en el sector sur, estas unidades se caracterizan por presentar sistemas de depósito marinos de mayor profundidad y espesor (Corbí, 2010; Corbí y Soria, 2016).

Trabajos recientes dividen el registro marino de esta cuenca en 5 unidades aloestratigráficas o sintemas: Tortoniense I, Tortoniense II, Tortoniense-Messiniense I, Messiniense II y Plioceno (Soria et al., 2008a, 2008b, Corbí y Soria, 2016). Estas unidades quedan, por tanto, separadas por discontinuidades representadas por superficies erosivas que corresponden con cambios paleogeográficos y tectonosedimentarios. Las figuras 3 y 4 ilustran el mapa geológico sintético, así como la organización estratigráfica para los sectores norte y sur de esta cuenca. En todo el registro sedimentario se ha documentado la completa sucesión de marcadores bioestratigráficos de foraminíferos planctónicos, desde la biozona MMi9 (Tortoniense inferior) al MPl3 (Plioceno inferior) (Corbí y Soria, 2016). Se refiere a Soria et al., 2008a, 2008b, Corbí et al. (2016), Corbí y Soria (2016) y Corbí (2010) para una descripción exhaustiva de las principales características estratigráficas, sedimentológicas y micropaleontológicas de las unidades representadas en esta cuenca. Es de destacar que, cada una de las unidades definidas contienen diferentes conjuntos de litofacies y biofacies y, por lo tanto, comprenden más de una unidad litoestratigráfica y bioestratigráfica. De acuerdo con estos autores, el análisis de los diferentes conjuntos de facies que integran estas unidades ha permitido diferenciar sistemas de depósito (sensu Fisher y McGowan, 1967) que cubren un amplio espectro de ambientes sedimentarios, desde batiales a lagunares (Figuras 3 y 4).

En este trabajo se ha realizado un análisis de correspondencia sin tendencia (DCA) en el conjunto de datos de distribución de foraminíferos representados en 139 muestras de 6 secciones estratigráficas. Para llevar a cabo este procedimiento estadístico se empleo el programa informático PAST© (Hammer et al., 2001). Una primera caracterización tafonómica de la asociación registrada pone de manifiesto que las conchas de la mayor parte de los foraminíferos analizadas presentan un buen estado de conservación (Corbí, 2010). Por tanto, en líneas generales las asociaciones pueden considerarse como acumuladas (sensu Fernández-López, 2000). No obstante, en algunos casos puntuales se han detectado foraminíferos que presentan signos de alteración como fracturas, perforaciones, perdida de transparencia y formación de concreciones lo cual podría ser indicativo de resedimentación y/o reelaboración.

Se eligió el análisis de correspondencia sin tendencia, ya que es el habitualmente utilizado cuando se tiene presencia de varias especies en las asociaciones de foraminíferos caracterizadas por bajas frecuencias. Además, permite representar en un mismo gráfico muestras y grupos ecológicos de foraminíferos bentónicos y planctónicos, lo cual facilita un análisis integral de toda la cuenca. Para ello, se han considerado los datos cuantitativos de distribución de asociaciones de foraminíferos bentónicos y planctónicos representados en los sistemas de depósito marinos y de medios de transición, principalmente lagunares, de edad Messiniense y Plioceno de las secciones representativas del sector norte (Elche, Torrellano y Colmenar) y del sector sur (Garruchal, La Pedrera y El Canal) de la Cuenca del Bajo Segura (véase Figuras 5 y 6 para localización de las muestras analizadas). Las muestras fueron previamente levigadas de acuerdo al procedimiento convencional, tras lo cual se analizó semicuantitativamente la fracción contenida entre los tamices de 125 y 500 μm, contabilizando y determinando a nivel de especie (siempre que fuera posible) un total de 300 ejemplares de foraminíferos por muestra (Corbí, 2010).

Con el propósito de analizar de forma integral, completa y sintética el registro marino y lagunar litoral representado, se han establecido grupos ecológicos de foraminíferos bentónicos (GEFB) y planctónicos (GEFP), tomando como punto de partida los datos cuantitativos de distribuciones de foraminíferos (porcentajes de representación) presentados en Corbí (2010). Los grupos ecológicos establecidos corresponden a conjuntos de especies de foraminíferos característicos de un determinado ambiente marino. Por tanto, en la matriz de datos generada para el análisis DCA, las asociaciones establecidas (grupos ecológicos) representan los “casos”, y las muestras, las “variables”. Esta agrupación ha permitido visualizar e integrar los datos a distintas escalas, como pueden ser los sectores de la cuenca, sintemas, secciones estratigráficas y sistemas de depósito (Figura 7).

Los grupos ecológicos de foraminíferos bentónicos (GEFB) considerados, expresados en la Tabla 1, están relacionados tanto con el ambiente marino como con el microhábitat (modo de vida). En esta tabla se muestran las especies representadas en la Cuenca del Bajo Segura agrupadas según el morfogrupo y microhábitat al que pertenecen. Respecto al microhábitat conviene señalar que los dos grupos establecidos, epifaunal-infaunal somero e infaunal profundo, son los propuestos por Jorissen et al. (1995). No obstante, y con el propósito de sintetizar los resultados, los grupos ecológicos de foraminíferos bentónicos (tabla 1) han sido agrupados y sintetizados en representativos de ambientes de tipo lagunar (L), infralitoral (I), circalitoral (C) y batial (B). El ambiente sedimentario característico de cada especie ha sido obtenido de la información autoecológica bibliográfica disponible, fundamentalmente de Murray (1991 y 2006). Ahora bien, cabe señalar que, para ciertos autores, algunas de estas especies poseen un rango posible de distribución más amplio del que se ha considerado en este trabajo, en el que sólo se ha tenido en cuenta su ambiente de más habitual. Así pues, las distribuciones de cada uno de estos grupos representan, en términos generales, el ambiente predominante en cada una de las muestras analizadas.

Por otra parte, el contenido en oxígeno es uno de los factores fundamentales que controlan la distribución de los foraminíferos bentónicos, además de encontrarse relacionado directamente con la salinidad, temperatura y circulación de la masa de agua (Van der Zwaan et al., 1999). Numerosos estudios ponen de manifiesto una relación directa entre el grado de oxigenación de las aguas y el modo de vida de los foraminíferos bentónicos, de tal forma que, generalmente, en condiciones de baja oxigenación predominan, en mayor o menor grado, los taxones de carácter infaunal profundo, denominados también como marcadores de estrés ambiental (Hinsbergen et al., 2005). Así pues, de acuerdo con este autor, el conjunto de los foraminíferos bentónicos puede subdividirse en taxones de aguas marinas abiertas “normales” y marcadores o indicativos de estrés ambiental (stress markers). Este último grupo, de gran importancia en la reconstrucción de las concentraciones de oxígeno, está constituido por las especies siguientes: Bolivina spp. (excepto B. plicatella Cushman, 1930 y B. pseudoplicata Heron- Allen & Earland, 1930), Bulimina spp. (excluidas las especies costuladas); Cancris spp., Chistomella spp., Fursenkoina spp., Globobulimina spp., Rectuvigerina spp., Stainforthia spp., Uvigerina spp. (se excluye U. semiornata) y Valvulineria spp. En definitiva, en el presente trabajo, con propósito de conocer el grado de estrés ambiental sobre el fondo de cada unidad, se ha considerado como variable el porcentaje de distribución de los foraminíferos bentónicos marcadores de estrés ambiental S, en sentido de Hinsbergen et al. (2005) que representa el porcentaje de marcadores de estrés ambiental sobre el fondo con respecto al total de la asociación de foraminíferos bentónicos.

Los grupos ecológicos de foraminíferos planctónicos (GEFP) considerados están relacionados directamente con la profundidad. Según Bè (1977) se pueden establecer tres grupos principales de foraminíferos planctónicos de acuerdo con su distribución batimétrica (aguas someras, intermedias y profundas). Atendiendo a la clasificación propuesta por este autor, en este trabajo se han considerado los siguientes grupos: a) especies de aguas someras (<50 metros de profundidad), grupo que incluye todas las especies de los géneros Globigerinoides y Turborotalita; b) especies de aguas intermedias (entre 50 y 100 metros de profundidad). Formado por las especies Globigerina bulloides d’Orbigny, 1826, Globigerinella spp., Orbulina universa d’Orbigny, 1839 y Globigerinita naparimaensis Brönnimann, 1951; y c) especies de aguas profundas (>100 metros de profundidad) representado por todas las especies de Globorotalia, Sphaeroidinellopsis spp., y Neogloboquadrina acostaensis (Blow, 1959).

El análisis multivariante realizado (DCA) ha permitido distribuir las muestras (casos) y los grupos ecológicos (variables: bentónicos, GEFB, planctónicos, GEFP) en un mismo gráfico bidimensional, expresado en la figura 7. Estos resultados estadísticos permiten establecer similitudes y diferencias entre el conjunto de datos estudiados, atendiendo principalmente a las siguientes cuestiones: a) distribución y dominancia de grupos ecológicos en los sintemas representados en la cuenca: Tortoniense-Messiniense, T-MI, en figura 7 en rojo, Messiniense II, MII, en amarillo y Plioceno, P, en verde; b) comparación sectores norte y sur de la cuenca (secciones Elche, Colmenar y Torrellano en el sector norte; secciones de Garruchal, La Pedrera y Canal en sector sur); y c) rasgos particulares (distribución de grupos ecológicos) de las secciones analizadas y sistemas de depósito representados.

La figura 7 muestra la ordenación bidimensional del total de datos analizados, de modo que se puede realizar una serie de afirmaciones estadísticas. A continuación se exponen y analizan las principales cuestiones representadas en el grafico, discutiendo los resultados en el contexto de los trabajos publicados hasta la fecha que abordan el estudio estratigráfico y micropaleontológico de la Cuenca del Bajo Segura.

Este estudio pretende demostrar como el análisis multivariante, en particular el análisis de corres-pondencia sin tendencia, representa un instrumento estadístico eficaz para el análisis integral de toda una cuenca sedimentaria basado en asociaciones de foraminíferos. Este trabajo supone un primer punto de partida para la exploración futura de otros pro-cedimientos estadísticos de detalle que contemplen otros aspectos relevantes en la interpretación de paleoambientes según las asociaciones de foraminíferos registradas, como, por ejemplo, las variaciones en el estado de conservación (gradientes tafonómicos) o la variedad intraespecífica entre distintas especies.

Este análisis estadístico, implementado aquí por primera vez para el registro marino de la Cuenca del Bajo Segura, permite extraer las siguientes conclusiones, siendo todas ellas contrastadas con la información bibliográfica de carácter estratigráfico, sedimentológico y micropaleontológico disponible: