Post-glacial sea-level changes around the Australian margin : a review

It has been known since Rhodes Fairbridge’s first attempt to establish a global pattern of Holocene sea-level change by combining evidence from Western Australia and from sites in the northern hemisphere that the details of sea-level history since the Last Glacial Maximum vary considerably across the globe. The Australian region is relatively stable tectonically and is situated in the ‘far-field’ of former ice sheets. It therefore preserves important records of post-glacial sea levels that are less complicated by neotectonics or glacio-isostatic adjustments. Accordingly, the relative sea-level record of this region is dominantly one of glacio-eustatic (ice equivalent) sea-level changes. The broader Australasian region has provided critical information on the nature of post-glacial sea level, including the termination of the Last Glacial Maximum when sea level was approximately 125 m lower than present around 21,000–19,000 years BP, and insights into meltwater pulse 1A between 14,600 and 14,300 cal. yr BP. Although most parts of the Australian continent reveals a high degree of tectonic stability, research conducted since the 1970s has shown that the timing and elevation of a Holocene highstand varies systematically around its margin. This is attributed primarily to variations in the timing of the response of the ocean basins and shallow continental shelves to the increased ocean volumes following ice-melt, including a process known as ocean siphoning (i.e. glacio-hydro-isostatic adjustment processes). Several seminal studies in the early 1980s produced important data sets from the Australasian region that have provided a solid foundation for more recent palaeo-sea-level research. This review revisits these key studies emphasising their continuing influence on Quaternary research and incorporates relatively recent investigations to interpret the nature of post-glacial sea-level change around Australia. These include a synthesis of research from the Northern Territory, Queensland, New South Wales, South Australia and Western Australia. A focus of these more recent studies has been the re-examination of: (1) the accuracy and reliability of different proxy sea-level indicators; (2) the rate and nature of post-glacial sea-level rise; (3) the evidence for timing, elevation, and duration of mid-Holocene highstands; and, (4) the notion of mid- to late Holocene sea-level oscillations, and their basis. Based on this synthesis of previous research, it is clear that estimates of past sea-surface elevation are a function of eustatic factors as well as morphodynamics of individual sites, the wide variety of proxy sea-level indicators used, their wide geographical range, and their indicative meaning. Some progress has been made in understanding the variability of the accuracy of proxy indicators in relation to their contemporary sea level, the inter-comparison of the variety of dating techniques used and the nuances of calibration of radiocarbon ages to sidereal years. These issues need to be thoroughly understood before proxy sea-level indicators can be incorporated into credible reconstructions of relative sea-level change at individual locations. Many of the issues, which challenged sea-level researchers in the latter part of the twentieth century, remain contentious today. Divergent opinions remain about: (1) exactly when sea level attained present levels following the most recent post-glacial marine transgression (PMT); (2) the elevation that sea-level reached during the Holocene sea-level highstand; (3) whether sea-level fell smoothly from a metre or more above its present level following the PMT; (4) whether sea level remained at these highstand levels for a considerable period before falling to its present position; or (5) whether it underwent a series of moderate oscillations during the Holocene highstand.

How yak farmers can hold back a glacial tsunami

As Himalayan glaciers melt, the natural dams formed beneath them become a dangerous threat to villages below. However, local yak farmers could soon have a simple solution.

RS and GIS in mapping of geomorphic records and understanding the local controls of glacial retreat from the Baspa Valley, Himachal Pradesh, India

