Quantitative sequence stratigraphy applied to the Barremian–lower Aptian Urgonian carbonates deposited in Provence (southeastern France) Open Access

Dynamic carbonate reservoir models are traditionally based on static geological frameworks derived from sequence stratigraphic well correlations, sedimentological interpretations, and an understanding of reservoir heterogeneities. Constructing carbonate stratigraphic architectures is a key step in carbonate sequence stratigraphy, especially for subsurface and reservoir modeling. Since the 1980s, sequence stratigraphy has been extensively used in both academia and the petroleum industry. It has proven invaluable for analyzing shallow-water carbonate sedimentary systems (e.g., Sarg, 1988; Handford and Loucks, 1993; Schlager, 1993, 2005; Pomar and Ward, 1995; Eberli et al., 2004; Droste, 2010; Catuneanu et al., 2011; Maurer et al., 2013; Pomar and Haq, 2016). This approach facilitates stratigraphic correlations across sedimentary basins and helps estimate influential factors such as sea-level changes and subsidence patterns.

However, recent research by Borgomano et al. (2020) highlights limitations in classical sequence stratigraphic models. To address these challenges, a quantitative carbonate sequence stratigraphic approach was developed. This method refines stratigraphic correlations and reduces uncertainties. It distinguishes between apparent accommodation and real accommodation. Real accommodation reflects the available sedimentary space over time, influenced by sea-level changes and subsidence relative to a fixed datum (Homewood et al., 2000). Apparent accommodation, in contrast, is inferred from measurable geological parameters, including facies bathymetry, bed thickness, and stratigraphic discontinuities. These measurements provide indirect insights into accommodation variations preserved in sedimentary records (Borgomano et al., 2020).

However, recent research by Borgomano et al. (2020) has underscored the limitations of classical sequence stratigraphic models, prompting the development of a new quantitative carbonate sequence stratigraphic approach to refine stratigraphic correlations and reduce associated uncertainties. It distinguishes between apparent accommodation and real accommodation. Real accommodation reflects the available space for sedimentation over time due to sea-level and subsidence changes relative to a fixed datum (Homewood et al., 2000). But neither sea level nor the subsidence is directly preserved as time-continuous and measurable variables in stratigraphic records. Apparent accommodation, in contrast, is inferred from measurable geological parameters, including facies bathymetry, bed thickness, and stratigraphic discontinuities. These measurements provide indirect insights into accommodation variations preserved in sedimentary records (Borgomano et al., 2020).

The new method relies on calculating and correlating apparent accommodation trends across sedimentary sections to define stratigraphic architectures. These interpretations can be tested and refined through iterative forward modeling. Its utility is demonstrated by applying it to the well-documented Barremian to lower Aptian (Lower Cretaceous) Urgonian carbonates of Provence. This case study refines our understanding of stratigraphic architecture in these carbonates.

These carbonates are well exposed in southern France (Figure 1) and have been studied in detail during the last decades. Masse (1976) developed the first large-scale stratigraphic model, which was later updated (Masse, 1993). More recent work identified sedimentological and stratigraphical evidence of at least four Barremian to the lower Aptian drowning surfaces (Masse and Fenerci-Masse, 2011). Few existing high-resolution stratigraphic models mainly focus on horizontal and vertical dimensions of rudist beds and their relationships with coral bodies (e.g., Borgomano et al., 2002; Fenerci-Masse et al., 2005). Léonide et al. (2012) provided a first high-resolution stratigraphic architecture highlighting the spatial distribution of rudist-corals facies and adjacent outer-platform bioclastic sediments. Finally, Tendil et al. (2018) published a recent updated stratigraphic model based on platform to basin correlations of ammonite-calibrated timelines (Frau et al., 2018). This study provides a well-constrained chronostratigraphic framework for the definition of high-resolution stratigraphic architectures of the Barremian-lower Aptian Urgonian carbonates of Provence.

The current western Provence domain in southern France covers approximately 9000 km². It is bounded by the Durance River, the Sainte-Victoire structure, and the Arc syncline to the east; by the Mediterranean Sea to the south; by the Carpentras Basin and the Crau Plain to the west; and, finally, by the Ventoux-Lure structure to the north (Figure 1). The study area comprises, from north to south by the Plateau d’Albion (altitude ∼850 m) and a series of east–west to northeast–southwest-trending folds and thrusts (La Nerthe, La Fare, Lubéron, and the Mont-Ventoux Lure structures). These are interrupted by north-northeast-trending strike-slip fault systems including the Nîmes, Salon-Cavaillon, Aix-en-Provence, and Mid-Durance faults (Tempier and Durand, 1981; Espurt et al., 2012).

During the Barremian to early Aptian, the study area was located at a paleolatitude of approximately 28°N. It constituted the southern margin of the Vocontian Basin, where shallow-water carbonate platforms developed (Masse et al., 2000; Masse and Fenerci-Masse, 2013). This period is considered part of the Cretaceous greenhouse, during which sediments recorded significant fluctuations in the δ¹³C curve, reflecting crises in primary productivity (Weissert and Erba, 2004). The mid-Barremian event, the Taxy event, and the early Aptian oceanic anoxic event (OAE1a) were reported to affect the carbonate production and sedimentation of the study area from the Barremian to early Aptian (Mutterlose et al., 2009; Huck et al., 2011, 2013; Tendil et al., 2018).

The Barremian to lower Aptian Urgonian carbonates in Provence were previously studied by d’Orbigny (1850), Leenhardt (1883), Masse (1976), Léonide et al. (2012), and Tendil et al. (2018). The stratigraphic units defined by the above studies depend on the geographic location.

In southern Provence, the Urgonian carbonates are divided into three lithostratigraphic units (Figure 2). The Urgonian member is a 400- to 500-m thick limestone formation dated from the lower to upper Barremian (Masse, 1976; Frau et al., 2018), characterized by massive rudist-bearing deposits. The overlying lithostratigraphic unit labeled the calcareous member by Moullade et al. (2000) is made up of alternations of upper Barremian to lower Aptian marls and limestones. Finally, the uppermost marly-calcareous member dated from the lower Aptian (Frau et al., 2015, 2018) is considered time equivalent to the OAE1a (Tendil et al., 2019).

In northern Provence, three Urgonian units were also defined (Leenhardt, 1883; Masse, 1976; Léonide et al., 2012; Tendil et al., 2018). The U1 unit is a 300- to 500-m-thick lower Barremian to basal upper Barremian unit, composed of cherty and bioclastic limestones. The U2 unit consists of a 50- to 250 m-thick rudist limestone formation dated to the upper Barremian. These last two units are time equivalent to the Urgonian and calcareous member in southern Provence. Finally, the U3 unit, dated from the lower Aptian, is faciologically equivalent to U1 and time-equivalent to the upper calcareous member in southern Provence.

