On the ability of the LODYC GCM to simulate the thermocline depth in the equatorial Atlantic

Using a multivariate model testing procedure that distinguishes between model inadequacies and data uncertainties, we investigate the ability of the LODYC GCM to simulate the evolution of the 20°C isotherm depth during the 1982–1984 FOCAL/SEQUAL experiment in the equatorial Atlantic. Two different versions of the model are considered: the Ri version which has a Richardson number dependent parameterization of vertical mixing and the new TKE version which uses a local estimation of the turbulent kinetic energy to parameterize vertical mixing. Some effects of the forcing uncertainties are considered by forcing the TKE version with three equally plausible wind stress fields whose differences are consistent with the measurement and sampling errors, and the drag coefficient indeterminacy. The resulting uncertainties in the model response are substantial and can be as large as the differences between simulations with the two GCM versions, which stresses the need to take the forcing uncertainties into account. Although only one Ri run is available, it is shown that the TKE parameterization significantly improves the representation of the equatorial upwelling and the simulation of the depth of the thermocline in the eastern Atlantic. However, there remain significant differences with the observations which cannot be explained by the forcing uncertainties that were considered. The two model versions perform better in the equatorial wave guide than in the 12°N-12°S domain, and they are better distinguished over large domains than along sections, which shows that a global multivariate view point must be used in model-reality comparisons. Finally, a comparison with a linear multimode model emphasizes the need for greater model complexity to properly simulate the equatorial upwelling and the thermocline variability in the tropical Atlantic.This paper was presented at the Second International Conference on Modelling of Global Climate Variability, held in Hamburg 7–11 September 1992 under the auspices of the Max-Planck-Institute for Meteorology. Guest Editor for these papers is L. Dümenil     

Numerical Modeling of Wave-Enhanced Turbulence in the Oceanic Upper Layer

A coupled model of air-wave-sea interaction is modified based on a new roughness formulation and the latest data. The model parameters for aerodynamic roughness from below (ARB) and wave-dependent roughness from above (ARA, z _0a ) are assumed equal. The combined roughness is assumed to be a function of friction velocity, gravity, air and seawater densities, and wave age (c _w ). The model is used in a study of wave-enhanced turbulence under breaking waves to predict turbulent dissipation (), ARA, and drag coefficient (C _d ). Both waves and shear production are considered as sources of ocean turbulent energy. The atmospheric part of the model is used only to specify a correct condition at the interface. Numerical experiments are performed to study the -distribution, z _0a and C _d , and to compare with data. The major achievement is model verification using all available data. The first full application of this model is in conjunction with an ocean circulation model in a coupled circulation-wave system. Simulations show that the -distribution is strongly dependent on local wind-forced wave heights. For each wind and wave state there is a particular wave-dependent depth that is verified by data. The comparison shows that the model predicted agrees well with the observed of the z ^–4 law distribution of Gargett (1989). Simulations also show that waves have an important role in causing to differ from the classical wall-layer theory and z _0a , with a value of 0.30 for the empirical constant a _a . The model-predicted , z _0a , C _d and C _gd agree well with data.     

Horizontal Pressure Gradient Errors of the Monterey Bay Sigma Coordinate Ocean Model with Various Grids

A coastal ocean -coordinate model of Monterey Bay (MOB) with realistic bottom topography and coastlines is developed using the Princeton Ocean Model (POM) and grid generation technique (GGT) to study the horizontal pressure gradient errors associated with the MOB steep topography. The submarine canyon in MOB features some of the steepest topography encountered anywhere in the world oceans. The MOB grids are designed using the EAGEAL View and GENIE⁺⁺ grid generation systems. A grid package developed by Ly and Luong (1993) is used in this study to couple grids to the model. The MOB model is tested with both orthogonal and curvilinear nearly-orthogonal (CNO) grids. The CNO grid has horizontal resolution which varies from 300 m to 2 km, while the resolution of the orthogonal grid is uniform with x = 1.25 km and y = 1.38 km. These grids cover a domain of 180 × 160 km with the same number of grid points of 131 × 131. Vertical resolutions of 25, 35 and 45 vertical sigma levels are tested. The error in the MOB are evaluated in terms of mean kinetic energy and velocity against various grids, vertical, horizontal resolution and distributions, and bottom topography smoothing. Simulations with various grids show that GGT can be used as another tool in reducing -coordinate errors in coastal ocean modeling besides increasing resolution and smoothing bottom topography. Topographical smoothing not only reduces topographic slope, but changes realistic topography. A CNO grid with a high grid density packed along steep slopes and Monterey Submarine Canyon reduces the errors by 40% compared to a rectangular grid with the same number of grid points. The CNO grid is more efficient than the rectangular grid, since it has most of its grids over water. The simulations show that the presented MOB -coordinate model can be used with a confidence regarding horizontal pressure gradient error.     

