Migração com amplitude verdadeira em meios com gradiente constante de velocidade

Este trabalho tem por objetivo a aplicação de um método de migração com amplitudes verdadeiras, considerando-se um meio acústico onde a velocidade de propagação varia linearmente com a profundidade. O método de migração é baseado na teoria dos raios e na integral de migração de Kirchhoff, procurando posicionar de forma correta os refletores e recuperar os respetivos coeficientes de reflexão. No processo de recuperação dos coeficientes de reflexão, busca-se corrigir o fator de espalhamento geométrico de reflexões sísmicas primárias, sem o conhecimento a priori dos refletores procurados. Ao considerar-se configurações fonte-receptor arbitrárias, as reflexões primárias podem ser imageadas no tempo ou profundidade, sendo as amplitudes do campo de ondas migrado uma medida dos coeficientes de reflexão (função do ângulo de incidência). Anteriormente têm sido propostos alguns algoritmos baseados na aproximação de Born ou Kirchhoff. Todos são dados em forma de um operador integral de empilhamento de difrações, que são aplicados à entrada dos dados sísmicos. O resultado é uma seção sísmica migrada, onde cada ponto de reflexão é imageado com uma amplitude proporcional ao coeficiente de reflexão no ponto. No presente caso, o processo de migração faz uso de um modelo com velocidade que apresenta uma distribuição que varia linearmente com a profundidade, conhecido também como gradiente constante de velocidade. O esquema de migração corresponde a uma versão modificada da migração de empilhamento por difração e faz uso explícito da teoria do raio, por exemplo, na descrição de tempos de trânsito e amplitudes das reflexões primárias, com as quais a operação de empilhamento e suas propriedades podem ser entendidas geometricamente. Efeitos como o espalhamento geométrico devido à trajetória do raio levam a distorção das amplitudes. Estes efeitos têm que ser corregidos durante o processamento dos dados sísmicos. Baseados na integral de migração de Kirchhoff e na teoria paraxial dos raios, foi derivada a função peso e o operador da integral por empilhamento de difrações para um modelo sísmico 2,5-D, e aplicado a uma serie de dados sintéticos em ambientes com ruído e livre de ruído. O resultado mostra a precisão e estabilidade do método de migração em um meio 2,5-D como ferramenta para obter informação sobre as propriedades de refletividade da subsuperfície da terra. Neste método não são levados em consideração a existência de caústicas nem a atenuação devido a fricção interna.

Migração em profundidade usando a solução numérica da equação da eiconal

Nos últimos anos tem-se verificado um interesse crescente no desenvolvimento de algoritmos de imageamento sísmico com a finalidade de obter uma imagem da subsuperfície da terra. A migração pelo método de Kirchhoff, por exemplo, é um método de imageamento muito eficiente empregado na busca da localização de refletores na subsuperficie, quando dispomos do cálculo dos tempos de trânsito necessários para a etapa de empilhamento, sendo estes obtidos neste trabalho através da solução da equação eiconal. Primeiramente, é apresentada a teoria da migração de Kirchhoff em profundidade baseada na teoria do raio, sendo em seguida introduzida a equação eiconal, através da qual são obtidos os tempos de trânsitos empregados no empilhamento das curvas de difrações. Em seguida é desenvolvido um algoritmo de migração em profundidade fazendo uso dos tempos de trânsito obtidos através da equação eiconal. Finalmente, aplicamos este algoritmo a dados sintéticos contendo ruído aditivo e múltiplas e obtemos como resultado uma seção sísmica na profundidade. Através dos experimentos feitos neste trabalho observou-se que o algoritmo de migração desenvolvido mostrou-se bastante eficiente e eficaz na reconstrução da imagem dos refletores.

Estudo da influência do lineamento Tibagi e Sinclinal de Torres na estruturação da porção norte da bacia de Pelotas

