Models of the rapid post‐rift subsidence in the eastern Qiongdongnan Basin,South China Sea: implications for the development of the deep thermal anomaly

The Qiongdongnan Basin is one of the largest Cenozoic rifted basins on the northern passive margin of the South China Sea. It is well known that since the Late Miocene, approximately 10 Ma after the end of the syn‐rift phase, this basin has exhibited rapid thermal subsidence. However, detailed analysis reveals a two‐stage anomalous subsidence feature of the syn‐rift subsidence deficit and the well‐known rapid post‐rift subsidence after 10.5 Ma. Heat‐flow data show that heat flow in the central depression zone is 70–105 mW m^?2, considerably higher than the heat flow (<70 mW m^?2) on the northern shelf. In particular, there is a NE‐trending high heat‐flow zone of >85 mW m^?2 in the eastern basin. We used a numerical model of coupled geothermal processes, lithosphere thinning and depositional processes to analyse the origin of the anomalous subsidence pattern. Numerical analysis of different cases shows that the stretching factor β_s based on syn‐rift sequences is less than the observed crustal stretching factor β_c, and if the lithosphere is thinned with β_c during the syn‐rift phase (before 21 Ma), the present basement depth can be predicted fairly accurately. Further analysis does not support crustal thinning after 21 Ma, which indicates that the syn‐rift subsidence is in deficit compared with the predicted subsidence with the crustal stretching factor β_c. The observed high heat flow in the central depression zone is caused by the heating of magmatic injection equivalently at approximately 3–5 Ma, which affected the eastern basin more than the western basin, and the Neogene magmatism might be fed by the deep thermal anomaly. Our results suggest that the causes of the syn‐rift subsidence deficit and rapid post‐rift subsidence might be related. The syn‐rift subsidence deficit might be caused by the dynamic support of the influx of warmer asthenosphere material and a small‐scale thermal upwelling beneath the study area, which might have been persisting for about 10 Ma during the early post‐rift phase, and the post‐rift rapid subsidence might be the result of losing the dynamic support with the decaying or moving away of the deep thermal source, and the rapid cooling of the asthenosphere. We concluded that the excess post‐rift subsidence occurs to compensate for the syn‐rift subsidence deficit, and the deep thermal anomaly might have affected the eastern Qiongdongnan Basin since the Late Oligocene.     

Isotopic and thermochronologic evidence of extremely cold lithosphere associated with a slab flattening in the Central Andes of Argentina

We present mineralogic, isotopic and thermochronologic analyses on psammopelitic and tuffaceous levels from the Bermejo and Vinchina basins – both foreland depocentres of the Central Andes of Argentina – that define a low‐temperature regime for the crust akin to a slab shallowing and flattening process. The contents of illite in illite/smectite interstratified (I/S) show a progressive illitization into the deeper parts of both basins. The distribution of I/S is compatible with theoretical simulations and predicted heat flow values of ca. 26 mW m^?2 in the 8–3.4 Ma interval for the Vinchina Basin and ca. 42 mW m^?2 since 9 Ma for the Bermejo Basin. The latter shows heat flow values that are comparable to those reported by magnetotelluric analysis (36–40 mW m^?2) in agreement with previously published heat flow calculations along the modern Andean foreland. The Rb–Sr isochrones in psammopelites (<2 μm fractions) show ages between 125 and 165 Ma, whereas the K–Ar ages decrease as the grain size is smaller (136–224 Ma for 1–2 μm, 112–159 Ma for 0.2–1 μm, 76–116 Ma for <0.2 μ and 39.3–42 Ma for <0.1 μm). These ages are significantly older than the sedimentation in the basins (ca. 16 Ma for the Vinchina Basin; U–Pb age), and can be explained by the presence of a significant amount of detrital components, mainly illite, even in the finer fractions. The preservation of detrital ages is consistent with the shallow diagenesis related to a low‐temperature regime, proposed here for the basins. Younger K–Ar ages (21.3–12 Ma) were obtained for a basal tuffaceous level. Clay mineralogy and R0 ordering in the deepest part of the Vinchina Basin, together with the evolution model of I/S with depth, suggest that the burial temperatures would have not exceeded ca. 100°C in agreement with (U–Th)/He analyses performed on apatite extracted from two tuffaceous units. Thermal indicators from both studied basins confirm the existence of a low‐temperature regime during flat subduction.     

