En este trabajo se estudian dos secciones, la de Arroyo Grande (Figs. 1 y 2A) y la de Arroyo Valfondo (Figs. 1 y 2B), descritas por Murelaga et al. (2012), en las que se identifican las unidades morfosedimentarias del Holoceno final (UM 2, 3 y 4) definidas por Sancho et al. (2008). La unidad UM 2 está formada por arenas y limos (Fig. 2), con intercalaciones esporádicas de conglomerados en la sección de Arroyo Grande (Murelaga et al., 2012). Esta unidad se depositó cubriendo una superficie de erosión irregular desarrollada sobre los materiales previos, y puede ser interpretada como un depósito fluvial en una secuencia de incisión y relleno de canales. La unidad UM 3 está constituida por arenas con yesos, alternantes con lutitas y limos laminados, que representan secuencias de incisión y relleno de canales con intercalaciones de depósitos de sedimentos finos posiblemente transportados por eventos esporádicos de lluvias torrenciales (Murelaga et al., 2012). Finalmente, la UM 4 aparece en pequeños afloramientos, estando fuertemente relacionada con los procesos actuales de erosión de barrancos y el relleno de canales incisos previos. En la sección de Arroyo Grande está compuesta por arenas intercaladas con limos (Fig. 2A), mientras que en la sección de Arroyo Valfondo está conformada por limos masivos y laminados (Fig. 2B) (Murelaga et al., 2012). Esta unidad representa el depósito de materiales transportados por paleoinundaciones.

La edad de estas secciones fue estimada por Murelaga et al. (2012) y está basada en nueve dataciones absolutas con ¹⁴C realizadas en restos de carbón (Fig. 2; Tabla 1; ver Murelaga et al., 2012 para detalles sobre el procesado de las muestras y la calibración de las edades obtenidas). De acuerdo a estas dataciones, las muestras analizadas abarcarían el intervalo temporal comprendido entre 4763±87 años cal. BP y 127±82 años cal. BP (Fig. 2; Tabla 1), es decir, desde el final del Holoceno medio hasta el Reciente.

En este trabajo se han estudiado 41 muestras: 30 de ellas obtenidas de la sección de Arroyo Grande (Fig. 2A) y las 11 restantes en la sección de Arroyo Valfondo (Fig. 2B).

En cada nivel seleccionado se han recogido aproximadamente 500 gramos de sedimento, con el objetivo de disponer de un número de ejemplares de ostrácodos que sea representativo de las asociaciones preservadas, dada la pobreza microfaunística del sedimento. Este sedimento ha sido lavado y tamizado, analizándose la fracción comprendida entre 150 µm y 1 mm de diámetro, de la que se ha extraído la totalidad de los individuos adultos y del último estadio juvenil (A-1) presentes en cada muestra. La fracción superior a 1 mm no se ha considerado ya que no contenía ejemplares clasificables de ostrácodos. Tanto las valvas desarticuladas como los ejemplares con caparazón completo se consideran como un único individuo. Debido al bajo número de ejemplares obtenidos en las muestras estudiadas, aquellas especies que aparecen representadas por más de 4 individuos por muestra se consideran “abundantes” en este trabajo.

Para efectuar la identificación taxonómica, a nivel de especie, de los ejemplares extraídos, se han seguido las clasificaciones definidas por Hartman & Puri (1974) y Horne et al. (2002), completadas con los trabajos de Kempf (1980–1997), Meisch (2000) y Fuhrmann (2012). Para las especies del género Ilyocypris se han tenido en consideración las descripciones de Janz (1994), Karanovic & Lee (2013) y Mazzini et al. (2014), mientras que para el género Paralimnocythere se ha seguido el trabajo de Martens ( 1992).

En los sedimentos fluvio-aluviales, como los que se estudian en este trabajo, es común el depósito de abundantes fragmentos de valvas de ostrácodos, fruto del transporte post-morten de las mismas por los flujos tractivos canalizados (e.g. Poquet & Mesquita-Joanes, 2011). Estos fragmentos pueden ser transportados largas distancias, produciendo acumulaciones mezcladas con las asociaciones in-situ una vez que estos flujos se desbordan en la llanura de inundación y depositan el sedimento acarreado. De esta manera, la asociación estudiada no permitiría efectuar una reconstrucción paleoecológica fiable del ambiente de depósito. Esta situación se ha observado en el presente trabajo, puesto que en las 41 muestras estudiadas (Fig. 2) son abundantes los fragmentos de valvas de ostrácodos, inidentificables siquiera a nivel genérico.