The area of this work is located along the Northern portion of the Pelotas´s basin. The onshore region extends from the Southern portion São Paulo State to the North of Rio Grande do Sul State, with a total area of 440 530 km ². The objective of the research is to characterize the main geological structures of the outcropping basement in the Paraná Basin and its correlation with the offshore portion of the Pelotas Basin. The characterization was performed through the integration of surface (SRTM images and outcrops) and subsurface data (seismic data and well logs), the method includes the use of remote sensing techniques and seismic interpretation. The main structures to be analyzed are the Tibagi Lineament and the “Torres Syncline in order to verify its influence on the geological evolution and deformation of the basin. As result, it was identified the continuity of the both structures inside the offshore basin (Pelotas Basin); but only the Tibagi Lineament represents an important structural feature in the offshore basin. In this sense, the Tibagi lineament shows a NW-SE direction along the both basins (Pelotas and Paraná) and it is defined by drainages and relief on the continental portion. It was classified as a graben structure limited by normal faults and it also was verified a trend of faults inside the lineament that cuts all the sedimentary package of Itararé Group from Paraná Basin. The Torres syncline is a structure that connects the arches of Ponta Grossa and Rio Grande, showing NW-SE orientation and a low angle dip in the NW direction. Its projection into the Pelotas Basin was not recognized by the methods applied in this work, because a preliminary analyze indicates that its projection corresponds to a structural high inside the Pelotas Basin

Caracterização dos reservatórios do grupo Macaé no Campo de Garoupa - Bacia de Campos - por meio de interpretação sísmica e perfis de poço

The study area is the Garoupa Field, located 75 km from the coast of Rio de Janeiro in bathymetric elevation of 125 m and it is the first important offshore discovery of the Campos Basin with hydrocarbon production in the Albian carbonates. The high significance of carbonate rock reservoirs for the global oil industry and its degree of complexity, combined with the fact that the Garoupa Field is a mature field, make this research project interesting to develop a re-exploratory vision of the study area. The objective of this paper is to characterize the physical properties of the Macae Group reservoirs in the Garoupa Field through well correlation techniques, petrophysical analysis and seismic interpretation (3D), in order to evaluate qualitatively and quantitatively the response of geophysical data to the presence of carbonate rock reservoirs. From the analysis of well profiles (gamma ray, resistivity, density and sonic) it was possible to identify six reservoir levels in the Quissamã Formation. As a result, the reservoirs have good correlation between the wells in the Garoupa Field and show lateral facies variation. In the seismic interpretation, structural contour maps were generated in time (ms) from the top of the Macae Group and two horizons that correspond to the reservoir levels. These maps show that the oil producing area of this field is a structural high located at the west of the interpreted seismic cube, and it was classified as a rollover limited by lateral faults forming a horst. The seismic attribute maps show a good response to the reservoirs of the Garoupa Field because they are associated with amplitude anomalies, verified by the correlation between the physical properties of the reservoir rock and seismic data through cross plots, with emphasis in a linear correlation between the resistivity profile and the values of Maximum Absolute Amplitude and RMS Amplitude attributes

Geologia das coquinas-reservatório do Campo de Pampo, Bacia de Campos

The Pampo oil field is located in the southwest of the trend Badejo-Linguado-Pampo in an accumulation of single - mixed trap: structural-stratigraphic-diagenetic. Its main reservoir is a coquina shell of bivalves (the lowest) in the Lagoa Feia Group. During the rift phase, the Badejo-Linguado-Pampo trend´s accumulations evolve according to three tectono-stratigraphic cycles. The first two cycles are formed by siliciclastic rocks with fining up sequences and carbonates coquinas bivalves. The youngest cycle related to Alagoas age is a transgressive event represented by the presence of an evaporitic layer in the top (anhydrite). This study aims to characterize the reservoir Coqueiros Formation based on the analysis of 2D and 3D seismic data and well data-correlation profiles. The structural map of the top of coquinas reservoir indicates a curvilinear contour of Pampo Fault as described on the literature. This fault was interpreted on seismic data as a basement´s high, and it doesn´t show influence on the horizons above the top of Lagoa Feia Group. The Pampo fault is responsible for the division of the field into two blocks: the hanginwall in the West and the footwall to the East. This division is well marked on the reservoir´s isopach map where a greater thickness of reservoir is observed on the lower block. In the Southeast extreme of Badejo-Linguado-Pampo trend, on Pampo Field, the thick siliciclastic´s interval ends laterally to the basement, and its lower´s cycle forms a wedge, as consequence the carbonate-coquina overlaps directly the basement. Another implication of the higher and distal position of Pampo field is that the third cycle is absent, truncated by the unconformity pre-Macaé Group (Albian)

Modelo petrofísico 3D do Arenito Namorado no Campo de Namorado - Bacia de Campos, RJ