Insulating effect of coals and organic rich shales: implications for topography-driven fluid flow, heat transport, and genesis of ore deposits in the Arkoma Basin and Ozark Plateau

ABSTRACT Sedimentary rocks rich in organic matter, such as coal and carbonaceous shales, are characterized by remarkably low thermal conductivities in the range of 0.2–1.0 W m^?1 °C^?1, lower by a factor of 2 or more than other common rock types. As a result of this natural insulating effect, temperature gradients in organic rich, fine‐grained sediments may become elevated even with a typical continental basal heat flow of 60 mW m^?2. Underlying rocks will attain higher temperatures and higher thermal maturities than would otherwise occur. A two‐dimensional finite element model of fluid flow and heat transport has been used to study the insulating effect of low thermal conductivity carbonaceous sediments in an uplifted foreland basin. Topography‐driven recharge is assumed to be the major driving force for regional groundwater flow. Our model section cuts through the Arkoma Basin to Ozark Plateau and terminates near the Missouri River, west of St. Louis. Fluid inclusions, organic maturation, and fission track evidence show that large areas of upper Cambrian rocks in southern Missouri have experienced high temperatures (100–140 °C) at shallow depths (< 1.5 km). Low thermal conductivity sediments, such as coal and organic rich mudstone were deposited over the Arkoma Basin and Ozark Plateau, as well as most of the mid‐continent of North America, during the Late Palaeozoic. Much of these Late Palaeozoic sediments were subsequently removed by erosion. Our model results are consistent with high temperatures (100–130 °C) in the groundwater discharge region at shallow depths (< 1.5 km) even with a typical continental basal heat flow of 60 mW m^?2. Higher heat energy retention in basin sediments and underlying basement rocks prior to basin‐scale fluid flow and higher rates of advective heat transport along basal aquifers owing to lower fluid viscosity (more efficient heat transport) contribute to higher temperatures in the discharge region. Thermal insulation by organic rich sediments which traps heat transported by upward fluid advection is the dominant mechanism for elevated temperatures in the discharge region. This suggests localized formation of ore deposits within a basin‐scale fluid flow system may be caused by the juxtaposition of upward fluid discharge with overlying areas of insulating organic rich sediments. The additional temperature increment contributed to underlying rocks by this insulating effect may help to explain anomalous thermal maturity of the Arkoma Basin and Ozark Plateau, reducing the need to call upon excessive burial or high basal heat flow (80–100 mW m^?2) in the past. After subsequent uplift and erosion remove the insulating carbonaceous layer, the model slowly returns to a normal geothermal gradient of about 30 °C km^?1.     

A model for temperature variations in sedimentary basins due to random radiogenic heat sources

The thermal structure of a sedimentary basin is controlled by its thermal conductivity, its boundary conditions, water flow, rate of sedimentation and erosion and radiogenic heat sources. The radiogenic heat production in the sediments is known to vary over several orders of magnitude, with the lowest values in evaporites and carbonates and the highest values in black shales. Due to a paucity of information available on the existing heat sources, this parameter can be represented with a known mean value and a Gaussian correlation structure rather than a deterministic function. In this paper, the 1-D steady-state thermal structure in a sedimentary basin has been modelled in a stochastic framework with a random radiogenic heat source, and analytical expressions for the first two moments of the temperature field have been obtained. A synthetic example has been examined to quantify the error bounds on the temperature field due to uncertainties in the radiogenic heat sources.     

Basin modelling of the Limón Back-arc Basin (Costa Rica): burial history and temperature evolution of an island arc-related basin-system