The Urgonian sedimentological model has been the topic of many studies (Masse, 1976; Fenerci-Masse et al., 2005; Léonide et al., 2012; Masse and Fenerci-Masse, 2013). Fourteen facies associations have been defined and updated using macro- and micro-observations gathered from field, cores, and thin sections (Michel et al., 2023). This paper does not aim to provide a detailed description of each facies. We encourage the reader to refer to the above-cited papers for this purpose. However, all sedimentological criteria of the facies are described in detail in the supplementary material available as AAPG Datashare 202 at www.aapg.org/datashare. The facies model includes five main environments: inner platform, platform margin, outer platform, slope, and basin (Figure 3). According to Masse et al. (2011, 2013) and Tendil et al. (2018), the carbonate platform transitioned from a ramp to an extensive flat platform, spanning 100–140 km in length. By the late Barremian, it evolved into an isolated platform, restricted to northern Provence, with a reduced length of 30–60 km.

Beach-fenestrate grainstone (FA1), miliolid rich wackestone to grainstone (FA2), and rudist facies (FA3) are interpreted to deposit in the inner platform. The platform margin is characterized by coarse rudist and coral-bearing limestones (FA4) and kilometers-wide patches of coral wackestones to grainstones (FA5) (Masse, 1976; Masse et al., 2011; Léonide et al., 2012). Oolitic (FA6) and heterozoan-rich bioclastic (FA7) facies associations deposited in more open and higher-energy environments in with sediments transportation by wave currents and storm is a dominant process. The slope environment was characterized by reduced hydrodynamic energy and low light levels. It favors the deposition of fine to very fine bioclastic packstones to grainstones, grouped into FA8 and FA9. These facies contain abundant annelids, well-rounded echinoderm fragments, bryozoan debris derived from the outer platform, and sponge spicules.

The basin environment encompasses hemipelagic and pelagic facies. Hemipelagic deposits (FA10) are associated with fine-grained micritic limestones, whereas pelagic deposits (FA11) consist predominantly of marl-dominated lithologies (Masse, 1976; Léonide et al., 2012; Tendil et al., 2018). Finally, heterotrophic conditions are represented by FA13 and FA14. These correspond to gastropod-rich and Palorbitolina-rich wackestones to packstones, deposited on the platform and slope, respectively.

The aim of this study is to test an alternative method for defining the stratigraphic architecture of carbonate sedimentary systems. The workflow is based on the spatial and temporal evolution of the facies water depth and quantified apparent accommodation changes combined with the classical analysis of stratigraphic surfaces.

Apparent accommodation is distinguished from real accommodation (Borgomano et al., 2020). The real accommodation represents in time the accommodation available created or destroyed from the evolution of both sea level (eustasy) and subsidence relative to a permanent datum (Borgomano et al., 2020). The apparent accommodation corresponds in depth to the preserved thickness of sediments added to the paleobathymetry for a given stratigraphic interval. This approach is rooted in the simple fact that neither sea level nor the subsidence is directly measurable in the stratigraphic record as a time variable. Facies paleobathymetry, bed thickness, porosity, and stratigraphic discontinuities are, on the contrary, directly measurable in the stratigraphic records and can thus be analyzed as indirect indications of real accommodation variations (Figure 4).

Consequently, values of apparent accommodation represent the evolution of the available space for sedimentation compared to the initial point, here referring to the base of a section. The finer the sedimentary descriptions are, the finer the vertical resolution of apparent accommodation is. Once the bathymetric and apparent accommodation curves have been built, increasing or decreasing trends emerge (Figure 4). A marker can then be placed at the interface between two different trends. Markers can be ranked according to the magnitude of the apparent changes.

In classical sequence stratigraphic practices, deposition sequences are identified from the vertical stacking of sediments and correlated laterally between the different sections in a sedimentary basin. These sequences define two main fundamental criteria, which are the sediment supply and real accommodation, commonly represented by triangles. These triangles can have several meanings, according to the authors; however, in theory, they represent regressive and transgressive trends when the accommodation/sediment supply ratio (AAC/S) is less than one and AAC/S is greater than one, respectively.

In this work, we propose to correlate apparent accommodation trends laterally from one section to another, but with additional criteria such as the bathymetric variation over thickness ratio (ΔB/T) to reinforce the correlation hypothesis. The principle lies in the incremental analysis of variables used to calculate the apparent accommodation and to define the markers: the bathymetric variation between two markers divided by the thickness of the interval defined by the markers.

This ratio is powerful, as it can be interpreted in terms of sedimentary processes (regression, transgression, and sediment supply). Five cases can be encountered (Figure 5).

1. A ΔB/T ratio < −1 suggests a water-depth decrease higher than the interval thickness. This case involves a reduction of accommodation that could correspond to a forced regression and be followed by subaerially exposure. Such an interval may be capped by a type-1 sequence boundary sensu Schlager (2005).

2. A ΔB/T ratio = −1 corresponds to a water-depth decrease equivalent to the interval thickness. For instance, a 5-m-thick interval recording a bathymetric decrease of 5 m would have a ΔB/T = −1. In that case, no reduction of accommodation space is needed; the latter is filled only by the sediment supply. This condition describes progradation.

3. A ΔB/T ratio between −1 and 0 would suggest a water-depth decrease lower than the interval thickness. For instance, a 100-m-thick interval recording a water-depth decrease of −65 m would belong to this scenario. Accommodation space is filled faster than it is created, which consequently leads to a progressive decrease in water depths, defining a possible progradational pattern.

4. A ΔB/T = 0 ratio indicates an interval characterized by a constant water depth from base to top, independently of the thickness. The sediment supply is equivalent to the accommodation, indicating an aggrading pattern.

    The last three configurations would result in a type-2 sequence boundary sensu Schlager (2005).

5. A ΔB/T > 0 ratio would, finally, indicate a water-depth increase independently of the thickness. The higher the ΔB/T > 0, the higher the water-depth increase is recorded. It would clearly indicate an accommodation gain and a possible drowning process if the gain is important. This kind of ratio is associated with a vertical retrogradation pattern. Such an interval may be capped by a type-3 boundary sensu Schlager (2005).

Therefore, this ratio can help to quantify processes of apparent accommodation creation or destruction through time and space and is used together with the accommodation curve trends for correlating the sections through the study area.

In this study, the water depths of sedimentary facies are considered as input data for the calculation of apparent accommodation. These depths have been attributed to each facies based on interpretations from Purdy (1961), Hughes (2000), Hillgärtner et al. (2003), Masse et al. (2003), Immenhauser et al. (2004), Burla et al. (2008), Bover-Arnal et al. (2012), Bastide (2014), and Michel et al. (2023), who studied analogous sedimentary systems (Figure 6). Three water-depth categories were defined for each facies (see Table 1), representing the most probable depth, the minimum depth, and the maximum depth. The most probable depth corresponds to the most common bathymetry found for each facies, based on literature and field observations. It represents the depth at which the facies are most frequently encountered in analogous depositional systems. The minimum and maximum depths represent the uncertainty ranges for the bathymetric deposition of the facies. The minimum depth corresponds to the shallowest deposition depth recorded, whereas the maximum depth represents the deepest depth still consistent with the facies characteristics. These ranges are derived from both the geological literature and the known water-depth ranges of the producers and elements associated with each facies.

