Comparing various approaches for assessing groundwater recharge at a regional scale in the Canadian Shield

ABSTRACT

Estimating groundwater recharge is crucial to ensuring the proper management of aquifers. In this study, net regional recharge and spatial potential recharge are estimated at four watersheds within the Charlevoix–Haute-Côte-Nord (CHCN) regions, Quebec Province, Canada. Four methods are applied based on available data. The first two approaches are regional water budget methods. These two methods differ in their estimation of vertical inflow (VI), which is estimated from two hydrological models: GR4J and HYDROTEL. The third method estimates potential recharge spatially over the study area. Finally, the streamflow data are analysed using the Eckhardt baseflow separation method to obtain an estimation of recharge, assuming that discharge is equal to recharge. According to the results of all investigated methods, the mean annual recharge for the CHCN region is approximately 183 mm, which is 18% of the total annual precipitation (P). The discussion section highlights uncertainties due to the assumptions of each method and the reliability of the data.     

The saturation of calcite and aragonite in the Arctic Ocean

We report on the chemical saturation of CaCO₃ in the waters of the Arctic Ocean calculated from total alkalinity (A_T) and total dissolved inorganic carbon (C_T). Data based on four different expeditions are presented: International Arctic Ocean Expedition (IAOE-91), Arctic Ocean Section 94 (AOS94), Polarstern Arctic '96 expedition (ACSYS 96), and Joint Ocean Ice Study 97 (JOIS 97). The results show a lysocline at around 3500 m for aragonite and that most of the Arctic Ocean sea floor lies above the lysocline for calcite. The only anomaly is the low degree of saturation at the shelf break depth in the Canadian Basin seen in the sections of the AOS94 and JOIS 97 cruises, correlated with nutrient maxima and very low O₂ concentration, suggesting decomposition of organic matter. The insignificant variability in degree of saturation between the deep waters of the different basins in the Arctic Ocean indicates a very low sedimentation/remineralisation of organic soft matter.     

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.     

Low 210Pb in the upper thermocline in the Canadian Basin: scavenge process over the Chukchi Sea

210Pb was measured during the 3rd Chinese National Arctic Research Expedition cruise to investigate its spatial pattern in the western Arctic Ocean, as well as its relation with the thermocline in the Canadian Basin. The specific activities varied from 0.04 to 2.72, 0.013 to 4.37, and 0.1 to 4.85 Bq/m3 for dissolved, particulate, and bulk 210Pb, respectively, corresponding to respective averages of 0.65, 0.43, and 1.08 Bq/m3. In the Canadian Basin, the minimum 210Pb activities occurred in the thermocline, which was characterized by low temperature of-1.52°C and salinity of 33.1. Combining the spatial distribution of 210Pb and hydrographical characteristics in the western Arctic Ocean, this scenario was ascribed to the effective scavenging of 210Pb when the Pacific water flowed across the Chukchi Shelf. Quantitatively, this interpretation was supported by both the shorter residence times and higher scavenging efficiencies(SE) of dissolved 210Pb over the Chukchi Shelf. The highest SE values were observed in the Herald Shoal and bottom waters over the slope.     

Quantifying sediment transport across an undisturbed prairie landscape using cesium-137 and high resolution topography

