Scalar advection schemes for ocean modelling on unstructured triangular grids

Several schemes for scalar advection on unstructured triangular grids are assessed for possible use in ocean modelling applications. Finite element, finite volume and finite volume–element approaches are evaluated. A series of tests, including a numerical order of convergence analysis, idealized rotating cone and cylinder experiments, and transport of a tracer through the Stommel Gyre representation of ocean basin-scale circulation, are carried out. Volume element Eulerian–Lagrangian and third-order Runge-Kutta discontinuous Galerkin schemes are recommended for use in tracer studies. Taylor–Galerkin and second-order Runge–Kutta discontinuous Galerkin are found to be robust and accurate second-order schemes. When positivity is required, a fluctuation redistribution scheme was found to be an easily implemented, accurate, and computationally efficient approach. Responsible editor: Phil Dyke     

Fatty acid profiles of benthic environment associated with artificial reefs in subtropical Hong Kong

Artificial reefs can enhance habitat heterogeneity, especially in seabed degraded by bottom-dredging and trawling. However, the trophodynamics of such reef systems are not well understood. This study provided baseline data on trophic relationships in the benthic environment associated with artificial reefs in late spring and mid summer of subtropical Hong Kong, using fatty acid profiles as an indicator. Data from sediments collected at the reef base, materials from sediment traps deployed on top and bottom of the reefs, total particulate matter from the water column and oyster tissues from reef surface were subjected to principal component analysis. Results showed variations of fatty acid profiles in the total particulate matter, upper sediment trap and oyster tissue samples collected in the two samplings, indicating seasonal, trophodynamic changes within the reef system. The wastes produced by fish aggregating at the reefs can also contribute a source of biodeposits to the nearby benthic environment.     

The influence of dunes on mixing in a migrating salt‐wedge: Fraser River estuary,Canada

Dune bedforms and salt‐wedge intrusions are common features in many estuaries with sand beds, and yet little is known about the interactions between the two. Flow visualization with an echosounder and velocity measurements with an acoustic Doppler current profiler over areas of flat‐bed and sand dunes in the highly‐stratified Fraser River estuary, Canada, were used to examine the effect of dunes on interfacial mixing. As the salt‐wedge migrates upstream over the flat‐bed, mixing is restricted to the lower portion of the water column. However, as the salt‐wedge migrates into the dune field from the flat bed, there is a dramatic change in the flow, and large internal in‐phase waves develop over each of the larger dunes, with water from the salt‐wedge reaching the surface of the estuary. The friction Richardson number shows that bed friction is more important in interfacial mixing over the dunes than over the flat‐bed, and a plot of internal Froude Number versus obstacle (dune) height shows that the salt‐wedge flow over the dunes is mainly supercritical. Such bedforms can be expected to cause similar effects in interfacial mixing in other estuaries and sediment‐laden density currents, and may thus be influential in fluid mixing and sediment transport. Copyright © 2010 John Wiley & Sons, Ltd.     

Geological controls on submarine groundwater discharge in Long Bay,South Carolina (USA)

A combination of geophysical methods including continuous electrical resistivity and high-resolution Chirp sub-bottom profiling were utilized to characterize geologic controls on pore fluid salinity in the nearshore of Long Bay, SC. Resistivity values ranged from less than 1 Ω m to greater than 40 Ω m throughout the bay. Areas of elevated electrical resistivity suggest the influence of relatively fresher water on pore water composition. Geophysical evidence alone does not eliminate all ambiguity associated with lithological and porosity variations that may also contribute to electrical structure of shallow marine sediments. The anomalous field is of sufficient magnitude that lithological variation alone does not control the spatial distribution of elevated electrical resistivity zones. Geographical distribution of electrical anomalies and structures interpreted from nearby sub-bottom profiles indicates abrupt changes in shallow geologic units control preferential pathways for discharge of fresh water into the marine environment. Shore parallel resistivity profiles show dramatic decreases in magnitude with increasing distance from shore, suggesting a significant portion of the terrestrially driven fresh SGD in Long Bay is occurring via the surficial aquifer within a few hundred meters of shore.     

Long‐term landscape evolution: linking tectonics and surface processes

