Numerical Analysis of Space–Time Finite Element Formulations for Miscible Displacements

A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates.     

Carbon budget for a subtropical seagrass dominated coastal lagoon: How important are seagrasses to total ecosystem net primary production?

It has been assumed that because seagrasses dominate macrophyte biomass in many estuaries they also dominate primary production. We tested this assumption by developing three carbon budgets to examine the contribution of autotrophic components to the total ecosystem net primary production (TENPP) of Lower Laguna Madre, Texas. The first budget coupled average photosynthetic parameters with average daily irradiance to calculate daily production. The second budget used average photosynthetic parameters and hourly in situ irradiance to estimate productivity. The third budget integrated temperature-adjusted photosynthetic parameters (using Q₁₀=2) and hourly in situ irradiance to estimate productivity. For each budget TENPP was calculated by integrating production from each autotroph based on the producers’ areal distribution within the entire Lower Laguna Madre. All budgets indicated that macroalgae account for 33–42% of TENPP and seagrasses consistently accounted for about 33–38%. The contribution by phytoplankton was consistently about 15–20%, and the contribution from the benthic microalgae varied between 8% and 36% of TENPP, although this may have been underestimated due to our exclusion of the within bed microphytobenthos component. The water column over the seagrass beds was net heterotrophic and consequently was a carbon sink consuming between 5% and 22% of TENPP, TENPP ranged between 5.41×10¹⁰ and 2.53×10¹¹ g C yr⁻¹, depending on which budget was used. The simplest, most idealized budget predicted the highest TENPP, while the more realistic budgets predicted lower values. Annual production rates estimated using the third budget forHalodule urightii andThalassia testudinum compare well with field data. Macroalgae and microalgae contribute 50–60% of TENPP, and seagrass may be more important as three-dimensional habitat (i.e., structure) than as a source of organic carbon to the water column in Lower Laguna Madre.     

Sub-wavebase cross-bedded grainstones on a distally steepened carbonate ramp, Upper Miocene, Menorca, Spain

Cross‐bedded grainstones on carbonate ramps and shelves are commonly related to the locus of major wave energy absorption such as shorelines, shoals or shelf breaks. In contrast, on the Early Tortonian carbonate platform of Menorca (Balearic Islands), coarse‐grained, cross‐bedded grainstones are found at a distance from the palaeoshoreline where they were deposited below the wavebase. Excellent exposures along continuous outcrops on the sea cliffs of Menorca reveal the depositional profile and three‐dimensional distribution of the different facies belts of the Tortonian ramp depositional system. Basinward from the palaeoshoreline, fan deltas and beach deposits pass into 5‐km‐wide gently dipping bioturbated dolopackstone (inner and middle ramp), then into 12–20°‐dipping dolograinstone/rudstone clinobeds (ramp slope) and, finally, into subhorizontal fine‐grained basinal dolowackestone to dolopackstone (outer ramp). In this Miocene example, coarse‐grained grainstones exist in five different settings other than beach deposits: (1) on the middle ramp, where cross‐bedded grainstones were deposited by currents roughly parallel to the shoreline at 40–70 m estimated water depth and are interbedded with gently dipping bioturbated dolomitized packstones; (2) on the upper slope, where clinobeds are composed mostly of in situ rhodoliths and red‐algae fragments; (3) on the lower slope, as small‐scale bedforms (small three‐dimensional subaqueous dunes) migrating parallel to the slope; (4) at the transition between the lower slope and the outer ramp, where mollusc‐rich and rhodolithic rudstones and grainstones, interbedded in dolomitized laminated wackestones containing abundant planktonic foraminifera, infill slide/slump scars as upslope‐backstepping bodies (backsets); (5) at the toe of the slope, where coarse skeletal grainstones indicate bedform migration parallel to the platform margin, induced by currents at more than 150 m estimated water depth. This Late Miocene example also illustrates how changes in intrabasinal environmental conditions (nutrients and/or temperature) may produce changes in stratal patterns and facies architecture if they affect the biological system. Two depositional sequences compose the Miocene platform on Menorca, where a reef‐rimmed platform prograded onto an earlier distally steepened ramp. The transition from the ramp to the reef‐rimmed platform was effected by an increase in accommodation space caused by ecological changes, promoting a shift from a grain‐ to a framework‐producing biota.     

Cretaceous and Tertiary terrane accretion in the Cordillera Occidental of the Andes of Ecuador

New field, geochronological, geochemical and biostratigraphical data indicate that the central and northern parts of the Cordillera Occidental of the Andes of Ecuador comprise two terranes. The older (Pallatanga) terrane consists of an early to late (?) Cretaceous oceanic plateau suite, late Cretaceous marine turbidites derived from an unknown basaltic to andesitic volcanic source, and a tectonic mélange of probable late Cretaceous age. The younger (Macuchi) terrane consists of a volcanosedimentary island arc sequence, derived from a basaltic to andesitic source. A previously unidentified, regionally important dextral shear zone named the Chimbo-Toachi shear zone separates the two terranes. Regional evidence suggests that the Pallatanga terrane was accreted to the continental margin (the already accreted Cordillera Real) in Campanian times, producing a tectonic mélange in the suture zone. The Macuchi terrane was accreted to the Pallatanga terrane along the Chimbo-Toachi shear zone during the late Eocene, probably in a dextral shear regime. The correlation of Cretaceous rocks and accretionary events in the Cordillera Occidental of Ecuador and Colombia remains problematical, but the late Eocene event is recognised along the northern Andean margin.     

