Downward hydrocarbon migration predicted from numerical modeling of fluid overpressure in the Paleozoic Anticosti Basin,eastern Canada

The Anticosti Basin is a large Paleozoic basin in eastern Canada where potential source and reservoir rocks have been identified but no economic hydrocarbon reservoirs have been found. Potential source rocks of the Upper Ordovician Macasty Formation overlie carbonates of the Middle Ordovician Mingan Formation, which are underlain by dolostones of the Lower Ordovician Romaine Formation. These carbonates have been subjected to dissolution and dolomitization and are potential hydrocarbon reservoirs. Numerical simulations of fluid‐overpressure development related to sediment compaction and hydrocarbon generation were carried out to investigate whether hydrocarbons generated in the Macasty Formation could migrate downward into the underlying Mingan and Romaine formations. The modeling results indicate that, in the central part of the basin, maximum fluid overpressures developed above the Macasty Formation due to rapid sedimentation. This overpressured core dissipated gradually with time, but the overpressure pattern (i.e. maximum overpressure above source rock) was maintained during the generation of oil and gas. The downward impelling force associated with fluid‐overpressure gradients in the central part of the basin was stronger than the buoyancy force for oil, whereas the buoyancy force for gas and for oil generated in the later stage of the basin is stronger than the overpressure‐related force. Based on these results, it is proposed that oil generated from the Macasty Formation in the central part of the basin first moved downward into the Mingan and Romaine formations, and then migrated laterally up‐dip toward the basin margin, whereas gas throughout the basin and oil generated in the northern part of the basin generally moved upward. Consequently, gas reservoirs are predicted to occur in the upper part of the basin, whereas oil reservoirs are more likely to be found in the strata below the source rocks. Geofluids (2010) 10 , 334–350     

Red abalone collecting and marine water temperature during the Middle Holocene occupation of Santa Cruz Island,California

Oxygen isotope values derived from prehistoric mussel (Mytilus californianus) shell calcite are used to determine whether sea-surface temperatures in the vicinity of the Punta Arena site (CA-SCRI-109) on Santa Cruz Island, California, were cooler than present between 6300 and 5300 cal BP. This site and others in the western sector of the island dating to this period are distinctive because of the presence of large red abalone (Haliotis rufescens) shells, a species that historically has been subtidal around Santa Cruz and the other northern Channel Islands. Comparison of temperature values derived from archaeological shells with those from modern mussel shells, along with water temperatures derived from satellites, indicates that waters were cooler during the 6300 to 5300 cal BP period. These results are consistent with a previous study and support the interpretation that collection of red abalone during the 6300–5300 cal BP period was partly the result of cooler water temperatures that made this large mollusk more available to foragers in the intertidal or shallow subtidal zone.     

Hydrothermal convection in and around mineralized fault zones: insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana

Among hydrogeological processes, free convection in faults has been cited as a possible cause of gold mineralization along major fault zones. Here, we investigate the effects of free convection to determine whether it can cause giant orogenic gold deposits and their regular spatial distribution along major fault/shear zones. The approach comprises: (i) coupled two- and three-dimensional numerical heat- and fluid-flow simulations of simplified geological models; and (ii) calculation of the rock alteration index (RAI) to delineate regions where precipitation/dissolution can occur. Then, comparing the deduced alteration patterns with temperature distribution, potential areas of gold mineralization, defined by T  > 200°C and RAI < 0, are predicted. The models are based on the orogenic Paleoproterozoic ore deposits of the Ashanti belt in western Africa. These deposits occur in the most permeable parts of the fault zone, where the lateral permeability contrast is the highest. For a simple geometry, with a fault zone adjacent to a sedimentary basin half as permeable, we note a transition from three-dimensional circulation within the fault to a two-dimensional convective pattern in the basin far from the fault. Moreover, whereas two-dimensional undulated isotherms dominate in the basin, three-dimensional corrugated isotherms result from the preferred convective pattern within the fault, thus enhancing a periodic distribution of thermal highs and lows. In our most elaborate three-dimensional model with an imposed lateral permeability gradient, the RAI distribution indicates that fluid circulation in fault zones gives rise to a spatial periodicity of alteration patterns consistent with field data.     

