El sondeo A (Fig. 1, A-B; longitud: 36º56´22´´N; latitud: 6º19´16´´W) se extrajo mediante un testificador manual de 20 mm de diámetro en las marismas fluvio-mareales situadas al sur del chenier de Las Nuevas, dentro del Parque Nacional de Doñana. Su litología y color (escala Munsell) fueron revisados directamente sobre el terreno (Ruiz et al., 2004). No se han tomado muestras de los últimos 47 cm del testigo estudiado por estos autores.

Se obtuvieron un total de 23 muestras de sedimento (15 g) para su análisis faunístico. Las muestras fueron levigadas a través de un tamiz de 63 μm de diámetro de malla, con un tratamiento previo mediante adición de pirofosfato sódico. El residuo fue secado a una temperatura constante de 60ºC.

Para el estudio de los foraminíferos bentónicos, se ha utilizado un microsplitter para dividir el residuo resultante en partes iguales. Cuando fue posible, al menos 300 individuos fueron extraídos de cada muestra e identificados de acuerdo con los criterios de Loeblich & Tappan (1987), Martins & Gomes (2004), Milker & Schmiedl (2012) y diversos trabajos referenciados en ellos, así como con el World Register of Marine Species. Las especies más abundantes (Fig. 2) fueron fotografiadas con un microscopio electrónico de barrido en el Departamento de Microscopía Electrónica y Rayos X de la Universidad de Huelva.

En cada muestra, se calcularon el número de especies y el número de individuos por gramo de arena. Además, se han obtenido diversos índices de diversidad (Shannon -H-, equitability -J-, α-Fisher), así como la relación entre formas planctónicas y bentónicas (P/B, en %).

Este análisis constó de los siguientes pasos: a) cálculo de los porcentajes de las especies halladas, aproximados al entero más cercano en el texto; b) determinación de los coeficientes de correlación entre las nueve especies más abundantes, que representan a más del 99% de todos los especímenes, para poder determinar las principales asociaciones de estas especies; y c) aplicación de un análisis clúster en modo Q (método de unión completa) para obtener agrupaciones de muestras según la distribución de las nueve especies citadas. En este análisis, los grupos se han diferenciado a partir de una distancia euclídea transformada de 10 unidades. Métodos semejantes y distancias similares o levemente superiores (10-20 unidades) son seleccionados usualmente para la identificación de grupos de muestras en estudios de foraminíferos bentónicos (p.e. Polovodova et al., 2009; Al-Dubai et al., 2017).

Se ha incluido una datación de Ruiz et al. (2004) obtenida en el nivel arenoso bioclástico del sondeo A (muestra A-10), así como otras dataciones de edad similar realizadas en niveles bioclásticos depositados por eventos de alta energía en los estuarios de los ríos Guadalquivir y Tinto-Odiel (Ruiz et al., 2004; González-Regalado et al., 2019a). Los resultados se presentan como edades calibradas (2σ), con la aplicación del efecto reservorio (-108 ± 31 ¹⁴C yr) calculado por Martins & Soares (2013) para este área. De forma indirecta, la edad de los tramos arcillosos inferior y superior de este testigo puede ser calculada aproximadamente si se extrapola la tasa media de sedimentación similar (~1,5 mm/año), obtenida para el Holoceno medio y superior en este sector por Lario et al. (2002) en dos testigos situados en las proximidades del lucio de Mari López y por Pozo et al. (2010) en un testigo situado en el sector suroccidental del chenier de Las Nuevas (Fig. 1, A).

El tramo basal del sondeo A (Fig. 1, B) está compuesto por unos 50 cm de arcillas limosas gris-azuladas (10YR 5/3) con un escaso contenido arenoso y frecuentes restos vegetales. Contiene un abundante registro de bivalvos, representados principalmente por valvas desarticuladas de Ruditapes decussatus (Linnaeus) en los primeros 35 cm del sondeo y esta especie y Cerastoderma edule (Linnaeus) en los 15 cm superiores del tramo basal.

A techo de estos sedimentos finos se sitúa un nivel fino de arenas limosas bioclásticas (8 cm; 10Y 7/10). En su base erosiva, aparecen frecuentes valvas desarticuladas de bivalvos, con Magallana angulata (Lamarck) y Ruditapes decussatus (Linnaeus) como especies más representativas. En sus 4 cm superiores, la densidad de los bivalvos disminuye, con menor representación de estas dos especies y una presencia esporádica de Acanthocardia tuberculata (Linnaeus).