For instance, in Provence (France), Masse et al. (2003) and Fenerci-Masse et al. (2005) noted that rudist size increases with water depth. Their work suggested that muddy facies with small rudists were deposited in water depths of less than 1 m, whereas larger rudist-bearing grainy facies were found in the higher-energy subtidal zones of the platform interior, where bathymetric depths ranged from 1 to 5 m. This interpretation is applied to FA3 (Table 1), where the most probable depth is assigned as 1–5 m. The assignment of water depths to all other facies follows a similar approach, based on depositional settings observed in the literature and the facies characteristics, including organism size, sediment texture, and hydrodynamic conditions. Although some uncertainty remains in the depth ranges, they are consistent with general trends observed in analogous environments and are supported by the available geological data.

Apparent accommodation was calculated along each sedimentary section. The main evolutions obtained for Fontaine de Vaucluse (7 in Figure 1) and Martigues (3 in Figure 1) sections are shown below.

The Fontaine de Vaucluse section is 868 m thick and represents one of the most complete sections in the study area (lower Barremian to lower Aptian; Léonide et al., 2012; Frau et al., 2018). It required a positive cumulative apparent accommodation of +818 m to be formed, considering the section thickness and the facies water-depth difference (−40 m) between the base (70 m) and the top of the section (30 m) (Figure 7). The 300 first meters of the section are characterized by major changes of water depths, highlighting positive (P1/p2: uppercase = major event, lowercase = minor event) and negative apparent accommodation events (N1, n2, N3). For instance, the second interval formed in basin conditions during the early Barremian resulted from a positive apparent accommodation estimated at +234.5 m and, thus, a ΔB/T > 0 (P2, Figure 7). The base of the upper Barremian (Toxaster vandeheckii ammonite subzone) displays a shallowing upward trend with a decrease of water depth approximately −50 m recorded along a 200-m-thick interval, which resulted thus from a positive accommodation pattern (−1 < ΔB/T < 0). Above, the trend is different. Four cycles of both water depth and apparent accommodation (from p5 to N7 in Figure 7) are clearly observed. Accommodation decrease during N7 is estimated at −57 m. On the contrary, the first overlying rudist unit deposited in an invariable positive accommodation pattern characterizing a ΔB/T = 0. This trend changes suddenly during the Imerites giraudi ammonite subzone, which was characterized by a drastic positive accommodation event (+88.5 m P9 in Figure 7), subsequently followed by a negative one, the overall being recorded along a 30-m-thick interval only (n8, Figure 7). Like previously, the second rudist unit formed during a constant positive accommodation context (ΔB/T = 0) during the end of the Barremian. Finally, the top of the Fontaine de Vaucluse section displays a last long-term bathymetric rise and fall, emphasizing a positive and then negative evolution of apparent accommodation in the early Aptian.

The Martigues section consists of 415-m-thick Hauterivian to lower Aptian carbonates (Masse, 1976; Frau et al., 2018). It required a positive accommodation of +445 m to be formed that is almost half of the apparent accommodation recorded at Fontaine de Vaucluse (Figure 8). The Hauterivian–Barremian transition was marked by a drastic bathymetric decrease, emphasizing an intense negative accommodation (ΔB/T < 0, N1 in Figure 8). Two cycles of water depth are recognized within the lower Barremian, but contrary to Fontaine de Vaucluse, the apparent accommodation is constantly positive (ΔB/T either < 0 or > 0 depending on the bathymetric evolution). The end of the early Barremian was marked by the exposure of the inner platform (Figure 8). At the beginning of the late Barremian (T. vandeheckii ammonite zone), aggradation of the inner platform occurred, resulting thus from a positive accommodation. Above, like at Fontaine de Vaucluse, the trend is different. Here, three cycles (instead of four at Fontaine de Vaucluse) of apparent accommodation are observed and involved three open-marine incursions onto the inner platform (from p5 to N7 in Figure 8). The top of the last negative apparent accommodation (−30 m) event (N7) was marked by exposure of the platform and the development of a paleosol. The overlying deposits characterized once again the aggradation of the inner platform in a positive accommodation context (ΔB/T = 0) (Figure 8). A drastic increase of apparent accommodation (+90 m, P9) leading to bathymetric rise was recorded during the end of the late Barremian, whereas shallower settings returned to northern Provence during the same time interval. The subsequent deep basin settings prevailed in the early Aptian.

Twelve major markers of apparent accommodation variations have been defined and then correlated between the different sections of study area (Figure 9). This correlation has been carried out within chronostratigraphic framework defined by Masse (1976) and Frau et al. (2018) and modified in this study by updated descriptions.

A first negative apparent accommodation event (N1, ΔB/T < −1) dated from the base of lower Barremian recognized at Martigues (3 in Figure 1) and Fontaine de Vaucluse (FdV; 7 in Figure 1) can be easily correlated at Orgon (5 in Figure 1; Figure 10). Its top defines the M0 marker. At Cassis (2 in Figure 1), M0 is placed at the uppermost part of an aggrading FA2 rich interval showcasing evidence of subaerial exposure (Masse, 1976). Above M0, a cyclic evolution of water depths and apparent accommodations (from P1 to N3 events) is discernible in southern and northern Provence, marking a change in the accommodation trends.

The M1 marker marks the end of the cyclic evolution. In northern Provence, M1 corresponds to the top of a negative apparent accommodation event (N3), above which apparent accommodation consistently increases (Figure 10). At Mont Ventoux, the same evolution is observed, and M1 corresponds here to a firmground surface (Tendil et al., 2018). In southern Provence, the M1 marker is defined by an irregular stratigraphic surface showing subaerial exposure criteria. This surface separates two adjacent intervals formed in a positive accommodation context but defined by different bathymetry evolution and ΔB/T ratio. The cumulative apparent accommodation of the M0–M1 interval varies geographically with values increasing northward.

Whatever the geographic location, a change in the accommodation trend has been reported from the top of the upper Barremian T. vandenheckii ammonite subzone, defining the M2 marker. In southern Provence, M2 marks the culmination of an aggrading inner-platform interval (ΔB/T = 0) characterized by subaerial exposure evidence at Cassis (2 in Figure 1; Masse, 1976). From Orgon (5 in Figure 1) to Saint-Christol (16 in Figure 1; including Fontaine de Vaucluse), it is placed at the top of a shallowing upward interval (ΔB/T ratio < 0) and above which the accommodation pattern is increasing suddenly (ΔB/T > 0) (Figure 10). In northward sections such as at Banon (17 in Figure 1) and Mont-Ventoux (23 in Figure 1), the M2 marker would correspond to the top of a negative apparent accommodation event marked by a firmground surface (Tendil et al., 2018). The cumulative apparent accommodation of the M1–M2 interval reveals distinct patterns between southern Provence (+30 m) and northern Provence (from +130 to +198 m).