The Namorado Oil Field represents the beginning of the oil exploration in Brazil, in the 70s, and it is still a subject of researches because the importance of this turbidite sandstone in the brazilian oil production. The Namorado’s production level was denominated “Namorado sandstone”, it is composed by turbidite sandstone deposited during the Albian-Cenomanian. In order to define the structural geometry of the main reservoir, geological and geophysical tools like RECON and Geographix (Prizm – Seisvision) softwares were used, and its application was focused on geological facies analysis, for that propose well logs, seismic interpretation and petrophysical calculations were applied. Along this work 15 vertical wells were used and the facies reservoirs were mapped of along the oil field; it is important to mentioned that the all the facies were calibrated by the correlation rock vs log profile, and 12 reservoir-levels (NA-1, NA-2, NA-3, NA-4, NA-5, NA-6, NA-7, NA-8, NA-9, NA-10, NA-11 e NA-12) were recognized and interpreted. Stratigraphic sections (NE-SW and NW-SE) were also built based on stratigraphic well correlation of each interpreted level, and seismic interpretation (pseudo-3D seismic data) on the southeastern portion of the oil field. As results it was interpreted on two- and three-dimensional maps that the deposition reservoir’s levels are hight controlled by normal faults systems. This research also shows attribute maps interpretation and its relationship with the selection of the reservoir attribute represented on it. Finally the data integration of stratigraphic, geophysical and petrophysical calculations lets us the possibility of obtain a detail geological/petrophysical 3D model of the main reservoir levels of “Namorado sandstone” inside the oil/gás field

Caracterização da distribuição física dos dolomitos do membro Búzios no Campo de Badejo - Bacia de Campos

In the region of Badejo oil Ffeld (Campos Basin) lies an Lower Albian dolomitic layer that shows reservoir rock and seal conditions, among others and, it is occasionally karstified. This work applies geostatistical techniques of ordinary and indicatior kriging, in an attempt to determine a method that assists the analysis of different scenarios offered for petroleum drilling from a quantitative point of view this fact is justified, because there are different visions and strategies that would be adopted under different dolomite conditions (stable - low porosity and high density; reservoir - high porosity and medium density; instable - high porosity and low dentisty). The main objective is a methodological introduction that has not been tested in dolomites before aiming to characterize the distribution of the three conditions described above by ordinary and indicatior kriging, which was based on the obtained data from the dolomite layer identification through well logs interpretation and correlation, followed by seismic interpretation. In addition, it was generated structural contour maps, based on 2D and 3D seismic data interpretation, and then, seismic attributes maps were calculated, in order to transform them into pseudo-density maps, i.e., maps that correlate the density values with the attribute values. As primary results, structural contour maps and seismic attributes were obtained and ordinary and indicatior kriging maps were done, on which it is possible to interpret the distribution of the main reservoir and risk probability for drilling of exploration wells obtaining trends N35E and N10W direction for areas of stability or reservoir levels, while the central part of the map presents a higher risk for loss of drilling fluid. The cut-off values levels were based on the values of first and third quartiles of cumulative histogram (instable and stable zones, respectively), as well as the reservoir level was set as the interval...

Tuzla City (BiH): an example of geohazard induced by salt extraction