The Limón back‐arc basin belongs to the southern Central American arc‐trench system and is situated at the east coast of Costa Rica. The basin‐fill consists of Late Cretaceous to Pleistocene sedimentary rocks. A northern and a southern sub‐basin can be defined, separated by the E–W‐trending Trans Isthmic Fault System. The North Limón Basin is nearly undeformed, whereas the South Limón Basin is characterized by a fold‐and‐thrust belt. Both sub‐basins have a very similar sedimentary fill and can act as a natural laboratory for reconstructing controlling factors of arc‐related sedimentary basins as well as the influence of deformation on a basin system. Modelling focused on burial history and temperature evolution. Two‐dimensional simulations were carried out with the software PetroMod^®. The geohistory curve of the North Limón Basin is overall linear, indicating continuous subsidence. The South Limón Basin is also characterized by continuous subsidence, but rates strongly increased at the beginning of the Neogene. Despite a rapid Plio‐Pleistocene deformation of the fold‐and‐thrust belt, the present‐day temperature field is not disturbed in that area. The modelling results indicate a mean heat flow of 60 mW m^?2 for the North Limón Basin and 41 mW m^?2 for the South Limón Basin. These values are low compared with other back‐arc basins. The lower values are attributed to the following effects: (1) underlying basaltic crust, (2) the lack of an initial rift phase, (3) the low extension rates, (4) absence of volcanic activity and (5) insulation effects of a thick sediment pile. The reasons for the locally lower heat flow in the southern sub‐basin can be found in the low‐angle subduction of the Cocos Ridge. Owing to the low subduction angle, the cool fore‐arc mantle‐wedge below the island‐arc is pushed backwards increasing the cooled area.     

Geothermal Characterization of the St. Lawrence Lowlands Sedimentary Basin,Québec,Canada