Para evitar esta eventualidad, en este estudio únicamente se han considerado los ejemplares que no presentan evidentes signos de transporte (e.g. fracturas, roturas, abrasiones o superficies de erosión) para poder así realizar una interpretación adecuada de las condiciones ambientales en el momento de depósito. Además, la presencia de valvas y caparazones de ejemplares adultos y juveniles en la misma asociación, junto con una proporción similar de mudas en cada muestra, apoya la suposición de la ausencia de transporte en estos individuos, validando la reconstrucción paleoecológica realizada en este trabajo.

De las 30 muestras estudiadas, estratigráficamente sucesivas (Fig. 2A), siete de las mismas no presentaban ostrácodos clasificables (niveles AG–7, 8, 13, 14, 16, 18 y 30). En el resto de los niveles analizados se han extraído un total de 131 ejemplares, pertenecientes a 17 especies (Tabla 2). Tanto la abundancia de ostrácodos como las especies más comunes identificadas en las tres Unidades Morfosedimentarias (UM 2, UM 3, UM 4) descritas en esta sección, son diferentes (Fig. 4).

Unidad Morfosedimentaria UM 2 (4763±87 a 2848±55 años cal. BP). Se identifica en la base de la sección e incluye desde el nivel AG–1 al nivel AG–12, donde la especie más abundante es I. bradyi, estando acompañada en algunos niveles por P. albicans (Fig. 4; Tabla 2).

Unidad Morfosedimentaria UM 3 (2848±55 a 836±65 años cal. BP). Se reconoce en la parte central de la sección y engloba las muestras desde AG–13 a AG–19 (Fig. 4). Esta unidad es muy pobre en ostrácodos, encontrándose ejemplares clasificables únicamente en tres de las muestras estudiadas (AG–15, 17 y 19), donde se registra una mezcla de escasos individuos pertenecientes a I. bradyi o P. albicans junto con otros del género Cyclocypris (Fig. 4; Tabla 2).

Unidad Morfosedimentaria UM 4 (191±97 a 127±82 años cal. BP). Se define hacia el techo de la sección (muestras desde AG–20 hasta AG–30) (Fig. 4). En la base de este intervalo aparecen algunos ejemplares de Candona sp. 1 (nivel AG–20) y Heterocypris incongruens (Ramdohr) (nivel AG–21), mientras que en el resto de las muestras estudiadas que componen esta UM 4 es más abundante I. bradyi, estando acompañada por Candona sp. 2 (Fig. 4; Tabla 2).

Se han estudiado 11 muestras recogidas a lo largo de esta sección (Fig. 2B), obteniéndose 228 ejemplares de ostrácodos agrupados en 9 especies (Tabla 2). De nuevo se observan diferencias en la abundancia relativa de las especies más comunes en las tres Unidades Morfosedimentarias identificadas (UM 2, UM 3, UM 4) (Fig. 5).

Unidad Morfosedimentaria UM 2 (4763±87 a 2848±55 años cal. BP). Se localiza en la base de la columna, englobando únicamente dos de los niveles muestreados (AV–1 y 2), donde aparece como especie más abundante I. bradyi, junto con Paralimnocythere aff. psammophila y P. albicans (Fig. 5; Tabla 2).

Unidad Morfosedimentaria UM 3 (2848±55 a 836±65 años cal. BP). Es la que presenta un mayor espesor en la sección de Arroyo Valfondo (Fig. 5) e incluye los niveles de AV–3 a AV–8. La especie más abundante en todo este intervalo es Paralimnocythere aff. psammophila, aunque también se encuentran ejemplares de P. albicans e I. bradyi a lo largo de esta unidad (Fig. 5; Tabla 2).

Unidad Morfosedimentaria UM 4 (191±97 a 127±82 años cal. BP). Se identifica a techo de la sección y engloba los niveles AV–9 a AV–11 (Fig. 5). En este último intervalo disminuye el número de ostrácodos, estando presentes las especies I. bradyi y P. albicans junto con Paralimnocythere aff. psammophila, todas ellas con escasos ejemplares.