The origin of overpressure in 'old' sedimentary basins: an example from the Cooper Basin, Australia

Overpressure in ‘old’ sedimentary basins that have not undergone rapid, recent sedimentation cannot be easily explained using traditional burial‐driven mechanisms. The last significant burial event in the Cooper Basin, Australia, was the Late Cretaceous deposition of the Winton Formation (98.5–90 Ma). Maximum temperature in the basin was attained during the Late Cretaceous, with cooling beginning prior to 75 Ma. Hence, overpressure related to rapid burial or palaeomaximum temperatures (e.g. hydrocarbon generation) must have developed prior to 75 Ma. Retaining overpressure for 75 Ma in ‘old’ basins such as the Cooper Basin requires extremely low seal permeabilities. An alternative explanation is that overpressure in the Cooper Basin has been generated because of an increase in mean stress associated with an increase in horizontal compressive stress since Late Cretaceous times. Structural observations and contemporary stress data indicate that there has been an increase in mean stress of approximately 50 MPa between Late Cretaceous times to that presently measured at 3780 m. The largest measured overpressure in the Cooper Basin is 14.5 MPa at 3780 m in the Kirby 1 well. Hence, disequilibrium compaction driven by increasing mean stress can explain the magnitude of the observed overpressure in the Cooper Basin. Increases in mean stress (tectonic loading) may be a feasible mechanism for overpressure generation in other ‘old’ basins that have undergone a recent increase in horizontal stress (e.g. Anadarko Basin).     

Overpressure and petroleum generation and accumulation in the Dongying Depression of the Bohaiwan Basin, China

The occurrence of abnormally high formation pressures in the Dongying Depression of the Bohaiwan Basin, a prolific oil‐producing province in China, is controlled by rapid sedimentation and the distribution of centres of active petroleum generation. Abnormally high pressures, demonstrated by drill stem test (DST) and well log data, occur in the third and fourth members (Es3 and Es4) of the Eocene Shahejie Formation. Pressure gradients in these members commonly fall in the range 0.012–0.016 MPa m^?1, although gradients as high as 0.018 MPa m^?1 have been encountered. The zone of strongest overpressuring coincides with the areas in the central basin where the principal lacustrine source rocks, which comprise types I and II kerogen and have a high organic carbon content (>2%, ranging to 7.3%), are actively generating petroleum at the present day. The magnitude of overpressuring is related not only to the burial depth of the source rocks, but to the types of kerogen they contain. In the central basin, the pressure gradient within submember Es32, which contains predominantly type II kerogen, falls in the range 0.013–0.014 MPa m^?1. Larger gradients of 0.014–0.016 MPa m^?1 occur in submember Es33 and member Es4, which contain mixed type I and II kerogen. Numerical modelling indicates that, although overpressures are influenced by hydrocarbon generation, the primary control on overpressure in the basin comes from the effects of sediment compaction disequilibrium. A large number of oil pools have been discovered in the domes and faulted anticlines of the normally pressured strata overlying the overpressured sediments; the results of this study suggest that isolated sandstone reservoirs within the overpressured zone itself offer significant hydrocarbon potential.     

Chemical and mineralogical characterization of historic mortars from the Santa Eulalia de Bóveda temple,NW Spain

Bulk samples of mortars from the Santa Eulalia de Bóveda temple (NW Spain) have been characterized from a chemical and mineralogical point of view, with the aims of contributing technological and provenance inferences, and evaluating possible correlations with the various identified stages along historical periods of construction.     

The lower jaw and dentition of the traversodontid Exaeretodon riograndensis Abdala,Barberena & Dornelles,from the Brazilian Triassic (Santa Maria 2 Sequence,Hyperodapedon Assemblage Zone)