The M3 marker punctuates a stratigraphic interval characterized by several cycles of apparent accommodation across the study area. The M3 is identified at the top of the last negative apparent accommodation event (N7) (Figure 10). This event is associated with the development and major subaerial exposure of an inner platform from southern to northern Provence (Gorges de la Nesque section 20–22 in Figure 1) (Figures 9, 10). A drop of apparent accommodation ranging from −1.5 to −57 m and a bathymetric decrease from −10 to −65 m has been calculated. Above M3, apparent accommodation returns to positive and stable values over the whole study area.

The M4 marker is placed at the top of the aggrading inner-platform interval deposited in positive AAC context (ΔB/T = 0) and above which stand regional major changes of both lithostratigraphy and accommodation evolution (P9, Figure 10). Whatever the geographic location, a sharp apparent accommodation gain (ΔB/T > 0) is observed above M4: +90 to +130 m in southern and middle Provence, +10 to +90 m in northern Provence (Figure 9). For the first time, apparent accommodation does not increase northward; an east–west trend is, rather, reported.

The M5 marker is supposed to date from the basal part of the upper Barremian Martelites sarasini ammonite subzone. It separates an underlying negative apparent accommodation interval (n8, ΔB/T < −1) from an overlying positive one (Figure 9). The M4–M5 interval recorded a cumulated positive accommodation ranging from +4.5 to +81.5 m from northern Provence to southern Provence. In northern Provence, the eastern part still exhibits lower apparent accommodation values (+10 m) than the western part (+25 m).

The M6 marker has been identified in middle and northern Provence only where it represents the top of an aggrading inner platform formed in a positive accommodation pattern during the late Barremian (M. sarasini ammonite subzone). In northern Provence, it corresponds to a major stratigraphic surface (D2 sensu Masse et al., 2011), resulting from subaerial exposure and correlated at Banon and Mont-Ventoux with basinal an aggrading interval (17 and 23, respectively, in Figure 1).

The M7 marker is identified at the Banon and Mont-Ventoux sections (Figure 9), marking the top of a major negative apparent accommodation event (N9). This event is associated with the deposition of platform-derived facies associations (FA9) onto basinal facies (FA11/FA12). The estimated drop in apparent accommodation would range from −40 to −70 m. Above, apparent accommodation drastically increased again.

In northern Provence, a new increase of accommodation (ΔB/T > 0) led to the drowning of the inner platform (P11 in Figures 7–12). But in the eastern part of northern Provence only, this new pulse of positive apparent accommodation dated here as early Aptian (base of Deshayesites oglanlensis ammonite subzone) is quickly followed by a negative AAC interval, whose top defines the M8 marker. The latter were identified from Font Jouvale to Banon sections (10 and 17 in Figure 1, respectively), which suggests a possible local tectonic event. Apparent accommodation from the M7–M8 interval ranges from 0 to +59 m and increases westward and northward from Oppedette, where no records have been observed (Figure 12).

A turnover of the apparent accommodation evolution has been reported at the top of early Aptian D. oglanlensis ammonite subzone (Tendil et al., 2018) from the stratigraphy in northern Provence. This change marks the M9 marker, which corresponds here to a hardground surface (Tendil et al., 2018). The M9 punctuates an interval defined either by a ΔB/T ratio = 0 or < 0 (i.e., Gorges de la Nesque) and above which a regional rise of apparent accommodation has been characterized (Figure 9). The apparent accommodation of the M7–M9 interval in northern Provence displays once again a positive trend westward. It may suggest the influence of an active fault system (such as the Salon-Cavaillon system fault) contributing to the progressive structuring of the western part of northern Provence and the deepening of the surrounding depositional areas (see further discussion).

Like the previous marker, the M10 marker corresponds to a hardground surface dated from the top of the lower Aptian Deshayesites furcata ammonite subzone (Tendil et al., 2018) in northern Provence. It separates an underlying interval consisting of outer-platform facies associations, deposited in a slightly positive AAC context with values ranging from +2 to +31 m. Above this interval, marly basinal facies dominate, indicating a significant increase in AAC values (Figure 9). The apparent accommodation calculated for the lower Aptian M9–M10 interval emphasizes high lateral variations in northern Provence. Here, apparent accommodations increase from +19 m at Oppedette to +31 m at Simiane la Rotonde (Figure 12), both sections being separated by only 4.5 km (13 and 14 in Figure 1). The same pattern is observed in the western part of northern Provence.

The M11 marker is dated from the top of the lower Aptian D. furcata ammonite subzone by Frau et al. (2018) and has been recognized in several sections only due to post-Durancian erosion. It corresponds to the base of the global increase of apparent accommodation (+100 m) in northern Provence (Figure 9).

The definition of the stratigraphic architectures of the Urgonian system of Provence is based on the interpretation of the spatial and temporal evolutions of the apparent accommodation, ΔB/T ratio, facies proportions, and nature of stratigraphic surfaces. For each interval, a broad environmental trend followed by more detailed descriptions is given.

The first M0–M1 interval mainly characterizes an inner platform in southern Provence while bioclastic sand shoals developed in open marine influences until Orgon in middle Provence (Figure 11). Northward, slope facies associations and basin conditions developed (Figures 11, 12). The basal P1 positive accommodation event marks a southward retrogradation of the facies belts. The latter is subsequently followed by a general progradation phase interpreted from negative ΔB/T ratios. During the same time interval, aggradation of a restricted inner platform is inferred at Cassis (Figures 11, 12). Finally, the second prograding interval would be associated with the N3 major negative apparent accommodation event, which led to the decrease of the water depth and the subaerial exposure of the inner platform from southern Provence to Saint Chamas. Spatial facies evolution from southern to northern Provence indicates a ramp profile sensu Tucker and Wright (1990). The angle of this ramp could range from 0.02° to 0.10° from Cassis to Orgon and northward from Orgon, respectively.

Restricted inner-platform conditions extended from southern to middle Provence during the M1–M2 interval. Northward open marine influences still dominated, and bioclastic sand shoal complex developed until Saint Christol and Gorges de la Nesque in northern Provence (Figures 11, 12). This pattern results in the flattening of the ramp profile with mean angle ranging from 0.01° to 0.04°. No subaerial exposure criteria have been observed at the top of the interval. The bathymetric evolutions and ΔB/T ratios calculated at the base of the M1–M2 interval in northern Provence would suggest a possible diachronous sedimentation pattern compared to southern and middle Provence. The first pulse of accommodation only recorded in northern Provence results in the development of margin and slope environments extending to middle Provence, while southern Provence remained exposed (Figures 11, 12). Onlapping geometries would be expected from bathymetric and accommodation patterns only. The uppermost part of the M1–M2 unit is characterized by a positive accommodation pattern recorded throughout the study area. The ΔB/T ratio = 0 in southern Provence still points to the aggradation of the inner platform. The ΔB/T ratios < 0 in middle Provence and northern Provence emphasize the efficiency of the carbonate production factory to fill the created available space during deposition of the M1–M2 unit (Figures 11, 12).

