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.     

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.     

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.     

Assessment of the Enhanced Geothermal System Resource Base of the United States

This paper describes an assessment of the enhanced geothermal system (EGS) resource base of the conterminous United States, using constructed temperature at depth maps. The temperature at depth maps were computed from 3 to 10 km, for every km. The methodology is described. Factors included are sediment thickness, thermal conductivity variations, distribution of the radioactive heat generation and surface temperature based on several geologic models of the upper 10 km of the crust. EGS systems are extended in this paper to include coproduced geothermal energy, and geopressured resources. A table is provided that summarizes the resource base estimates for all components of the EGS geothermal resource. By far, the conduction-dominated components of EGS represent the largest component of the U.S. resource. Nonetheless, the coproduced resources and geopressured resources are large and significant targets for short and intermediate term development. There is a huge resource base between the depths of 3 and 8 km, where the temperature reaches 150–250°C. Even if only 2% of the conventional EGS resource is developed, the energy recovered would be equivalent to roughly 2,500 times the annual consumption of primary energy in the U.S. in 2006. Temperatures above 150°C at those depths are more common in the active tectonic regions of the western conterminous U.S., but are not confined to those areas. In the central and eastern U.S. there are identified areas of moderate size that are of reasonable grade and probably small areas of much higher grade than predicted by this analyses. However because of the regional (the grid size is 5′ × 5′) scale of this study such potentially promising sites remain to be identified. Several possible scenarios for EGS development are discussed. The most promising and least costly may to be developments in abandoned or shut-in oil and gas fields, where the temperatures are high enough. Because thousands of wells are already drilled in those locations, the cost of producing energy from such fields could be significantly lowered. In addition many hydrocarbon fields are producing large amounts of co-produced water, which is necessary for geothermal development. Although sustainability is not addressed in this study, the resource is so large that in at least some scenarios of development the geothermal resource is sustainable for long periods of time.     

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.     

High Potential Regions for Enhanced Geothermal Systems in Canada

Previous estimates of geothermal energy potential in Canada give an indication of available heat to be ‘farmed’ at depth. This article examines in more detail depth–temperature relationships near large population centers in western Canada, as well as remote communities in northern Canada, in order to provide a first order assessment of Enhanced Geothermal Systems (EGS) potential for electrical generation. Quantities of EGS thermal power output and electrical generation are dependent on output temperature and flow rate. We relate these potential power rates as a whole to drilling and installation cost for the doublet systems and triplet system. Results show areas with significant EGS potential in northern Alberta, northeastern British Columbia, and southern Northwest Territories related to high heat flow and thermal blanketing of thick sedimentary cover. Estimated installation costs in 2008 dollars are under 2 mln$/MW_e. We also estimate significant reductions in CO₂ emissions by conversion to geothermal electric production.     

Anomalous sea-ice reduction in the Eurasian Basin of the Arctic Ocean during summer 2010

During the summer of 2010 ice concentration in the Eurasian Basin, Arctic Ocean was unusually low. This study examines the sea-ice reduction in the Eurasian Basin using ice-based autonomous buoy systems that collect temperature and salinity of seawater under the ice along the course of buoy drift. An array of GPS drifters was deployed with 10 miles radius around an ice-based profiler, enabling the quantitative discussion for mechanical ice divergence/convergence and its contribution to the sea-ice reduction. Oceanic heat fluxes to the ice estimated using buoy motion and mixed-layer (ML) temperature suggest significant spatial difference between fluxes under first-year and multi-year ice. In the former, the ML temperature reached 0.6 K above freezing temperature, providing >60–70 W m^?2 of heat flux to the overlying ice, equivalent to about 1.5 m of ice melt over three months. In contrast, the multiyear ice region indicates nearly 40 W m^?2 at most and cumulatively produced 0.8 m ice melt. The ice concentration was found to be reduced in association with an extensive low pressure system that persisted over the central Eurasian Basin. SSM/I indicates that ice concentration was reduced by 30–40% while the low pressure persisted. The low ice concentration persisted for 30 days even after the low dissipated. It appears that the wind-forced ice divergence led to enhanced absorption of incident solar energy in the expanded areas of open water and thus to increased ice melt.     

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.     

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.     