Use of Chemotaxonomy To Study the Influence of Benzalkonium Chloride on Bacterial Populations in Biodegradation Testing

Contradictory results are reported for the behaviour of quaternary ammonium compounds (QACs) in sewage treatment plants (STPs). QACs may sorb onto activated sludge. Only little information is available with respect to effects of QACs against bacteria in STPs. Only 5 to 15 % of bacteria present in sewage sludge can be detected by means of culture dependent microbiological methods. The shift of the bacterial populations due to effects of test compounds have not been studied up to now with culture independent methods. The microbial populations shift was studied in situ using culture independent chemotaxonomy profiling ubiquinones and polyamines. Additionally, toxic effects of QACs against bacteria present in the test vessels of the Zahn‐Wellens test (OECD 302 B) were assessed with a toxicity control in the test. The ubiquinone profiles representing changes in Gram‐negative populations mainly showed that the activated sludge was affected only in test vessels containing benzalkonium chloride. According to chemotaxonomy Acinetobacter or/and some members of Pseudomonas spp. have been selected by benzalkonium chloride after some adaptation period (8 to 12 days).     

Sediment yield and storage variations in the Négron River catchment (south western Parisian Basin,France) during the Holocene period

In the Négron River catchment area (162 km²), surface‐sediment stores are composed of periglacial calcareous ‘grèze’ (5 × 10⁶ t) and loess (21 × 10⁶ t), and Holocene alluvium (12·6 × 10⁶ t), peat (0·6 × 10⁶ t) and colluvium (18·5 × 10⁶ t). Seventy‐five per cent of the Holocene sediments is stored along the thalwegs. Present net sediment yield, calculated from solid discharge at the Négron outlet, is low (0·6 t km^?2 a^?1) due to the dominance of carbonate rocks in the catchment. Mean sediment yield during the Holocene period is 7·0 t km^?2 a^?1 from alluvium stores and 7·6 t km^?2 a^?1 from colluvium stores. Thus, the gross sediment yield during the Holocene period is about 18·7 t km^?2 a^?1 and the sediment delivery ratio 3 per cent. The yield considerably varies from one sub‐basin to another (3·9 to 24·5 t km^?2 a^?1) according to lithology: about 25 per cent and 50 per cent of initial stores of periglacial grèze and loess respectively were reworked during the Holocene period. Sediment yield has increased by a factor of 6 in the last 1000 years, due to the development of agriculture. The very high rate of sediment storage on the slope during that period (88 per cent of the yield) can be accounted for by the formation of cultivation steps (‘rideaux’). It is predicted that the current destruction of these steps will result in a sediment wave reaching the valley floors in the coming decades. Subboreal and Subatlantic sediments and pollen assemblages in the Taligny marsh, where one‐third of the alluvium is stored, show the predominant influence of human activity during these periods in the Négron catchment. Copyright © 2002 John Wiley & Sons, Ltd.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 John Wiley & Sons, Ltd.     

Inelastic displacement ratios for evaluation of existing structures

Results of a detailed statistical study of constant relative strength inelastic displacement ratios to estimate maximum lateral inelastic displacement demands on existing structures from maximum lateral elastic displacement demands are presented. These ratios were computed for single‐degree‐of‐freedom systems with different levels of lateral strength normalized to the strength required to remain elastic when subjected to a relatively large ensemble of recorded earthquake ground motions. Three groups of soil conditions with shear wave velocities higher than 180m/s are considered. The influence of period of vibration, level of lateral yielding strength, site conditions, earthquake magnitude, distance to the source, and strain‐hardening ratio are evaluated and discussed. Mean inelastic displacement ratios and those associated with various percentiles are presented. A special emphasis is given to the dispersion of these ratios. It is concluded that distance to the source has a negligible influence on constant relative strength inelastic displacement ratios. However, for periods smaller than 1s earthquake magnitude and soil conditions have a moderate influence on these ratios. Strain hardening decreases maximum inelastic displacement at a fairly constant rate depending on the level of relative strength for periods of vibration longer than about 1.0s while it decreases maximum inelastic displacement non‐linearly as the period of vibration shortens and as the relative‐strength ratio increases for periods of vibration shorter than 1.0s. Finally, results from non‐linear regression analyses are presented that provide a simplified expression to be used to approximate mean inelastic displacement ratios during the evaluation of existing structures built on firm sites. Copyright © 2003 John Wiley & Sons, Ltd.     

Landslide susceptibility analysis using GIS and artificial neural network

The purpose of this study is to develop landslide susceptibility analysis techniques using an arti?cial neural network and to apply the newly developed techniques to the study area of Yongin in Korea. Landslide locations were identi?ed in the study area from interpretation of aerial photographs, ?eld survey data, and a spatial database of the topography, soil type and timber cover. The landslide‐related factors (slope, curvature, soil texture, soil drainage, soil effective thickness, timber age, and timber diameter) were extracted from the spatial database. Using those factors, landslide susceptibility was analysed by arti?cial neural network methods. The landslide susceptibility index was calculated by the back‐propagation method, which is a type of arti?cial neural network method, and the susceptibility map was made with a geographic information system (GIS) program. The results of the landslide susceptibility analysis were veri?ed using landslide location data. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide location. A GIS was used to ef?ciently analyse the vast amount of data, and an arti?cial neural network to be an effective tool to maintain precision and accuracy. The results can be used to reduce hazards associated with landslides and to plan land use and construction. Copyright © 2003 John Wiley & Sons, Ltd.