The M2–M3 interval is characterized by the continued development of inner-platform environments in southern and middle Provence, while open marine conditions prevailed in northern Provence. Slope environments in the northernmost sections, such as Mont-Ventoux and Banon, transition into basinal settings (Figures 11, 12). As described in the Apparent Accommodation Calculation and Marker Identification section, this interval is marked by three to four cycles of positive and negative apparent accommodation. In northern Provence, four distinct cycles, with thicknesses ranging from 15 to 100 m, have been identified. These cycles reflect ∼10-km-scale retrogradations and progradations of facies belts, inferred from the alternations between slope and outer-platform environments (Figures 11, 12). In contrast, in southern Provence, three cycles of apparent accommodation are recorded, with facies alternations between inner-platform and outer-platform environments over accumulations of 15–25 m.

The second cycle within this interval exhibits maximum retrogradation of facies belts, corresponding to a general bathymetric increase. This retrogradation is evidenced in both southern and northern Provence. The uppermost N7 negative apparent accommodation episode recorded in most of the sections from the study area results in (1) a drastic seaward migration of the facies belts, (2) the overall extension of inner-platform environments from southern to northern Provence, and (3) a regional subaerial exposure till the Gorges de la Nesque and Rustrel sections. In southern Provence, this episode is coeval with the last occurrence of Agriopleura rudist fauna.

In middle and northern Provence, the base of the M3–M4 interval, formed in a positive apparent accommodation framework, is marked by the accumulation of Agriopleura-rich rudist facies (FA3). In southern Provence, the latter rudist fauna disappeared with the underlying negative apparent accommodation event and exposure of the platform. This consequently reflects the development of an inner platform from middle Provence to northern Provence onlapping the subaerially exposed southern Provence (Figures 11, 12). As apparent accommodation continues to increase, inner-platform environments developed from southern to northern Provence, indicating a very flat topography over more than 120 km with some structural irregularities such as at Saint-Chamas. The inner platform was protected from the high hydrodynamic open-marine conditions by an ∼3-km-wide coral complex at the Gorges de la Nesque section and at the Simiane la Rotonde section eastward (20–22 and 14 in Figure 1; Figures 11, 12). Finally, lower hydrodynamic open-marine environments are recorded northward. No subaerial exposure criteria have been observed at the top of the interval, except some local beach-rocks (FA1) at Gorges de la Nesque section (Léonide et al., 2012).

The regionally observed positive P9 event points to the drowning of the underlying inner platform, but the latter has been differently expressed, depending on the geographic position. In northern Provence, high hydrodynamic shallow outer platform to slope environments developed in its eastern and western parts, respectively. The first Palorbitolina unit is deposited in the western part of northern Provence exclusively. Basin environments still characterize northern sections from Banon and Mont-Ventoux (Figures 11, 12). In southern Provence, a higher pulse of apparent accommodation than northward boosts the development of slope to basin conditions. Consequently, at this time, the shallowest environments are restricted to an isolated zone located in the eastern part of northern Provence. The overlying n8 negative apparent accommodation event is associated with the decrease of water depths and the return of shallow subtidal inner-platform conditions in only northern Provence, while the southern Provence remains under basin settings.

In northern Provence, this interval is characterized by the aggradation of inner-platform environments. No coral complex has been observed on the platform margin to protect the inner restricted realm as for the underlying intervals. Besides, the accumulation of coral gravel and coarse rudist facies associations indicates a higher hydrodynamic inner platform than in previous intervals. These environments extended until La Flachère and Gorges de la Nesque (15 and 20 in Figure 1) through 10- to 15-km northward progradation, which rapidly evolved to the outer platform and basin conditions. Similarly, a 5- to 10-km southward progradation of the inner-platform environments is inferred at Rustrel and Lagnes. Consequently, the M5–M6 interval characterizes the general development of the northern Provence isolated carbonate platform, on either side of which slope and basin environments developed in less than ∼10 km, suggesting high slope angle from 0.4° to 1.15°. The top of this interval exhibits a very sharp planar surface to irregular epikarstic surface interpreted by Tendil et al. (2018) as emersion surface. A second high-hydrodynamic inner carbonate platform is recorded at Saint-Chamas (4 in Figure 1) in middle Provence, while southern Provence remains under the influence of deep basin conditions. Recent erosion and lack of lateral continuity of outcrops prevent knowing if such environments developed until Orgon (5 in Figure 1) and if lateral connectivity between this platform and the northern Provence inner platform exists.

The major N7 negative apparent accommodation event (ΔB/T < −1) during the M6–M7 interval caused a significant decrease in water depth, resulting in the development of slope environments that replaced the deeper conditions in the northernmost part of the platform (Figures 10, 11). During this interval, the supposed subaerial exposure of the northern Provence inner platform thus evolved to slope environments in less than 10 km, suggesting slope angle from 0.2° to 1.0°. Basin conditions still prevailed in middle and southern Provence.

This interval, dated to the Barremian–Aptian transition, characterizes an apparent accommodation cycle recorded only in the eastern part of northern Provence. The base of the cycle marks the drowning of the previously exposed inner carbonate platform in northern Provence, while southern and middle Provence remained influenced by basin conditions (Figures 10, 11). The post-Aptian Durancian erosion does not allow for study of the relationships between the middle to southern Provence areas. The westward increase of apparent accommodation from Oppedette (13 in Figure 1), associated with the progressive deepening of depositional environments, points to a structuring of the sedimentary profile, evolving to a 60-km-wide gentle slope from Oppedette toward Fontaine de Vaucluse and Lagnes sections (Figure 12). A similar environmental evolution is described from Oppedette toward Banon northward in only 7 km, suggesting a more abrupt ramp profile. The M7–M8 would thus highlight possible structural influences on the sedimentation pattern.

Although the M8–M9 interval is characterized by an aggrading trend in the western part of northern Provence, a different pattern is observed in its eastern part. Indeed, bioclastic sand shoals finally covered the previous exposed Oppedette area (13 in Figure 1). These conditions extended to Banon at the top of the interval, resulting from an uppermost negative apparent accommodation event locally recorded (ΔB/T ratio < −1). This points once again to the activity of fault systems, which would have led to the uplift of the Banon area and the appearance of shallower and higher hydrodynamic conditions northward (Figure 12). The top of this interval is marked by a hardground regionally observed (Tendil et al., 2018). In southern and middle Provence, deep and low energy basin to slope environments still dominated.