Thermal properties of sediments in southern Israel: a comprehensive data set for heat flow and geothermal energy studies

To decipher the thermal structure of the sedimentary veneer in southern Israel, new values of thermal conductivity and porosity as well as of the radiogenic heat budget are provided. Thermal conductivity is measured for lithotypes and scaled up for geological formations. The new data are higher than most of the previously measured values, in particular for sandstones and siltstones, whose mean values are 5.0 and 2.9 W m^?1 K^?1. Mean values of the most abundant lithotypes, which are dolomites and limestones, are on the order of 4.1 and 2.7 W m^?1 K^?1, respectively. The total radiogenic heat production of the sedimentary cover varies slightly over southern Israel, due to variable lithology and total sediment thickness, yielding a maximum heat flow on the order of 4 mW m^?2 where the sedimentary section is thickest (ca. 7 km). A temperature prognosis was made by calculating temperature profiles to the top of the crystalline basement at five well locations applying the new thermal‐conductivity data set and three scenarios of surface heat flow (50, 55 and 60 mW m^?2). The calculated temperatures best match with measured drillstem‐test temperatures by using heat‐flow values close to the upper bound of range. Surface heat flow on the order of 55–60 mW m^?2 is supported by a reevaluation of an existing temperature log and the application of thermal conductivity from this study. The temperature prediction for southern Israel shows values of 100–120°C at 3500–4500 m depth, indicating a geothermal potential that can be used for heating as well as electricity production.     

Soil hydrological and erosional responses in areas of woody encroachment, pasture and woodland in semi-arid Australia

In arid and semi-arid areas, woody encroachment is the increase in density, cover, extent and/or biomass of woody plants. Woody encroachment is often associated with increased runoff and soil erosion. Hydrological and erosional responses of woody encroachment and of pastures established after management of encroachment in semi-arid Australia are not well understood. This study compared the hydrological and erosional responses across vegetation states comprising woody plant encroachment (>1200 stems ha⁻¹), recently established pastures (<23 years of age), long-established pasture (50-100 years of age) and open woodland (<330 stems ha⁻¹) in semi-arid eastern Australia. Responses were measured using rainfall simulation with intensity of 35 mm h⁻¹ for 30 min applied on 1 -m² plots. Runoff and sediment production did not differ significantly between vegetation states. Average runoff in woody encroachment was 9.0 mm h⁻¹, followed by recent pasture (8.2 mm h⁻¹), long-established pasture (5.9 mm h⁻¹) and open woodland (4.2 mm h⁻¹). Total sediment production in recent pasture was 11.6 g m⁻², followed by woody encroachment (9.0 g m⁻²), long-established pasture (7.3 g m⁻²) and open woodland (4.3 g m⁻²). Runoff and sediment production were significantly lower at one pasture site (0.9 mm h⁻¹ and 1.3 g m⁻²) where rotational grazing and minimum tillage had been implemented than in the adjacent paired woody encroachment site (10.3 mm h⁻¹and 6.5 g m⁻², respectively). This example of a pasture that had been managed to increase ground cover illustrated the effect of pasture management on reducing runoff and sediment production. Across all vegetation states, small scale runoff and sediment production were minimal or zero where total ground cover was 73% or higher.     

Investigation of the thermal history of the Delaware Basin (West Texas,USA) using carbonate clumped isotope thermometry

