Glacial isostatic adjustment at the Laurentide ice sheet margin: Models and observations in the Great Lakes region

Glacial Isostatic Adjustment (GIA) modelling in North America relies on relative sea level information which is primarily obtained from areas far away from the uplift region. The lack of accurate geodetic observations in the Great Lakes region, which is located in the transition zone between uplift and subsidence due to the deglaciation of the Laurentide ice sheet, has prevented more detailed studies of this former margin of the ice sheet. Recently, observations of vertical crustal motion from improved GPS network solutions and combined tide gauge and satellite altimetry solutions have become available. This study compares these vertical motion observations with predictions obtained from 70 different GIA models. The ice sheet margin is distinct from the centre and far field of the uplift because the sensitivity of the GIA process towards Earth parameters such as mantle viscosity is very different. Specifically, the margin area is most sensitive to the uppermost mantle viscosity and allows for better constraints of this parameter. The 70 GIA models compared herein have different ice loading histories (ICE-3/4/5G) and Earth parameters including lateral heterogeneities. The root-mean-square differences between the 6 best models and the two sets of observations (tide gauge/altimetry and GPS) are 0.66 and 1.57 mm/yr, respectively. Both sets of independent observations are highly correlated and show a very similar fit to the models, which indicates their consistent quality. Therefore, both data sets can be considered as a means for constraining and assessing the quality of GIA models in the Great Lakes region and the former margin of the Laurentide ice sheet.     

Glacial isostatic adjustment: New developments from advanced observing systems and modeling

The measurement of glacial isostatic adjustment (GIA) is one of the key ways in which geophysicists probe the long-term mantle rheology and Pleistocene ice history. GIA models are also tied to global and regional relative sea-level (RSL) histories, to 20th century tide-gauge (TG) data and to space and terrestrial geodetic measurements. Two new types of observation are related to the high-resolution space–gravity data recovered from the Gravity and Climate Experiment (GRACE) satellite pair and the soon-to-be launched Gravity and Ocean Circulation Experiment (GOCE) with on-board three-component gradiometer. Gravity mapping has the unique capability of isolating those regions that lack isostatic equilibrium. When coupled with other space and terrestrial geodetic measurements, such as those of the Global Positioning System (GPS) networks and with multi-decade terrestrial gravity data, new constraints on GIA are in the offing and should soon illuminate new interpretations of ice-sheet history and mantle response. GIA studies also incorporate space-based altimetry data, which now provide multi-decadal coverage over continents, oceans and lakes. As we are approaching 72 monthly solutions of GRACE gravity coefficients for determining the Earth's secular component of gravity change over the continents, a new issue has surfaced: the problem of relying on interannual hydrological modeling to determine the hydrological contribution to the linear trend in the gravity field. Correctly extracting this contribution is germane to using the GIA-driven component for modeling solid-Earth and paleo-climatic parameters.Seismic and heat-flux-based models of the Earth's interior are emerging with ever higher levels of sophistication regarding material strength (or viscosity). A basic question raised is: how good are traditional Newtonian and non-Newtonian viscosity models that only allow radial variations of Earth parameters? In other words: under what circumstances must this assumption be abandoned for joint interpretations of new and traditional data sets. In this short review we summarize the issues raised in the papers forming this special issue (SI) dedicated to GIA.     

Ground surface temperature history in southern Canada: Temperatures at the base of the Laurentide ice sheet and during the Holocene

We use temperature profiles from 7 deep (≈ 2000 m) boreholes located in southern Canada to infer ground surface temperature histories (GSTH) during the Last Glacial Maximum (LGM) and the Holocene. Visual inspection of the heat flow and of the reduced temperature depth profiles reveals significant regional differences with some sites showing conspicuous signs of post glacial warming, and other indicating only very small changes in ground surface temperature. These differences are confirmed by the inversions of the temperature profiles. The most prominent variations in GST are found at the Sudbury, Ontario, sites where the present ground surface temperature is high. With the exception of Sept-Iles, Quebec, the other sites only show moderate or no variation in GST. For all the sites, except possibly Sept-Iles, temperatures at the base of the ice sheet during the LGM were at or slightly below the melting point of ice. Temperatures might have been lower, a few degrees below 0 °C, at Sept-Iles. These results are consistent with field observations and model predictions suggesting high velocity basal flows in the ice sheet above the studied regions. These new data on basal temperatures will provide better quantitative constraints on glacier flow dynamics. The inversions give a chronology for the retreat of the ice sheet comparable to other proxies. Inversion and direct modeling show that, following the ice retreat, there was a warm period between 2 and 5 ka with temperatures 1–2 K higher than present. The inversion yields a time for this episode 1–2 kyr more recent than that inferred by other proxies for the Holocene climate optimum (HCO).     