The M9–M10 interval characterizes the broad extension of a high hydrodynamic outer-platform environment across northern Provence surrounded by slope and basin environments both to the southern and northward from the Gorges de la Nesque area. In more detail, it appears that kilometer-wide areas upon which only outer-platform environments developed (such as at Joucas and Oppedette) are separated by deeper zones where slope environments prevailed (Figures 11, 12). The progradation of bioclastic sand shoal into these depressions led to the decrease of water depth and the general extension of outer-platform conditions. The development of the deeper depositional conditions coincides with the geographic distribution of the highest apparent accommodation values and ΔB/T ratio > 0 at the base of the interval, suggesting the impact of local fault activity. According to Kandel (1992), kilometer-wide half graben systems may have been created throughout northern Provence in the early Aptian, affecting the distribution of the depositional environments.

The M10–M11 interval, dated from the lower Aptian Deshayesites forbesi ammonite subzone, characterizes an overall increase of apparent accommodation (+100 m), associated with the widespread development of basin conditions from southern to northern Provence.

Stratigraphic architecture of the Urgonian platform is interpreted with respect to the paleoclimatic, oceanographic, eustatic, and tectonic evolution in the Early Cretaceous.

The Barremian and early Aptian are commonly considered greenhouse periods, with the study area located approximately 28–30°N from the paleoequator, creating favorable conditions for intense carbonate accumulation on the southern margin of the Vocontian Basin.

During the Hauterivian to early Barremian, mesotrophic to eutrophic conditions, inferred from high nutrient levels and a humid climate (Godet et al., 2006; Huck et al., 2011, 2013), may have (1) limited the development of photosynthetic organisms and (2) restricted inner platform growth to southern Provence during the M0–M1 interval. The mid-Barremian event affected the western Tethys at the end of the early Barremian (Mutterlose et al., 2009; Huck et al., 2011, 2013), but its impact is not clearly visible in the stratigraphy of southern and northern Provence.

From the late Barremian (T. vandehenckii to I. giraudi ammonite subzones), bioclastic sand shoal complexes developed widely (M1–M2 interval), followed by the expansion of inner-platform environments into northern Provence (M2–M4 interval). This environmental change, the extensive growth of photozoan organisms (e.g., rudists and green algae), and the evolution from a ramp to a carbonate platform suggest increased carbonate production, probably influenced by a return to a colder, arid climate (Bodin et al., 2009; Mutterlose et al., 2009). The M4–M5 interval (upper Barremian I. Giraudi subzone) records the decline of photozoan organisms (Mutti and Hallock, 2003) and a resurgence of heterozoan groups, likely due to mesotrophic to eutrophic conditions associated with the warmer, short-lived, humid Taxy event, which led to higher nutrient levels and organic-rich layers in the western Tethys (Huck et al., 2011, 2013). The end of the late Barremian was characterized by a return to oligotrophic conditions, fostering the growth of a carbonate platform in northern Provence. This pattern, also observed in the Jura and Helvetic regions, would be concomitant with the global return to arid and cold conditions (Bodin et al., 2009; Mutterlose et al., 2009).

The early Aptian was marked by extreme temperature increases, eutrophic conditions, humidity, and rising volcanic activity. In northern Provence, the drowning of the late Barremian inner platform and the expansion of outer-platform environments dominated by heterozoan organisms during the M6–M9 interval indicate the return of mesotrophic conditions and deteriorating climatic and environmental conditions. The final drowning of the northern Provence outer platform during the M10–M11 interval indicates a decrease in carbonate production, coinciding with the OAE1a (Erba et al., 2015).

Stratigraphic architectures can be good indicators of eustatic variations under particular geodynamic conditions. In a nonsubsiding area, sea-level falls can result in a seaward migration of the facies belt by sedimentary progradations and a common (but not systematic) subaerial exposure of the inner part of the sedimentary profile. On the contrary, sea-level rises can be associated with the progressive deepening of the water depths if sediment supply is lower than the accommodation, which results in a landward migration of the facies belt and a retrogradation. The quantitative sequence stratigraphic approach (Borgomano et al., 2020) can help to reconstruct the eustatic sea-level changes that are discussed in the following paragraphs.

The consistency of the eustatic curve from Haq (2014) is discussed according to the apparent accommodation and stratigraphic architectures of Barremian to lower Aptian sequences. Four regional negative and positive apparent accommodation events associated with progradational and retrogradational patterns, respectively, have been identified in the Urgonian carbonates of Provence (until the Deshayesites deshayesi ammonite subzone). During the same time interval, Haq (2014) reported four third-order eustatic cycles, each of them lasting approximately 1.4 to 1.7 m.y.

The first negative event of apparent accommodation (below the M0 marker) in the early Barremian would be time-equivalent to the worldwide drastic eustatic fall (–95 m) dated at 130.4 Ma by Haq (2014) (Figure 13). In Provence, the apparent accommodation was approximately −50 m for M0, suggesting that sea-level fall would have been compensated by regional subsidence. Following M0, a regional increase in apparent accommodation and facies retrogradation occurred during the early Barremian (Nicklesia nicklesi to N. pulchella ammonite zones; 130.4 to 129.7 Ma; Emericiceras sp. ammonite found in Fontaine de Vaucluse section by Masse, 1976). This event aligns with the sea-level rise reported by Haq (2014) in the early Barremian, associated with a +110-m rise. The latter, however, is lower than the apparent accommodation reported in middle and northern Provence (+250 to +360 m), suggesting once again regional subsidence influences.

A subsequent worldwide eustatic drop was reported by Haq (2014) before the M1 geochron at the base of the upper Barremian. However, criteria such as the northward facies progradation, decreasing water depths, and inner-platform subaerial exposure in southern Provence suggest an accommodation space reduction, related to eustacy and/or tectonics. As the M1 marker separates two adjacent units deposited under different climatic and environmental conditions, whose transition is dated from the top of the Moutoniceras moutonianum and Paraceratites elegans ammonite subzones in the Tethyan and Boreal regions, respectively (McArthur et al., 2004), M1 would be consequently coeval with the end of the eustatic fall reported by Haq (2014). Nonetheless, the fact that no such event would be recognized in the adjacent Tethyan realm (Frau et al., 2020) suggests either a limited fall of sea level and/or a sea-level fall greatly compensated by an important regional subsidence. Indeed, the high variability of interval thickness of the sections points also to an important differential subsidence effect.