Los 49 cm superiores del sondeo A están constituidos por arcillas limosas verdosas (5Y 4/2), masivas o con una leve laminación paralela. La malacofauna (C. edule, M. angulata) solo está presente en los primeros 20 cm, en tanto que esporádicos caparazones del gasterópodo Rissoa sp. aparecen en los 24 cm finales del sondeo.

El sondeo A muestra una relativa homogeneidad en los porcentajes de sus principales especies. H. germanica y A. tepida, las dos especies dominantes (Tab. 1: 92,9% del total entre ambas) son especies eurihalinas que soportan grandes variaciones de salinidad (<32⁰/₀₀) y condiciones de estrés ambiental (ver revisión en Murray, 2006). Coexisten en medios litorales y parálicos, con especial abundancia en lagos y lagunas costeras, marismas, estuarios, bahías, etc. En estos ambientes, H. germanica es más resistente a cambios ambientales y prefiere en general zonas más restringidas, de sustrato limo-arcilloso y con una mayor elevación topográfica, como las zonas intermareales. A. tepida es más abundante en áreas más abiertas, más profundas y con mayor renovación mareal, como los canales mareales, áreas más próximas a la bocana de estas lagunas o incluso zonas infralitorales someras (Jorissen, 1988; Debenay & Guillou, 2002; Blázquez, 2005; Frontalini et al., 2011, 2013).

El notable predominio de H. germanica, los bajos índices de diversidad observados, la escasa presencia y progresiva desaparición de formas planctónicas, así como la paulatina disminución del número de individuos, apuntan a la presencia de un medio lagunar intermareal situado en el sector marino de un estuario (Cearreta, 1998; Murray, 2006). Estos datos, unidos a los importantes porcentajes de A. tepida, así como a la ausencia de especies típicas de las marismas actuales de los estuarios suratlánticos de España (González-Regalado et al., 2001, 2019; Ruiz et al., 2005), como Trochammina inflata (Montagu) o Entzia macrescens (Brady), o de sus zonas más externas, como Ammonia beccarii (Linnaeus), situarían al sondeo A en una llanura mareal interna o en una zona poco profunda de una laguna sometida a un importante estrés medioambiental (p.e., variaciones de salinidad, periodos prolongados de emersión) durante el periodo romano. Esta reconstrucción paleogeográfica del sector meridional del estuario se ha observado en lagunas atlánticas y mediterráneas actuales (p.e. Ruiz et al., 2006; Cearreta et al., 2007; Coccioni et al., 2009; Calvo-Marcilese & Langer, 2010; Alday et al., 2013). Esta reconstrucción paleogeográfica coincide con los textos históricos y con la inferida en otros estudios pluridisciplinares (Ruiz et al., 2010; Borja Barrera, 2013; Borja Barrera et al., 2018).

La presencia del resto de especies minoritarias también está en concordancia con esta interpretación paleogeográfica. Las especies más abundantes del género Elphidium (E. cuvillieri -como Cribroelphidium cuvillieri-, E. oceanense -como E. oceanensis-), B. ordinaria y P. granosum (como Elphidium granosum o Cribroelphidium granosum) son comunes en lagunas costeras restringidas, pero con un aún apreciable intercambio mareal (Cearreta et al., 2002; Frontalini et al., 2011; Alves et al., 2015). Estas especies también pueden estar presentes en ambientes marinos someros (Mendes, 2010).

Los elevados porcentajes de H. germanica y A. tepida proporcionan datos paleoambientales adicionales, en conjunción con la presencia de numerosos restos vegetales. H. germanica se asocia a una moderada disponibilidad de microfitobentos (p.e. diatomeas) en el medio, del que extrae los cloroplastos que permanecen funcionales durante al menos dos semanas mediante el fenómeno conocido como cleptoplastia (Jauffrais et al., 2016; Wukovits et al., 2018). Las preferencias alimenticias de A. tepida son más amplias, ya que esta especie se alimenta de detritos, cianobacterias o materia orgánica refractaria (Goldstein & Corliss, 1994; Moodley et al., 2000; De Nooijer, 2007). En consecuencia, H. germanica prefiere sedimentos más oxigenados, con alto Eh, en tanto que A. tepida es más tolerante a cierto nivel de hipoxia, derivado del enriquecimiento en materia orgánica (Bouchet et al., 2007; Martins et al., 2015).