Glacial isostatic adjustment in Fennoscandia from GRACE data and comparison with geodynamical models

The Earth’s gravity field observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission shows variations due to the integral effect of mass variations in the atmosphere, hydrosphere and geosphere. Several institutions, such as the GeoForschungsZentrum (GFZ) Potsdam, the University of Texas at Austin, Center for Space Research (CSR) and the Jet Propulsion Laboratory (JPL), Pasadena, provide GRACE monthly solutions, which differ slightly due to the application of different reduction models and centre-specific processing schemes. The GRACE data are used to investigate the mass variations in Fennoscandia, an area which is strongly influenced by glacial isostatic adjustment (GIA). Hence the focus is set on the computation of secular trends. Different filters (e.g. isotropic and non-isotropic filters) are discussed for the removal of high frequency noise to permit the extraction of the GIA signal. The resulting GRACE based mass variations are compared to global hydrology models (WGHM, LaDWorld) in order to (a) separate possible hydrological signals and (b) validate the hydrology models with regard to long period and secular components. In addition, a pattern matching algorithm is applied to localise the uplift centre, and finally the GRACE signal is compared with the results from a geodynamical modelling. The GRACE data clearly show temporal gravity variations in Fennoscandia. The secular variations are in good agreement with former studies and other independent data. The uplift centre is located over the Bothnian Bay, and the whole uplift area comprises the Scandinavian Peninsula and Finland. The secular variations derived from the GFZ, CSR and JPL monthly solutions differ up to 20%, which is not statistically significant, and the largest signal of about 1.2 Gal/year is obtained from the GFZ solution. Besides the GIA signal, two peaks with positive trend values of about 0.8 Gal/year exist in central eastern Europe, which are not GIA-induced, and also not explainable by the hydrology models. This may indicate that the recent global hydrology models have to be revised with respect to long period and secular components. Finally, the GRACE uplift signal is also in quite good agreement with the results from a simple geodynamical modelling.     

Geothermal control on flow patterns in the Last Glacial Maximum ice sheet of Iceland

Because it is located both on the Mid‐Atlantic Ridge and on a mantle plume, Iceland is a region of intense tectonics and volcanism. During the last glaciation, the island was covered by an ice sheet approximately 1000 m thick. A reconstruction of the ice flow lines, based on glacial directional features, shows that the ice sheet was partly drained through fast‐flowing streams. Fast flow of the ice streams has been recorded in megascale lineations and flutes visible on the currently deglaciated bedrock, and is confirmed by simple mass balance considerations. Locations of the major drainage routes correlate with locations of geothermal anomalies, suggesting that ice stream activity was favoured by lubrication of the bed by meltwater produced in regions of high geothermal heat flux. Similar control of ice flow by geothermal activity is expected in ice sheets currently covering tectonically and volcanically active area such as the West Antarctic ice sheet. Copyright © 2000 John Wiley & Sons, Ltd.     

冰川均衡调整对南极冰质量平衡监测的影响及其不确定性