Marine and Brackish Ecosystems of South Brittany (Lorient and Vilaine Bays) with Particular Reference to the Effect of the Turbidity Maxima

The south Brittany coast of France offers considerable opportunities to study the transition between estuaries and the sea. The sea coast is semi-enclosed, due to the local geomorphology, with sheltered conditions. The infaunal macrobenthic communities have been identified and characterized by species richness (S), abundance (A) and biomass (B), related to environmental factors. The spatial succession of infauna along the estuarine gradient shows the presence of an area with a much reduced fauna, and an oligobiotic area with low values of S, A and B. In these strongly stratified estuaries, concentration of organic matter and sediment can result in oxygen depletion, which is not favourable for survival of the infauna. The oligobiotic area seems to be related to the presence of the turbidity maximum, the location and importance of which are functions of river discharge and hydrodynamic mixing. The turbidity maximum which causes these conditions is shown to be at a higher salinity than has previously been described for this phenomenon, and is shown to break the spatial gradient of marine to brackish communities which has been described in larger estuaries.     

Improving fisheries management in New Zealand: Developing dialogue between fisheries science and management (FSM) and ecosystem science and management (ESM)

This paper attempts to overcome the dichotomy between the broadly different and largely separate fisheries science and management (FSM) and ecosystem science and management (ESM) knowledge systems that characterise the international literature and are found in fisheries management practice in different countries. The paper argues that the construction of a heuristic we term the fisheries problematic, around issues and contexts, reveals the breadth of international fisheries management concerns and the variety of contexts in which these concerns are being faced. Adopting a political economy informed nature-society approach the paper considers ecological and socio-economic processes in their institutional settings in an attempt to shift from the either/or arguments around fish or ecosystems found in the FSM or ESM literatures to investigation that is grounded in understandings of the historically and geographically specific trajectories of fisheries related interactions and understandings of how knowledge about the trajectories and their interactions is fashioned. Drawing on recent conceptual innovations in the field, the paper develops a matrix-centred approach to explore ecological, industry, community and policy domains in New Zealand’s Quota Management System (QMS) and Individual Transferable Quota (ITQ) fisheries management regime. The extended framework prioritises scrutiny of the interaction amongst the four domains, as a strategy to help develop institutional frameworks that facilitate behaviours that are societally inclusive. The paper offers three conclusions. First, the landscape of New Zealand fisheries issues is very much a product of the contingent interaction of the QMS, a management regime designed around the principles of a FSM approach and laid down in a neo-liberal political environment and Maori aspirations encompassing the fisheries sector. Second, the conceptual mapping of FSM and ESM perspectives over New Zealand’s fisheries management experience highlights that a number of management issues have been down played by the commitment to FSM, a situation that has led to on-going tensions between commercial, recreational and customary stakeholders regarding fisheries management. Put another way, there is more to running a sustainable fishery (as defined in the Fisheries Act 1996) than QMS and other tools and dialogue about the development of these should be a priority. Third and more generally, improved dialogue on fisheries questions is likely to be most expeditiously advanced by studies that explicitly conceptualise and contextualise ecological and socio-economic processes and their institutional arrangements.     

Rift propagation at craton margin.: Distribution of faulting and volcanism in the North Tanzanian Divergence (East Africa) during Neogene times

A revised kinematic model is proposed for the Neogene tectono-magmatic development of the North Tanzanian Divergence where the axial valley in S Kenya splits southwards into a wide diverging pattern of block faulting in association with the disappearance of volcanism. Propagation of rifting along the S Kenya proto-rift during the last 8 Ma is first assumed to have operated by linkage of discrete magmatic cells as far S as the Ngorongoro–Kilimanjaro transverse volcanic belt that follows the margin of cratonic blocks in N Tanzania. Strain is believed to have nucleated throughout the thermally-weakened lithosphere in the transverse volcanic belt that might have later linked the S Kenya and N Tanzania rift segments with marked structural changes along-strike. The North Tanzanian Divergence is now regarded as a two-armed rift pattern involving: (1) a wide domain of tilted fault blocks to the W (Mbulu) that encompasses the Eyasi and Manyara fault systems, in direct continuation with the Natron northern trough. The reactivation of basement fabrics in the cold and intact Precambrian lithosphere in the Mbulu domain resulted in an oblique rift pattern that contrasts with the orthogonal extension that prevailed in the Magadi–Natron trough above a more attenuated lithosphere. (2) To the E, the Pangani horst-like range is thought to be a younger (< 1 Ma) structure that formed in response to the relocation of extension S of the Kilimanjaro magmatic center. A significant contrast in the mechanical behaviour of the stretched lithosphere in the North Tanzanian diverging rift is assumed to have occurred on both sides of the Masai cratonic block with a mid-crustal decoupling level to the W where asymmetrical fault-basin patterns are dominant (Magadi–Natron and Mbulu), whereas a component of dynamical uplift is suspected to have caused the topographic elevation of the Pangani range in relation with possible far-travelled mantle melts produced at depth further N.     

