A depth-integrated 2D coastal and estuarine model with conformal boundary-fitted mesh generation

Details are given of the development and application of a 2D depth-integrated, conformal boundary-fitted, curvilinear model for predicting the depth-mean velocity field and the spatial concentration distribution in estuarine and coastal waters. A numerical method for conformal mesh generation, based on a boundary integral equation formulation, has been developed. By this method a general polygonal region with curved edges can be mapped onto a regular polygonal region with the same number of horizontal and vertical straight edges and a multiply connected region can be mapped onto a regular region with the same connectivity. A stretching transformation on the conformally generated mesh has also been used to provide greater detail where it is needed close to the coast, with larger mesh sizes further offshore, thereby minimizing the computing effort whilst maximizing accuracy. The curvilinear hydrodynamic and solute model has been developed based on a robust rectilinear model. The hydrodynamic equations are approximated using the ADI finite difference scheme with a staggered grid and the solute transport equation is approximated using a modified QUICK scheme. Three numerical examples have been chosen to test the curvilinear model, with an emphasis placed on complex practical applications

Inferring animal husbandry strategies in coastal zones through stable isotope analysis: new evidence from the Flemish coastal plain (Belgium, 1st-15th century AD)

In a proof-of-concept study, Britton et al. (2008) demonstrated that the isotopic composition of halophytic plants can be traced in the skeletal tissues of their animal consumers. Here we apply the method to domestic herbivore remains (n = 303) from nine archaeological sites in or near the Flemish coastal plain (Belgium), where, prior to embankments, salt-marshes offered extensive pasture grounds for domestic herbivores. The sites span a period of ∼1500 years (Roman to late medieval period), during which the coastal landscape was progressively transformed from little managed wetlands to a fully embanked polder area. The bulk collagen data show variations between sites and over time, which are consistent with this historical framework and are interpreted as reflecting environmental change and differences in animal management in the coastal plain throughout the late Holocene. The study demonstrates the immense value of faunal stable isotope analysis for characterising coastal husbandry strategies beyond the means of traditional zooarchaeological techniques.

Generation of a buoyancy-driven coastal current by an Antarctic Polynya

Descent and spreading of high salinity water generated by salt rejection during sea ice formation in an Antarctic coastal polynya is studied using a hydrostatic, primitive equation three-dimensional ocean model called the Proudman Oceanographic Laboratory Coastal Ocean Modeling System (POLCOMS). The shape of the polynya is assumed to be a rectangle 100 km long and 30 km wide, and the salinity flux into the polynya at its surface is constant. The model has been run at high horizontal spatial resolution (500 m), and numerical simulations reveal a buoyancy-driven coastal current. The coastal current is a robust feature and appears in a range of simulations designed to investigate the influence of a sloping bottom, variable bottom drag, variable vertical turbulent diffusivities, higher salinity flux, and an offshore position of the polynya. It is shown that bottom drag is the main factor determining the current width. This coastal current has not been produced with other numerical models of polynyas, which may be because these models were run at coarser resolutions. The coastal current becomes unstable upstream of its front when the polynya is adjacent to the coast. When the polynya is situated offshore, an unstable current is produced from its outset owing to the capture of cyclonic eddies. The effect of a coastal protrusion and a canyon on the current motion is investigated. In particular, due to the convex shape of the coastal protrusion, the current sheds a dipolar eddy.

Flood inundation dynamics and socioeconomic vulnerability under environmental change

Floods are a major threat to human existence and historically have both caused the collapse of civilizations and forced the emergence of new cultures. The physical processes of flooding are complex. Increased population, climate variability, change in catchment and channel management, modified landuse and land cover, and natural change of floodplains and river channels all lead to changes in flood dynamics, and as a direct or indirect consequence, social welfare of humans. Section 5.16.1 explores the risks and benefits brought about by floods and reviews the responses of floods and floodplains to climate and landuse change. Section 5.08.2 reviews the existing modeling tools, and the top–down and bottom–up modeling frameworks that are used to assess impacts on future floods. Section 5.08.3 discusses changing flood risk and socioeconomic vulnerability based on current trends in emerging or developing countries and presents an alternative paradigm as a pathway to resilience. Section 5.08.4 concludes the chapter by stating a portfolio of integrated concepts, measures, and avant-garde thinking that would be required to sustainably manage future flood risk.

Reducing inconsistencies in point observations of maximum flood inundation level

