Modelação cronológica de sequências sedimentares do Estuário da Ribeira de Bensafrim com recurso à estatística Bayesiana

MODELAÇÃO CRONOLÓGICA DE SEQUÊNCIAS SEDIMENTARES DO ESTUÁRIO DA RIBEIRA DE BENSAFRIM COM RECURSO À ESTATÍSTICA BAYESIANA Paulo J. Cesário Portela1 José M. Matos Martins1 antónio M. Monge soares1 ana ramos-Pereira2 João araújo-Gomes3 andré torres3 RESUMO – O estuário da ribeira de Bensafrim, no Barlavento algarvio, foi um dos sistemas flúviomarinhos utilizado para avaliar a evolução de ambientes costeiros ao longo dos últimos 5000 anos. Datou-se pelo radiocarbono a matéria orgânica sedimentar presente nos testemunhos verticais colhidos neste estuário. tendo em conta o registo estratigráfico associado aos dados cronológicos obtidos, foi utilizada uma abordagem estatística bayesiana, fazendo uso do software de calibração OxCal, para obter um enquadramento geocronológico robusto para a sequência sedimentar. a análise efetuada, baseada na construção de um modelo de deposição para a sequência sedimentar, permitiu identificar a presença de outliers, bem como estabelecer um enquadramento geocronológico para os últimos 5000 anos no Barlavento algarvio. Para além deste enquadramento, verifica-se que o modelo geocronológico idealizado, associado à caracterização sedimentológica de um dos testemunhos, aponta para a presença de eventos de alta energia, englobando, eventualmente, o tsunami que se seguiu ao terramoto de Lisboa de 1755. Palavras-chave: radiocarbono; estatística bayesiana; estuário da ribeira de Bensafrim; sedimentos; tsunami de 1755. recebido: janeiro 2019. aceite: 1 Centro de Ciências e tecnologias nucleares (C2tn), instituto superior técnico, Universidade de Lisboa, estrada nacional 10 (km 139,7), 2695-066, Bobadela, Portugal. e-mail: paulo[email protected]; [email protected]; amsoar[email protected] 2 Professora Catedrática e investigadora efetiva do Centro de estudos Geográficos, instituto de Geografia e Ordenamento do território, Universidade de Lisboa, Lisboa, Portugal. e-mail: [email protected] 3 Centro de estudos Geográficos, instituto de Geografia e Ordenamento do território, Universidade de Lisboa, Lisboa, Portugal. e-mail: jpag[email protected]m; andre.roch.torr[email protected] Finisterra, LIV(112), 2019, pp. 113 ‑130 issn: 0430-5027 doi: 10.18055/finis16084 artigo 114 Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 ABSTRACT – aGe MODeLLinG Of seDiMentarY seQUenCes UsinG a BaYesian aPPrOaCH – tHe BensafriM estUarY seDiMentarY infiLLinG seQUenCe as a Case stUDY. The Bensafrim estuary, located in southern Portugal (algarve coast), was selected to assess the evolution of the fluvial-marine interaction along the Portuguese coast during the last 5000 years. sedimentary sequences were collected in this estuary and radiocarbon dated, providing the necessary chronological framework to establish a reliable palaeo-environmental reconstruction. a Bayesian approach was employed in an attempt to overcome some difficulties related with the chronological data, which was performed using the OxCal calibration software, which allows the incorporation of stratigraphic information and an outlier analysis of the dataset. The analysis enabled the identification of outliers and the assessment of a reliable deposition model for the sedimentary sequence, establishing a robust chronological framework for the last 5000 years in this region of the Portuguese coast. also, this model associated with sedimentological characterization points to the presence of high-energy events, including eventually the tsunami that followed the Lisbon earthquake of 1755. Keywords: radiocarbon; Bayesian statistics; Bensafrim estuary; sediments; Lisbon tsunami 1755. RÉSUMÉ – MODÉLisatiOn CHrOnOLOGiQUe sÉDiMentaire À L’aiDe De La statistiQUe BaYÉsienne – Une ÉtUDe De Cas: La sÉQUenCe De reMPLisaGe De L’estUaire De riBeira De BensafriM. L’estuaire de la ribeira de Bensafrim, en algarve auprès de Lagos, a été utilisé comme un exemple de système fluvio-marin pour évaluer l’évolution des environnements côtiers au cours des 5 000 dernières années. La matière organique sédimentaire présente dans les échantillons verticaux collectés dans cet estuaire a été datée au radiocarbone. Compte tenu de l’enregistrement stratigraphique associé aux données chronologiques obtenues, une approche statistique bayésienne a été utilisée avec le logiciel de calibrage OxCal, pour obtenir un cadre géochronologique robuste pour la séquence sédimentaire. L’analyse, basée sur la construction d’un modèle de dépôt pour la séquence sédimentaire, a permis d’identifier la présence de valeurs aberrantes (outliers), ainsi que d’établir un cadre géochronologique pour les 5000 dernières années sur l’Ouest de la côte méridional. Par ailleurs à ce cadre, le modèle géochronologique idéalisé, associé à la caractérisation sédimentologique, montre plusieurs évènements de haute énergie responsable par l’ensablement de l’estuaire, y compris éventuellement le tsunami qui a suivi le tremblement de terre de 1755 à Lisbonne. Mots clés: radiocarbone; statistique bayésienne; estuaire de Bensafrim; sédiments; tsunami de 1755. RESUMEN – MODeLaCiÓn CrOnOLÓGiCa De seCUenCias seDiMentaLes DeL estUariO De La riBera De BensafriM COn Base en La estaDÍstiCa BaYesiana el estuario de la ribera Bensafrim, situado al oeste de algarve, fue uno de los sistemas fluviales y marinos utilizado para evaluar la evolución de ambientes costeros en los últimos 5000 años. a través del método de radiocarbono, se dató la materia orgánica sedimentaria presente en los núcleos verticales recogidos en este estuario. Considerando el registro estratigráfico asociado con los datos cronológicos obtenidos, se utilizó un enfoque estadístico bayesiano, utilizando el software de calibración OxCal, para obtener un marco geocronológico robusto para la secuencia sedimentaria. el análisis, basado en la construc- 115Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 ción de un modelo de deposición para la secuencia sedimentaria, permitió identificar la presencia de valores atípicos (outliers), así como establecer un marco geocronológico para los últimos 5000 años en el Barlavento algarvio. Más allá de este marco, se verifica que el modelo geocronológico idealizado, asociado con la caracterización sedimentológica de uno de los núcleos, apunta a la presencia de eventos de alta energía, posiblemente abarcando el tsunami luego del terremoto de Lisboa de 1755. Palabras clave: radiocarbono; estadisticas bayesianas; estuario de ribeira de Bensafrim; sedimentos; tsunami de 1755. Na parte vestibular dos vales portugueses (…) seria necessário precisar melhor as modalidades e o ritmo da colmatagem sofrida pelos diversos estuários, distinguindo a sedimentação fluvial da marinha e da propriamente estuarina e, datando, na medida do possível, os episódios principais. S. Daveau (1987. p.101) i. intrODUÇÃO Um enquadramento cronológico fiável é essencial para a construção de modelos paleoambientais e paleoclimáticos baseados em análises multi-proxy de registos sedimentares. Devido ao limite de idade (ca. 50000 anos) e à variedade de amostras susceptíveis de serem datadas pelo radiocarbono (14C), este é um dos métodos mais utilizados para a construção de esqueletos cronológicos onde ancorar a variabilidade dos diversos parâmetros identificados no registo sedimentar. em ambientes estuarinos, onde a interação flúvio-marinha é mais intensa, a remobilização de material orgânico afeta frequentemente as datas de 14C, o que dificulta a construção de um enquadramento cronológico robusto. a estatística bayesiana oferece as ferramentas necessárias para interpretar com fiabilidade conjuntos de datas, permitindo identificar inversões estratigráficas e possíveis outliers. Uma vez que, entre nós, a estatística bayesiana tem sido muito pouco utilizada neste domínio, julgamos útil que o objetivo deste trabalho seja o de exemplificar este tipo de análise com um estudo de caso, o do enchimento sedimentar do estuário da ribeira de Bensafrim, para o qual foram obtidos dois conjuntos de datas de 14C, correspondentes a duas sondagens, Bf3 (Portela, 2013) e LagosConv (araújo-Gomes, 2010), efetuadas nesse estuário. ii. O estUÁriO Da riBeira De BensafriM e a reCOLHa Da infOrMaÇÃO no âmbito da investigação da evolução de sistemas estuarinos nos últimos 5000 anos, selecionaram-se diversos sistemas, entre os quais o estuário da ribeira de Bensafrim, 116 com o objetivo de individualizar o balanço dos elementos forçadores, sejam de natureza climática, marinha ou antrópica (ramos-Pereira et al., 2011). a ribeira de Bensafrim é um pequeno curso de água com 16,2km de comprimento e a sua bacia alongada tem uma área de 85,1km2 (fig. 1). Desenvolve-se em duas unidades morfoestruturais, de litologias variadas: turbiditos paleozóicos a norte, na serra algarvia, e complexos detríticos e carbonatados mesozóicos e cenozóicos no troço médio e jusante. também do ponto de vista geomorfológico se distingue: um vale muito encaixado na serra, que se alarga quando atravessa a depressão marginal desenvolvida no Grés de silves, à saída no Maciço antigo (Maciço Hespérico). a influência litológica é marcante na morfologia e, em especial, na largura do fundo do vale, que apenas se estreita quando corta as formações calcárias. O estuário desta ribeira tem, assim, uma forma em funil, cuja foz original tinha a largura de 550m, mas está hoje completamente artificializado e o canal de saída do rio tem apenas 62m de largura (fig.1). fig. 1 – Bacia de drenagem da ribeira de Bensafrim e vista aérea do estuário e da cidade de Lagos. O quadrado branco indica o local de amostragem. Fig. 1 – The Bensafrim drainage basin and aerial view of estuary and the town of Lagos. The white square stands for the coring site. 1. As sondagens e o tratamento dos sedimentos Para alcançar os objetivos propostos foram realizadas diversas sondagens na planície aluvial do vale da ribeira de Bensafrim, sendo aqui utilizadas apenas duas. estas foram realizadas na área do estuário não artificializado, cerca de 1900m a montante da foz (fig. 1). a localização destas duas sondagens, ambas na margem direita, não foi arbitrária, mas Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 117 sim exatamente em frente da foz atual e original da ribeira, de forma que o enchimento sedimentar pudesse registar não só os eventos terrestres, mas também os marinhos, especialmente os mais energéticos, designadamente galgamentos oceânicos ou ondas de tsunamis, cujos efeitos na cidade de Lagos estão amplamente documentados nos documentos históricos e arqueológicos. a sondagem LagosConv, exploratória, foi realizada 220m a montante da ponte velha que liga Lagos à Meia Praia, com um trado manual Edelman, atingido uma profundidade de 1,95m e com a boca da sondagem a 37º 06’ 45,30’’ n e 08º 40’ 43,87’’ W (araújo-Gomes, 2010). a outra sondagem, a Bf3, dista 20cm da LagosConv (boca da sondagem a 37º 06’ 45,50’’ n e 08º 40’ 43,86’’ W) e foi realizada com um amostrador contínuo (Shelby, com 80mm de diâmetro interior) com linner em PVC transparente, para recolha de amostras intactas, tendo alcançado a profundidade de 5,68m (sondagem realizada por uma empresa especializada). a boca de ambas as sondagens estava a 1,29m acima do nível médio do mar (Datum de Cascais). na sondagem Bf3 e apesar de todas as precauções, nomeadamente o linner de PVC transparente capaz de proteger as paredes e impedindo a queda de sedimento, não foi possível obter a amostragem completa de todo o sedimento devido à presença de um aquicludo. Os sedimentos desta sondagem foram subamostrados com 1-2cm de intervalo, a fim de permitir determinar pequenas variações na sua composição fundamental. foram obtidas 249 amostras para análise sedimentológica. na sondagem LagosConv a amostragem dos sedimentos foi feita por unidades litoestratigráficas e, para cada uma delas, selecionadas as amostras para datação por radiocarbono (araújo-Gomes, 2010). O tratamento das amostras, após prévia destruição da matéria orgânica, seguiu o protocolo estabelecido no quadro do projeto, que incluiu a) secagem das amostras a 40ºC, em estufa, b) pesagem e tratamento com peróxido de hidrogénio, c) lavagem e secagem para determinação da percentagem de matéria orgânica, d) crivagem e atenuação de raios X (Sedigraph) para obtenção dos dados da análise dimensional. O método de crivagem foi utilizado quando a dimensão do sedimento predominante era superior a 0,063mm e o sedigraph para os restantes. este tratamento foi realizado no laboratório do departamento de Geologia da Universidade do Minho. Os parâmetros estatísticos para a caraterização elementar dos sedimentos, quer de uma, quer de outra sondagem, foram obtidos através do software seDPC (Henriques, 1998; 2003; 2004). Utilizou-se o método de folk e Ward para o cálculo da média (Mz (ϕ)), desvio padrão (σ(ϕ)), assimetria (skew (ϕ)) e kurtose (Kurt (ϕ)) para descrever a distribuição do tamanho de grão. 2. O enchimento sedimentar do estuário O testemunho da sondagem Bf3, representado na figura 2, evidencia os três sub- -conjuntos (tarolos) que foi possível amostrar e que ilustram não só a homogeneidade da quase totalidade da sondagem, mas também a diversidade do conjunto superior. Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 118 fig. 2 – resultados da análise granulométrica (da gravilha, à direita na coluna, da argila, à esquerda) e parâmetros estatísticos da distribuição do tamanho de grão para o testemunho Bf3. Ls – secção inferior; Us – secção superior. Fig. 2 – Results of size analysis (from gravel, on the right, to clay, on the left) and statistical distribution parameters of grain size for core BF3. LS – lower section; US – upper section. este testemunho Bf3 apresenta, assim, 2 secções distintas: i) a secção inferior (Ls), desde 558cm até 109cm, que comporta dois pacotes sedimentares, é monótona e dominada por areia fina, pouco calibrada, com assimetria e kurtose variáveis, onde os siltes e argilas não ultrapassam 5%, (com exceção da camada inferior); ii) a secção superior (Us), de 109cm ao topo da sondagem, divide-se em três sub-secções. a subsecção inferior, de 109cm a 81cm, é uma camada argilosa (argila sempre superior a 80%), com silte (entre 13 e 22%) e areia (abaixo de 8%). a subsecção intermédia, entre 81cm e 43cm, é arenosa (60 a 80%) com conteúdo variável em argila, silte e areias finas. a subsecção superior, acima da anterior e até ao topo da sondagem, é semelhante à subsecção inferior, com excepção de uma intrusão arenosa entre 27 e 22cm, com um máximo aos 24cm, onde a areia atinge 90%. Os sedimentos da sondagem Bf3 representam, por conseguinte, dois ambientes sedimentares distintos: Ls possui características de um ambiente mais energético, enquanto Us a de um ambiente mais protegido, interrompido episodicamente por eventos de elevada energia. a análise do core LagosConv, muito próximo de Bf3 (a cerca de 20cm), mostrou variações idênticas (araújo-Gomes, 2010; ramos-Pereira et al., 2010), pelo que nos dispensamos de também a apresentar aqui. 3. As datações de radiocarbono além da análise sedimentológica atrás referida resumidamente, as amostras do testemunho Bf3 foram também sujeitas a uma análise elementar e isotópica (Portela, 2013), cujos resultados não serão tomados em conta neste trabalho, uma vez que não representam qualquer mais-valia para a elaboração do modelo bayesiano que irá ser proposto. Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 119 além destas análises, algumas amostras de Bf3 foram também, como já atrás mencionado, datadas pelo radiocarbono (quadro i). De igual modo, algumas amostras do testemunho LagosConv foram objecto de datação pelo mesmo método, apresentando-se os resultados no quadro ii. Quadro i – Datas de radiocarbono obtidas a partir de amostras do testemunho Bf3. Table I – Radiocarbon dates from core BF3 samples. Ref. Laboratório Amostra # Prof. (cm) Espessura (cm) δ13C (‰) 14C age (BP) Beta-355636 1,14 31 2 -24,4 490 ± 30 Beta-355637 1,32 56 2 -23,3 1750 ± 30 Beta-355691 2,1 60 3 -20,9 3220 ± 30 Beta-355692 2,9 72 1 -25,0 1440 ± 30 Beta-355693 4,3 180 3 -20,7 2810 ± 30 Beta-355638 6,33 568 3 -22,3 6680 ± 40 Quadro ii – Datas de radiocarbono obtidas a partir de amostras do testemunho LagosConv. Table II – Radiocarbon dates from core LagosConv samples. Ref. Laboratório Amostra # Prof. Média (cm) Espessura da Amostra (cm) δ13C (‰) 14C age (yr BP) Beta-282768 B11-13 95 2 -22,2 2780±40 Beta-282769 B16-18 102 2 -21,3 2700±40 Beta-282770 B22-33 110 11 -18,8 2880±40 Beta-282771 B34-44 125 10 -19,3 3000±40 Beta-282772 C1-26 160 25 -20,5 3330±40 Beta-282773 C30-40 174 10 -20,5 3360±40 Beta-282774 C41-48 186 7 -19,0 3610±40 De maneira a definir as fronteiras das variações identificadas no ambiente sedimentar estuarino, bem como dos episódios de alta energia detetados, foram testadas diferentes aproximações com um modelo bayesiano, aproximações essas que se apresentam e se discutem a seguir. iii. O MODeLO BaYesianO 1. A Estatística Bayesiana. Algumas noções Com base nos dados cronológicos obtidos dos sedimentos procurou-se construir um modelo estatístico, que levasse em conta a sequência estratigráfica, fazendo uso de uma análise bayesiana, o que iria permitir a construção de modelos de deposição sedimentar, a identificação de eventuais inversões estratigráficas e outliers e, consequentemente, a criação de um enquadramento cronológico robusto. Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 120 a estatística Bayesiana, na qual a modelação de datas de 14C utilizada se baseia, é descrita em detalhe em vários artigos (Buck et al., 1991; 1992; Bronk ramsey, & allen, 1995, por exemplo), assim como os modelos matemáticos usados pelo programa OxCal, o software que se utilizou para a modelação estatística e calibração das datas convencionais de radiocarbono (Bronk ramsey, 2001; 2008; 2009). a estatística Bayesiana baseia- -se no teorema de Bayes (1763), sendo os cálculos estatísticos baseados na Markov Chain Monte Carlo (MCMC), incorporada no software OxCal utilizado. a sua aplicação à calibração de datas de radiocarbono permite a redução da incerteza associada ao processo de conversão das datas convencionais de radiocarbono em datas de calendário solar, devido à incorporação de informações relacionadas com constrangimentos variados, designadamente relações estratigráficas entre as amostras ou, mais especificamente, no nosso caso, a profundidade das amostras datadas. Uma vez que se proceda à incorporação dessas informações no modelo, obtém-se uma distribuição de incerteza posterior para cada data e, consequentemente, novos limites para o intervalo de tempo da data calibrada (Bronk ramsey, 2001; 2008, 2009). Deve ser realçado que a estatística Bayesiana pode ser aplicada a qualquer método de datação, isto é, permite incorporar informação cronológica com origem noutros métodos de datação absoluta para além do radiocarbono (OsL, tL ou datação histórica, por exemplo) na arquitetura de um modelo simples ou complexo (Bronk ramsey, 2001; 2008; 2009; rhodes et al., 2003). a modelação bayesiana é usada para construir, entre outros, modelos cronológicos baseados em sequências sedimentares. nesta modelação, os dados cronológicos são associados à sequência estratigráfica e, consequentemente, à profundidade. a tais modelos bayesianos dá-se o nome de modelos de deposição (Bronk ramsey, 2008), dos quais o modelo de deposição P_sequence é largamente utilizado (Lienkaemper & Bronk ramsey, 2009). este tipo de modelo permite levar em conta a aleatoriedade do processo de sedimentação ao longo do tempo, integrando os vários constrangimentos que afetam a sedimentação. Um fator adicional (fator k) é, normalmente, introduzido na arquitetura de um modelo P_sequence. este fator k define-se como o número de eventos por unidade de comprimento. Valores elevados de k refletem pequenas variações das taxas de sedimentação, enquanto valores baixos refletem flutuações maiores neste parâmetro. Quando a informação cronológica é reduzida, a taxa de deposição pode ser obtida usando um k elevado, o que permitirá que os índices de concordância com o modelo (individuais e global) sejam iguais ou superiores a 60%. esta percentagem encontra-se definida, pela modelação bayesiana, como o critério estatístico acima do qual se considera que as datas individuais estão em concordância com o modelo proposto. O mesmo limite (60%) é utilizado na análise global do modelo. inclui a análise estatística dos constrangimentos impostos ao modelo de deposição ou, por outras palavras, a ordem e a posição (profundidade) das datas absolutas obtidas para a sequência sedimentar, bem como a interpolação usada (k) (Bronk ramsey, 2008). a utilização de um índice de concordância individual (a) e outro de concordância global do modelo (amodel) permite verificar até que ponto o modelo conceptual definido pelo utilizador está de acordo com os dados (datações absolutas, estratigrafia), Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 121 dando uma indicação do grau de confiança do modelo. a sua utilização é de grande importância e quando não é levada em consideração provoca, em muitos casos, discordância entre os modelos cronológicos e de deposição criados e os dados experimentais obtidos (Bronk ramsey, 2012). todas as datas de 14C cuja concordância individual (a) seja inferior a 60%, isto é, que não respeitem o critério imposto, podem ser rejeitadas e um novo modelo deverá ser construído. este tipo de modelos de deposição permite determinar taxas de sedimentação (Dr) para qualquer secção da sondagem sedimentar, com um intervalo de confiança de 1σ ou de 2σ. adicionalmente, pode também ser calculada a data calibrada para qualquer profundidade, da qual não tivesse sido datada qualquer amostra. 2. O tipo de modelo proposto Optou-se por um modelo P_Sequence, uma vez que este, como referido atrás, permite a incorporação da natureza aleatória do processo sedimentar ao longo do tempo, num dado local, resultando num modelo estatístico que leva em conta uma multiplicidade de condições que possam provocar alterações sedimentares (Bronk ramsey, 2008; 2009). neste tipo de modelo bayesiano é necessária a utilização de fronteiras para estabelecer um modelo cronológico correto. Dada a ausência de informação cronológica (de radiocarbono) para o topo da sondagem, a fronteira inicial para o modelo, considerou-se que a fronteira inicial (para o modelo) corresponderia ao ano de amostragem, 2009. a fronteira final para o modelo foi definida como se situando 1cm abaixo da profundidade da última amostra datada (corresponderia, por conseguinte, a 569cm). relativamente ao processo sedimentar, ele foi condicionado por um fator adicional (k) introduzido na arquitetura do modelo P_Sequence. Utilizou-se k = 0,4, de que resulta uma interpolação de 2,5cm. este valor representa o valor máximo de k que pode ser introduzido no modelo de maneira a obter índices de concordância superiores a 60% para os dados cronológicos utilizados. além disto, a modelação bayesiana proposta seria realizada com o programa de calibração OxCal 4.2.3 (Bronk ramsey, 2008; 2013), fazendo uso da curva de calibração intCal13 (reimer et al., 2013). 3. A construção do modelo. Discussão numa primeira aproximação, considerou-se que todos os dados do testemunho Bf3 eram válidos para a sequência sedimentar, tendo sido, por isso, utilizadas todas as seis datas de radiocarbono de Bf3 (quadro i), obtidas por aMs a partir da matéria orgânica presente nas amostras sedimentares datadas. Como resultado dessa modelação, obteve-se uma distribuição nULa, com concordâncias individuais e global muito pobres (a < 60,0%), provocadas pela inversão observada na figura 3. Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 128 fig. 6 – representação gráfica final da idade dos sedimentos do core Bf3 em função da profundidade. Fig. 6 – Final plot of age versus depth for core BF3. iV. COnCLUsÕes a análise bayesiana aplicada aos dados obtidos na sequência sedimentar colhida no estuário da ribeira de Bensafrim provou ser uma ferramenta valiosa no estabelecimento de modelos de deposição, facilitando a criação de um enquadramento cronológico robusto para uma posterior reconstrução paleoambiental e paleoclimática fiável desta região costeira. a subida do nível médio do mar durante o Holocénico Médio e superior induziu o enchimento sedimentar dos estuários, com alterações no regime de sedimentação e, consequentemente, do tipo de sedimentos e das taxas de deposição. as variações sedimentares registadas no estuário em causa representam uma evidência deste facto, mostrando diferentes ambientes energéticos. a primeira fase de sedimentação (Ls), essencialmente constituída por sedimentos arenosos, ocorreu em ambiente de estuário aberto, uma ria em assoreamento como foi designada em ramos-Pereira et al. (2010) e araújo-Gomes (2010), com taxas de deposiPortela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 129 ção de c. de 1mm/ano-1, que duraram de 8000 a 3000 cal BP (6000 a 1000 cal BC). estes valores são concordantes com os obtidos para outros estuários localizados na costa do algarve. a segunda fase (Us) corresponde a uma maior variação das condições de sedimentação, com taxas de sedimentação menores do que na primeira fase. essa situação é consistente com um estuário protegido, provavelmente parcialmente fechado por um cordão ou barra arenosa, como existia ainda em meados do século XX, com diversos episódios de overwash, sendo um bem definido a 24cm. este evento poderá corresponder ao tsunami que se seguiu ao terramoto de Lisboa de 1755, se outros indicadores também apontarem nessa direcção. a ser aceitável, essa correspondência seria totalmente compatível com o modelo bayesiano que se construiu. aGraDeCiMentOs esta investigação foi financiada pela fundação para a Ciência e a tecnologia através do projecto PTDC/CTEGIX/ 104035/2008 – FMI 5000: Environmental changes: Fluvio-marine interactions over the last 5000 yrs. P. Portela agradece a BiC do mesmo projecto e J. Martins a bolsa PhD sfrH/BD/45528/2008, financiada pela mesma instituição. referÊnCias BiBLiOGrÁfiCas Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130 araújo-Gomes, J. (2010). Estuário da Ribeira de Bensafrim: Leitura Geo-arqueossismológica [ribeira de Bensafrim estuary: Geo-archaeismological reading]. (Dissertação de Mestrado). instituto de Geografia e Ordenamento do território, Universidade de Lisboa, Lisboa. Bayes, t. r. (1763). an essay towards solving a problem in the doctrine of chances. Philosophical Transactions of the Royal Society, 53, 370-418. Boski, t., Camacho, s., Moura, D., fletcher, W., Wilamowski, a., Veiga-Pires, C… santana, P. (2008). Chronology of the sedimentary processes during the postglacial sea level rise in two estuaries of the algarve coast, southern Portugal. Estuarine, Coastal and Shelf Science, 77, 230-244. Bronk ramsey, C. (2013). Oxcal 4.2.3. interface build: 78, Updated: 9/10/2013. Bronk ramsey, C. (2012). Oxcal 4.1 Manual. Web interface build number: 69. retrieved from http://c14.arch.ox.ac.uk/oxcalhelp/hlp_contents.html Bronk ramsey, C. (2001). Development of the radiocarbon calibration program OxCal. Radiocarbon, 43(2a), 355-363. Bronk ramsey, C. (2009). Bayesian analysis of radiocarbon dates. Radiocarbon, 51(1), 337-360. Bronk ramsey, C. (2008). Deposition models for chronological records. Quaternary Science Reviews, 27(1-2), 42-60. Bronk ramsey, C., & allen, M. J. (1995). analysis of the radiocarbon dates and their archaeological significance. in r. M. J. Cleal, K. e. Walker & r. Montague (eds.), Stonehenge in its landscape: twentieth century excavations (pp. 526-535). London: english Heritage Buck, C. e., Kenworthy, J. B., Litton, C. D., & smith, a. f. M. (1991). Combining archaeological and radiocarbon information – a Bayesian-approach to Calibration. Antiquity, 65(249), 808-821. Buck, C. e., Litton, C. D., & smith, a. f. M. (1992). Calibration of radiocarbon results Pertaining to related archaeological events. Journal of Archaeological Science, 19(5), 497-512. 130 Daveau, s. (1987). Comentários e actualizações [Comments and updates]. in O. ribeiro, H. Lautensach & s. Daveau (eds.), Geografia de Portugal. Vol. 1. A posição Geográfica e o Território [Geography of Portugal. Vol. 1. Geographic Position and territory] (pp. 83-117). Lisboa: edições João sá da Costa, 83-117. Delgado, J., Boski, t., nieto, J. M., Pereira, L., Moura, D., Gomes, a… García-tenorio, r. (2012). sea-level rise and anthropogenic activities recorded in the late Pleistocene/Holocene sedimentary infill of the Guadiana estuary (sW iberia). Quaternary Science Reviews, 33, 121-141. font, e., nascimento, C., Omira, r., Baptista, M. a., & silva, P. f. (2010). identification of tsunamiinduced deposits using numerical modeling and rock magnetism techniques: a study case of the 1755 Lisbon tsunami algarve, Portugal. Physics of the Earth and Planetary Interiors, 182(3-4), 187-198. font, e., Veiga-Pires, C., Pozo, M., nave, s., Costas, s., ruiz Muñoz, f… rodríguez-Vidal, J. (2013). Benchmarks and sediment source(s) of the 1755Lisbon tsunamideposit at Boca do rio estuary. Marine Geology,343(1),1-14. Henriques, r. (1998). Propostas Metodológicas para a Monitorizaçãodas Zonas Costeiras – Aspectos Sedimentológicos [Methodological Proposals for Coastal Zone Monitoring – sedimentological aspects]. (Dissertação de Mestrado). faculdade de Ciências, Universidade do Porto, Porto. Henriques, r. (2004). seDMaC/seDPC: an application to support particle size analysis of unconsolidated sediments. 32nd International Geological Congress, florence-italy, abs. 154-6, p. 726. Henriques, r. (2003). seDMaC/seDPC: programa informático de apoio à análise dimensional de populações detríticas [seDMaC / seDPC: computer program to support dimensional analysis of detrital populations]. Ciências da Terra, 5 (special issue), Vi Congresso nacional de Geologia, Lisbon, 40. Lienkaemper, J. J., & Bronk ramsey, C. (2009). OxCal: Versatile tool for Developing Paleoearthquake Chronologies – a Primer. Seismological Research Letters, 80, 431-434. Kortekaas, s., & Dawson, a. G. (2007). Distinguishing tsunami and storm deposits: an example from Martinhal, sW Portugal. Sedimentary Geology, 200, 208-221. Portela, P. J. C. (2013). Análise Elementar e Isotópica de Sedimentos Estuarinos do Rio Alcabrichel (Torres Vedras) e da Ribeira de Bensafrim (Lagos) [elemental and isotopic analysis of estuarine sediments of the alcabrichel river (torres Vedras) and the ribeira de Bensafrim (Lagos)]. (Dissertação de Mestrado). faculdade de Ciências, Universidade de Lisboa, Lisboa. ramos-Pereira, a., araújo-Gomes, J., trindade, J., & ramos, C. (2010). Da ria à marina de Lagos. a evolução de um ambiente estuarino no Holocénico superior [from the ria to the Marina of Lagos. The evolution of an upper Holocene estuarine environment]. Actas do V Congresso Nacional de Geomorfologia, associação Portuguesa de Geomorfólogos, Porto, 49-52. ramos-Pereira, a., ramos, C, trindade, J., araújo- -Gomes, J., rocha, J., Granja, H… Martins, J. M. (2011). fMi 5000: um projecto sobre mudanças ambientais holocénicas [fMi 5000: a project on holocene environmental changes]. Finisterra – Revista Portuguesa de Geografia, XLVI(91), 99-106. reimer, P. J., Bard, e., Bayliss, a., Beck, J. W., Blackwell, P. G., Bronk ramsey, C… van der Plicht, J. (2013). intCal13 and Marine13 radiocarbon age Calibration Curves 0-50,000 Years cal BP. Radiocarbon, 55(4), 1869–1887. rhodes, e. J., Bronk ramsey, C., Outram, Z., Batt, C., Willis, L… Bond, J. (2003). Bayesian method applied to the interpretation of multiple OsL dates: high precision sediment ages from Old scatness Broch excavations, shetland isles. Quaternary Science Reviews, 22, 1231-1244. Portela, P. J. C., Martins, J. M. M., Monge soares, a. M., ramos-Pereira, a., araújo-Gomes, J., torres, a. Finisterra, LIV(112), 2019, pp. 113-130