Exaeretodon riograndensis Abdala, Barberena, & Dornelles, 2002 Abdala, F., Barberena, M.C. and Dornelles, J.E.F. 2002. A new species of the traversodontid cynodont Exaeretodon from the Santa Maria Formation (Middle/Late Triassic) of Southern Brazil. Journal of Vertebrate Paleontology, 22: 313–325. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] is the most abundant traversodontid preserved in the basal Santa Maria 2 Sequence (Hyperodapedon Assemblage Zone), southern Brazil and is closely related to Exaeretodon argentinus Cabrera, 1943 Cabrera, A. 1943. El primer hallazgo de terápsidos en la Argentina. Notas del Museo de La Plata, Paleontología, 8: 317–331.  [Google Scholar] from the Ischigualasto Formation, Argentina. Previous anatomical studies of E. riograndensis have focused mainly on cranial material and little is known about the morphology of its lower dentition or postdentary bones. We describe the first fairly complete postdentary series of the mandible of E. riograndensis and provide additional information on its lower dentition. The postdentary bones of E. riograndensis include a complex coronoid, an angular with a delicate reflected lamina and a stout retroarticular process of the articular, contrasting with the morphology reconstructed for Argentinean specimens, which possess a small retroarticular process. Apart from that, the postdentary bones do not differ significantly from those known for E. argentinus, a fact expected due to the great similarity between other skeletal features of these species. Furthermore, the lower postcanines of E. riograndensis have virtually the same structure as those of E. argentinus, with an approximate quadrangular shape in occlusal view. Moreover, the transverse cusp row is placed anteriorly and comprises a lingual and a buccal cusp, and the occlusal basin delimited by the four main cusps is relatively deep. The new material does not add any taxonomically diagnostic features to E. riograndensis. However, the fossils greatly improve our understanding of the anatomy of the Brazilian species.     

Linked thermal convection of the basement and basinal fluids in formation of the unconformity‐related uranium deposits in the Athabasca Basin,Saskatchewan, Canada

World‐class unconformity‐related U deposits in the Athabasca Basin (Saskatchewan, Canada) are generally located within or near fault zones that intersect the unconformity between the Athabasca Group sedimentary basin rocks and underlying metamorphic basement rocks. Two distinct subtypes of unconformity‐related uranium deposits have been identified: those hosted primarily in the Athabasca Group sandstones (sediment‐hosted) and those hosted primarily in the underlying basement rocks (basement‐hosted). Although significant research on these deposits has been carried out, certain aspects of their formation are still under discussion, one of the main issues being the fluid flow mechanisms responsible for uranium mineralization. The intriguing feature of this problem is that sediment‐hosted and basement‐hosted deposits are characterized by oppositely directed vectors of fluid flow via associated fault zones. Sediment‐hosted deposits formed via upward flow of basement fluids, basement‐hosted deposits via downward flow of basinal fluids. We have hypothesized that such flow patterns are indicative of the fluid flow self‐organization in fault‐bounded thermal convection (Transport in Porous Media, 110, 2015, 25). To explore this hypothesis, we constructed a simplified hydrogeologic model with fault‐bounded thermal convection of fluids in the faulted basement linked with fluid circulation in the overlying fault‐free sandstone horizon. Based on this model, a series of numerical experiments was carried out to simulate the hypothesized fluid flow patterns. The results obtained are in reasonable agreement with the concept of fault‐bounded convection cells as an explanation of focused upflow and downflow across the basement/sandstone unconformity. We then discuss application of the model to another debated problem, the uranium source for the ore‐forming basinal brines.     

Evolution of deformation and fault‐related fluid flow within an ancient methane seep system (Eocene,Varna, Bulgaria)