Salt deposits characterize the subsurface of Tuzla (BiH) and made it famous since the ancient times. Archeological discoveries demonstrate the presence of a Neolithic pile-dwelling settlement related to the existence of saltwater springs that contributed to make the most of the area a swampy ground. Since the Roman times, the town is reported as “the City of Salt deposits and Springs”; "tuz" is the Turkish word for salt, as the Ottomans renamed the settlement in the 15th century following their conquest of the medieval Bosnia (Donia and Fine, 1994). Natural brine springs were located everywhere and salt has been evaporated by means of hot charcoals since pre-Roman times. The ancient use of salt was just a small exploitation compared to the massive salt production carried out during the 20th century by means of classical mine methodologies and especially wild brine pumping. In the past salt extraction was practised tapping natural brine springs, while the modern technique consists in about 100 boreholes with pumps tapped to the natural underground brine runs, at an average depth of 400-500 m. The mining operation changed the hydrogeological conditions enabling the downward flow of fresh water causing additional salt dissolution. This process induced severe ground subsidence during the last 60 years reaching up to 10 meters of sinking in the most affected area. Stress and strain of the overlying rocks induced the formation of numerous fractures over a conspicuous area (3 Km2). Consequently serious damages occurred to buildings and infrastructures such as water supply system, sewage networks and power lines. Downtown urban life was compromised by the destruction of more than 2000 buildings that collapsed or needed to be demolished causing the resettlement of about 15000 inhabitants (Tatić, 1979). Recently salt extraction activities have been strongly reduced, but the underground water system is returning to his natural conditions, threatening the flooding of the most collapsed area. During the last 60 years local government developed a monitoring system of the phenomenon, collecting several data about geodetic measurements, amount of brine pumped, piezometry, lithostratigraphy, extension of the salt body and geotechnical parameters. A database was created within a scientific cooperation between the municipality of Tuzla and the city of Rotterdam (D.O.O. Mining Institute Tuzla, 2000). The scientific investigation presented in this dissertation has been financially supported by a cooperation project between the Municipality of Tuzla, The University of Bologna (CIRSA) and the Province of Ravenna. The University of Tuzla (RGGF) gave an important scientific support in particular about the geological and hydrogeological features. Subsidence damage resulting from evaporite dissolution generates substantial losses throughout the world, but the causes are only well understood in a few areas (Gutierrez et al., 2008). The subject of this study is the collapsing phenomenon occurring in Tuzla area with the aim to identify and quantify the several factors involved in the system and their correlations. Tuzla subsidence phenomenon can be defined as geohazard, which represents the consequence of an adverse combination of geological processes and ground conditions precipitated by human activity with the potential to cause harm (Rosenbaum and Culshaw, 2003). Where an hazard induces a risk to a vulnerable element, a risk management process is required. The single factors involved in the subsidence of Tuzla can be considered as hazards. The final objective of this dissertation represents a preliminary risk assessment procedure and guidelines, developed in order to quantify the buildings vulnerability in relation to the overall geohazard that affect the town. The historical available database, never fully processed, have been analyzed by means of geographic information systems and mathematical interpolators (PART I). Modern geomatic applications have been implemented to deeply investigate the most relevant hazards (PART II). In order to monitor and quantify the actual subsidence rates, geodetic GPS technologies have been implemented and 4 survey campaigns have been carried out once a year. Subsidence related fractures system has been identified by means of field surveys and mathematical interpretations of the sinking surface, called curvature analysis. The comparison of mapped and predicted fractures leaded to a better comprehension of the problem. Results confirmed the reliability of fractures identification using curvature analysis applied to sinking data instead of topographic or seismic data. Urban changes evolution has been reconstructed analyzing topographic maps and satellite imageries, identifying the most damaged areas. This part of the investigation was very important for the quantification of buildings vulnerability.

Geological and Structural evolution of the Eurasia Africa plate boundary in the Gulf of Cadiz Central Eastern Atlantic Sea.

Iberia Africa plate boundary, cross, roughly W-E, connecting the eastern Atlantic Ocean from Azores triple junction to the Continental margin of Morocco. Relative movement between the two plate change along the boundary, from transtensive near the Azores archipelago, through trascurrent movement in the middle at the Gloria Fracture Zone, to transpressive in the Gulf of Cadiz area. This study presents the results of geophysical and geological analysis on the plate boundary area offshore Gibraltar. The main topic is to clarify the geodynamic evolution of this area from Oligocene to Quaternary. Recent studies have shown that the new plate boundary is represented by a 600 km long set of aligned, dextral trascurrent faults (the SWIM lineaments) connecting the Gloria fault to the Riff orogene. The western termination of these lineaments crosscuts the Gibraltar accretionary prism and seems to reach the Moroccan continental shelf. In the past two years newly acquired bathymetric data collected in the Moroccan offshore permit to enlighten the present position of the eastern portion of the plate boundary, previously thought to be a diffuse plate boundary. The plate boundary evolution, from the onset of compression in the Oligocene to the Late Pliocene activation of trascurrent structures, is not yet well constrained. The review of available seismics lines, gravity and bathymetric data, together with the analysis of new acquired bathymetric and high resolution seismic data offshore Morocco, allows to understand how the deformation acted at lithospheric scale under the compressive regime. Lithospheric folding in the area is suggested, and a new conceptual model is proposed for the propagation of the deformation acting in the brittle crust during this process. Our results show that lithospheric folding, both in oceanic and thinned continental crust, produced large wavelength synclines bounded by short wavelength, top thrust, anticlines. Two of these anticlines are located in the Gulf of Cadiz, and are represented by the Gorringe Ridge and Coral Patch seamounts. Lithospheric folding probably interacted with the Monchique – Madeira hotspot during the 72 Ma to Recent, NNE – SSW transit. Plume related volcanism is for the first time described on top of the Coral Patch seamount, where nine volcanoes are found by means of bathymetric data. 40Ar-39Ar age of 31.4±1.98 Ma are measured from one rock sample of one of these volcanoes. Analysis on biogenic samples show how the Coral Patch act as a starved offshore seamount since the Chattian. We proposed that compression stress formed lithospheric scale structures playing as a reserved lane for the upwelling of mantle material during the hotspot transit. The interaction between lithospheric folding and the hotspot emplacement can be also responsible for the irregularly spacing, and anomalous alignments, of individual islands and seamounts belonging to the Monchique - Madeira hotspot.