Soil erosion is a global environmental problem, and anthropogenic fallout radionuclides offer a promising tool for describing and quantifying soil redistribution on decadal time scales. To date, applications of radioactive fallout to trace upland sediment transport have been developed primarily on lands disturbed by agriculture, grazing, and logging. Here we use ¹³⁷Cs to characterize and quantify soil erosion at the Konza Prairie Long-Term Ecological Research (LTER) site, an undisturbed grassland in northeastern Kansas. We report on the small scale (< 10 m) and landscape scale (10 to 1000 m) distribution of fallout ¹³⁷Cs, and show significant variability in the concentrations and amounts of ¹³⁷Cs in soils at our site. ¹³⁷Cs soil concentrations and amounts typically vary by 10% to 30% on small scales, which most likely represents the spatial heterogeneity of the depositional processes. Landscape scale variability of soil ¹³⁷Cs was significantly higher than small scale variability. Most notably, soils collected on convex (divergent) landforms had ¹³⁷Cs inventories of 2500 to 3000 Bq m^− 2, which is consistent with the expected atmospheric inputs to the study area during the 1950s and 1960s. Concave landforms, however, had statistically lower inventories of 1800 to 2300 Bq m^− 2. The distribution of ¹³⁷Cs on this undisturbed landscape contrasts significantly with distributions observed across disturbed sites, which generally have accumulations of radioactive fallout in valley bottoms. Because the upslope contributing area at each sampling point had a significant negative correlation with the soil inventory of ¹³⁷Cs, we suggest that overland flow in convergent areas dominates soil erosion at Konza on time scales of decades. Very few points on our landscape had ¹³⁷Cs inventories significantly above that which would be predicted from direct deposition of ¹³⁷Cs on the soil surface; we conclude therefore that there is little net sediment storage on this undisturbed landscape.     

Assessment Methodology and Results From the Canadian Gas Potential Committee 2001 Report

Since 1991 volunteers from the Canadian Gas Potential Committee (CGPC) have conducted assessments of undiscovered gas potential in Canada. Reports were published in 1997 and 2001. The 2001 CGPC report assessed all established and some conceptual exploration plays in Canada and incorporated data from about 29,000 discovered gas pools and gas fields. Mainly year-end 1998 data were used in the analysis of 107 established exploration plays. The CGPC assessed gas in place without using economic cut offs. Estimates of nominal marketable gas were made, based on the ratio between gas in place and marketable gas in discovered pools. Only part of the estimated nominal marketable gas actually will be available, primarily because of restrictions on access to exploration and the small size of many accumulations. Most plays were assessed using the Petrimes program where it could be applied. Arps-Roberts assessments were made on plays where too many discovered pools were present to use the Petrimes program. Arps-Roberts assessments were corrected for economic truncation of the discovered pool sample. Several methods for making such corrections were tried and examples of the results are shown and compared with results from Petrimes. In addition to assessments of established plays, 12 conceptual plays, where no discoveries have been made, were assessed using Petrimes subjective methodology. An additional 65 conceptual plays were recognized, discussed, and ranked without making a quantitative assessment. No nominal marketable gas was attributed to conceptual plays because of the high risk of failure in such plays. Nonconventional gas in the form of coalbed methane, gas hydrates, tight gas, and shale gas are discussed, but no nominal marketable gas is attributed to those sources pending successful completion of pilot study projects designed to demonstrate commercially viable production. Conventional gas resources in Canada include 340 Tcf of gas in place in discovered pools and fields and 252 Tcf of undiscovered gas in place. Remaining nominal marketable gas includes 96 Tcf in discovered pools and fields and 138 Tcf of undiscovered nominal marketable gas. The Western Canada Sedimentary Basin holds 61% of the remaining nominal marketable gas. Future discoveries from that area will be mainly in pools smaller than 2.5 Bcf of marketable gas and increasing levels of exploratory drilling will be required to harvest this undiscovered resource. A pragmatic, geologically focussed approach to the assessment of undiscovered gas potential by the CGPC provides a sound basis for future exploration and development planning. Peer reviewed assessment on a play-by-play basis for entire basins provides both detailed play information and the ability to evaluate new exploration results and their impact on overall potential.     

Relict nebkhas (pimple mounds) record prolonged late Holocene drought in the forested region of south-central United States

