Using the pseudospectral method on curved grids for 2D elastic forward modelling1

When applying the conventional Fourier pseudospectral method (FSM) on a Cartesian grid that has a sufficient size to propagate a pulse, spurious diffractions from the staircase representation of the curved interfaces appear in the wavefield. It is demonstrated that these non-physical diffractions can be eliminated by using curved grids that conform to all the interfaces of the subsurface. Methods for solving the 2D acoustic wave equation using such curved grids have been published previously by the authors. Here the extensions to the full 2D elastic wave equations are presented. The curved grids are generated by using the so-called multiblock strategy which is a well-known concept in computational fluid dynamics. In principle the sub-surface is divided into a number of contiguous subdomains. A separate grid is generated for each subdomain patching the grid lines across domain boundaries to obtain a globally continuous grid. Using this approach, even configurations with pinch outs can be handled. The curved grid is taken to constitute a generalized curvilinear coordinate system. Thus, the elastic equations have to be written in a curvilinear frame before applying the numerical scheme. The method implies that twice the number of spatial derivatives have to be evaluated compared to the conventional FSM on a Cartesian grid. However, it is demonstrated that the extra terms are more than compensated for by the fewer grid points needed in the curved approach.     

PCBs and chlorinated pesticides in the atmosphere and aquatic organisms of Ross Island, Antarctica

PCBs, p,p′-DDT, p,p′-DDE and lindane (γ-hexachlorocyclohexane) were monitored in the lower atmosphere of Ross Island, in Antarctica for 2 yr. Geometrical means were 15.2 pg m⁻³ for PCBs, 2.0 pg m⁻³ for p,p′-DDT, 1.0 pg m⁻³ for p,p′-DDE and 25.8 pg m⁻³ for lindane. Atmospheric levels of lindane were positively correlated with temperature, and a significant difference was found between spring-summer and summer-winter concentrations. No season related differences were found for the other chlorinated hydrocarbons, possibly owing to their lower vapour pressure and the cold climate. Periods with increased atmospheric levels of PCBs and DDT compounds were recorded. Lindane, p,p′-DDE and PCBs were present in fish and zooplankton sampled close to Ross Island. Pollutant levels in the zooplankton (on an extractable fat basis) were highest during the Antarctic spring and autumn and were inversely correlated to their fat content.     

Handling of waste in ports

A present major problem for ports is arranging adequate reception facilities for ship-generated waste, as the lack of such facilities is a problem to many shipowners. To preserve the environment in and around the ports it is necessary also to deal efficiently with the port-generated waste.

Although international conventions on the retention of waste on board ships for subsequent discharge to shore reception facilities have been ratified by most seafaring nations, a considerable part of the ship-generated waste still goes into the sea. Much of the waste generated in ports will also be dumped in quantities that at least affect the local marine environment. The regulations governing the handling of port-generated waste are often national and/or local legislation, whereas the handling of ship-generated waste is governed by the MARPOL Convention in most parts of the world.

The handling of waste consists of two main phases—collection and treatment. Waste has to be collected in every port and on board every ship, whereas generally only some wastes are treated and to a certain degree in ports and on board ships.

This paper considers the different kinds of waste generated in both ports and on board ships, where and how it is generated, how it could be collected and treated. The two sources are treated together to show how some ship-generated waste may be treated in port installations primarily constructed for the treatment of the port-generated waste, making integrated use of the available treatment facilities.

The different types of wastes can be divided into different categories. A world-wide accepted categorization of ship-generated waste has been established in the MARPOL 73/78 Convention. This will also be used for grouping the port-generated waste.     

Neodymium isotopes in seawater from the Barents Sea and Fram Strait Arctic-Atlantic gateways