Spatial and Temporal Evolution of the Volcanics and Sediments of the Kenya Rift

New volumetric and mass flux estimates have been calculated for the Kenya Rift. Spatial and temporal histories for volcanic eruptions, lacustrine deposition, and hominin fossil sites are presented, aided by the compilation of a new digital geologic map. Distribution of volcanism over time indicates several periods of southward expansion followed by relative positional stasis. Volcanism occurs throughout the activated rift length, with no obvious abandonment as the rift system migrated. The main exception is a period of volcanic concentration around 10 Ma, when activity was constrained within 2° of the equator. Volumes derived from seismic data indicate a total volume of c. 310,000 km3 (2.47 x 1010 kg/yr ), which is significantly more than the map-derived volumes found here or published previously. Map-based estimates are likely affected by a bias against recognizing small volume events in the older record. Such events are, however, the main driver of erupted volume over the last 5 Ma. A technique developed here to counter this bias results in convergence of the two volume estimation techniques. Relative erupted composition over time is variable. Overall, the erupted material has a mafic to silicic ratio of 0.9:1. Basalts are distinctly more common in the Turkana region, which previously experienced Mesozoic rifting. Despite the near equal ratio of mafic to silicic products, the Kenya Rift otherwise fits the definition of a SLIP. It is proposed that the compositions would better fit the published definition if the Turkana region was not twice-rifted. Lacustrine sedimentation post-dates initial volcanism by about 5 million years, and follows the same volcanic trends, showing south and eastward migration over time. This sedimentation delay is likely related to timing of fault displacements. Evidence of hominin habitation is distinctly abundant in the northern and southern sections of the Kenya Rift, but there is an observed gap in the equatorial rift between 4 and 0.5 million years ago. After 0.5 Ma, sites appear to progress towards the equator. The pattern and timing of hominid site distributions suggests that the equatorial gap in habitation may be the result of active volcanic avoidance.

Possible earthquake trigger for 6th century mass wasting deposit at Lake Ohrid (Macedonia/Albania)

Abstract. Lake Ohrid shared by the Republics of Albania and Macedonia is formed by a tectonically active graben within the south Balkans and suggested to be the oldest lake in Europe. Several studies have shown that the lake provides a valuable record of climatic and environmental changes and a distal tephrostratigraphic record of volcanic eruptions from Italy. Fault structures identified in seismic data demonstrate that sediments have also the potential to record tectonic activity in the region. Here, we provide an example of linking seismic and sedimentological information with tectonic activity and historical documents. Historical documents indicate that a major earthquake destroyed the city of Lychnidus (today: city of Ohrid) in the early 6th century AD. Multichannel seismic profiles, parametric sediment echosounder profiles, and a 10.08m long sediment record from the western part of the lake indicate a 2m thick mass wasting deposit, which is tentatively correlated with this earthquake. The mass wasting deposit is chronologically well constrained, as it directly overlays the AD472/AD 512 tephra. Moreover, radiocarbon dates and cross correlation with other sediment sequences with similar geochemical characteristics of the Holocene indicate that the mass wasting event took place prior to the onset of the Medieval Warm Period, and is attributed it to one of the known earthquakes in the region in the early 6th century AD.

Reconstructed glaciation limits in Eurasia (digitized from map Fig.1)

