Heat flow and ground water flow in the Great Plains of the United States

Regional groundwater flow in deep aquifers adds advective components to the surface heat flow over extensive areas within the Great Plains province. The regional groundwater flow is driven by topographically controlled piezometric surfaces for confined aquifers that recharge either at high elevations on the western edge of the province or from subcrop contacts. The aquifers discharge at lower elevations to the east. The assymetrical geometry for the Denver and Kennedy Basins is such that the surface areas of aquifer recharge are small compared to the areas of discharge. Consequently, positive advective heat flow occurs over most of the province. The advective component of heat flow in the Denver Basin is on the order of 15 mW m⁻² along a zone about 50 km wide that parallels the structure contours of the Dakota aquifer on the eastern margin of the Basin. The advective component of heat flow in the Kennedy Basin is on the order of 20 mW m⁻² and occurs over an extensive area that coincides with the discharge areas of the Madison (Mississippian) and Dakota (Cretaceous) aquifers. Groundwater flow in Paleozoic and Mesozoic aquifers in the Williston Basin causes thermal anomalies that are seen in geothermal gradient data and in oil well temperature data. The pervasive nature of advective heat flow components in the Great Plains tends to mask the heat flow structure of the crust, and only heat flow data from holes drilled into the crystalline basement can be used for tectonic heat flow studies.     

Comparisons of borehole temperature–depth profiles and surface air temperatures in the northern plains of the USA

Temperature–depth profiles measured in boreholes contain a record of temperature changes at the Earth's surface. The degree to which these profiles and surface air temperature records track each other is quantitatively assessed for the northern plains of the USA. Surface air temperature records are used as a forcing function to generate synthetic transient temperature profiles which are compared with transient temperatures derived from borehole temperature–depth data. These comparisons indicate that surface air and ground temperatures are correlated. Furthermore, these comparisons yield a long-term mean temperature tied to the meteorological record which provides a context for interpreting contemporary warming trends. Our results indicate that warming recorded in surface air temperature time series represents a positive departure above baseline temperature estimates.     

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.