Faults are often important in fuelling methane seep systems; however, little is known on how different components in fault zones control subsurface fluid circulation paths and how they evolve through time. This study provides insight into fault‐related fluid flow systems that operated in the shallow subsurface of an ancient methane seep system. The Pobiti Kamani area (NE Bulgaria) encloses a well‐exposed, fault‐related seep system in unconsolidated Lower Eocene sandy deposits of the Dikilitash Formation. The Beloslav quarry and Beloslav N faults displace the Dikilitash Formation and are typified by broad, up to 80 m wide, preferentially lithified hanging wall damage zones, crosscut by deformation bands and deformation band zones, smaller slip planes and fault‐related joints. The formation of a shallow plumbing system and chimney‐like concretions in the Dikilitash Formation was followed by at least two phases of fault‐related methane fluid migration. Widespread fluid circulation through the Dikilitash sands caused massive cementation of the entire damage zones in the fault hanging walls. During this phase, paths of ascending methane fluids were locally obstructed by decimetre‐thick, continuous deformation band zones that developed in the partly lithified sands upon the onset of deformation. Once the entire damage zone was pervasively cemented, deformation proceeded through the formation of slip planes and joints. This created a new network of more localized conduits in close vicinity to the main fault plane and around through‐going slip planes. ¹³C‐depleted crustiform calcite cements in several joints record the last phase of focused methane fluid ascent. Their formation predated Neogene uplift and later meteoric water infiltration along the joint network. This illustrates how fault‐related fluid pathways evolved, over time, from ‘plumes’ in unconsolidated sediments above damage zones, leading to chimney fields, over widespread fluid paths, deflected by early deformation structures, to localized paths along fracture networks near the main fault.     

Basin‐Wide Sediment Grain‐Size Numerical Analysis and Paleo‐Climate Interpretation in the Shiyang River Drainage Basin

Basin‐wide sediment transport affects estimates of basin sediment yield, which is a fundamental scientific issue in drainage basin studies. Many studies have been conducted to examine erosion and deposition rates in drainage networks. In this study, we proposed a new approach using grain‐size standard deviation model of sedimentary samples from different geomorphological units for numerical analysis and paleo‐climate interpretation in the Shiyang River drainage basin, arid China. 1043 sedimentary samples were obtained from the upper reaches, the midstream alluvial plain and the terminal lake area; chronological frames were established based on 58 radiocarbon ages. Grain‐size standard deviation model was introduced to examine sediment components according to grain‐size and transport forces. In addition, transient paleo‐climate simulations, including the Community Climate System Model version 3 and the Kiel models, were synthesized, as well as the results from PMIP 3.0 project, to detect the long‐term climate backgrounds. Totally, we found four major common components, including fine particulates (<2 μm), fine silt (2–20 μm), sandy silt (20–200 μm), coarse sand (>200 μm), from basin‐wide sedimentary samples. The fine particulates and fine silt components exist in all the sedimentary facies, showing long‐term airborne aerosol changes and its transport by suspended load. There are some differences in ranges of sandy silt and coarse sand components, due to lake and river hydrodynamics, as well as the distance with the Gobi Desert. Paleo‐climate simulations have shown that the strong Asian summer monsoon during the transition of the Last Deglaciation and Holocene was conducive to erosion and transport of basin‐wide suspended load, also enhancing sediment sorting effects due to strong lake hydrodynamics. Our findings provide a new approach in research of long‐term basin‐wide sediment transport processes.     

Impacts of Pleistocene glaciation on large‐scale groundwater flow and salinity in the Michigan Basin