We utilized carbonate clumped isotope thermometry to explore the thermal history of the Delaware Basin, West Texas, USA. Carbonate wellbore cuttings from five oil/gas wells across the basin yielded clumped isotope temperatures (T(Δ₄₇)) ranging from 27°C to 307°C, interpreted to reflect a combination of initial precipitation/recrystallization temperature and solid-state C-O bond reordering during burial. Dolomite samples generally record lower apparent T(Δ₄₇)s than calcite, reflecting greater resistance to reordering in dolomite. In all five wells, clumped isotope temperatures exceed modern downhole temperature measurements, indicating higher heat flow in the past. Using modelled burial curves based on sedimentological history, we created unique time-temperature histories by linearly applying a geothermal gradient. Applying two different thermal history reordering models, we modelled the extent of solid-state C–O bond reordering to iteratively find the time-averaged best-fit geothermal gradients for each of the five wells. Results of this modelling suggest that the shallower, southwestern portion of the study area experienced higher geothermal gradients throughout the sediment history (~45°C/km) than did the deeper, southeastern portion (~32°C/km), with the northern portion experiencing intermediate geothermal gradients (~35–38°C/km). This trend is in agreement with the observed gas/oil ratios of the Delaware Basin, increasing from east to west. Furthermore, our clumped isotope temperatures agree well with previously published vitrinite reflectance data, confirming previous observations and demonstrating the utility of carbonate clumped isotope thermometry to reconstruct basin thermal histories.     

Assessing vegetation cover and biomass in restored erosion areas in Iceland using SPOT satellite data

Due to highly erodible volcanic soils and a harsh climate, livestock grazing in Iceland has led to serious soil erosion on about 40% of the country's surface. Over the last 100 years, various revegetation and restoration measures were taken on large areas distributed all over Iceland in an attempt to counteract this problem. The present research aimed to develop models for estimating percent vegetation cover (VC) and aboveground biomass (AGB) based on satellite data, as this would make it possible to assess and monitor the effectiveness of restoration measures over large areas at a fairly low cost. Models were developed based on 203 vegetation cover samples and 114 aboveground biomass samples distributed over five SPOT satellite datasets. All satellite datasets were atmospherically corrected, and digital numbers were converted into ground reflectance. Then a selection of vegetation indices (VIs) was calculated, followed by simple and multiple linear regression analysis of the relations between the field data and the calculated VIs.Best results were achieved using multiple linear regression models for both %VC and AGB. The model calibration and validation results showed that R² and RMSE values for most VIs do not vary very much. For percent VC, R² values range between 0.789 and 0.822, leading to RMSEs ranging between 15.89% and 16.72%. For AGB, R² values for low-biomass areas (AGB < 800 g/m²) range between 0.607 and 0.650, leading to RMSEs ranging between 126.08 g/m² and 136.38 g/m². The AGB model developed for all areas, including those with high biomass coverage (AGB > 800 g/m²), achieved R² values between 0.487 and 0.510, resulting in RMSEs ranging from 234 g/m² to 259.20 g/m². The models predicting percent VC generally overestimate observed low percent VC and slightly underestimate observed high percent VC. The estimation models for AGB behave in a similar way, but over- and underestimation are much more pronounced.These results show that it is possible to estimate percent VC with high accuracy based on various VIs derived from SPOT satellite data. AGB of restoration areas with low-biomass values of up to 800 g/m² can likewise be estimated with high accuracy based on various VIs derived from SPOT satellite data, whereas in the case of high biomass coverage, estimation accuracy decreases with increasing biomass values. Accordingly, percent VC can be estimated with high accuracy anywhere in Iceland, whereas AGB is much more difficult to estimate, particularly for areas with high-AGB variability.     

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.     

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.     

Arsenic contamination of lake sediments in Florida: evidence of herbicide mobility from watershed soils