本文研究了新的全球冰川均衡调整(GIA)模型对南极冰盖质量平衡监测的影响,考虑现有冰川负荷模型和地幔黏滞度模型的差异,完整评估了结果的不确定性,最后结合GRACE和卫星测高的结果进行了对比分析.结果表明,GIA对GRACE监测的等效水柱变化有重大影响,较大的GIA影响出现在西南极,沿罗斯冰架-卡姆布冰流-罗尼冰架-南极半岛分布,最大值在卡姆布冰流,达到29.8 mm/a;GIA对南极整体冰质量平衡的影响达到134±28 Gt/a.在不确定性的方差中,西南极和东南极分别以冰负荷模型差异和地幔黏滞度差异影响为主,对整个南极,冰模型差异影响占88.4%;在一些典型地区,GRACE监测的等效水柱在扣除GIA前后分别是,卡姆布冰流~32.8 mm/a和~6.3 mm/a,阿蒙森海湾~-95.3 mm/a和~-102.5 mm/a,Enderby Land~13.6 mm/a和~8.1 mm/a.整个南极冰盖总质量变化在扣除GIA贡献后为-82±29 Gt/a,该估计与卫星测高结果较吻合.此外,GIA对卫星测高监测的冰面高程变化的影响一般不超过8%.本研究为空间大地测量监测南极冰质量平衡提供了新的改正模型.     

Glacial isostatic adjustment in 3-D earth models: Implications for the analysis of tide gauge records along the U.S. east coast

Tide gauge records of recent sea-level change along the U.S. east coast have received significant attention within the literature of glacial isostatic adjustment (GIA). Geographic trends in these tide gauge rates are not reduced by a GIA correction based on a commonly adopted radial viscosity profile (characterized, in particular, by a lower mantle viscosity 1−2×10²¹ Pa s), and this has led to speculation that the residual trends reflect contributions from neotectonics or oceanographic processes. While the trends can be significantly reduced by adopting an Earth model with a stiffer lower mantle, such a model appears to be incompatible with independent constraints from post-glacial decay times in Hudson Bay. We use a finite-element model of the GIA process to investigate whether 3-D viscosity variations superimposed onto the “common” radial viscosity profile may provide a route to reconciling the east coast sea-level trends. We find that the specific 3-D structure we impose has little impact on the geographic trends in the GIA-corrected rates. However, we do find that the imposed lateral variations in lower mantle viscosity introduce a nearly uniform upward shift of 0.5 mm/yr in GIA-induced sea-level rates along the U.S. east coast. Thus, inferences of regional (U.S. east coast) sea-level rise due to modern melting of ice reservoirs, based on tide gauge rates corrected using 1-D GIA models, may be significantly biased by this simplifying assumption.     

Glacial curvilineations found along the southern sector of the Laurentide Ice sheet and a hypothesis of formation involving subglacial slope failure in tunnel valleys and subglacial lakes

First discovered in Poland, glacial curvilineations (GCLs) are enigmatic landforms comprising parallel sets of sinuous ridges and troughs of metres amplitude and around 150 m wavelength, found within kilometres‐wide valleys interpreted as being produced by meltwater flowing subglacially. Their morphological and sedimentary characteristics and association with tunnel valleys has been described for some prominent Polish examples. From these observations the existing hypothesis is that they form as a consequence of erosion by longitudinal vortices that develop in subglacial floods. Here we report, for the first time, GCLs found along the southern sector of the Laurentide Ice Sheet in three northern states of the USA. Using mapping and topographic analysis from high resolution digital elevation models we report observations on their morphological properties and landform associations. We find aspects of their context and morphology difficult to explain using the existing hypothesis. We instead suggest that these glacial curvilineations are produced by subglacial bank and slope failures that locally widen tunnel valleys, or that occur near subglacial lake shorelines. Further investigation is required to test this hypothesis and to ascertain the mechanisms of proposed mass movements, which may have occurred by rotational or translational slope failure or by creep deformation. Our preferred mechanism is that such movements occurred where subglacial water was emplaced over previously perma‐frozen ground. Under such circumstances, sediment blocks thawed by the water may then easily glide over a frozen décollement at low slope angles; analogous to subaerial active‐layer glides in permafrost environments. Permafrost spring sapping may have provided lines of weakness for slope failure. If the requirement for permafrost is found to hold, then GCLs may become an important indicator of the palaeo‐distribution of permafrost. Copyright © 2017 John Wiley & Sons, Ltd.     

Glacial geomorphology and ice ages in Tibet and the surrounding mountains