Evolution of Mesozoic fluvial systems along the SE flank of the West Siberian Basin, Russia

The Mesozoic stratigraphy in the subsurface of the West Siberian Basin contains prolific hydrocarbon accumulations, and thus the depositional environments of marine and marginal marine Jurassic and Cretaceous age sediments are well-established. However, no information is currently available on strata of equivalent age that crop out along the SE basin margin in the Mariinsk–Krasnoyarsk region, despite the potential of these exposures to supply important information on the sediment supply routes into the main basin. Detailed sedimentological analysis of Jurassic–Cretaceous clastic sediments, in conjunction with palaeo-botanical data, reveals five facies associations that reflect deposition in a range of continental environments. These include sediments that were deposited in braided river systems, which were best developed in the Early Jurassic. These early river systems infilled the relics of a topography that was possibly inherited from earlier Triassic rifting. More mature fluvial land systems evolved in the Mid to Late Jurassic. By the Mid Jurassic, well-defined overbank areas had become established, channel abandonment was commonplace, and mudrocks were deposited on floodplains. Coal deposition occurred in mires, which were subject to periodic incursions by crevasse splay processes. Cretaceous sedimentation saw a renewed influx of sand-grade sediment into the region. It is proposed that landscape evolution throughout the Jurassic was driven simply by peneplanation rather than tectonic processes. By contrast, the influx of sandstones in the Cretaceous is tentatively linked to hinterland rejuvenation/ tectonic uplift, possibly coeval with the growth of large deltaic clinoform complexes of the Neocomian in the basin subsurface.     

Segment,Linkage,and Extensional Fault-Related Fold in Western Liaodong Bay Subbasin,Northeastern Bohai Sea,China

The western Liaodong (辽东) Bay subbasin displays examples of segment,linkage of extensional fault,and fault-related folds.The Liaoxi (辽西) extensional fault system consists of a series of NNE- and NE-trending segments that were linked through relay ramps.The fault hanging walls are characterized by a series of en echelon synclines with axial traces sub-parallel to the faults.The synclines are doubly plunging located on the hanging wall of normal faults,with the strata dip sub-parallel to the fault.These folds result from along-strike displacement variations of the individual fault segments,as well as from extensional fault-related folding.In the study area,the synclines are separated by transverse intra-basin highs and relay ramps that formed where segment linkage occurred.These hanging wall synclines and their relation to fault displacement variations indicate that they are formed by extensional fault-related fold.     

Feasibility Analysis of the Use of Binary Genetic Algorithms as Importance Samplers Application to a 1-D DC Resistivity Inverse Problem

DC resistivity inverse problems are ill-posed. For the Vertical Electrical Sounding method the acceptable solutions lie in very narrow elongated-shape regions in the model space. To characterize this ensemble of solutions is a central question. In a Bayesian framework this issue is solved adopting as solution the so-called model posterior probability distribution. However, due to the nonlinearity of the problem, this distribution is not explicitly known, or it is difficult to calculate. Therefore, algorithms that efficiently sample the model space according to it (importance sampling) are very desirable. The main goal of this paper is to numerically explore the performance of binary genetic algorithms as posterior importance sampling strategies. Their behavior will be firstly analyzed using 2D synthetic posterior test functions bearing the relevant properties of the real geo-electrical inverse problem. The conclusions will be again checked through the histogram reconstruction of parameters in a synthetic VES case, and eventually, in a real, higher dimensional, sea-water coastal intrusion problem, by comparing the results with those obtained with a theoretically correct Metropolis-Hasting importance sampler (simulated annealing without cooling). Percentile curves are introduced as an effective tool for risk assessment. We show that binary genetic algorithms perform well under very general assumptions. When the roulette wheel is the selection method used, mutation is over 10%, and the algorithm does not incorporate elitism. The results do not depend on the values of the remaining tuning parameters. Finally, to improve the efficiency of the sampling strategy, we introduce a binary genetic algorithm with oriented search space. This is done with the help of linearization of the forward operator and singular value decomposition around the maximum posterior estimate. It is shown, also, that the logarithmic model parameterization is adequate for this task.