Defining temperature at depth to identify geothermal resources relies on the evaluation of the Earth heat flow based on equilibrium temperature measurements as well as thermal conductivity and heat generation rate assessment. Such high-quality geothermal data can be sparse over the region of interest. This is the case of the St. Lawrence Lowlands sedimentary basin covering 20,000 km² to the south of Québec, Canada, and enclosing only three wells up to a depth of 500 m with equilibrium heat flow measurements. However, more than 250 oil and gas exploration wells have been drilled in this area, providing for this study (parce que c'est 93 sinon) 81 locations with bottom-hole temperature up to a depth of 4300 m, however, not at equilibrium. Analyzing these data with respect to the deep geothermal resource potential of this sedimentary basin requires evaluating the thermal conductivity and heat generation rate of its geological units to properly extrapolate temperature downward. This was done by compiling literature and recent thermal conductivity measurements in outcrop and core samples as well as new heat generation rate estimates from spectral gamma ray logs to establish a first thermal assessment of geological units deep down into the basin. The mean thermal conductivity of the thermal units varies from 2.5 to 6.3 W/m·K, with peak values in the basal sandstones, while the heat generation rate varies from 1.6 to 0.3 µW/m³, decreasing from the upper caprocks toward the base of the sequence. After correcting the bottom-hole temperatures for drilling disturbance with the Harrison correction and subsequently for paleoclimate variations, results indicate a mean geothermal gradient of 23.1 °C/km, varying from 14 to 40 °C/km. Evaluating the basin thermal state from oil and gas data is a significant challenge facilitated by an understanding of its thermal properties.     

Deep vs. shallow controlling factors of the crustal thermal field – insights from 3D modelling of the Beaufort‐Mackenzie Basin (Arctic Canada)

Significant lateral variations in observed temperatures in the Beaufort‐Mackenzie Basin raise the question on the temperature‐controlling factors. Based on the structural configuration of the sediments and underlying crust in the area, we calculate the steady‐state 3D conductive thermal field. Integrated data include the base of the relic permafrost layer representing the 0 °C‐isotherm, public‐domain temperature data (from 227 wells) and thermal conductivity data. For >75% of the wells the predicted temperatures deviate by <10 K from the observed temperatures, which validates the overall model setup and adopted thermal properties. One important trend reproduced by the model is a decrease in temperatures from the western to the eastern basin. While in the west, a maximum temperature of 185 °C is reached at 5000 m below sea level, in the east the maximum temperature is 138 °C. The main cause for this pattern lies in lateral variations in thermal conductivity indicating differences in the shale and sand contents of the different juxtaposed sedimentary units. North‐to‐south temperature trends reveal the superposition of deep and shallow effects. At the southern margin, where the insulating effect of the low‐conductive sediments is missing, temperatures are lowest. Farther north, where the sub‐sedimentary continental crust is thick enough to produce considerable heat and a thick pile of sediments efficiently stores heat, temperatures tend to be highest. Temperatures decrease again towards the northernmost distal parts of the basin, where thinned continental and oceanic crust produce less radiogenic heat. Wells with larger deviations of the purely conductive model from the temperature observations (>15 K at 10% of the wells) and their basin‐wide pattern of misfit tendency (too cold vs. too warm temperature predictions) point to a locally restricted coupling of heat transport to groundwater flow.     

Terrestrial heat flow and lithospheric thermal structure in the Chagan Depression of the Yingen-Ejinaqi Basin,north central China

The Chagan Depression in the Yingen-Ejinaqi Basin, located at the intersection of the Paleo-Asian Ocean and the Tethys Ocean domains is an important region to gain insights on terrestrial heat flow, lithospheric thermal structure and deep geodynamic processes. Here, we compute terrestrial heat flow values in the Chagan Depression using a large set of system steady-state temperature data from four representative wells and rock thermal conductivity. We also estimate the “thermal” lithospheric thickness, mantle heat flow, ratio of mantle heat flow to surface heat flow and Moho temperature to evaluate the regional tectonic framework and deep dynamics. The results show that the heat flow in the Chagan Depression ranges from 66.5 to 69.8 mW/m², with an average value of 68.3 ± 1.2 mW/m². The Chagan Depression is characterized by a thin “thermal” lithosphere, high mantle heat flow, and high Moho temperature, corresponding to the lithospheric thermal structure of “cold mantle and hot crust” type. We correlate the formation of the Yingen-Ejinaqi Basin to the Early Cretaceous and Cenozoic subduction of the western Pacific Plate and the Cenozoic multiple extrusions. Our results provide new insights into the thermal structure and dynamics of the lithospheric evolution in central China.     

The Characteristics and Formation of A High-Arctic Proglacial Icing

Well‐known from permafrost hydrology, icings (naled or Aufeis) are also frequently encountered at the margins of high‐latitude glaciers. The morphology of a proglacial icing at Scott Turnerbreen in the Norwegian Arctic archipelago of Svalbard is described, and the process of formation is considered in detail. Ground thermal‐regime modelling indicates an equilibrium permafrost depth of at least 200 m in the studied catchment, and it appears unlikely that groundwater contributes to icing formation. Meltwater flow through ice‐marginal drainage channels is accompanied by estimated heat fluxes of up to about 190 W m^?2, suggesting that stored meltwater may continue to percolate through thawed sub‐channel sediments when surface runoff is absent during winter. A hydraulic conductivity of 6.9 × 10^?3 m s^?1 is implied, which is consistent with other studies of glacier drainage systems. The long residence time of winter‐draining meltwater, and solute rejection by refreezing water, account for high observed concentrations of solute in interstitial water in the icing. It has often been asserted that the presence of a proglacial icing indicates that a glacier is polythermal. However, as Scott Turnerbeen is entirely non‐temperate, the presence of an icing cannot always be treated as a reliable guide to the thermal regime of a glacier.     

A model for underpressure development in a glacial valley,an example from Adventdalen,Svalbard

The underpressure observed in the glacial valley Adventdalen at Svalbard is studied numerically with a basin model and analytically with a compartment model. The pressure equation used in the basin model, which accounts for underpressure generation, is derived from mass conservation of pore fluid and solid, in addition to constitutive equations. The compartment model is derived as a similar pressure equation, which is based on a simplified representation of the basin geometry. It is used to derive analytical expressions for the underpressure (overpressure) from a series of unloading (loading) intervals. The compartment model gives a characteristic time for underpressure generation of each interval, which tells when the pressure state is transient or stationary. The transient pressure is linear in time for short‐time spans compared to the characteristic time, and then it is proportional to the weight removed from the surface. We compare different contributions to the underpressure generation and find that porosity rebound from unloading is more important than the decompression of the pore fluid during unloading and the thermal contraction of the pore fluid during cooling of the subsurface. Our modelling shows that the unloading from the last deglaciation can explain the present day underpressure. The basin model simulates the subsurface pressure resulting from erosion and unloading in addition to the fluid flow driven by the topography. Basin modelling indicates that the mountains surrounding the valley are more important for the topographic‐driven flow in the aquifer than the recharging in the neighbour valley. The compartment model turns out to be useful to estimate the orders of magnitude for system properties like seal and aquifer permeabilities and decompaction coefficients, despite its geometric simplicity. We estimate that the DeGeerdalen aquifer cannot have a permeability that is higher than 1 · 10^?18 m², as otherwise, the fluid flow in the aquifer becomes dominated by topographic‐driven flow. The upper value for the seal permeability is estimated to be 1 · 10^?20 m², as higher values preclude the generation and preservation of underpressure. The porosity rebound is estimated to be <0.1% during the last deglaciation using a decompaction coefficient α_r = 1 · 10^?9 Pa^?1.     

A New Model for Heat Flow in Extensional Basins: Estimating Radiogenic Heat Production

Radiogenic heat production (RHP) represents a significant fraction of surface heat flow, both on cratons and in sedimentary basins. RHP within continental crust—especially the upper crust—is high. RHP at any depth within the crust can be estimated as a function of crustal age. Mantle RHP, in contrast, is always low, contributing at most 1 to 2 mW/m² to total heat flow. Radiogenic heat from any noncrystalline basement that may be present also contributes to total heat flow. RHP from metamorphic rocks is similar to or slightly lower than that from their precursor sedimentary rocks. When extension of the lithosphere occurs—as for example during rifting—the radiogenic contribution of each layer of the lithosphere and noncrystalline basement diminishes in direct proportion to the degree of extension of that layer. Lithospheric RHP today is somewhat less than in the distant past, as a result of radioactive decay. In modeling, RHP can be varied through time by considering the half lives of uranium, thorium, and potassium, and the proportional contribution of each of those elements to total RHP from basement. RHP from sedimentary rocks ranges from low for most evaporites to high for some shales, especially those rich in organic matter. The contribution to total heat flow of radiogenic heat from sediments depends strongly on total sediment thickness, and thus differs through time as subsidence and basin filling occur. RHP can be high for thick clastic sections. RHP in sediments can be calculated using ordinary or spectral gamma-ray logs, or it can be estimated from the lithology.     

Thermal state of the Roer Valley Graben,part of the European Cenozoic Rift System

We performed a detailed analysis of the thermal state of the Cenozoic Roer Valley Graben, the north–western branch of the European Cenozoic Rift System, based on a new set of temperature data. We developed a numerical technique for correcting bottom hole temperatures, including an evaluation of the uncertainty of thermal parameters. Comparison with drill stem test temperatures indicated that the uncertainty in corrected bottom hole temperatures using a two‐component numerical model is approximately ± 4 °C, which is much more accurate than the up to 15 °C errors encountered in often‐used line‐source or Horner correction methods. The subsurface temperatures and the derived regional heat flow estimates of 53 ± 6 to 63 ± 6 mW m^?2 show no significant difference between the central rift and the adjacent structural highs. The absence of an elevated heat flow is attributed to the low amount of lithospheric thinning during the Cenozoic rifting phase (β=1.06–1.15). A local thermal anomaly exceeding +10 °C was found in five wells in the north–western part of the rift basin at depths of 1000–1500 m, and is most likely caused by the upward flow of fluids along faults, whereas lower temperatures in the upper 1500 m in the southern part of the rift basin could indicate cooling by topography‐driven groundwater flow. Conflicting ideas exist on the active or passive rifting mechanisms responsible for the formation of the different rift basins of European Cenozoic Rift System. The low spatial variation in heat flow found in this study suggests that the mechanism responsible for forming the Roer Valley Graben is passive rifting.     

The Astrakhan Arch of the Pricaspian basin: Geothermal analysis and modelling

The Astrakhan Arch (ASAR) region contains one of the largest sub‐salt carbonate structures of the Pricaspian salt basin (located to the northwest of the Caspian Sea), where prospects for hydrocarbon generation and accumulation in the Devonian to Carboniferous deposits are considered to be high. We evaluate the regional vertical temperature gradient within stratigraphic units based on the analysis of 34 boreholes drilled in the region. To show that the thermal gradient is altered in the vicinity of salt diapirs, we study measured temperatures in six deep boreholes. We develop a three‐dimensional geothermal model of the ASAR region constrained by temperature measurements, seismic stratigraphic and lithological data. The temperatures of the sub‐salt sediments predicted by the geothermal model range from about 100 °C to 200 °C and are consistent with the temperatures obtained from the analysis of vitrinite reflectivity and from previous two‐dimensional geothermal models. Temperature anomalies are positive in the uppermost portions of salt diapirs as well as within the salt‐withdrawal basins at the depth of 3.5 km depth and are negative beneath the diapirs. Two areas of positive temperature anomalies in the sub‐salt sediments are likely to be associated with the deep withdrawal basins above and with the general uplift of salt/sub‐salt interface in the southern part of the study region. This implies an enhancement of thermal maturity of any organically rich source rocks within these areas. The surface heat flux in the model varies laterally from about 40 to 55 mW m^?2. These variations in the heat flux are likely to be associated with structural heterogeneities of the sedimentary rocks and with the presence of salt diapirs. The results of our modelling support the hypothesis of oil and gas condensate generation in the Upper Carboniferous to Middle Devonian sediments of the ASAR region.     

The potential and actual primary productivity of southern Ontario's agro-ecosystem

There is a generally accepted assumption that the primary productivity of the natural vegetation cover of an area is the most appropriate indicator of the level of primary productivity that may be attained by agricultural systems in the same area. However, this relationship is not clearly understood. This paper assesses the potential rate of primary productivity for southern Ontario and compares these results with the equivalent values for the agro-ecosystems. Analysis is conducted at the level of the county. Climatically determined potential (natural) primary productivity values are adjusted by a soil performance index yielding a range of values from 986 g m^?2 y^?1 to 44 g m^?2 y^?1. The actual primary productivity for agriculture ranges from 3593 g m^?2 y^?1 to 515 g m^?2 y^?1. Over southern Ontario actual (agricultural) productivity exceeds potential (natural) productivity by an average of 2.6. The variations about this mean are examined. The results indicate where, and by which agricultural processes, the humanized landscape can improve the organic resource base.     

Unravelling the multi-stage burial history of the Swiss Molasse Basin: integration of apatite fission track, vitrinite reflectance and biomarker isomerisation analysis

A complex basin evolution was studied using various methods, including thermal constraints based on apatite fission‐track (AFT) analysis, vitrinite reflectance (VR) and biomarker isomerisation, in addition to a detailed analysis of the regional stratigraphic record and of the lithological properties. The study indicates that (1) given the substantial amount of data, the distinction and characterisation of successive stages of heating and burial in the same area are feasible, and (2) the three thermal indicators (AFT, VR and biomarkers) yield internally consistent thermal histories, which supports the validity of the underlying kinetic algorithms and their applicability to natural basins. All data pertaining to burial and thermal evolution were integrated in a basin model, which provides constraints on the thickness of eroded sections and on heat flow over geologic time. Three stages of basin evolution occurred in northern Switzerland. The Permo‐Carboniferous strike–slip basin was characterised by high geothermal gradients (80–100°C km^?1) and maximum temperature up to 160°C. After the erosion of a few hundreds of metres in the Permian, the post‐orogenic, epicontinental Mesozoic basin developed in Central Europe, with subsidence triggered by several stages of rifting. Geothermal gradients in northern Switzerland during Cretaceous burial were relatively high (35–40°C km^?1), and maximum temperature typically reached 75°C (top middle Jurassic) to 100°C (base Mesozoic). At least in the early Cretaceous, a stage of increased heat flow is needed to explain the observed maturity level. After erosion of 600–700 m of Cretaceous and late Jurassic strata during the Paleocene, the wedge‐shaped Molasse Foreland Basin developed. Geothermal gradients were low at this time (≤20°C km^?1). Maximum temperature of Miocene burial exceeded that of Cretaceous burial in proximal parts (<35 km from the Alpine front), but was lower in more distal parts (>45 km). Thus, maximum temperature as well as maximum burial depth ever reached in Mesozoic strata occurred at different times in different regions. Since the Miocene, 750–1050 m were eroded, a process that still continues in the proximal parts of the basin. Current average geothermal gradients in the uppermost 2500 m are elevated (32–47°C km^?1). They are due to a Quaternary increase of heat flow, most probably triggered by limited advective heat transport along Paleozoic faults in the crystalline basement.     

Stratabound Geothermal Resources in North Dakota and South Dakota

Geothermal energy resources in North Dakota and South Dakota occur as low (T < 90°C) and intermediate (T < 150°C) temperature geothermal waters in regional-scale aquifers within the Williston and Kennedy Basins. The accessible resource base is approximately 21.25 exajoules (10¹⁸ J = 1 exajoule, 10¹⁸ J ~ 10¹⁵ Btu = 1 quad) in North Dakota and 12.25 exajoules in South Dakota. Resource temperatures range from 40°C at depths of about 700 m to 150°C at 4500 m in the Williston Basin in North Dakota. In South Dakota, resource temperatures range from 44°C at a depth of 550 m near Pierre to 100°C at a depth of 2500 m in the northwestern corner. This resource assessment raises the identified accessible resource base by 31% above the previous assessments and by 310% over an earlier assessment. The large increases in the identified accessible resource bases reported in this study result from including all potential geothermal aquifers and better understanding of the thermal regime of the region. These results imply that a reassessment of stratabound geothermal resources in the United States that includes all geothermal aquifers would increase significantly the identified accessible resource base. The Williston Basin in North Dakota is characterized by conductive heat flows ranging from 43 to 68 mW m^–2 and averaging 55 mW m^–2. Comparisons of calculated and bottomhole temperatures measured in oil fields over the Nesson Anticline and the Billings Nose show temperature differences which suggest that upward groundwater flow in fractures on the westward sides of the structures slightly perturbs the otherwise conductive thermal field. The maximum heat-flow disturbance is estimated to be of the order of 10 to 20 mW m^–2. These thermal anomalies do not alter significantly the accessible geothermal resource base. Anomalous heat flow in south-central South Dakota is caused by heat advection in gravity-driven groundwater flow in regional aquifers. Heat flow is anomalously high (Q > 130 mW m^–2) in the discharge area in south-central South Dakota and anomalously low (30 mW m²) in the recharge area near the Black Hills and along the western limb of the Kennedy Basin in western South Dakota. Heat-flow disturbances are the result of vertical groundwater flow through fractures in the discharge area of the regional flow system in South Dakota are minor and may be significant only in deeply incised stream valleys. An important factor that controls the temperature of the resource in both North Dakota and South Dakota is the insulating effect of a thick (500–2000 m) layer of low thermal-conductivity shales that overlie the region. The effective thermal conductivity of the shale layer is approximately 1.2 W m^–1 K^–1 in contrast to sandstones and carbonates, which have conductivities of 2.5 to 3.5 W m^–1 K^–1. This low conductivity leads to high geothermal gradients (dT/dz > 50°C km^–1), even where heat flow has normal continental values, that is 40–60 mW m^–2. Engineering studies show that geothermal space heating using even the lowest temperature geothermal aquifers (T 40 °C) in North Dakota and South Dakota is cost effective at present economic conditions. The Inyan Kara Formation of the Dakota Group (Cretaceous) is the preferred geothermal aquifer in terms of water quality and productivity. Total dissolved solids in the Inyan Kara Formation ranges from 3,000 to more than 20,000 mg L^–1. Porosities normally are higher than 20%, and the optimum producing zones generally are thicker than 30 m. The estimated water productivity index of a productive well in the Inyan Kara Formation is 0.254179 l s^–1 Mpa^–1. Deeper formations have warmer waters, but, in general, are less permeable and have poorer water quality than the Inyan Kara.     

Anomalous thermal maturities caused by carbonaceous sediments

Sedimentary rocks such as coal and carbonaceous mudstone which contain abundant carbonaceous matter are characterized by thermal conductivity much lower than that exhibited by other common rock types, by a factor of 5–10. As a result, temperature gradients in such sediments can range up to 0.25 °Cm^-1 even under conditions of average heat flow. When such steep gradients extend over a significant sedimentary thickness, temperatures of underlying rock units are elevated, causing both organic and inorganic phases to record what seem to be anomalously high levels of thermal maturity. This carbonaceous blanket insulating effect may help to explain unusual levels of maturity observed at shallow depths in the Appalachian Basin, Michigan Basin and other regions of the world with significant carbonaceous strata.     

The interplay between deformation,erosion and sedimentation in the deep‐water Mexican Ridges foldbelt,western Gulf of Mexico basin

We analyse a regional 2D seismic section of the Mexican Ridges foldbelt (MRFB), western Gulf of Mexico, and construct excess‐area diagrams for each of the structures comprising the foldbelt to estimate shortening, the onset of folding and the degradation of the folded seafloor. From the chronostratigraphy, we derive rates of tectonic and superficial mass transport and illustrate how they change across the MRFB. The resulting tectonic transport in the MRFB is 11.8 km forming a train of twelve buckle folds above a detachment at a depth of ca. 6 s of two‐way travel time, with an average strain of ca. 10%. The fold train grew at a mean uplift rate of ca. 0.21 mm year^?1. Cross‐sectional balancing demonstrates that shortening balances the down‐slip motion of the Quetzalcoatl extensional system (QES), suggesting that horizontal compaction, volume loss and other penetrative deformation mechanisms are negligible. By assuming steady‐state denudation, we are able to distinguish sediments derived locally from sediments transported from distant sources. The constant of mass diffusivity, a parameter controlling the degradation rate, is ca. 0.42 m² year ^?1, which is characteristic of rapid, episodic, superficial mass movements. The combined sedimentation rate from both, local and distal sources is ca. 0.23 mm year ^?1. Those values are not constant; structures proximal to the continental shelf are rising rapidly and are being degraded more intensely than those in the distal part of the MRFB, where sedimentation outweighs tectonic uplift. Our results indicate deformation initiated up to 3 Myr earlier than estimated from stacking patterns. Moreover, we find deformation started synchronously during the Late Miocene throughout the MRFB and not in two episodes as the stacking relations suggest. The discrepancy can be explained by a delay in the sedimentary response to folding. During early fold growth, nearly constant thickness strata are deposited before a progressive unconformity and other converging geometries develop. The development of growth strata is fast in the folds near the QES, which are being uplifted rapidly and degraded vigorously. Under these conditions, the stratigraphic relations give only a broad estimate of the pretectonic/syntectonic limit when compared to the excess‐area method. On the other hand, the development of growth strata took twice as much time for folds near the abyssal plain, which are being uplifted at a slower rate and where degradation is less intense. Consequently, the delay takes more time, and the use of stratigraphic relations introduces an even more pronounced bias towards younger ages in the identification of the onset of folding.     

Late Archaean to Palaeoproterozoic geotherms in the Kaapvaal craton, South Africa: constraints on the thermal evolution of the Witwatersrand Basin

The c. 2.97–2.71 Ga Witwatersrand Basin located in the Kaapvaal craton of South Africa represents a remnant of a large Late Archaean sedimentary basin that hosts the world's premier gold deposit within a series of conglomerate horizons. Evidence of postdepositional gold mobility within these conglomerates associated with hydrothermal–metamorphic activity has led to speculation about the Late Archaean to Palaeoproterozoic geothermal gradients in the basin. We use surface heat flow and heat production data from rocks in the basin and its environs in order to calculate detailed temperature profiles for the central Kaapvaal craton that show that the steady state crustal geotherm during the Late Archaean and Palaeoproterozoic was relatively cool at 15–20 K km^?1. The geotherm in the upper crustal strata is also largely unaffected by substantial increases in the heat flow into the base of the crust. Consequently, regional greenschist facies metamorphism of the basin sediments could only have been achieved during a transient thermal event that advected heat into the upper crust. The most likely candidate for this is the Bushveld magmatic event at 2.06 Ga.     

Heat refraction and heat production in and around granite plutons in north-east England

Summary. Existing surface heat flow determinations in north-east England indicate a difference in thermal structure between the Alston block of the North Pennines and the coastal regions. New heat flow determinations in the Askrigg block and on the coast support the original indications of higher heat flow in the positive blocks, which are underlain by 400 Myr old granites, than in their marginal sedimentary troughs along the present coastline. Although the two blocks are geophysically and geologically similar, surface heat flow in the Alston block is 30mW m^-2 higher than it is in the Askrigg block 50km to its south. Heat flow refraction around high conductivity plutons with high heat production may significantly affect the interpretation of heat flow measurements made on a regional scale but in the present case its effect, if any, is to emphasize the heat flow difference between the two otherwise similar blocks. The chemical differences between the Wensleydale and the Weardale granites – cupolas of the granites underlying the Askrigg and Alston may be representative of the granites as a whole and the blocks respectively – difference in surface heat flows may result from a contrast in radiogenic heat production between the two granites extending to mid-crustal depths.