Abstract   Since 1973 new data have been obtained on the maximum extent of glaciation in High Asia. Evidence for an ice sheet covering Tibet during the last glacial period means a radical rethinking about glaciation in the northern hemisphere. The ice sheet's subtropical latitude, vast size (2.4 million km²) and high elevation (6000 m a.s.l.) are supposed to have resulted in a substantial, albedo-induced cooling of the Earth's atmosphere and the disruption of summer monsoon circulation. Moraines were found to reach down to 460 m a.s.l. on the southern flank of the Himalayas and to 2300 m a.s.l. on the northern slope of the Tibetan Plateau, in the Qilian Shan region. On the northern slopes of the Karakoram, Aghil and Kuen-Lun Mountains, moraines occur as far down as 1900 m a.s.l. In southern Tibet, radiographic analyses of erratics suggest a former ice thickness of at least 1200 m. Glacial polish and roches moutonnées in the Himalayas and Karakoram suggest former glaciers as thick as 1200–2700 m. On the basis of this evidence, a 1100–1600-m lower equilibrium line has been reconstructed, resulting in an ice sheet of 2.4 million km², covering almost all of Tibet. Radiometric ages, obtained by different methods, classify this glaciation as isotope stage 3–2 in age (Würmian, the last glacial period, ca  60 000–18 000 years ago).     

Vulnerability of water resources under a changing climate and human activity in the lower Great Lakes region

Globally, the number of people experiencing water stress is expected to increase by millions by the end of the century. The Great Lakes region, representing 20% of the world's surface freshwater, is not immune to stresses on water supply due to uncertainties on the impacts of climate and land use change. It is imperative for researchers and policy makers to assess the changing state of water resources, even if the region is water rich. This research developed the integrated surface water-groundwater GSFLOW model and investigated the effects of climate change and anthropogenic activities on water resources in the lower Great Lakes region of Western New York. To capture a range of scenarios, two climate emission pathways and three land development projections were used, specifically RCP 4.5, RCP 8.5, increased urbanization by 50%, decreased urbanization by 50%, and current land cover, respectively. Model outputs of surface water and groundwater discharge into the Great Lakes and groundwater storage for mid- and late century were compared to historical to determine the direction and amplitude of changes. Both surface water and groundwater systems show no statistically significant changes under RCP 4.5 but substantial and worrisome losses with RCP 8.5 by mid-century and end of century. Under RCP 8.5, streamflow decreased by 22% for mid-century and 42% for late century. Adjusting impervious surfaces revealed complex land use effects, resulting in spatially varying groundwater head fluctuations. For instance, increasing impervious surfaces lowered groundwater levels from 0.5 to 3.8 m under Buffalo, the largest city in the model domain, due to reduced recharge in surrounding suburban areas. Ultimately, results of this study highlight the necessity of integrated modelling in assessing temporal changes to water resources. This research has implications for other water-rich areas, which may not be immune to effects of climate change and human activities.     

A system dynamics approach for urban water reuse planning: a case study from the Great Lakes region

Water reclamation and reuse practices are recently receiving growing attention due to increasing water scarcity, concerns about the effect of wastewater discharges on receiving water, and availability of high-performing and cost-effective water reuse technologies. However, incorporation of water reuse schemes into water/wastewater infrastructure systems is a complex decision making process, involving various economical, technological, and environmental criteria. System dynamics (SD) allows modeling of complex systems and provides information about the temporal and feedback behavior of the system. In this sense, a SD model of the existing water/wastewater system in Kalamazoo-Michigan, an urban area in the Great Lakes region, was created with the hypothetical incorporation of water reuse. The model simulates and optimizes the overall water system cost (including water, wastewater and water reuse components), accounting for future scenarios of population, economic growth and climate change. Results indicate significant levels of water reuse after an infrastructure build delay. The model also indicates that a decision to implement water reuse yields remarkably lower water withdrawals and lower water treatment costs even in a location with a relatively abundant water supply like Kalamazoo. This study emphasizes the fact that a true understanding of the practice of water reuse cannot be achieved without taking regional and climatic parameters into account.     

Glacial meltwater and sea ice meltwater in the Prydz Bay, Antarctica

