An eclectic morphostratigraphic model for the sedimentary response to Holocene sea-level rise in northwest Europe

The improved empirical understanding of silt facies in Holocene coastal sequences provided by such as diatom, foraminifera, ostracode and testate amoebae analysis, combined with insights from quantitative stratigraphic and hydraulic simulations, has led to an inclusive, integrated model for the palaeogeomorphology, stratigraphy, lithofacies and biofacies of northwest European Holocene coastal lowlands in relation to sea-level behaviour. The model covers two general circumstances and is empirically supported by a range of field studies in the Holocene deposits of a number of British estuaries, particularly, the Severn. Where deposition was continuous over periods of centuries to millennia, and sea level fluctuated about a rising trend, the succession consists of repeated cycles of silt and peat lithofacies and biofacies in which series of transgressive overlaps (submergence sequences) alternate with series of regressive overlaps (emergence sequences) in association with the waxing and waning of tidal creek networks. Environmental and sea-level change are closely coupled, and equilibrium and secular pattern is of the kind represented ideally by a closed limit cycle. In the second circumstance, characteristic of unstable wetland shores and generally affecting smaller areas, coastal erosion ensures that episodes of deposition in the high intertidal zone last no more than a few centuries. The typical response is a series of regressive overlaps (emergence sequence) in erosively based high mudflat and salt-marsh silts that record, commonly as annual banding, exceptionally high deposition rates and a state of strong disequilibrium. Environmental change, including creek development, and sea-level movement are uncoupled. Only if deposition proceeds for a sufficiently long period, so that marshes mature, are equilibrium and close coupling regained. (C) 2002 Elsevier Science B.V. All rights reserved.

An Assessment of the Recent Past and Future Climate Change, Glacier Retreat, and Runoff in the Caucasus Region Using Dynamical and Statistical Downscaling and HBV-ETH Hydrological Model

The paper discusses the observed and projected warming in the Caucasus region and its implications for glacier melt and runoff. A strong positive trend in summer air temperatures of 0.05 degrees C a(-1) is observed in the high-altitude areas providing for a strong glacier melt and continuous decline in glacier mass balance. A warming of 4-7 degrees C and 3-5 degrees C is projected for the summer months in 2071-2100 under the A2 and B2 emission scenarios respectively, suggesting that enhanced glacier melt can be expected. The expected changes in winter precipitation will not compensate for the summer melt and glacier retreat is likely to continue. However, a projected small increase in both winter and summer precipitation combined with the enhanced glacier melt will result in increased summer runoff in the currently glaciated region of the Caucasus (independent of whether the region is glaciated at the end of the twenty-first century) by more than 50% compared with the baseline period.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Link between DOC in near surface peat and stream water in an upland catchment

Hydrologic transport of dissolved organic carbon (DOC) from peat soils may differ to organo-mineral soils in how they responded to changes in flow, because of differences in soil profile and hydrology. In well-drained organo-mineral soils, low flow is through the lower mineral layer where DOC is absorbed and high flow is through the upper organic layer where DOC is produced. DOC concentrations in streams draining organo-mineral soils typically increase with flow. In saturated peat soils, both high and low flows are through an organic layer where DOC is produced. Therefore, DOC in stream water draining peat may not increase in response to changes in flow as there is no switch in flow path between a mineral and organic layer. To verify this, we conducted a high-resolution monitoring study of soil and stream water at an upland peat catchment in northern England. Our data showed a strong positive correlation between DOC concentrations at − 1 and − 5 cm depth and stream water, and weaker correlations between concentrations at − 20 to − 50 cm depth and stream water. Although near surface organic material appears to be the key source of stream water DOC in both peat and organo-mineral soils, we observed a negative correlation between stream flow and DOC concentrations instead of a positive correlation as DOC released from organic layers during low and high flow was diluted by rainfall. The differences in DOC transport processes between peat and organo-mineral soils have different implications for our understanding of long-term changes in DOC exports. While increased rainfall may cause an increase in DOC flux from peat due to an increase in water volume, it may cause a decrease in concentrations. This response is contrary to expected changes in DOC exports from organo-mineral soils, where increase rainfall is likely to result in an increase in flux and concentration.

The sea-level conundrum: case studies from palaeo-archives