In the study area, a second major positive accommodation event recognized at the base of the M2–M3 interval is associated with general bathymetric deepening and a regional retrogradation of facies belts. No climatic or environmental changes are reported to explain this evolution. The creation of accommodation space by a rise in sea level may be a possibility. Nevertheless, no eustatic rise has been reported on the reference Cretaceous curve during the late Barremian (G. sartousiana ammonite zone). Indeed, the eustatic rise reported by Haq (2014) would end earlier in the late Barremian (T. vandenheckii ammonite zone) (Figure 13). It would thus be naturally associated with the deposition of the underlying M1–M2 interval. Other processes such as tectonic activity might have influenced the accommodation pattern during the M2–M3 interval in Provence, leading to a regional increase of apparent accommodation.

The sudden seaward migration and subaerial exposure of inner-platform environments at the top of the M2–M3 interval indicate a regional drop in apparent accommodation (ΔB/T ratio < −1). This event, dated to the late Barremian (G. sartousiana ammonite zone), aligns with time-equivalent observations in adjacent Tethyan realms (Bas-Vivarais, Switzerland, Spain; Frau et al., 2020) and may point to a global fall in sea level. Haq (2014) reported such a worldwide eustatic fall at the beginning of the late Barremian, which can thus be considered to have also impacted southern France (Figure 13).

The drastic regional increase of apparent accommodation during the M4–M5 interval, leading to the drowning of the underlying inner platform and the reappearance of deeper depositional environments from southern to northern Provence can be interpreted in different ways. On the one hand, the return of climatic-driven mesotrophic conditions may have contributed to weakening the activity of carbonate factory, limiting the carbonate accumulations, and favoring the increase of water depths. Nevertheless, according to Haq (2014), a sea-level rise phase ended during the late Barremian I. giraudi ammonite subzone (Figure 13). Therefore, continuous carbonate production of the inner platform (M3–M4 interval) likely initially counteracted the sea-level rise, maintaining a relatively stable water depth and promoting aggradation. However, the subsequent drowning of the latter may thus reflect the cumulative effect of the eustatic rising and the weakening activity of the carbonate factory due to the mesotrophic climatic conditions. Finally, considering the nonuniform geographic distribution of the apparent accommodation from northern Provence to southern Provence and from east to west in northern Provence, tectonic influence may have contributed to the stratigraphic architectures during the upper Barremian M4–M5 interval.

During the late Barremian, the return to arid and cold climatic conditions led to the reappearance of photozoan organisms in the M5–M7 interval in northern Provence. These accumulations are marked by regional stratigraphic surfaces likely due to subaerial exposure (Tendil et al., 2018). Similar criteria observed in other Mediterranean Tethys platforms (Frau et al., 2020) suggest a global eustatic fall at the end of the Barremian, as reported by Haq (2014). In northern Provence, subaerial exposure of the inner platform coincided with a negative apparent accommodation event in basin sections (e.g., Banon, Mont-Ventoux). Borgomano et al. (2020) noted that proximal records are incomplete due to exposures and hiatuses, whereas distal sections preserve the complete accommodation signal. This supports the hypothesis of a global eustatic fall at the end of the Barremian (Haq, 2014). Conversely, the positive apparent accommodation in southern Provence indicates local structural processes causing high subsidence.

At the beginning of the early Aptian, the inner exposed platform in northern Provence drowned. The basal retrogradation of facies belts and increased apparent accommodation could have been induced by a eustatic rise, such as the one reported by Haq (2014) at the base of the lower Aptian, which would apparently be slightly postponed related to the drowning of the inner platform. Furthermore, except in Spain, where basin conditions developed onto the underlying inner exposed platforms, synsedimentary hiatus occurred in other platforms from the Mediterranean Tethys (Frau et al., 2020), pointing probably to a limited eustatic influence. Besides, the westward and northward nonuniform increase of apparent accommodation from the Oppedette area and the creation of kilometer-scale depocenters and topographic highs suggest a structuring of the depositional area by active fault systems. Whatever the eustatic and/or tectonic origins, drowning is coeval with the worldwide degradation of the climatic and environmental conditions, which marks the return of mesotrophic conditions. Such conditions could have weakened the activity of carbonate factory, contributing to limit the carbonate accumulation and consequently to increase the water depths. Finally, the M10–M11 interval is characterized by a general increase of apparent accommodation, and the development of deeper basin conditions may result from the cumulative effect of the sea-level rise reported on the Cretaceous reference curve (Haq 2014) during the OAE1a.

Carbonate sequence stratigraphic models typically assume that a single A/S signal represents the entire stratigraphic architecture, allowing correlation throughout the carbonate system. However, this assumption frequently overlooks the variability in accommodation and sediment supply over time. Additionally, A/S ratio changes are commonly interpreted as sea-level changes (Sharland et al., 2001), but such models are rarely validated through quantitative forward modeling.

The first point is the fact that this method offers the opportunity for challenging the interpretation of the facies association, their relationships, and, finally, the evolution of the deposition environment (i.e., water depth). Although this study uses water depths from previous research, they can also be inferred from sedimentological interpretations of facies. Incorporating uncertainty in water depths allows for the construction of multiple accommodation scenarios, which can be tested, validated, or adjusted using stratigraphic forward modeling tools. Unlike the assumption of uniform eustatic variations, focusing on accommodation reveals interactions with other factors, such as subsidence patterns, which can vary geographically. Consequently, considering a reference eustatic curve for a given time interval, it is possible to quantify the influence of the apparent accommodation resulting from differential subsidence in every location of a sedimentary basin.

The second point is that we have demonstrated that an explicit method is possible. Apparent accommodation curves can be constructed for all locations within a carbonate system using measurable data from outcrops or cores. A high-resolution stratigraphic framework can then be developed through (1) lateral correlation of apparent accommodation trends and ΔB/T ratios, and (2) interpretation of major stratigraphic surfaces (exposure or drowning surfaces). However, without a well-constrained chronostratigraphic framework, these correlations can lead to inconsistencies and diachronic issues. For example, correlating negative accommodation events from the inner platform to the outer shelf can be complex, but with accurate chronostratigraphic data, we can align inner-platform exposures (ΔB/T ratios −1 to 0) with outer-platform and outer-shelf negative accommodations (ΔB/T ratio < −1). In the inner platform, because of exposures and stratigraphic hiatuses, apparent accommodation does not match the real accommodation systematically, whereas the record of the latter is more continuous in the deeper environments due to higher accommodation space.

This method has improved the characterization of stratigraphic geometries in the studied area, identifying lowstand intervals in northern Provence, while southern Provence remained subaerially exposed. By analyzing accommodation below and above major stratigraphic surfaces, we achieve a more robust interpretation of stratigraphic architecture.

Nevertheless, this method has limitations. It relies heavily on prior estimation of facies depositional environments. Sharp changes in apparent accommodation smaller than paleo water depth ranges may not be detected. A well-constrained chronostratigraphic framework is also essential for accurate time-thickness functions and carbonate production rates. Without it, accommodation markers may intersect real chronostratigraphic timelines, distorting reservoir architectures. For subsurface data with uncertain chronostratigraphy, we suggest multiple scenarios of accommodation and sedimentation rates tested using stratigraphic forward modeling tools.