It has been well known by oceanographers that the World Ocean Circulation originates inNorth Atlantic near the Greenland, where in wintertime the cooled surface water descends to form North Atlantic Deep Waters (NADW). The NADW, when passing the Antarctic…     

Role of the Lakes in Groundwater Recharge and Discharge in the Young Glacial Area,Northern Poland

The aim of this research was to delineate characteristic hydrogeological lake types in the Young Glacial Area (YGA). The YGA is in the central and east part of the Kashubian Lake District (KLD) in Northern Poland, an area covered by deposits of Quaternary glaciation. All the bigger lakes were investigated in the area of about 1500 km² (39 lakes). The role of lakes in groundwater recharge and discharge was determined from total dissolved solids (TDS) in lake waters and also from groundwater flow simulation. The general trend was that gaining lakes, as determined by flow modeling, had higher values of TDS than losing lakes. In addition to typical gaining lakes (with TDS > 250 mg/l), there were losing lakes perched on glacial till deposits with very low TDS (<100 mg/l). Two groups of losing lakes were delineated: ones with very low TDS and another group with slightly higher TDS (due to local contact with groundwater). Flow‐through lakes with TDS of 170–200 mg/l were also delineated.     

Radioactivity in sediments of the Great Lakes: Post-depositional redistribution by deposit-feeding organisms

At two locations in southern Lake Huron (U.S.A.), twelve 35.5-cm diameter cores of fine-grained sediments were taken for comparison of the vertical distributions of²¹⁰Pb and falllout¹³⁷Cs with the distributions of benthic macroinvertebrates, mainly oligochaete worms (Tubificidae) and the amphipod,Pontoporeia affinis. Locations were selected on the basis of²¹⁰Pb distributions measured a year earlier which indicated contrasting depths of mixing of surface sediments. At one location the activity of²¹⁰Pb is uniform down to about 6 cm and 95% of total invertebrates occur within this zone; at the other location the zone of constant activity is only 3 cm deep but 90% of the invertebrates occur within it. Comparison of published tubificid defecation rates with sediment accumulation rates based on²¹⁰Pb shows that oligochaetes alone can account for mixing in one case while the effects of amphipods may be required in the case of shallower mixing. If mixing is represented as a diffusional process, eddy diffusion coefficients are at least 5.8 and 3.3 cm² yr^?1 at respective locations. In comparison with bioturbation, molecular diffusion is of minor importance in the post-depositional mobility of¹³⁷Cs. The necessity for introducing a diffusion coefficient varying continuously with depth is indicated by characteristics of the distribution of¹³⁷Cs. Biological reworking of near-surface sediments is an important process affecting radioactivity and chemical profiles in profundal deposits of this and probably other Great Lakes.     

Features of isostatic gravity anomaly and seismic activity in the Central Asian region

Ｆｅａｔｕｒｅｓｏｆｉｓｏｓｔａｔｉｃｇｒａｖｉｔｙａｎｏｍａｌｙａｎｄｓｅｉｓ┐ｍｉｃａｃｔｉｖｉｔｙｉｎｔｈｅＣｅｎｔｒａｌＡｓｉａｎｒｅｇｉｏｎＳＨＥＮＧ－ＭＩＮＧＦＡＮＧ１）（方盛明）ＲＵＩＦＥＮＧ２）（冯锐）ＣＨＡＮＧ－ＺＨＥＮＧＴＩＡＮ...     

Many-year level variations in the Great Lakes of Eurasia and North America

Regularities in the secular variations of water level in the Great Lakes of Eurasia and North America under natural climate and anthropogenic impacts are studied.     

Temperature variations in the mesopause region according to the hydroxyl-emission observations at midlatitudes

The seasonal temperature variations in the mesopause region and the inter-day and nighttime temperature variability, the measure of which is standard deviations, have been studied based on the hydroxyl emission spectral observations at the Zvenigorod station of the Obukhov Institute of Atmospheric Physics in 2000–2011 and Institute of Solar Terrestrial Physics geophysical station (Tory) in 2008–2011. The long-term variations in all temperature variability parameters have been analyzed.