Research in landscape evolution over millions to tens of millions of years slowed considerably in the mid‐20th century, when Davisian and other approaches to geomorphology were replaced by functional, morphometric and ultimately process‐based approaches. Hack's scheme of dynamic equilibrium in landscape evolution was perhaps the major theoretical contribution to long‐term landscape evolution between the 1950s and about 1990, but it essentially ‘looked back’ to Davis for its springboard to a viewpoint contrary to that of Davis, as did less widely known schemes, such as Crickmay's hypothesis of unequal activity. Since about 1990, the field of long‐term landscape evolution has blossomed again, stimulated by the plate tectonics revolution and its re‐forging of the link between tectonics and topography, and by the development of numerical models that explore the links between tectonic processes and surface processes. This numerical modelling of landscape evolution has been built around formulation of bedrock river processes and slope processes, and has mostly focused on high‐elevation passive continental margins and convergent zones; these models now routinely include flexural and denudational isostasy. Major breakthroughs in analytical and geochronological techniques have been of profound relevance to all of the above. Low‐temperature thermochronology, and in particular apatite fission track analysis and (U–Th)/He analysis in apatite, have enabled rates of rock uplift and denudational exhumation from relatively shallow crustal depths (up to about 4 km) to be determined directly from, in effect, rock hand specimens. In a few situations, (U–Th)/He analysis has been used to determine the antiquity of major, long‐wavelength topography. Cosmogenic isotope analysis has enabled the determination of the ‘ages’ of bedrock and sedimentary surfaces, and/or the rates of denudation of these surfaces. These latter advances represent in some ways a ‘holy grail’ in geomorphology in that they enable determination of ‘dates and rates’ of geomorphological processes directly from rock surfaces. The increasing availability of analytical techniques such as cosmogenic isotope analysis should mean that much larger data sets become possible and lead to more sophisticated analyses, such as probability density functions (PDFs) of cosmogenic ages and even of cosmogenic isotope concentrations (CICs). PDFs of isotope concentrations must be a function of catchment area geomorphology (including tectonics) and it is at least theoretically possible to infer aspects of source area geomorphology and geomorphological processes from PDFs of CICs in sediments (‘detrital CICs’). Thus it may be possible to use PDFs of detrital CICs in basin sediments as a tool to infer aspects of the sediments' source area geomorphology and tectonics, complementing the standard sedimentological textural and compositional approaches to such issues. One of the most stimulating of recent conceptual advances has followed the considerations of the relationships between tectonics, climate and surface processes and especially the recognition of the importance of denudational isostasy in driving rock uplift (i.e. in driving tectonics and crustal processes). Attention has been focused very directly on surface processes and on the ways in which they may ‘drive’ rock uplift and thus even influence sub‐surface crustal conditions, such as pressure and temperature. Consequently, the broader geoscience communities are looking to geomorphologists to provide more detailed information on rates and processes of bedrock channel incision, as well as on catchment responses to such bedrock channel processes. More sophisticated numerical models of processes in bedrock channels and on their flanking hillslopes are required. In current numerical models of long‐term evolution of hillslopes and interfluves, for example, the simple dependency on slope of both the fluvial and hillslope components of these models means that a Davisian‐type of landscape evolution characterized by slope lowering is inevitably ‘confirmed’ by the models. In numerical modelling, the next advances will require better parameterized algorithms for hillslope processes, and more sophisticated formulations of bedrock channel incision processes, incorporating, for example, the effects of sediment shielding of the bed. Such increasing sophistication must be matched by careful assessment and testing of model outputs using pre‐established criteria and tests. Confirmation by these more sophisticated Davisian‐type numerical models of slope lowering under conditions of tectonic stability (no active rock uplift), and of constant slope angle and steady‐state landscape under conditions of ongoing rock uplift, will indicate that the Davis and Hack models are not mutually exclusive. A Hack‐type model (or a variant of it, incorporating slope adjustment to rock strength rather than to regolith strength) will apply to active settings where there is sufficient stream power and/or sediment flux for channels to incise at the rate of rock uplift. Post‐orogenic settings of decreased (or zero) active rock uplift would be characterized by a Davisian scheme of declining slope angles and non‐steady‐state (or transient) landscapes. Such post‐orogenic landscapes deserve much more attention than they have received of late, not least because the intriguing questions they pose about the preservation of ancient landscapes were hinted at in passing in the 1960s and have recently re‐surfaced. As we begin to ask again some of the grand questions that lay at the heart of geomorphology in its earliest days, large‐scale geomorphology is on the threshold of another ‘golden’ era to match that of the first half of the 20th century, when cyclical approaches underpinned virtually all geomorphological work. Copyright © 2007 John Wiley & Sons, Ltd.     