Los depósitos sedimentarios correspondientes a este intervalo temporal (UM 2) presentan un registro más amplio en la sección de Arroyo Grande (Fig. 2A). La asociación de ostrácodos identificada en las muestras recogidas en este intervalo está dominada por I. bradyi (Fig. 4; Tabla 2). Esta especie es indicativa de ambientes con agua corriente, tales como pequeños arroyos, riachuelos o manantiales helocrenos y reocrenos (Henderson, 1990; Roca & Baltanás, 1993; Külköylüoğlu, 1999; Wu et al., 2001; Mezquita et al., 2005), aunque también puede aparecer en charcas temporales o ambientes lacustres de aguas estancadas (Henderson, 1990) asociada a la entrada ocasional de agua corriente desde arroyos cercanos (Curry, 1999). A lo largo de la UM 2, y en especial a techo de dicha unidad en Arroyo Valfondo (Figs. 4 y 5; Tabla 2), aparece P. albicans. Esta especie es típica de bordes de lagos con cierta corriente interna del agua (Meisch, 2000; Mezquita et al., 2005) (Tabla 3), por lo que su aparición confirmaría la presencia de aguas corrientes en esta área para este período.

Esta UM 2 está constituida por varias secuencias solapadas de canales arenosos y depósitos de desbordamiento (ver Fig. 2A) indicativos del desarrollo de un sistema fluvial de canales activos y áreas encharcadas (Murelaga et al., 2012). Estas secuencias sedimentarias identificadas en este intervalo temporal en Bardenas Reales de Navarra, podrían representar la respuesta de pequeños sistemas fluviales ante unas condiciones ambientales caracterizadas por una distribución relativamente regular de las precipitaciones, unas laderas fitoestabilizadas y escasez de eventos hidrológicos extremos, lo que se traduce en una cierta estabilidad geomorfológica y en el desarrollo de paleosuelos (Murelaga et al., 2012). Los estudios palinológicos efectuados en Bardenas Reales de Navarra, definen el desarrollo de un bosque ribereño relacionado con flujos de agua permanente en la depresión de Bardena Blanca en 3495±35 años BP (Iriarte & Meaza, 1996), en torno a la Época Fría del Hierro. Así mismo, los estudios faunísticos (gasterópodos) e isotópicos previamente realizados, permiten inferir una tendencia hacia condiciones ambientales húmedas y frías durante este periodo en esta zona (Murelaga et al., 2012). Por otro lado, Thorndycraft & Benito (2006) han descrito un descenso en las avenidas durante el final del Holoceno medio en la Península Ibérica que apunta a una fase de cierta estabilidad hidrológica.

La interpretación paleoambiental establecida a partir de la asociación de ostrácodos identificada en este trabajo es concordante con estas observaciones previas, puesto que define el desarrollo de un sistema de pequeños arroyos o riachuelos con agua corriente en una llanura de inundación entre 4763±87 y 2848±55 años cal. BP (UM 2).

A pesar de que ambas secciones estudiadas presentan una potencia similar de sedimentos depositados durante este intervalo temporal (UM 3), es en las muestras recogidas en Arroyo Valfondo (Fig. 2B) donde se ha encontrado un mayor número de ostrácodos clasificables.

En todos los niveles correspondientes a este intervalo, Paralimnocythere aff. psammophila es la especie dominante (Fig. 5; Tabla 2). Considerando unos requisitos ecológicos para esta especie similares a los descritos para P. psammophila, es posible efectuar una interpretación paleoambiental para este intervalo temporal. Esta especie es típica de cuerpos de agua dulce temporales o zonas litorales de ambientes lacustres sin apenas corrientes de agua (Meisch, 2000; Mezquita et al., 2005; Gascón et al., 2009) (Tabla 3). Por tanto, implicaría el establecimiento de un medio acuático efímero y de aguas estancadas. La aparición, tanto en la Sección de Arroyo Grande como en la de Arroyo Valfondo, de ejemplares de I. bradyi y de P. albicans, especies características de aguas corrientes, indicaría episodios esporádicos de precipitaciones irregulares, tal vez torrenciales, que aumentarían el hidrodinamismo en dichos cuerpos de agua estancada.

Los materiales que conforman la UM 3 se organizan en secuencias arenosas granodecrecientes, que presentan yeso, bioturbaciones y acumulaciones de gasterópodos y restos vegetales en la base de cada una de ellas (ver Fig. 2B). Esta unidad se interpreta como un periodo de elevada tasa de actividad aluvial, relacionada con lluvias esporádicas y torrenciales en condiciones climáticas más áridas que las acontecidas durante el intervalo temporal previo (Murelaga et al., 2012). En valles holocenos de la parte central de la Cuenca del Ebro se han descrito secuencias aluviales similares (Peña et al., 2004; Constante et al., 2010). Además, en materiales lacustres muestreados en las Sierras Exteriores de los Pirineos, Morellón et al . (2009) identifican una fase árida acontecida entre 870 y 750 años cal. BP.