Pleistocene melting of kilometer‐thick continental ice sheets significantly impacted regional‐scale groundwater flow in the low‐lying stable interiors of the North American and Eurasian cratons. Glacial meltwaters penetrated hundreds of meters into the underlying sedimentary basins and fractured crystalline bedrock, disrupting relatively stagnant saline fluids and creating a strong disequilibrium pattern in fluid salinity. To constrain the impact of continental glaciation on variable density fluid flow, heat and solute transport in the Michigan Basin, we constructed a transient two‐dimensional finite‐element model of the northern half of the basin and imposed modern versus Pleistocene recharge conditions. The sag‐type basin contains up to approximately 5 km of Paleozoic strata (carbonates, siliciclastics, and bedded evaporites) overlain by a thick veneer (up to 300 m) of glacial deposits. Formation water salinity increases exponentially from <0.5 g l^?1 total dissolved solids (TDS) near the surface to >350 g l^?1 TDS at over 800 m depth. Model simulations show that modern groundwater flow is primarily restricted to shallow glacial drift aquifers with discharge to the Great Lakes. During the Pleistocene, however, high hydraulic heads from melting of the Laurentide Ice Sheet reversed regional flow patterns and focused recharge into Paleozoic carbonate and siliciclastic aquifers. Dilute waters (<20 g l^?1 TDS) migrated approximately 110 km laterally into the Devonian carbonate aquifers, significantly depressing the freshwater‐saline water mixing zones. These results are consistent with ¹⁴C ages and oxygen isotope values of confined groundwaters in Devonian carbonates along the basin margin, which reflect past recharge beneath the Laurentide Ice Sheet (14–50 ka). Constraining the paleohydrology of glaciated sedimentary basins, such as the Michigan Basin, is important for determining the source and residence times of groundwater resources, in addition to resolving geologic forces responsible for basinal‐scale fluid and solute migration.     

Impacts of Pleistocene glacial loading on abnormal pore‐water pressure in the eastern Michigan Basin

The hydromechanical effects of Pleistocene glacial loading on the Michigan Basin are assessed using numerical analysis based on coupled stress and pore‐water pressure. The two‐dimensional model domain included the Basin cross section and extended 10 km into the Precambrian. In the analysis, we considered the effects of the number of glacial loading cycles, the presence and connectedness of a deep Cambrian aquifer, the direction of glacial advance, the effect of a wet versus dry glacier/soil interface, topographic effects, density‐driven flow effects, and lithosphere flexure on the development of anomalous pressures. The modeling results were compared with data collected from highly instrumented wells completed in the eastern margin of the Basin. The present‐day results define regions of significant underpressure in the upper Ordovician and lower Silurian formations characterized by very low hydraulic conductivity and regions of overpressure where hydraulic conductivity is higher. To achieve the degree of underpressure observed in the instrumented wells using the model, a specific loading cycle applied over 100 000 years was repeated four times. As the number of loading cycles increased, the Paleozoic formations reached a state where the underpressures remain constant. Our results also illustrate the difference in poroelastic modeling between the application of mechanical loads on the land surface and the application of an equivalent hydraulic head, where the latter developed overpressures rather than the observed underpressures. The modeling also shows that the overpressures observed in the Cambrian formations are most likely to be the result of density‐driven flow. Finally, the simulations show that the effects of lithosphere flexure in the hydromechanical models results in the development of lateral stresses that generate overpressures rather than underpressures in the southern half of the domain. As there are no suitable observation points, these results remain unconfirmed, and further study is warranted.     

Archaeological survey,paleogeography, and the search for Late Pleistocene Paleocoastal peoples of Santa Rosa Island,California

Abstract

The northern Pacific Coast is an important area for understanding human colonization of the Americas, but Late Pleistocene coastal sites are rare and interglacial sea level rise has inundated the continental shelf and the primary areas where Paleocoastal archaeological sites are likely to occur. Here we outline a terrestrial archaeological survey project designed to identify Paleocoastal sites on Santa Rosa Island, California. Using reconstructions of ancient shorelines and paleogeography, we predicted that sites might be found where lithic resources, freshwater springs, caves or rockshelters, and strategic vistas drew Paleocoastal peoples into the island interior. We identified nine new Paleocoastal sites, including four radiocarbon dated to >11,000 cal b.p. that are among the oldest sites on North America’s Pacific coast. Our targeted survey demonstrates an important technique for investigating island and coastal settings where sea level rise remains a significant challenge for locating early sites.     