Organic arsenical herbicides, which include monosodium methylarsonate (MSMA), have been applied to golf courses and lawns throughout Florida, USA, since the 1950s. These products convert rapidly to inorganic forms of arsenic (As) in soils and are mobilized readily. Leachates have been known to contaminate groundwater and surface waters, although past studies have not examined whether use of these products has led to significant As accumulation in lake sediments. We used paleolimnological methods to document the depositional history and inventories of total As in sediments and porewaters of Little Lake Jackson in Florida, which is adjacent to three golf courses. Six sediment cores, four of which were ²¹⁰Pb dated, showed porewater total As concentrations as high as 435 μg l⁻¹, and dry-sediment total As concentrations as high as 148 mg kg⁻¹. Approximately 537 kg of total As is present in >19,000 metric tons of sediment (dry mass), and an additional 18 kg of As is dissolved in 10.8 × 10⁴ m³ of porewaters. Total As content in surface sediments (mean = 47.3 mg kg⁻¹) exceeds the consensus-based sedimentary concentration for probable toxicity effects in freshwater benthic fauna. Surface and subsurface waters flow to the lake from topographically higher areas to the west, where golf courses and residential areas are located. Total As concentrations were elevated highly in monitoring wells and in a stream that flows between the golf courses and lake, but As was below detection limits in wells that were located at the distal perimeter of the golf courses. Subsurface and surface waters exit the lake towards topographically lower areas to the east. Nearly all As in sediments remains bound in the solid phase, indicating that As sedimentary profiles largely reflect depositional history. Sedimentary As concentrations are correlated strongly with aluminum and iron, which suggests that As was scavenged from lake waters during the past. Sedimentary As concentrations increased until the 1980s, then declined somewhat to the present time. Dissolved As was scavenged efficiently from the water column when hypolimnetic waters were oxygenated persistently, but after eutrophication led to a seasonally anoxic hypolimnion in the 1980s, apparently less As was co-precipitated, and more was lost to hydrological outflow. Arsenic accumulation in sediments might be common in areas where As derived from organic arsenical herbicide applications is directed by shallow water tables towards adjacent lakes.     

Late Ordovician,deep‐water gravity‐flow deposits,palaeogeography and tectonic setting,Tarim Basin,Northwest China

The Upper Ordovician in the Tarim Basin contains 5000–7000 m of siliciclastic and calciclastic deep‐water, gravity‐flow deposits. Their depositional architecture and palaeogeographical setting are documented in this investigation based on an integrated analysis of seismic, borehole and outcrop data. Six gravity‐flow depositional–palaeogeomorphological elements have been identified as follows: submarine canyon or deeply incised channels, broad and shallow erosional channels, erosional–depositional channel and levee–overbank complexes, frontal splays‐lobes and nonchannelized sheets, calciclastic lower slope fans and channel lobes or sheets, and debris‐flow complexes. Gravity‐flow deposits of the Sangtamu and Tierekeawati formations comprise a regional transgressive‐regressive megacycle, which can be further classified into six sequences bounded by unconformities and their correlative conformities. A series of incised valleys or canyons and erosional–depositional channels are identifiable along the major sequence boundaries which might have been formed as the result of global sea‐level falls. The depositional architecture of sequences varies from the upper slope to abyssal basin plain. Palaeogeographical patterns and distribution of the gravity‐flow deposits in the basin can be related to the change in tectonic setting from a passive continental margin in the Cambrian and Early to Middle Ordovician to a retroarc foreland setting in the Late Ordovician. More than 3000 m of siliciclastic submarine‐fan deposits accumulated in south‐eastern Tangguzibasi and north‐eastern Manjiaer depressions. Sedimentary units thin onto intrabasinal palaeotopographical highs of forebulge origin and thicken into backbulge depocentres. Sediments were sourced predominantly from arc terranes in the south‐east and the north‐east. Slide and mass‐transport complexes and a series of debris‐flow and turbidite deposits developed along the toes of unstable slopes on the margins of the deep‐water basins. Turbidite sandstones of channel‐fill and frontal‐splay origin and turbidite lobes comprise potential stratigraphic hydrocarbon reservoirs in the basin.     