El intervalo temporal de este depósito aluvial coincide con un incremento en la erosión de laderas ocurrida durante los periodos Romano y post-Romano, dentro de una tendencia general hacia la aridificación que culmina en el “Periodo Cálido Medieval” o “Anomalía Climática Medieval” (Moreno et al., 2012). Durante este periodo (datado aproximadamente en 800–1200 d.C.) las temperaturas de algunas zonas europeas alcanzaron los valores máximos de los últimos 4000 años (Barclay et al., 2009). Estas condiciones climáticas explicarían la disminución en el hidrodinamismo del sistema aluvial en este periodo con respecto al intervalo previo en esta área, tal y como han evidenciado las especies de ostrácodos identificadas. Por otro lado, en la parte central de la Cuenca del Ebro se ha constatado que estos procesos aluviales pueden acentuarse localmente debido a la actividad humana, ya sea por un aumento del pastoreo o por la deforestación (Iriarte & Meaza, 1996; Peña et al., 2004; González-Sampériz et al., 2009), por lo que no se puede despreciar la influencia antrópica como explicación de la variación en las condiciones del medio aluvial detectada en este estudio.

Este intervalo temporal (UM 4) presenta una mayor potencia de sedimentos depositados en la sección de Arroyo Grande (Fig. 2A). De manera general, la asociación específica de ostrácodos identificada en esta unidad está dominada por I. bradyi (Fig. 4; Tabla 2), lo que implicaría, nuevamente, el establecimiento de un medio aluvial con agua corriente en el Reciente. La presencia de algunos ejemplares de otras especies menos abundantes a lo largo de esta Unidad Morfosedimentaria definiría cambios de menor escala temporal dentro de este intervalo. Así, en la base de la sección de Arroyo Grande aparece H. incongruens (Tabla 2), lo que implica el desarrollo de pequeños cuerpos de agua temporales con cierta corriente interna (Meisch, 2000; Mezquita et al., 2005). Por su parte, a techo de Arroyo Valfondo abunda Paralimnocythere aff. psammophila (Tabla 2), indicando la presencia de cuerpos de agua estancada en la llanura de inundación aluvial.

A partir del análisis de las asociaciones de gasterópodos y la señal isotópica de sus conchas, Murelaga et al. (2012) identificaron el establecimiento de unas condiciones relativamente frías y con alternancias de fases húmedas y secas en el periodo Reciente en Bardenas Reales de Navarra. Estos autores aluden a la variabilidad climática acontecida durante la “Pequeña Edad de Hielo” (Little Ice Age, LIA) como posible causante del incremento en las condiciones de frío y los cambios hidrológicos extremos de menor escala deducidos durante este intervalo temporal en esta área. Este evento, datado aproximadamente en 1500–1850 d.C., consiste en un primer avance de los glaciares continentales entre 1180–1320 d.C. (Barclay et al., 2009), seguido por dos etapas frías principales, de aproximadamente un siglo de duración cada una de ellas, acontecidas en los siglos XVII y XIX de nuestra era, con dos fases relativamente cálidas ocurridas en los siglos XVI y XVIII (Bradley & Jones, 1992; Barclay et al., 2009). Las fases frías supusieron un nuevo avance de los glaciares en Europa, Asia y Norteamérica, con una expansión del hielo oceánico en el Atlántico Norte (Lamb, 1982). Sin embargo, algunos autores han puesto en duda que se trate de un evento global (e.g. Bradley & Jones, 1992). En el norte de la Península Ibérica, la LIA se caracteriza por un incremento en la humedad y un descenso en la temperatura (Valero-Garcés et al., 2008; Morellón et al., 2009), con una alta variabilidad climática de menor escala que incluye fases de inundaciones extremas y periodos de sequía (e.g. Álvarez Vázquez, 1986; Benito et al., 1996; Saz, 2003). Estas condiciones climáticas serían las causantes de la influencia de un agua corriente, posiblemente asociada a arroyos o riachuelos de pequeña entidad, alternante con momentos en los que se desarrollan cuerpos de agua efímeros y estancados durante el periodo Reciente en Bardenas Reales de Navarra, tal y como se ha constatado en la reconstrucción paleoambiental efectuada a partir de las asociaciones de ostrácodos. Sin embargo, tal y como sucedía en el intervalo temporal previo, no se puede descartar la posible influencia antrópica, consistente principalmente en prácticas de deforestación, durante la LIA en esta área (Murelaga et al., 2012).