En este escenario general relativamente permanente, el nivel bioclástico arenoso representa un aumento puntual de la energía del medio, cuyas consecuencias más relevantes son: a) una erosión de las arcillas infrayacentes; b) un incremento del tamaño de grano; c) la presencia de una abundante malacofauna (M. angulata, R. decussatus) típica de medios marinos (Ferreira & Vale, 1998; Chicharo et al., 2001; Gómez, 2017) con evidencias tafonómicas de transporte (fracturas, desarticulación), una de las características de depósitos holocenos de tormentas (Ruiz et al., 2007); d) una disminución del número de foraminíferos bentónicos, observada también en registros pliocénicos de tormentas (González-Regalado et al., 2009); e) un leve incremento relativo de formas planctónicas; f) una asociación dominante similar a las arcillas inferiores, aunque con un leve incremento de H. germanica (~65% a ~74%) y una disminución de A. tepida (23,5% a 17,3%); y g) una ausencia casi total de foraminíferos bentónicos marinos someros. Todo ello apunta a la acción de una tormenta que erosionó los sedimentos de fondo de las zonas más externas de la laguna y los depositó en los sectores más internos ya descritos. La edad más probable de este evento (finales del siglo I d.C.) coincide con un periodo de tormentas en este sector atlántico, así como en la zona mediterránea adyacente (Sabatier et al., 2012; Orme et al., 2017), que ocasionaron el depósito de abanicos de derrame (washover fans) sobre cheniers previos en la cercana Isla de Saltés (Fig. 1-Tab. 3; González-Regalado et al., 2019).

Los ostrácodos del sondeo A fueron estudiados dentro de un análisis multidisciplinar de varios sondeos situados en el sector meridional del Parque Nacional de Doñana (sondeo C de Ruiz et al., 2004). Entre las muestras A-1 y A-20, sus asociaciones están dominadas por Cyprideis torosa (Jones), Loxoconcha elliptica (Brady) y Leptocythere castanea (Sars), especies típicas de zonas internas muy someras de lagunas salobres (salinidad inferior a 25^o/oo en su mayoría) con conexión marina (Marocco et al.,1996; Montenegro & Pugliese, 1996; Ruiz et al., 2000). Además, la progresiva reducción de la influencia marina también queda reflejada en los ostrácodos, con una disminución del número de individuos hacia techo acompañada de una desaparición paulatina de especies marinas, que aparecen puntualmente en el tramo arcilloso inferior del sondeo (Carinocythereis whitei (Baird), Palmoconcha spp., Semicytherura sp.). También desaparecen las abundantes espinas de equinodermos observadas en la parte basal y media del sondeo A. En ostrácodos, el nivel bioclástico también queda diferenciado por una disminución en el número de individuos y un aumento en el número de adultos en relación con el tramo arcilloso inferior. En consecuencia, la reconstrucción paleoambiental general del sondeo A es similar en función de las asociaciones de foraminíferos y ostrácodos.

La principal diferencia entre ambos análisis se sitúa en las muestras A-21 a A-23. No se observan cambios sustanciales en las asociaciones de foraminíferos (Fig. 2 y Tab. 1) ni su análisis estadístico separa a estas muestras (Fig. 3), si bien disminuye su número de individuos. Por el contrario, la asociación de ostrácodos cambia paulatinamente hasta consolidarse un dominio de formas dulceacuícolas. En este tramo superior, es muy abundante Cyprinotus salinus (Brady) y también aparecen de forma puntual Cypridopsis vidua (Müller), Cypris bispinosa Lucas, Herpetocypris chevreuxi (Sars) e Ilyocypris gibba (Ramdohr), una asociación límnico-oligohalina típica de aguas estancadas, lagunas, lagos, marismas o incluso zonas muy internas de estuarios (Margalef, 1983; Bronsthein, 1988; Fernandes et al., 2010). Esta asociación está actualmente presente en las lagunas temporales interiores, lucios y caños del Parque Nacional de Doñana (Armengol, 1976; Ruiz et al., 1996; Alcorlo et al., 2014). De acuerdo con las tasas de sedimentación indicadas, esta transición hacia medios dulceacuícolas se produjo durante los siglos IV-V d.C.

En consecuencia, los ostrácodos detectan con mayor precisión los cambios paleoambientales en esta área, y especialmente las oscilaciones de la salinidad. Esta conclusión también ha sido obtenida en otros estudios comparativos entre ostrácodos y foraminíferos como indicadores (paleo-) ambientales de actuales y antiguas lagunas costeras (p.e. Ruiz et al., 2012; Barbieri & Vaiani, 2018).