We report here the role of remote sensing (RS) and geographical information system (GIS) in the identification of geomorphic records and understanding of the local controls on the retreat of glaciers of the Baspa Valley, Himachal Pradesh, India. The geomorphic records mapped are accumulation zone, exposed ablation zone, moraine-covered ablation zone, snout, deglaciated valley, lateral moraine, medial moraine, terminal moraine and hanging glacier. Details of these features and stages of deglaciation have been extracted from RS data and mapped in a GIS environment. Glacial geomorphic data have been generated for 22 glaciers of the Baspa Valley. The retreat of glaciers has been estimated using the glacial maxima observed on satellite images. On the basis of percentage of retreat and the critical analysis of glacial geomorphic data for 22 glaciers of the Baspa Valley, they are classified into seven categories of very low to very very high retreat. From the analysis of the above 22 glaciers, it has been found that other than global warming, the retreat of glaciers of the Baspa Valley is inversely proportional to the size of the accumulation zone and the ratio of the moraine covered ablation/exposed ablation zone.

Did southern Western Ghats of peninsular India serve as refugia for its endemic biota during the Cretaceous volcanism?

The Western Ghats (WG) of south India, a global biodiversity hotspot, has experienced complex geological history being part of Gondwana landmass and encountered extensive volcanic activity at the end of Cretaceous epoch. It also has a climatically and topographically heterogeneous landscape. Thus, the WG offer a unique setting to explore the influence of ecological and geological processes on the current diversity and distribution of its biota. To this end, three explicit biogeographical scenarios were hypothesized to evaluate the distribution and diversification of wet evergreen species of the WG - (1) southern WG was a refuge for the wet evergreen species during the Cretaceous volcanism, (2) phylogenetic breaks in the species phylogeny would correspond to geographic breaks (i.e., the Palghat gap) in the WG, and (3) species from each of the biogeographic subdivisions within the WG would form distinct clades. These hypotheses were tested on the centipede genus Digitipes from the WG which is known to be an ancient, endemic, and monophyletic group. The Digitipes molecular phylogeny was subjected to divergence date estimation using Bayesian approach, and ancestral areas were reconstructed using parsimony approach for each node in the phylogeny. Ancestral-area reconstruction suggested 13 independent dispersal events to explain the current distribution of the Digitipes species in the WG. Among these 13 dispersals, two dispersal events were at higher level in the Digitipes phylogeny and were from the southern WG to the central and northern WG independently in the Early Paleocene, after the Cretaceous Volcanism. The remaining 11 dispersal events explained the species' range expansions of which nine dispersals were from the southern WG to other biogeographic subdivisions in the Eocene-Miocene in the post-volcanic periods where species-level diversifications occurred. Taken together, these results suggest that southern WG might have served as a refuge for Digitipes species during Cretaceous volcanism.

The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India