Phase‐field modeling of epitaxial growth of polycrystalline quartz veins in hydrothermal experiments

Mineral precipitation in an open fracture plays a crucial role in the evolution of fracture permeability in rocks, and the microstructural development and precipitation rates are closely linked to fluid composition, the kind of host rock as well as temperature and pressure. In this study, we develop a continuum thermodynamic model to understand polycrystalline growth of quartz aggregates from the rock surface. The adapted multiphase‐field model takes into consideration both the absolute growth rate as a function of the driving force of the reaction (free energy differences between solid and liquid phases), and the equilibrium crystal shape (Wulff shape). In addition, we realize the anisotropic shape of the quartz crystal by introducing relative growth rates of the facets. The missing parameters of the model, including surface energy and relative growth rates, are determined by detailed analysis of the crystal shapes and crystallographic orientation of polycrystalline quartz aggregates in veins synthesized in previous hydrothermal experiments. The growth simulations were carried out for a single crystal and for grain aggregates from a rock surface. The single crystal simulation reveals the importance of crystal facetting on the growth rate; for example, growth velocity in the c‐axis direction drops by a factor of ~9 when the faceting is complete. The textures produced by the polycrystal simulations are similar to those observed in the hydrothermal experiments, including the number of surviving grains and crystallographic preferred orientations as a function of the distance from the rock wall. Our model and the methods to define its parameters provide a basis for further investigation of fracture sealing under varying conditions.     

Changes in chemical composition caused by water–rock interactions across a strike‐slip fault zone: case study of the Atera Fault,Central Japan

In some cases, water–rock interactions in fault zones can affect radionuclide migration. Here, we analyzed the chemical compositions of well‐exposed fault rocks from the strike‐slip Atera Fault, Central Japan, in order to understand the variability and behavior of major and selected trace elements. The fault zone has a 1.2‐m‐wide, smectite‐rich fault core and paired damage zones that developed within welded tuff on one side of the core and within granite on the other side. The 30‐cm‐wide, kaolinite‐rich fault gouge is developed in granite cataclasite, and it shows indications of the latest fault activity, while the 1.2‐m‐wide fault core appears to be older. Hydrogen and oxygen isotope ratios in the clay‐rich fault gouges, and carbon and oxygen isotope ratios in carbonates indicate that the two major clay‐rich zones formed in bedrock near the surface, consistent with observed deformation structures. Based on chemical analyses, we identified (1) a slight depletion in SiO₂, Na₂O, K₂O, and light rare earth elements at the edges of the 1.2‐m‐wide fault core, (2) a clear depletion in SiO₂, Na₂O, K₂O, and all rare earth elements except Eu in the 30‐cm‐wide fault gouge, and (3) an increase in CaO, MnO, and heavy rare earth elements across the entire 1.2‐m‐wide fault core. Findings (1) and (2) reflect water–rock interactions in the 1.2‐m‐wide fault core and in the 30‐cm‐wide fault gouge that resulted in the formation of smectite and kaolinite. Finding (3) reflects carbonate precipitation caused by the addition of basalt fragments from a nearby site to the 1.2‐m‐wide fault core during faulting, and subsequent sorption reactions of heavy rare earth elements via processes such as complexation with the carbonates.     

Fluid inclusion and isotopic constraints on the origin of ore‐forming fluid of the Jinman–Liancheng vein Cu deposit in the Lanping Basin,western Yunnan,China