Flood simulation models and hazard maps are only as good as the underlying data against which they are calibrated and tested. However, extreme flood events are by definition rare, so the observational data of flood inundation extent are limited in both quality and quantity. The relative importance of these observational uncertainties has increased now that computing power and accurate lidar scans make it possible to run high-resolution 2D models to simulate floods in urban areas. However, the value of these simulations is limited by the uncertainty in the true extent of the flood. This paper addresses that challenge by analyzing a point dataset of maximum water extent from a flood event on the River Eden at Carlisle, United Kingdom, in January 2005. The observation dataset is based on a collection of wrack and water marks from two postevent surveys. A smoothing algorithm for identifying, quantifying, and reducing localized inconsistencies in the dataset is proposed and evaluated showing positive results. The proposed smoothing algorithm can be applied in order to improve flood inundation modeling assessment and the determination of risk zones on the floodplain.

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.

New models of North West European Holocene palaeogeography and inundation

This paper presents new 500 year interval palaeogeographic models for Britain, Ireland and the North West French coast from 11000 cal. BP to present. These models are used to calculate the varying rates of inundation for different geographical zones over the study period. This allows for consideration of the differential impact that Holocene sea-level rise had across space and time, and on past societies. In turn, consideration of the limitations of the models helps to foreground profitable areas for future research.

Season of birth and sheep husbandry in late Roman and Medieval coastal Flanders: a pilot study using tooth enamel d18O analysis

From the early Roman period, there is archaeological evidence for the exploitation of the Flemish coastal plain (Belgium) for a range of activities, such as sheep herding on the then developing salt-marshes and salt-meadows for the production of wool. During the early Middle Ages, this culminated in the establishment of dedicated ‘sheep estates’. This phase of exploitation was followed by extensive drainage and land reclamation measures in the high Medieval period, transforming areas into grassland, suited for cattle breeding. As part of a larger project investigating the onset, intensification and final decline of sheep management in coastal Flanders in the historical period, this pilot study presents the results of sequential sampling and oxygen isotope analysis of a number of sheep teeth (M2, n = 8) from four late Roman and Medieval sites (dating from 4th to 15th century AD), in order to assess potential variations in season of birth between the different sites and through time. In comparison with published data from herds of known birth season, incremental enamel data from the Flemish sites are consistent with late winter/spring births, with the possibility of some instances of slightly earlier parturition. These findings suggest that manipulation of season of birth was not a feature of the sheep husbandry-based economies of early historic Flanders, further evidencing that wool production was the main purpose of contemporary sheep rearing in the region. Manipulation of season of birth is not likely to have afforded economic advantage in wool-centred economies, unlike in some milk- or meat-based regimes.

Anthropogenic nutrient sources and loads from a Mediterranean catchment into a coastal lagoon: Mar Menor, Spain

The Mar Menor is a coastal lagoon increasingly threatened by urban and agricultural pressures. The main watercourse draining into the lagoon is the Rambla del Albujón. A fortnightly campaign carried out over one annual cycle enabled us to characterize the treated urban sewage effluents and agricultural sources which contribute to the nutrient fluxes in the watercourse. Multivariate analysis provided information for establishing chemical signatures and for assessing the relative influence of the various sources on the water quality at the outlet. Mass balances were used to examine net gains and losses, and cross-correlations with rainfall to analyze climatic influence and control factors in the trends of the nutrient flux. The rainfall pattern was significantly cross-correlated with nitrate and phosphorus fluxes from agricultural sources, while fluctuations in the resident population explained the phosphorus flux trend in urban sources. 50% of dissolved inorganic nitrogen was from agricultural sources, while 70% of total phosphate and 91% of total organic carbon were from urban point sources. The net amounts of all the nutrients fell as a result of plant uptake and/or denitrification in the channel. The control of urban point sources (phosphorus-enriched) is suggested as a promptly action for improving the health of the coastal lagoon.

Low-level mixing height detection in coastal locations with a scanning Doppler lidar

Mixing layer height (MLH) is one of the key parameters in describing lower tropospheric dynamics and capturing its diurnal variability is crucial, especially for interpreting surface observations. In this paper we introduce a method for identifying MLH below the minimum range of a scanning Doppler lidar when operated at vertical. The method we propose is based on velocity variance in low-elevation-angle conical scanning and is applied to measurements in two very different coastal environments: Limassol, Cyprus, during summer and Loviisa, Finland, during winter. At both locations, the new method agrees well with MLH derived from turbulent kinetic energy dissipation rate profiles obtained from vertically pointing measurements. The low-level scanning routine frequently indicated non-zero MLH less than 100 m above the surface. Such low MLHs were more common in wintertime Loviisa on the Baltic Sea coast than during summertime in Mediterranean Limassol.