The neodymium concentration, C_Nd, and isotopic composition, ε_Nd, in seawater have been determined in the water column at five sites in the Barents Sea-Fram Strait area where most of the water exchange between the Arctic Ocean and the North Atlantic takes place. In the main Arctic Ocean inflow branch across the Barents Sea the concentration and isotopic composition (C_Nd = 15.5 pmol/kg and ε_Nd = −10.8) are similar to those reported for the northeastern Nordic Seas, which is consistent with this region being a source area for the Arctic inflow. Due to the addition of Nd from Svalbard shelf sediments, the C_Nd in the surface waters above 150 m, in the Fram Strait inflow branch is higher by a factor of 2 and the ε_Nd is shifted to lower values (−11.8).In the stratified Nansen Basin, where cold low salinity water overlies warmer Atlantic water the C_Nd and ε_Nd do not vary with the vertical temperature-salinity structure but are essentially constant and similar to those of the Atlantic inflow throughout the entire water column, down to 3700 m depth, which indicates that the Nd is to a large extent of Atlantic origin.Compared to the Atlantic inflow water, the Nd in the major Arctic Ocean outflow, the Fram Strait, show higher C_Nd in the surface waters above 150 m, and a higher ε_Nd (−9.8) throughout the entire water column down to 1300 m depth. Sources for the more radiogenic Nd isotopic composition in deep water of the Fram Strait outflow most likely involve boundary exchange with sediments on the shelf and slope as the water passes along the Canadian archipelago. River water is a possible source in the surface water but it also seems likely that Pacific water Nd, modified by interactions on the shelf, is an important component in the Fram Strait surface outflow. Changes in the relative proportions of inflow of river water and flow of Pacific water through the Arctic Ocean could thus influence the isotopic composition of Nd in the North Atlantic.     

Iron isotope variations in Holocene sediments of the Gotland Deep, Baltic Sea

Holocene sediments from the Gotland Deep basin in the Baltic Sea were investigated for their Fe isotopic composition in order to assess the impact of changes in redox conditions and a transition from freshwater to brackish water on the isotope signature of iron. The sediments display variations in δ⁵⁶Fe (differences in the ⁵⁶Fe/⁵⁴Fe ratio relative to the IRMM-14 standard) from −0.27 ± 0.09‰ to +0.21 ± 0.08‰. Samples deposited in a mainly limnic environment with oxygenated bottom water have a mean δ⁵⁶Fe of +0.08 ± 0.13‰, which is identical to the mean Fe isotopic composition of igneous rocks and oxic marine sediments. In contrast, sediments that formed in brackish water under periodically euxinic conditions display significantly lighter Fe isotope signatures with a mean δ⁵⁶Fe of −0.14 ± 0.19‰. Negative correlations of the δ⁵⁶Fe values with the Fe/Al ratio and S content of the samples suggest that the isotopically light Fe in the periodically euxinic samples is associated with reactive Fe enrichments and sulfides. This is supported by analyses of pyrite separates from this unit that have a mean Fe isotopic composition of −1.06 ± 0.20‰ for δ⁵⁶Fe. The supply of additional Fe with a light Fe isotopic signature can be explained with the shelf to basin Fe shuttle model. According to the Fe shuttle model, oxides and benthic ferrous Fe that is derived from dissimilatory iron reduction from shelves is transported and accumulated in euxinic basins. The data furthermore suggest that the euxinic water has a negative dissolved δ⁵⁶Fe value of about −1.4‰ to −0.9‰. If negative Fe isotopic signatures are characteristic for euxinic sediment formation, widespread euxinia in the past might have shifted the Fe isotopic composition of dissolved Fe in the ocean towards more positive δ⁵⁶Fe values.     

Depositional systems in multiphase rifts: seismic case study from the Lofoten margin,Norway

The evolution of depositional systems in multiphase rifts is influenced by the selective reactivation of faults between subsequent rift phases. The Middle Jurassic to Palaeocene tectonic history of the Lofoten margin, a segment of the North Atlantic rift system, is characterised by three distinct rift phases separated by long (>20 Myr) inter‐rift periods. The initial rift phase comprised a distinct fault initiation and linkage stage, whereas the later rift phases were characterised by selective reactivation of previously linked through‐going faults which resulted in immediate rift climax. Using 2‐D and 3‐D seismic reflection data in conjunction with shallow core data we present a 100 Myr record of shallow to deep marine depositional environments that includes deltaic clinoform packages, slope aprons and turbidite fans. The rapid re‐establishment of major faults during the later rift phases impacts on drainage systems and sediment supply. Firstly, the immediate localisation of strain and accumulation of displacement on few faults results in pronounced footwall uplift and possible fault block rotation along those faults, which makes it more likely for any antecedent fault‐transverse depositional systems to become reversed. Secondly, any antecedent axially‐sourced depositional systems that are inherited from the foregoing rift phase(s) are likely to be sustained after reactivation because such axial systems have already been directed around fault tips. Hence, the immediate localisation of strain through selective reactivation in the later rift phases restricts fault‐transverse sediment supply more than axial sediment supply, which is likely to be a key aspect of the tectono‐sedimentary evolution of multiphase rifts.     

Runoff modelling at two field slopes: use of in situ measurements of air permeability to characterize spatial variability of saturated hydraulic conductivity