Finally, the method’s suitability depends on the sedimentary system and case study. It is most effective when significant spatial and temporal variations in thickness and facies water-depth exist. Even in cases with uniform water depth, such as Middle Eastern carbonate reservoirs, spatial thickness variations can still reveal accommodation changes. These variations, influenced by local tectonics, can be quantified and tested through stratigraphic forward modeling.

When seismic images or continuous outcrops are unavailable, constructing three-dimensional (3-D) stratigraphic architectures of carbonate platforms typically relies on standard sequence stratigraphic correlations. However, this method, despite its common use, is implicit, qualitative, and based on strong assumptions.

This study introduces a novel approach to establish the stratigraphic architectures applied to the Barremian–lower Aptian Urgonian carbonate platform in Provence. It uses a recently published quantitative method that integrates the correlation of apparent accommodation trends across distant locations and introduces a new parameter: the ΔB/T ratio. Unlike the traditional A/S ratio, the ΔB/T ratio is directly measurable and enhances understanding of the processes governing the formation of specific stratigraphic intervals and their correlation.

Four successive major negative apparent accommodation events were identified, predominantly in outer-shelf domains, and were correlated with subaerially exposed units in shallower inner environments. Positive accommodation events were easier to correlate, as they were consistently recorded across all locations of the carbonate sedimentary system.

This new high-resolution quantified stratigraphic framework reveals a morphologic evolution of carbonate platform during the early Barremian to the beginning of the late Barremian. It highlights the cumulative influence of carbonate production and negative apparent accommodation events, driving multiple seaward progradations of the facies belt. Spatial variability in the cumulative apparent accommodation suggests significant tectonic influences on platform evolution. The end of the late Barremian was marked by an inversion of the accommodation polarity toward southern Provence, pointing to geodynamic structuring of the study area. In the early Aptian, successive accommodation pulses and a climatic factor would have led to the progressive drowning of the platform.

In conclusion, the method of quantitative sequence stratigraphy and the introduction of the ΔB/T ratio represent an interesting advancement. They provide an explicit and nuanced understanding of the sedimentological processes driving stratigraphic evolution. This method not only facilitates robust sedimentological interpretations, but also allows for the formulation of multiple geological scenarios of accommodation evolution and needs to then be tested using 3-D forward sedimentological reservoir simulations.

When seismic data or continuous outcrops are unavailable, constructing 3-D stratigraphic architectures of carbonate platforms frequently relies on standard sequence stratigraphic correlations. Although widely used, this traditional method is largely qualitative, implicit, and based on strong assumptions, which can limit its predictive capabilities.

This study introduces a quantitative approach to define stratigraphic architectures, applied to the Barremian–lower Aptian Urgonian carbonate platform in Provence. The approach integrates apparent accommodation trends across geographically distant locations with a new parameter, the ΔB/T ratio. Unlike the traditional A/S ratio, which is inferred indirectly, the ΔB/T ratio is directly measurable and provides additional insights into the processes driving the formation and correlation of specific stratigraphic intervals.

Four successive major negative apparent accommodation events were identified, predominantly in outer-platform environments, and were correlated with subaerial exposure surfaces in shallower inner-platform settings. In contrast, positive apparent accommodation events were consistently recorded across all locations, making them easier to correlate.

The results reveal a detailed morphologic evolution of the carbonate platform from the early Barremian to the beginning of the late Barremian. The platform evolved through multiple seaward progradations of facies belts, driven by the combined influence of carbonate production and episodic negative apparent accommodation events. Spatial variability in cumulative apparent accommodation highlights significant tectonic controls on platform evolution. By the end of the late Barremian, a shift in accommodation polarity toward southern Provence indicates regional geodynamic restructuring. During the early Aptian, successive accommodation pulses, coupled with climatic influences, likely drove the progressive drowning of the platform.

The method of quantitative sequence stratigraphy, combined with the ΔB/T ratio, represents a significant advancement in understanding carbonate platform dynamics. This approach provides a nuanced, explicit interpretation of the sedimentological processes governing stratigraphic evolution. Furthermore, it enables the construction of robust geological scenarios for accommodation evolution, which can be validated and refined through 3-D forward sedimentological and reservoir simulations.

The authors thank the reviewers for their invaluable contributions, which have significantly improved the quality of this paper. We also extend our gratitude to TotalEnergies for funding this R&D project and granting permission to publish this work.

Supplementary material is available in an electronic version on the AAPG website (www.aapg.org/datashare) as Datashare 202.

Mickaël Barbier earned his Ph.D. in geosciences in 2012 from Aix Marseille University, France. His career spans roles in industrial and research-integrated projects aimed at understanding multiscale distribution of sedimentological, diagenetic, and structural heterogeneities in carbonate reservoirs. Currently at Akkodis, his expertise focuses on the characterization and three-dimensional modeling of carbonate reservoirs to provide advanced predictive modeling techniques.

Jean Borgomano earned his Ph.D. in carbonate geology in 1987 from Aix-Marseille University, France. He worked for Shell International as a carbonate geologist in exploration–production (1988–2003). As a professor, he led the Carbonate Reservoir Laboratory at Aix-Marseille University. He worked as a carbonate expert at Total (2013–2015) and then returned as a professor at Aix-Marseille University. His research focuses on the characterization and geological–numerical modeling of carbonate reservoirs.

Philippe Léonide earned his Ph.D. in sedimentology in 2007 from Aix-Marseille University. He joined the sedimentology and marine geology group at Vrije Universiteit Amsterdam (2009–2011). He is currently an assistant professor in carbonate sedimentology at Aix-Marseille University. His research focuses on the evolution of carbonates through time, which have importance for the characterization of petrophysical properties in carbonate reservoirs.

Gérard Massonnat earned his Ph.D. in hydrogeology and then graduated from IFP School in Petroleum Engineering and Project Development. After several positions in field monitoring, multidisciplinary studies, and research and development (R&D), he is presently an International Expert for Reservoir Geology and Geomodeling and R&D Fellow for TotalEnergies. His work now focuses on the development of the next generation of modeling tools for both matrix and dual porosity reservoirs.

Charles Danquigny obtained his Ph.D. in fluid mechanics in porous media from the University of Strasbourg (France) in 2003. He has been an associate professor at the University of Avignon (France) since 2006 and seconded to TotalEnergies since 2016. His research focuses on characterization and modeling of geological reservoirs, with a particular interest in predictive modeling of heterogeneity and upscaling.

Jean-Louis Lesueur earned his Ph.D. in lacustrine carbonate geology in 1991 from Bordeaux 3 University, France. With 18 years of experience abroad, he has held various roles including reservoir geology manager and geoscience manager. At headquarters, he served as deputy for reservoir geology, department deputy for geological specialties, and senior advisor for the exploration director. He joined the R&D entity 4 years ago.