Extensive quartz–carbonate–Cu sulfide veins occur in clastic rocks and are spatially related to Paleocene granites in the western border of the Lanping Basin, western Yunnan, China. Abundant aqueous‐carbonic fluid inclusions occur in these veins but their origin is debated. In the Jinman–Liancheng deposit, individual primary inclusion groups contain either exclusively liquid‐rich inclusions (Gl), or coexisting liquid‐rich and vapor‐rich inclusions (Glv). Microthermometry and estimate of CO₂ content indicate that type Gl inclusions either have homogenization temperatures (Th) 238–263°C and contain c. 3.9–5.5 mole % CO₂, or have Th 178–222°C and contain c. 1.6–3.2 mole % CO₂. Type Glv inclusions are thought to represent samples of fluid unmixing that occurred at 183–218°C. At that time, the liquid phase in the unmixing fluid may contain c. 2.0–3.3 mole % CO₂. As such, the correlation of CO₂ content with Th for type Gl inclusions is thought to be caused by fluid unmixing with decreasing temperature and subsequent CO₂ escape. δ¹⁸O and δD values of the parent water mainly fall in the field below that of primary magmatic water, indicative of fluid derivation from degassed (in open system) magmatic water, with no contributions from basinal or meteoric water. Initial Sr isotopic compositions of hydrothermal carbonates suggest that the fluid was magmatic, probably derived from the Paleogene granites. δ¹³C_PDB values (?4‰ to ?7‰) of the hydrothermal carbonates and δ³⁴S_V‐CDT values of sulfides (mainly ?11‰ to +5‰) indicate that the carbon and sulfur can be derived from (degassed) magma and/or nonmagmatic sources. The CO₂‐rich and magmatic‐water‐derived fluid at Jinman–Liancheng differs from the CO₂‐poor and basinally derived fluid in sediment‐hosted stratiform Cu (SSC) deposits, which suggests that there are no genetic linkages between the vein Cu and SSC deposits in the Lanping Basin.     

Palaeo‐formation water evolution in the Latrobe aquifer,Gippsland Basin,south‐eastern Australia continental shelf

Offshore fresh or brackish groundwater has been observed around the globe and represents an interesting but unusual freshwater reserve. Formation waters in sedimentary basins evolve at geological time through fluid–rock interactions and water movements in aquifers. However, the mechanism and timing of freshwater displacing and mixing with pre‐existing formation water offshore under the seafloor has not been investigated in many cases. The growing need for developing freshwater resources in deeper parts of sedimentary basins that have not been economic or technically feasible in the past, may potentially lead to an increasing conflict with petroleum production or injection of carbon dioxide. For being able to assess and mitigate possible impacts of fluid production or injection on groundwater flow and quality, a better understanding of the natural history of the interaction between fresh meteoric water and deep basin formation water is necessary. A low‐salinity wedge of meteoric origin with less than 5000 ppm currently extends to about 20 km offshore in the confined Latrobe aquifer in the Gippsland Basin (Australia). The Latrobe aquifer is a freshwater resource in the onshore, hosts major petroleum reservoirs and has been considered for carbon dioxide storage in the offshore parts of the basin. The objective of this study is to constrain the evolution of formation water in the Latrobe aquifer by investigating the water naturally trapped in fluid inclusions during burial. The measured palaeo‐salinities from onshore and offshore rock samples have a minimum of about 12 500 ppm (NaCl equivalent) and a maximum of about 50 000 ppm. Most of the salinities are in the 32 000–35 000 ppm range. There is no evidence for freshwater in fluid inclusions and the variation in palaeo‐salinity across the basin is consistent with the palaeogeography of deposition of the sedimentary rocks. The current low‐salinity water wedge must have started to form recently after most of the diagenetic processes that led to the trapping of water in fluid inclusions happened. The minimum homogenisation temperatures (T_h) recorded are consistent with current formation temperature. However, they are generally higher than present day suggesting that hotter temperatures were attained in the past. The T_h and salinity data together suggest that the fluid inclusions record the diagenetic modification of connate water to higher salinities over a time period that was accompanied by an increase in temperature, consistent with a westward palaeo‐fluid flow from the deeper part of the basin through the aquifer. Subsequent pore‐water evolution from palaeo‐ to current day conditions is consistent with an influx of fresher and cooler meteoric water into the Latrobe Group. The meteoric recharge originates from the area of the Baragwanath anticline in the onshore part of the basin where the Latrobe Group subcrops at high elevations.