Using basin thermal history to evaluate the role of Miocene–Pliocene flat‐slab subduction in the southern Central Andes (27° S–30° S)

Studies in both modern and ancient Cordilleran‐type orogenic systems suggest that processes associated with flat‐slab subduction control the geological and thermal history of the upper plate; however, these effects prove difficult to deconvolve from processes associated with normal subduction in an active orogenic system. We present new geochronological and thermochronological data from four depositional areas in the western Sierras Pampeanas above the Central Andean flat‐slab subduction zone between 27° S and 30° S evaluating the spatial and temporal thermal conditions of the Miocene–Pliocene foreland basin. Our results show that a relatively high late Miocene–early Pliocene geothermal gradient of 25–35 °C km⁻¹ was typical of this region. The absence of along‐strike geothermal heterogeneities, as would be expected in the case of migrating flat‐slab subduction, suggests that either the response of the upper plate to refrigeration may be delayed by several millions of years or that subduction occurred normally throughout this region through the late Miocene. Exhumation of the foreland basin occurred nearly synchronously along strike from 27 to 30° S between ca. 7 Ma and 4 Ma. We propose that coincident flat‐slab subduction facilitated this wide‐spread exhumation event. Flexural modelling coupled with geohistory analysis show that dynamic subsidence and/or uplift associated with flat‐slab subduction is not required to explain the unique deep and narrow geometry of the foreland basin in the region implying that dynamic processes were a minor component in the creation of accommodation space during Miocene–Pliocene deposition.     

Spatial variability of soil hydrophobicity after wildfires in Montana and Colorado

Increases in runoff and erosion after wildfires are often attributed to the development of hydrophobic soils. The potential for increased overland flow depends on the spatial contiguity of the hydrophobicity as well as its overall strength, but there is limited information on the spatial variability of soil hydrophobicity. We conducted spatially intensive hydrophobicity measurements in 225 m² and 1 m² plots in forested areas of Montana and Colorado burned at moderate to high severity, and in unburned control plots. Both the burned and unburned 225 m² plots contained 10–23 hydrophobic soil patches in which hydrophobicity was strongest at the surface and declined rapidly with depth. The hydrophobic patches were closer together and up to 3 times larger in the burned plots. Consequently, 19% to 76% of the burned plots were hydrophobic compared to just 11% of the unburned plots. In five of the six burned plots, the patches were not laterally connected, suggesting that in most cases Hortonian overland flow generated from hydrophobic patches will infiltrate near its point of origin. The 1 m² plots were smaller than most of the hydrophobic patches, so they did not capture the spatial characteristics of soil hydrophobicity. Characterization of the spatial variability of soil hydrophobicity should be based on measurements conducted at  1 m intervals across areas of > 100 m². Due to the patchiness of soil hydrophobicity at the 10⁰ to 10¹ meter scale, overland flow measurements in small ( 1 m²) plots may overestimate the magnitude and variability of runoff from burned catchments.     

Regional variations of heat flow differences with depth in Alberta, Canada

Summary. Differences between estimated average heat flow values for the Mesozoic and Cenozoic formations ( Q ₁) and estimated average heat flow values for the Palaeozoic formations below the erosional unconformity ( Q ₂) are calculated for the Alberta part of the western Canadian sedimentary basin. Significant heat flow differences exist for these two intervals and the map of Δ Q = Q ₁– Q ₂ shows that Q ₂ is generally greater than Q ₁ in the western and south-western part of Alberta, while in the northern part of the province Q ₂ is generally less than Q ₁. The regional variations of Δ Q are large, with standard deviation of 26 mW m⁻² and average value –13.5 mW m⁻². A regional trend of Δ Q correlates with topographic relief and the hydraulic head variations in the basin. It is shown that there is a heat flow increase with depth in water recharge areas and a decrease in heat flow with depth in the low topographic elevation water discharge areas when comparing the average heat flow in Mesozoic + Cenozoic and Palaeozoic formations.