We investigated area changes in glaciers covering an area of similar to 200 km(2) in the Tista basin, Sikkim, Eastern Indian Himalaya, between similar to 1990 and 2010 using Landsat Thematic Mapper (TM) and Indian Remote-sensing Satellite (IRS) images and related the changes to debris cover, supraglacial lakes and moraine-dam lakes. The glaciers lost an area of 3.3 +/- 0.8% between 1989/90 and 2010. More detailed analysis revealed an area loss of 2.00 +/- 0.82, 2.56 +/- 0.61 and 2.28 +/- 2.01 km(2) for the periods 1989-97, 1997-2004/05 and 2004-2009/10, respectively. This indicates an accelerated retreat of glaciers after 1997. On further analysis, we observed (1) the formation and expansion of supraglacial lakes on many debris-covered glaciers and (2) the merging of these lakes over time, leading to the development of large moraine-dam lakes. We also observed that debris-covered glaciers with lakes lose a greater area than debris-covered glaciers without lakes and debris-free glaciers. The climatic data for 24 years (1987-2011), measured at the Gangtok meteorological station (1812 m a.s.l.), showed that the region experienced a 1.0 degrees C rise in the summer minimum temperature and a 2.0 degrees C rise in the winter minimum temperature, indicating hotter summers and warmer winters. There was no significant trend in the total annual precipitation. We find that glacier retreat is caused mainly by a temperature increase and that debris-covered glaciers can retreat at a faster rate than debris-free glaciers, if associated with lakes.

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6 +/- 1 Gta(-1), which decreases about sixfold to -35 +/- 2 Gta(-1) by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12 +/- 2 Gta(-1) by the 2080s. We also find that 10.6 and 27 % of the glaciers could face `eventual disappearance' by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.

The short-timescale behavior of glacial ice

Glaciers are often assumed to deform only at slow (i.e., glacial) rates. However, with the advent of high rate geodetic observations of ice motion, many of the intricacies of glacial deformation on hourly and daily timescales have been observed and quantified. This thesis explores two such short timescale processes: the tidal perturbation of ice stream motion and the catastrophic drainage of supraglacial meltwater lakes. Our investigation into the transmission length-scale of a tidal load represents the first study to explore the daily tidal influence on ice stream motion using three-dimensional models. Our results demonstrate both that the implicit assumptions made in the standard two-dimensional flow-line models are inherently incorrect for many ice streams, and that the anomalously large spatial extent of the tidal influence seen on the motion of some glaciers cannot be explained, as previously thought, through the elastic or viscoelastic transmission of tidal loads through the bulk of the ice stream. We then discuss how the phase delay between a tidal forcing and the ice stream’s displacement response can be used to constrain in situ viscoelastic properties of glacial ice. Lastly, for the problem of supraglacial lake drainage, we present a methodology for implementing linear viscoelasticity into an existing model for lake drainage. Our work finds that viscoelasticity is a second-order effect when trying to model the deformation of ice in response to a meltwater lake draining to a glacier’s bed. The research in this thesis demonstrates that the first-order understanding of the short-timescale behavior of naturally occurring ice is incomplete, and works towards improving our fundamental understanding of ice behavior over the range of hours to days.

The deep ocean density structure at the Last Glacial Maximum : What was it and why?

The search for reliable proxies of past deep ocean temperature and salinity has proved difficult, thereby limiting our ability to understand the coupling of ocean circulation and climate over glacial-interglacial timescales. Previous inferences of deep ocean temperature and salinity from sediment pore fluid oxygen isotopes and chlorinity indicate that the deep ocean density structure at the Last Glacial Maximum (LGM, approximately 20,000 years BP) was set by salinity, and that the density contrast between northern and southern sourced deep waters was markedly greater than in the modern ocean. High density stratification could help explain the marked contrast in carbon isotope distribution recorded in the LGM ocean relative to that we observe today, but what made the ocean's density structure so different at the LGM? How did it evolve from one state to another? Further, given the sparsity of the LGM temperature and salinity data set, what else can we learn by increasing the spatial density of proxy records?

We investigate the cause and feasibility of a highly and salinity stratified deep ocean at the LGM and we work to increase the amount of information we can glean about the past ocean from pore fluid profiles of oxygen isotopes and chloride. Using a coupled ocean--sea ice--ice shelf cavity model we test whether the deep ocean density structure at the LGM can be explained by ice--ocean interactions over the Antarctic continental shelves, and show that a large contribution of the LGM salinity stratification can be explained through lower ocean temperature. In order to extract the maximum information from pore fluid profiles of oxygen isotopes and chloride we evaluate several inverse methods for ill-posed problems and their ability to recover bottom water histories from sediment pore fluid profiles. We demonstrate that Bayesian Markov Chain Monte Carlo parameter estimation techniques enable us to robustly recover the full solution space of bottom water histories, not only at the LGM, but through the most recent deglaciation and the Holocene up to the present. Finally, we evaluate a non-destructive pore fluid sampling technique, Rhizon samplers, in comparison to traditional squeezing methods and show that despite their promise, Rhizons are unlikely to be a good sampling tool for pore fluid measurements of oxygen isotopes and chloride.

Long-term natural acidification of upland sites in Cumbria: Evidence from post-glacial lake sediments

Studies by the Freshwater Biological Association over the last 25 years have supplied data relevant to the levels of acidity in local soils and water before the onset of industrial pollution and current interest in acid rain. This article reviews published analysis from cores of lake sediments, in or near the catchment of the River Duddon. Electron spin resonance spectra of humic acids and iodine values confirm evidence from pollen analysis for a history of progressive acidification of the source material of lake sediments since before 5000 radiocarbon years, in upland catchments of the Lake District. Processes involved included: removal of basic ions from soils by rainfall, the effects of which were intensified by removal by man of deciduous forest; acidification of soils and waters by decomposition products of Calluna and further acidification of waters by Sphagnum species which colonized habitats where drainage became impeded by paludification processes.

Comparison of survey methods for estimating abundance of harbor seals (Phoca vitulina) in glacial fjords

The importance of glacial ice habitats to harbor seals (Phoca vitulina) in Alaska has become increasingly apparent. However, enumerating harbor seals hauled out on ice in glacial fjords has been difficult. At Johns Hopkins Inlet in Glacier Bay, Alaska, we compared a shore-based counting method to a large-format aerial photography method to estimate seal abundance. During each aerial survey, shore-based observers simultaneously counted seals from an observation post. Both survey methods incurred errors in double-counting and missing seals, especially when ice movements caused seals to drift between survey zones. Advantages of shore-based counts included the ability to obtain multiple counts for relatively little cost, distinguish pups from adults, and to distinguish mobile seals from shadows or glacial debris of similar size. Aerial photography provided a permanent record of each survey, allowing both a reconciliation of counts in overlapping zones and the documentation of the spatial distribution of seals and ice within the fjord.

A survey of ichthyofauna of Lake Kanyaboli and other small waterbodies in Kenya: alternative refugia for endangered fish species

In 1988, the World Conservation Union (WCU) Red Book of Endangered Species listed hundreds of endemic fishes of Lake Victoria under a single heading - "ENDANGERED". Most of the endemic native food fishes are either endangered or extinct. However, a survey of the fauna of Lake Kanyaboli, revealed that a few remaining samples of these native fishes are actually thriving. These include several unidentified Haplochromis spp., Oreochromis esculentus and Oreochromis variabilis. As a result, a stock rehabilitation and management strategy has been designed to use Lake Kanyaboli and other small waterbodies as conservation 'Refugia' for endangered fish species of the larger Lake Victoria.

A nearsurface, daytime occurrence of two mesopelagic fish species (Stenobrachius leucopsarus and Leuroglossus schmidti) in a glacial fjord

The northern lampfish (Stenobrachius leucopsarus, family Myctophidae) and northern smoothtongue (Leuroglossus schmidti, family Bathylagidae) are mesopelagic fishes, defined by their vertical distribution in the mesopelagic zone (200–1000 m) during daylight hours. Northern lampfish range from the Bering Sea to southern California (Shimada, 1948), where their abundance is highest along the continental slope and decreases over the continental shelf. They are the most abundant species in the mesopelagic zone of the Bering Sea (Pearcy et al., 1977; Sobolevsky et al., 1996), the Gulf of Alaska (Purcell, 1996), and the eastern North Pacific Ocean off Oregon (Pearcy, 1964; Pearcy et al., 1977). Northern smoothtongue also concentrate in areas bordering the continental slope and are widely distributed from southern British Columbia to the Bering Sea (Peden, 1981) and are very abundant in the Okhotsk Sea (Sobolevsky et al., 1996).

Centennial oscillation of the climatic crescendo during the last glacial maximum: evidence from Lake Estancia, central New Mexico

Sediments deposited in late Pleistocene Lake Estancia, central New Mexico, contain a paleoclimatic record that includes the last glacial maximum and deglacial episode. Stratigraphic reconstruction of an interval representing the highstand of the lake that occurred during the last glacial maximum reveals ~2000-, ~600-, and ~200-year oscillations in lake level and climate. Shifting position of the polar jetstream in response to expansion and contraction of the North American ice sheet may be partly responsible for the millenial-scale changes in Lake Estancia but probably does not explain the centennial-scale oscillations.