Uncertainties in sea-level projections for the 21st century have focused ice sheet modelling efforts to include the processes that are thought to be contributing to the recently observed rapid changes at ice sheet margins. This effort is still in its infancy, however, leaving us unable to make reliable predictions of ice sheet responses to a warming climate if such glacier accelerations were to increase in size and frequency. The geological record, however, has long identified examples of nonlinear ice sheet response to climate forcing (Shackleton NJ, Opdyke ND. 1973. Oxygen isotope and paleomagnetic stratigraphy of equatorial Pacific core V28–239, late Pliocene to latest Pleistocene. Geological Society of America Memoirs145: 449–464; Fairbanks RG. 1989. A 17,000 year glacio-eustatic sea level record: influence of glacial melting rates on the Younger Dryas event and deep ocean circulation. Nature342: 637–642; Bard E, Hamelin B, Arnold M, Montaggioni L, Cabioch G, Faure G, Rougerie F. 1996. Sea level record from Tahiti corals and the timing of deglacial meltwater discharge. Nature382: 241–244), thus suggesting an alternative strategy for constraining the rate and magnitude of sea-level change that we might expect by the end of this century. Copyright © 2009 John Wiley & Sons, Ltd.

Hyperresolution global land surface modeling: meeting a grand challenge for monitoring Earth's terrestrial water

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (∼10–100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface‐subsurface interactions due to fine‐scale topography and vegetation; improved representation of land‐atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a “grand challenge” to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Ultra-distal tephra deposits from super-eruptions: examples from Toba, Indonesia and Taupo Volcanic Zone, New Zealand

Voluminous rhyolitic eruptions from Toba, Indonesia, and Taupo Volcanic Zone (TVZ), New Zealand, have dispersed volcanic ash over vast areas in the late Quaternary. The ~74 ka Youngest Toba Tuff (YTT) eruption deposited ash over the Bay of Bengal and the Indian subcontinent to the west. The ~340 ka Whakamaru eruption (TVZ) deposited the widespread Rangitawa Tephra, dominantly to the southeast (in addition to occurrences northwest of vent), extending across the landmass of New Zealand, and the South Pacific Ocean and Tasman Sea, with distal terrestrial exposures on the Chatham Islands. These super-eruptions involved ~2500 km^3 and ~1500 km3 of magma (dense-rock equivalent; DRE), respectively. Ultra-distal terrestrial exposures of YTT at two localities in India, Middle Son Valley, Madhya Pradesh, and Jurreru River Valley, Andhra Pradesh, at distances of >2000 km from the source caldera, show a basal ‘primary’ ashfall unit ~4 cm thick, although deposits containing reworked ash are up to ~3 m in total thickness. Exposures of Rangitawa Tephra on the Chatham Islands, >900 km from the source caldera, are ~15-30 cm thick. At more proximal localities (~200 km from source), Rangitawa Tephra is ~55-70 cm thick and characterized by a crystal-rich basal layer and normal grading. Both distal tephra deposits are characterized by very-fine ash (with high PM10 fractions) and are crystal-poor. Glass chemistry, stratigraphy and grain-size data for these distal tephra deposits are presented with comparisons of their correlation, dispersal and preservation. Using field observations, ash transport and deposition were modeled for both eruptions using a semi-analytical model (HAZMAP), with assumptions concerning average wind direction and strength during eruption, column shape and vent size. Model outputs provide new insights into eruption dynamics and better estimates of eruption volumes associ- ated with tephra fallout. Modeling based on observed YTT distal tephra thicknesses indicate a relatively low (<40 km high), very turbulent eruption column, consistent with deposition from a co-ignimbrite cloud extending over a broad region. Similarly, the Whakamaru eruption was modeled as producing a predominantly Plinian column (~45 km high), with dispersal to the southeast by strong prevailing winds. Significant ash fallout of the main dispersal direction, to the northwest of source, cannot be replicated in this modeling. The widespread dispersal of large volumes of fine ash from both eruptions may have had global environmental consequences, acutely affecting areas up to thousands of kilometers from vent.

Reconciling observed and modelled phytoplankton dynamics in a major lowland UK river system, the Thames

This study aims to elucidate the key mechanisms controlling phytoplankton growth and decay within the Thames basin through the application of a modified version of an established river-algal model and comparison with observed stream water chlorophyll-a concentrations. The River Thames showed a distinct simulated phytoplankton seasonality and behaviour having high spring, moderate summer and low autumn chlorophyll-a concentrations. Three main sections were identified along the River Thames with different phytoplankton abundance and seasonality: (i) low chlorophyll-a concentrations from source to Newbridge; (ii) steep concentration increase between Newbridge and Sutton; and (iii) high concentrations with a moderate increase in concentration from Sutton to the end of the study area (Maidenhead). However, local hydrologic (e.g. locks) and other conditions (e.g. radiation, water depth, grazer dynamics, etc.) affected the simulated growth and losses. The model achieved good simulation results during both calibration and testing through a range of hydrological and nutrient conditions. Simulated phytoplankton growth was controlled predominantly by residence time, but during medium–low flow periods available light, water temperature and herbivorous grazing defined algal community development. These results challenge the perceived importance of in-stream nutrient concentrations as the perceived primary control on phytoplankton growth and death.

Scheduling satellite-based SAR acquisition for sequential assimilation of water level observations into flood modelling

Satellite-based Synthetic Aperture Radar (SAR) has proved useful for obtaining information on flood extent, which, when intersected with a Digital Elevation Model (DEM) of the floodplain, provides water level observations that can be assimilated into a hydrodynamic model to decrease forecast uncertainty. With an increasing number of operational satellites with SAR capability, information on the relationship between satellite first visit and revisit times and forecast performance is required to optimise the operational scheduling of satellite imagery. By using an Ensemble Transform Kalman Filter (ETKF) and a synthetic analysis with the 2D hydrodynamic model LISFLOOD-FP based on a real flooding case affecting an urban area (summer 2007,Tewkesbury, Southwest UK), we evaluate the sensitivity of the forecast performance to visit parameters. We emulate a generic hydrologic-hydrodynamic modelling cascade by imposing a bias and spatiotemporal correlations to the inflow error ensemble into the hydrodynamic domain. First, in agreement with previous research, estimation and correction for this bias leads to a clear improvement in keeping the forecast on track. Second, imagery obtained early in the flood is shown to have a large influence on forecast statistics. Revisit interval is most influential for early observations. The results are promising for the future of remote sensing-based water level observations for real-time flood forecasting in complex scenarios.

Geoarchaeological investigations of midden-formation processes in the early to late Neolithic levels at Catalhoyuk, Turkey ca. 8550-8370 cal BP.

Anthropogenic midden deposits are remarkably well preserved at the Neolithic settlement of atalhöyük and provide significant archaeological information on the types and nature of activities occurring at the site. To decipher their complex stratigraphy and to investigate formation processes, a combination of geoarchaeological techniques was used. Deposits were investigated from the early ceramic to late Neolithic levels, targeting continuous sequences to examine high resolution and broader scale changes in deposition. Thin-section micromorphology combined with targeted phytolith and geochemical analyses indicates they are composed of a diverse range of ashes and other charred and siliceous plant materials, with inputs of decayed plants and organic matter, fecal waste, and sedimentary aggregates, each with diverse depositional pathways. Activities identified include in situ burning, with a range of different fuel types that may be associated with different activities. The complexity and heterogeneity of the midden deposits, and thus the necessity of employing an integrated microstratigraphic approach is demonstrated, as a prerequisite for cultural and palaeoenvironmental reconstructions.

Global snow mass measurements and the effect of stratigraphic detail on inversion of microwave brightness temperatures

Snow provides large seasonal storage of freshwater, and information about the distribution of snow mass as Snow Water Equivalent (SWE) is important for hydrological planning and detecting climate change impacts. Large regional disagreements remain between estimates from reanalyses, remote sensing and modelling. Assimilating passive microwave information improves SWE estimates in many regions but the assimilation must account for how microwave scattering depends on snow stratigraphy. Physical snow models can estimate snow stratigraphy, but users must consider the computational expense of model complexity versus acceptable errors. Using data from the National Aeronautics and Space Administration Cold Land Processes Experiment (NASA CLPX) and the Helsinki University of Technology (HUT) microwave emission model of layered snowpacks, it is shown that simulations of the brightness temperature difference between 19 GHz and 37 GHz vertically polarised microwaves are consistent with Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) and Special Sensor Microwave Imager (SSM/I) retrievals once known stratigraphic information is used. Simulated brightness temperature differences for an individual snow profile depend on the provided stratigraphic detail. Relative to a profile defined at the 10 cm resolution of density and temperature measurements, the error introduced by simplification to a single layer of average properties increases approximately linearly with snow mass. If this brightness temperature error is converted into SWE using a traditional retrieval method then it is equivalent to ±13 mm SWE (7% of total) at a depth of 100 cm. This error is reduced to ±5.6 mm SWE (3 % of total) for a two-layer model.

PERSiST: a flexible rainfall-runoff modelling toolkit for use with the INCA family of models

Runoff generation processes and pathways vary widely between catchments. Credible simulations of solute and pollutant transport in surface waters are dependent on models which facilitate appropriate, catchment-specific representations of perceptual models of the runoff generation process. Here, we present a flexible, semi-distributed landscape-scale rainfall-runoff modelling toolkit suitable for simulating a broad range of user-specified perceptual models of runoff generation and stream flow occurring in different climatic regions and landscape types. PERSiST (the Precipitation, Evapotranspiration and Runoff Simulator for Solute Transport) is designed for simulating present-day hydrology; projecting possible future effects of climate or land use change on runoff and catchment water storage; and generating hydrologic inputs for the Integrated Catchments (INCA) family of models. PERSiST has limited data requirements and is calibrated using observed time series of precipitation, air temperature and runoff at one or more points in a river network. Here, we apply PERSiST to the river Thames in the UK and describe a Monte Carlo tool for model calibration, sensitivity and uncertainty analysis