Based on field investigations in northern Russia and interpretation of offshore seismic data, we have made a preliminary reconstruction of the maximum ice-sheet extent in the Barents and Kara Sea region during the Early/Middle Weichselian and the Late Weichselian. Our investigations indicate that the Barents and Kara ice sheets attained their maximum Weichselian positions in northern Russia prior to 50 000 yr BP, whereas the northeastern flank of the Scandinavian Ice Sheet advanced to a maximum position shortly after 17000 calendar years ago. During the Late Weichselian (25 000-10000 yr BP), much of the Russian Arctic remained ice-free. According to our reconstruction, the extent of the ice sheets in the Barents and Kara Sea region during the Late Weichselian glacial maximum was less than half that of the maximum model which, up to now, has been widely used as a boundary condition for testing and refining General Circulation Models (GCMs). Preliminary numerical-modelling experiments predict Late Weichselian ice sheets which are larger than the ice extent implied for the Kara Sea region from dated geological evidence, suggesting very low precipitation.

Strike-slip faults mediate the rise of crustal-derived fluids and mud volcanism in the deep sea

We report on newly discovered mud volcanoes located at about 4500 m water depth 90 km west of the deformation front of the accretionary wedge of the Gulf of Cadiz, and thus outside of their typical geotectonic environment. Seismic data suggest that fluid flow is mediated by a >400-km-long strike-slip fault marking the transcurrent plate boundary between Africa and Eurasia. Geochemical data (Cl, B, Sr, 87Sr/86Sr, Delta18O, DeltaD) reveal that fluids originate in oceanic crust older than 140 Ma. On their rise to the surface, these fluids receive strong geochemical signals from recrystallization of Upper Jurassic carbonates and clay-mineral dehydration in younger terrigeneous units. At present, reports of mud volcanoes in similar deep-sea settings are rare, but given that the large area of transform-type plate boundaries has been barely investigated, such pathways of fluid discharge may provide an important, yet unappreciated link between the deeply buried oceanic crust and the deep ocean.

Physical properties of 4 cores from the Bengal Fan

New stratigraphic and high-resolution seismic data from the Bengal Fan indicate that the world's largest fan shows active growth during the most recent sea-level rise and the recent highstand. This unique phenomenon contradicts common sequence-stratigraphic models, and the sediment preserved provides new insight into the sedimentological response of a fan system to sea-level rise, climatic terminations, and monsoon intensity during the past climatic cycle. We present a detailed dated sequence of turbidite sedimentation based on a core transect perpendicular to the active channel-levee system in the upper mid-fan area. Between the two major terminations 1a (12 800 14C yr B.P.) and 1b (9700 14C yr B.P.), and especially at the end of the Younger Dryas, a 13-km-wide channel built up levees 50 m high. With decreasing sediment supply, continued sea-level rise, and increasing monsoon intensity during the early Holocene, turbidity currents were confined to the channel and gradually filled it. The canyon "Swatch of No Ground," a shelf depocenter that serves as the source for frequent turbidity currents, and the channel-levee system provide the unique opportunity for studying an active highstand system. Many fans showed this behavior only during lowered sea-level.

Turbidite sedimentology and history at DSDP Site 89-585

At Site 585 in the East Mariana Basin, a 900-m section of Aptian-Albian to Recent sediments was recovered. The upper 590 m are pelagic components (carbonate, siliceous, and clay); small-scale graded sequences and laminations are common. The underlying sediments are volcaniclastic sandstones with a large proportion of shallow-water carbonate debris; sedimentary structures including complete Bouma sequences, cross-laminae, and scouring are common. These structures indicate that the entire section was deposited by turbidity currents. The change in lithology upward in the section reflects the evolution of the surrounding seamounts, from their growth stages during the middle of the Cretaceous to the later subsidence phases. Several black layers containing pyritized organic debris and associated turbidite structures were cored near the Cenomanian/Turonian boundary; this material has been transported from the flanks of the seamounts where it was deposited within a shallow anoxic zone. Seismic data extends the stratigraphy across the entire Basin, showing the reflectors onlapping the seamounts, and indicating at least 1200 m of sediment at Site 585. The crust is placed at 6900 m after correcting for sediment loading, and the subsidence curve indicates that the Basin has been deeper than 5500 m since before the Aptian.