Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.     

Isolation of protein from Chlorella sorokiniana CY1 using liquid biphasic flotation assisted with sonication through sugaring-out effect

Microalgae, a sustainable source of multi beneficial components has been discovered and could be utilised in pharmaceutical, bioenergy and food applications. This study aims to investigate the sugaring-out effect on the recovery of protein from wet green microalga, Chlorella sorokiniana CY1 which was assisted with sonication. A comparison of monosaccharides and disaccharides as one of the phaseforming constituents shows that the monosaccharides, glucose was the most suitable sugar in forming the phases with acetonitrile to enhance the production of protein(52% of protein). The protein productivity of microalgae was found to be significantly influenced by the volume ratio of both phases, as the yield of protein increased to 77%. The interval time between the sonication as well as the sonication modes were influencing the protein productivity as well. The optimum protein productivity was obtained with 10 s of resting time in between sonication. Pulse mode of sonication was suitable to break down the cell wall of microalgae compared to continuous mode as a lower protein yield was obtained with the application of continuous mode. The optimum condition for protein extraction were found as followed: 200 g/L glucose as bottom phase with volume ratio of I:1.25, 10 s of resting time for ultrasonication, 5 s of ultrasonication in pulse mode and 0.25 g of biomass weight. The high yield of protein about 81% could be obtained from microalgae which demonstrates the potential of this source and expected to play an important role in the future.     

Astronomical calibration of the Matuyama-Brunhes boundary: Consequences for magnetic remanence acquisition in marine carbonates and the Asian loess sequences

We have compiled 19 records from marine carbonate cores in which the Matuyama-Brunhes boundary (MBB) has been reasonably well constrained within the astronomically forced stratigraphic framework using oxygen isotopes. By correlation of the δ¹⁸O data to a timescale based on astronomical forcing, we estimate astronomical ages for each of the MBB horizons. In all but one record the MBB occurs within Stage 19.

Most magnetostratigraphic sections in Asian Loess place the MBB within a loess interval. Since loess deposition is presumed to be associated with glacial intervals, loess horizons should correspond to even-numbered oxygen isotope stages. A glacial age for the MBB is at odds with the results presented here, which firmly place the MBB within interglacial Stage 19. Inconsistency among the many loess sections and between the loess and the marine records suggests that the magnetic interpretation of loess sections may be more complicated than hitherto supposed.

The mean of the Stage 19 age estimates for the MBB is 777.9 ± 1.8 (N = 18). Inclusion of the single Stage 20 age results in a mean of 778.8 ± 2.5 (N = 19). The astronomical age estimate of the MBB compares favorably with an (unweighted) mean of 778.2 ± 3.5 (N = 10) from a compilation of ⁴⁰Ar/39^Ar results of transitional lava flows. Combining the two independent data sets yields a grand mean of 778.0 ± 1.7 (N = 28).

The new compilation shows virtually no trend in placement of the MBB within isotope Stage 19 as a function of sediment accumulation rate. We interpret this to mean that the average depth of remanence acquisition is within a few centimeters of the sediment-water interface.

Separating the cores into two geographic regions (an Indo-Pacific-Caribbean [IPC] Group and an Atlantic Group) results in a significant difference in the position of the mid-point of the reversal with respect to the astronomical time scale. The data presented here suggest a difference of several thousand years between the two regions. This observation could be caused by systematic differences between the two regions in sedimentation rate within the interval of interest, systematic differences in remanence acquisition, or by genuine differences in the timing of the directional changes between the two regions.     

Rapid least-squares inversion of apparent resistivity pseudosections by a quasi-Newton method1

A fast inversion technique for the interpretation of data from resistivity tomography surveys has been developed for operation on a microcomputer. This technique is based on the smoothness-constrained least-squares method and it produces a two-dimensional subsurface model from the apparent resistivity pseudosection. In the first iteration, a homogeneous earth model is used as the starting model for which the apparent resistivity partial derivative values can be calculated analytically. For subsequent iterations, a quasi-Newton method is used to estimate the partial derivatives which reduces the computer time and memory space required by about eight and twelve times, respectively, compared to the conventional least-squares method. Tests with a variety of computer models and data from field surveys show that this technique is insensitive to random noise and converges rapidly. This technique takes about one minute to invert a single data set on an 80486DX microcomputer.     

Laboratory sulfurisation of the marine microalga Nannochloropsis salina