MT+, Integrating Magnetotellurics to Determine Earth Structure, Physical State, and Processes

As one of the few deep-earth imaging techniques, magnetotellurics provides information on both the structure and physical state of the crust and upper mantle. Magnetotellurics is sensitive to electrical conductivity, which varies within the earth by many orders of magnitude and is modified by a range of earth processes. As with all geophysical techniques, magnetotellurics has a non-unique inverse problem and has limitations in resolution and sensitivity. As such, an integrated approach, either via the joint interpretation of independent geophysical models, or through the simultaneous inversion of independent data sets is valuable, and at times essential to an accurate interpretation. Magnetotelluric data and models are increasingly integrated with geological, geophysical and geochemical information. This review considers recent studies that illustrate the ways in which such information is combined, from qualitative comparisons to statistical correlation studies to multi-property inversions. Also emphasized are the range of problems addressed by these integrated approaches, and their value in elucidating earth structure, physical state, and processes.     

High-order kernels for Riemannian wavefield extrapolation

Riemannian wavefield extrapolation is a technique for one‐way extrapolation of acoustic waves. Riemannian wavefield extrapolation generalizes wavefield extrapolation by downward continuation by considering coordinate systems different from conventional Cartesian ones. Coordinate systems can conform with the extrapolated wavefield, with the velocity model or with the acquisition geometry. When coordinate systems conform with the propagated wavefield, extrapolation can be done accurately using low‐order kernels. However, in complex media or in cases where the coordinate systems do not conform with the propagating wavefields, low order kernels are not accurate enough and need to be replaced by more accurate, higher‐order kernels. Since Riemannian wavefield extrapolation is based on factorization of an acoustic wave‐equation, higher‐order kernels can be constructed using methods analogous to the one employed for factorization of the acoustic wave‐equation in Cartesian coordinates. Thus, we can construct space‐domain finite‐differences as well as mixed‐domain techniques for extrapolation. High‐order Riemannian wavefield extrapolation kernels improve the accuracy of extrapolation, particularly when the Riemannian coordinate systems does not closely match the general direction of wave propagation.     

Thermal history of the Rio Muni (West Africa)–NE Brazil margins during continental breakup

We document the thermal record of breakup of the conjugate Rio Muni (West Africa) and NE Brazil margins using apatite fission track analysis, vitrinite reflectance data and stratigraphic observations from both margins. These results permit determination of the timing of four cooling episodes, and the temperature of samples at the onset of each episode. All samples are interpreted to have experienced higher temperatures in the geological past due to i) elevated basal heatflow (palaeogeothermal gradient in Rio Muni-1 well decaying from 58 °C/km during the Mid Cretaceous to 21.5 °C/km in the Late Cenozoic) and ii) progressive exhumation from formerly greater burial depth. A well constrained history of changing palaeogeothermal gradient allows for much more precise quantification of the thickness of eroded section (exhumation) than if a constant heatflow is assumed. Cooling episodes identified from the palaeotemperature data at 110–95 Ma (both margins) and 85–70 Ma (Rio Muni only) coincide with major unconformities signifying, respectively, the cessation of rifting (breakup) and compressional shortening that affected the African continent following the establishment of post-rift sedimentation (drift). The interval between these separate unconformities is occupied by allochthonous rafts of shallow-water carbonates recording gravitational collapse of a marginal platform. The rift shoulder uplift that triggered this collapse was enhanced by local transpression associated with the obliquely divergent Ascension Fracture Zone, and thermal doming due to the coeval St Helena and Ascension Plumes. The data also reveal a c.45–35 Ma cooling episode, attributed to deep sea erosion at the onset of Eo-Oligocene ice growth, and a c.15–10 Ma episode interpreted as the record of Miocene exhumation of the West African continental margin related to continent-wide plume development. Integration of thermal history methods with traditional seismic- and stratigraphy-based observations yields a dynamic picture of kilometre-scale fluctuations in base level through the breakup and early drift phases of development of these margins. Major unconformities at ocean margins are likely to represent composite surfaces recording not only eustasy, but also regional plate margin-generated deformation, local ‘intra-basinal’ reorganization, and the amplifying effect of negative feedbacks between these processes.