Mercury and methylmercury bioaccumulation by polychaete worms is governed by both feeding ecology and mercury bioavailability in coastal mudflats

Polychaete worms are abundant in many mudflats but their importance to coastal food web Hg biomagnification is not known. We sampled sediments and polychaete worms from mudflats in the Bay of Fundy to investigate the bioaccumulation of mercury (Hg) and methylmercury (MeHg) in the coastal invertebrate food web. Hg concentrations in the sediments were low (<20 μg kg−1). Labile Hg (methanol/KOH sediment extraction) in surface sediments (0–1 cm) was positively correlated with Hg bioaccumulation by surface sediment-ingesting polychaetes but, surprisingly, there was a negative correlation between δ15N (i.e. trophic level) and THg bioaccumulation factors in polychaete worms. Worms feeding on deeper sediments contained the greatest MeHg concentrations (69.6 μg kg−1). Polychaetes are an important vector for Hg biomagnification to the coastal avian food web. This research demonstrates that feeding depth and method of feeding are more important than trophic position or sediment Hg concentrations for predicting Hg bioaccumulation.

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

Shifting sand, shifting livelihoods? Reflections on a coastal gold rush in Ghana

Artisanal and small-scale mining (ASM) is often characterised as rush-type activity undertaken by people looking to ‘get-rich quick’. An alternative view posits ASM as poverty-driven activity which provides a valuable source of employment to a rural population where agriculture has failed to provide an adequate income. Against this dichotomy, this article critically reviews the short-lived, but high profile, gold rush which occurred on the beaches of Elmina, a fishing town in Ghana’s Central Region. At its peak, more than 1000 people, including local fishermen and career miners from the Western Region and Eastern Region of the country, were mining on the beaches. Drawing on interviews with miners and local business owners, the analysis explores the dynamics of the rush and its implications for understanding the ASM sector in Ghana. The discussion highlights the challenges associated with managing resources in dynamic and heterogeneous contexts, reconciling separate and conflicting sectoral priorities and the need to go beyond simplified and dualistic representations of the sector.

Understanding 2H/1H systematics of leaf wax n-alkanes in coastal plants at Stiffkey saltmarsh, Norfolk, UK

Interpretation of sedimentary n-alkyl lipid d2H data is complicated by a limited understanding of factors controlling interspecies variation in biomarker 2H/1H composition. To distinguish between the effects of interrelated environmental, physical and biochemical controls on the hydrogen isotope composition of n-alkyl lipids, we conducted linked d2H analyses of soil water, xylem water, leaf water and n-alkanes from a range of C3 and C4 plants growing at a UK saltmarsh (i) across multiple sampling sites, (ii) throughout the 2012 growing season, and (iii) at different times of the day. Soil waters varied isotopically by up to 35& depending on marsh sub-environment, and exhibited site-specific seasonal shifts in d2H up to a maximum of 31 per mil. Maximum interspecies variation in xylem water was 38 per mil, while leaf waters differed seasonally by a maximum of 29 per mil. Leaf wax n-alkane 2H/1H, however, consistently varied by over 100 per mil throughout the 2012 growing season, resulting in an interspecies range in the ewax/leaf water values of -79 per mil to –227 per mil. From the discrepancy in the magnitude of these isotopic differences, we conclude that mechanisms driving variation in the 2H/1H composition of leaf water, including (i) spatial changes in soil water 2H/1H, (ii) temporal changes in soil water 2H/1H, (iii) differences in xylem water 2H/1H, and (iv) differences in leaf water evaporative 2H-enrichment due to varied plant life forms, cannot explain the range of n-alkane d2H values we observed. Results from this study suggests that accurate reconstructions of palaeoclimate regimes from sedimentary n-alkane d2H require further research to constrain those biological mechanisms influencing species-specific differences in 2H/1H fractionation during lipid biosynthesis, in particular where plants have developed biochemical adaptations to water-stressed conditions. Understanding how these mechanisms interact with environmental conditions will be crucial to ensure accurate interpretation of hydrogen isotope signals from the geological record.