To understand more fully the mode of preservation of organic matter in marine sediments, laboratory sulfurisation of intact cells of the cultured microalga Nannochloropsis salina was performed using inorganic polysulfides in seawater at 50°C. Solvent extractable and non-extractable materials were analysed by GC–MS and Py–GC–MS, respectively, to study the incorporation of sulfur into the microalgal organic matter. No GC-amenable sulfur-containing compounds were found in the extracts apart from some minor thiophenes with a phytanyl carbon skeleton. The residue after extraction and hydrolysis contained abundant macromolecular sulfur-containing moieties as revealed by the presence of dominant C₂₈–C₃₂ thiols, thiophenes, thianes and thiolanes in the flash pyrolysates. These products are thought to be formed from moieties derived from sulfurisation of C₂₈–C₃₂ diols and alkenols, characteristic lipids of N. salina. C₁–C₂ alkylated thiophenes were also found in the pyrolysates and probably result from moieties formed upon sulfurisation of carbohydrates. The highly resistant biomacromolecule (algaenan) synthesised by N. salina remains unaffected by sulfurisation. The non-hydrolysable residue isolated from the sulfurised N. salina thus comprises algaenan and (poly)sulfide-bound long alkyl chains. The sulfurisation experiments show that both selective preservation of algaenans and lipid and carbohydrate “vulcanisation” can be involved in the preservation of algal organic matter in marine environments.     

Analysis of Cassini/CIRS limb spectra of Titan acquired during the nominal mission II: Aerosol extinction profiles in the 600–1420 cm spectral range

We have analyzed the continuum emission of limb spectra acquired by the Cassini/CIRS infrared spectrometer in order to derive information on haze extinction in the 3–0.02 mbar range (∼150–350 km). We focused on the 600–1420 cm⁻¹ spectral range and studied nine different limb observations acquired during the Cassini nominal mission at 55°S, 20°S, 5°N, 30°N, 40°N, 45°N, 55°N, 70°N and 80°N. By means of an inversion algorithm solving the radiative transfer equation, we derived the vertical profiles of haze extinction coefficients from 17 spectral ranges of 20-cm⁻¹ wide at each of the nine latitudes. At a given latitude, all extinction vertical profiles retrieved from various spectral intervals between 600 and 1120 cm⁻¹ display similar vertical slopes implying similar spectral characteristics of the material at all altitudes. We calculated a mean vertical extinction profile for each latitude and derived the ratio of the haze scale height (H_haze) to the pressure scale height (H_gas) as a function of altitude. We inferred H_haze/H_gas values varying from 0.8 to 2.4. The aerosol scale height varies with altitude and also with latitude. Overall, the haze extinction does not show strong latitudinal variations but, at 1 mbar, an increase by a factor of 1.5 is observed at the north pole compared to high southern latitudes. The vertical optical depths at 0.5 and 1.7 mbar increase from 55°S to 5°N, remain constant between 5°N and 30°N and display little variation at higher latitudes, except the presence of a slight local maximum at 45°N. The spectral dependence of the haze vertical optical depth is uniform with latitude and displays three main spectral features centered at 630 cm⁻¹, 745 cm⁻¹ and 1390 cm⁻¹, the latter showing a wide tail extending down to ∼1000 cm⁻¹. From 600 to 750 cm⁻¹, the optical depth increases by a factor of 3 in contrast with the absorbance of laboratory tholins, which is generally constant. We derived the mass mixing ratio profiles of haze at the nine latitudes. Below the 0.4-mbar level all mass mixing ratio profiles increase with height. Above this pressure level, the profiles at 40°N, 45°N, 55°N, at the edge of the polar vortex, display a decrease-with-height whereas the other profiles increase. The global increase with height of the haze mass mixing ratio suggest a source at high altitudes and a sink at low altitudes. An enrichment of haze is observed at 0.1 mbar around the equator, which could be due to a more efficient photochemistry because of the strongest insolation there or an accumulation of haze due to a balance between sedimentation and upward vertical drag.     

Enhancing the relevance of Shared Socioeconomic Pathways for climate change impacts,adaptation and vulnerability research

This paper discusses the role and relevance of the shared socioeconomic pathways (SSPs) and the new scenarios that combine SSPs with representative concentration pathways (RCPs) for climate change impacts, adaptation, and vulnerability (IAV) research. It first provides an overview of uses of social–environmental scenarios in IAV studies and identifies the main shortcomings of earlier such scenarios. Second, the paper elaborates on two aspects of the SSPs and new scenarios that would improve their usefulness for IAV studies compared to earlier scenario sets: (i) enhancing their applicability while retaining coherence across spatial scales, and (ii) adding indicators of importance for projecting vulnerability. The paper therefore presents an agenda for future research, recommending that SSPs incorporate not only the standard variables of population and gross domestic product, but also indicators such as income distribution, spatial population, human health and governance.