The time required at a field site to obtain a few measurements of saturated hydraulic conductivity (K_s) will allow for many measurements of soil air permeability (k_a). This study investigates if k_a measured in situ (k_{a, in situ}) can be a substitute for measurement of K_s in relation to infiltration and surface runoff modelling. Measurements of k_{a, in situ} were carried out in two small agricultural catchments. A spatial correlation of the log‐transformed values existed having a range of approximately 100 m. A predictive relationship between K_s and k_a measured on 100‐cm³ soil samples in the laboratory was derived for one of the field slopes and showed good agreement with an earlier suggested predictive K_s–k_a relationship. In situ measurements of K_s and k_a suggested that the predictive relationships also could be used at larger scale. The K_s–k_a relationships together with the k_{a, in situ} data were applied in a distributed surface runoff (DSR) model, simulating a high‐intensity rainfall event. The DSR simulation results were highly dependent on whether the geometric average of k_{a, in situ} or kriged values of k_{a, in situ} was used as model input. When increasing the resolution of K_s in the DSR model, a limit of 30–40 m was found for both field slopes. Below this limit, the simulated runoff and hydrograph peaks were independent of resolution scale. If only a few randomly chosen values of K_s were used to represent the spatial variation within the field slope, very large deviations in repeated DSR simulation results were obtained, both with respect to peak height and hydrograph shape. In contrast, when using many predicted K_s values based on a K_s–k_a relationship and measured k_{a, in situ} data, the DSR model generally captured the correct hydrograph shape although simulations were sensitive to the chosen K_s–k_a relationship. As massive measurement efforts normally will be required to obtain a satisfactory representation of the spatial variability in K_s, the use of k_{a, in situ} to assess spatial variability in K_s appears a promising alternative. Copyright © 2004 John Wiley & Sons, Ltd.     

Stochastic Modeling of Subseismic Faults Conditioned on Displacement and Orientation Maps

Mathematical Geosciences - Subseismic faults are small faults or fractures that may be difficult to determine but can have large consequences for fluid flow and pressure communication in the...     

Palaeogeographic,climatic and tectonic change in southeastern Australia: the Late Neogene evolution of the Murray Basin

The Murray Basin is a low-lying but extensive intracratonic depocentre in southeastern Australia, preserving an extraordinary record of Late Neogene sedimentation. New stratigraphic and sedimentologic data allow the long-term evolution of the basin to be re-evaluated and suggest a significant role for: (1) tectonism in controlling basin evolution, and (2) progressive and step-wise climatic change beginning in the early Pleistocene. Tectonic change is associated with regional uplift, occurring at approximately the same rate from the early Pliocene until the present day, and possibly associated with changing mantle circulation patterns or plate boundary processes. This uplift led to the defeat and re-routing of the Murray River, Australia’s major continental drainage system. Key to our interpretation is recognition of timing relationships between four prominent palaeogeographic features – the Loxton-Parilla Sands strandplain, the Gambier coastal plain, palaeo megalake Bungunnia and the Kanawinka Escarpment. Geomorphic and stratigraphic evidence suggest that during the Early Pliocene the ancestral Murray River was located in western Victoria, flowing south along the Douglas Depression. Relatively small amounts of regional uplift (<200 m) defeated this drainage system, dramatically changing the palaeogeography of southeastern Australia and forming Plio-Pleistocene megalake Bungunnia. At its maximum extent Lake Bungunnia covered more than 50,000 km², making it one of the largest known palaeo- or modern-lakes in an intracontinental setting. Magnetostratigraphic constraints suggest lake formation c. 2.4 Ma. The formation of Lake Bungunnia influenced the Pliocene coastal dynamics, depriving the coastline of a sediment source and changing the coastal system from a prograding strandline system to an erosional one. Erosion during this period formed the Kanawinka Escarpment, a palaeo sea-cliff and one of the most prominent and laterally extensive geomorphic features in southeastern Australia. Marine sediments c. 800 ka to c. 1.16 Ma represent the time of re-establishment of depositional coastal dynamics and of a permanent outlet for the Murray River. This age range is consistent with our best estimate of the age of the youngest Lake Bungunnia sediments and points towards an early Pleistocene age for the demise of the lake system. The youngest Lake Bungunnia sediment, present on a number of distinct terraces, suggests that progressive, step-wise climatic change played a role in the demise of the lake. However, in order for the ancestral Murray River system to have been able to breach the pre-existing tectonic dam, it is likely that tectonic change and/or temporarily enhanced discharge was also significant. This scenario indicates that the modern Murray River has only been in existence for at most 700 ka.