The origin and significance of pimple mounds (low, elliptical to circular dune-like features found across much of the south-central United States) have been debated for nearly two centuries. We cored pimple mounds at four sites spanning the Ozark Plateau, Arkansas River Valley, and Gulf of Mexico Coastal Plain and found that these mounds have a regionally consistent textural asymmetry such that there is a significant excess of coarse-grained sediment within their northwest flanks. We interpret this asymmetry as evidence of an eolian depositional origin of these mounds and conclude they are relict nebkhas (coppice dunes) deposited during protracted middle to late Holocene droughts. These four mounds yield optically stimulated luminescence ages between 2400 and 700 yr that correlate with well-documented periods of eolian activity and droughts on the southern Great Plains, including the Medieval Climate Anomaly. We conclude vegetation loss during extended droughts led to local eolian deflation and pimple mound deposition. These mounds reflect landscape response to multi-decadal droughts for the south-central U.S. The spatial extent of pimple mounds across this region further underscores the severity and duration of late Holocene droughts, which were significantly greater than historic droughts.     

A densely feathered ornithomimid (Dinosauria: Theropoda) from the Upper Cretaceous Dinosaur Park Formation,Alberta, Canada

A recently discovered articulated partial skeleton of Ornithomimus from the Upper Cretaceous Dinosaur Park Formation of Alberta, Canada is remarkable in the extent and quality of preservation of integumentary structures including feathers. It is the first ornithomimid to preserve a tail bearing extensive plumaceous feathers that are slightly more elongate in comparison to those present on the remainder of the body. However, the underside of the tail and the hind limb distal to the middle of the femur appear devoid of plumage. Overall, the plumage pattern in Ornithomimus is similar to that of Struthio camelus (ostrich) and other large palaeognaths, indicating a probable function in thermoregulation. The specimen also preserves the body outline around the legs, including a skin contour anterior to the femur, analogous to skin webs in extant birds. Whereas the knee web of birds bridges the knee to the abdomen, in Ornithomimus it spans from the mid-femoral shaft to the abdomen, and is herein referred to as an anterior femoral web. This is the first report of such soft tissue structures in non-avian theropods. It may indicate that the resting position of the femur was positioned more anteroventrally in ornithomimids than in most theropods, and in that sense may have been transitional to the situation in modern birds.     

The thermal structure and thickness of continental roots

We compare heat flow data from the Precambrian shields in North America and in South Africa. We also review data available in other less well-sampled Shield regions. Variations in crustal heat production account for most of the variability of the heat flow. Because of this variability, it is difficult to define a single average crustal model representative of a whole tectonic province. The average heat flow values of different Archean provinces in Canada, South Africa, Australia and India differ by significant amounts. This is also true for Proterozoic provinces. For example, the heat flow is significantly higher in the Proterozoic Namaqua–Natal Belt of South Africa than in the Grenville Province of the Canadian Shield (61 vs. 41 mW m⁻² on average). These observations indicate that it is not possible to define single value of the average heat flow for all provinces of the same crustal age. Large amplitude short wavelength variations of the heat flow suggest that most of the difference between Proterozoic and Archean heat flow is of crustal origin. In eastern Canada, there is no good correlation between the local values of heat flow and heat production. In the Archean, Proterozoic and Paleozoic provinces of eastern Canada, heat flow values through rocks with the same heat production are not significantly different. There is therefore no evidence for variations of the mantle heat flow beneath these different provinces. After removing the local crustal heat production from the surface heat flow, the mantle (Moho) heat flow was estimated to be between 10–15 mW m⁻² in the Archean, Proterozoic and Paleozoic provinces of eastern Canada. Estimates of the mantle heat flow in the Kaapvaal craton of South Africa may be slightly higher (≈17 mW m⁻²). Large-scale variations of bulk crustal heat production are well-documented in Canada and imply significant differences of deep lithospheric thermal structure. In thick lithosphere, surficial heat flow measurements record a time average of heat production in the lithospheric mantle and are not in equilibrium with the instantaneous heat production. The low mantle heat flow and current estimates of heat production in the lithospheric mantle do not support a mechanical (conductive) lithosphere thinner than 200 km and thicker than 330 km. Temperature anomalies with surrounding oceanic mantle extend to the convective boundary layer below the conductive layer, and hence to depths greater than these estimates. Mechanical and thermal stability of the lithosphere require the mantle part of the lithosphere to be chemically buoyant and depleted in radiogenic elements. Both characteristics are achieved simultaneously by partial melting and melt extraction.