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.     

Invariance and scaling properties in the distributions of contributing area and energy in drainage basins

The cumulative probability distributions for stream order, stream length, contributing area, and energy dissipation per unit length of channel are derived, for an ordered drainage system, from Horton's laws of network composition. It is shown how these distributions can be related to the fractal nature of single rivers and river networks. Finally, it is shown that the structure proposed here for these probability distributions is able to fit the observed frequency distributions, and their deviations from straight lines in a log-log plot.     

The solubility of iron sulphides in synthetic and natural waters at ambient temperature

A critical evaluation of literature values for the solubility products, K _sp ^NBS = [Fe²⁺][HS^–] Fe²⁺ HS^– (H _NBS ⁺ )^–1, of various iron sulphide phases results in consensus values for the pKs of 2.95 ± 0.1 for amorphous ferrous sulphide, 3.6 ± 0.2 for mackinawite, 4.4 ± 0.1 for greigite, 5.1 ± 0.1 for pyrrhotite, 5.25 ± 0.2 for troilite and 16.4 ± 1.2 for pyrite.Where the analogous ion activity products have been measured in anoxic freshwaters in which there is evidence for the presence of solid phase FeS, the values lie within the range of 2.6–3.22, indicating that amorphous iron sulphide is the controlling phase. The single value for a groundwater of 2.65 (2.98 considering carbonate complexation) agrees. In seawater four values range between 3.85 to 4.2, indicating that mackinawite or greigite may be the controlling phase. The single low value of 2.94 is in a situation where particularly high fluxes of Fe (II) and S (–II) may result in the preferential precipitation of amorphous iron sulphide. Formation of framboidal pyrite in these sulphidic environments may occur in micro-niches and does not appear to influence bulk concentrations. Calculations show that the formation of Fe₂S₂ species probably accounts for very little of the iron or sulphide in most natural waters. Previously reported stability constants for the formation of Fe (HS)₂ and (Fe (HS)₃)^– are shown to be suspect, and these species are also thought to be negligible in natural waters. In completely anoxic pore waters polysulphides also have a negligible effect on speciation, but in tidal sediments they may reach appreciable concentrations and lead to the direct formation of pyrite. Concentrations of iron and sulphide in pore waters can be controlled by the more soluble iron sulphide phase. The change in the IAP with depth within the sediment may reflect ageing of the solid phase or a greater flux of Fe (II) and S (–II) nearer the sediment surface. This possible kinetic influence on the value of IAPs has implications for their use in geochemical studies involving phase formation.     

Quaternary fault segmentation and interaction in the epicentral area of the 1561 earthquake (Mw = 6.4), Vallo di Diano, southern Apennines, Italy

The main structural characteristics of the Caggiano and Polla faults, exposed in the epicentral area of the 1561 earthquake (Mw = 6.4), southern Italy, have been investigated in detail to assess their spatial and temporal properties, and to evaluate their seismogenic potential. These right stepping normal faults show an overlap of about 7 km and an across strike separation of about 4 km. The geometric relationships between the Caggiano and Polla faults, but also the displacement distribution along each fault, demonstrate that they have been strongly interacting throughout the Pleistocene. Nevertheless, geological evidence of Holocene tectonic activity was mainly recognized along the Caggiano Fault (faulted late glacial deposits) and in the southernmost part of the Polla Fault (faulted deposits of probably Late Pleistocene age). This suggests that the Caggiano Fault can be considered as the most tectonically active fault in the Vallo di Diano Fault System. By calculating Coulomb stress changes, we have constrained modes of mechanical interactions between the two faults in a scenario compatible with the 1561 earthquake. This approach allows us to argue that both the Caggiano and the Polla Faults are probably linked at depth, and part of the same seismogenic structure which may be potentially responsible for composite ruptures with magnitude ≥ 6.5.     

Revised age of proximal deposits in the Zagros foreland basin and implications for Cenozoic evolution of the High Zagros

The regionally extensive, coarse-grained Bakhtiyari Formation represents the youngest synorogenic fill in the Zagros foreland basin of Iran. The Bakhtiyari is present throughout the Zagros fold-thrust belt and consists of conglomerate with subordinate sandstone and marl. The formation is up to 3000 m thick and was deposited in foredeep and wedge-top depocenters flanked by fold-thrust structures. Although the Bakhtiyari concordantly overlies Miocene deposits in foreland regions, an angular unconformity above tilted Paleozoic to Miocene rocks is expressed in the hinterland (High Zagros).

The Bakhtiyari Formation has been widely considered to be a regional sheet of Pliocene–Pleistocene conglomerate deposited during and after major late Miocene–Pliocene shortening. It is further believed that rapid fold growth and Bakhtiyari deposition commenced simultaneously across the fold-thrust belt, with limited migration from hinterland (NE) to foreland (SW). Thus, the Bakhtiyari is generally interpreted as an unmistakable time indicator for shortening and surface uplift across the Zagros. However, new structural and stratigraphic data show that the most-proximal Bakhtiyari exposures, in the High Zagros south of Shahr-kord, were deposited during the early Miocene and probably Oligocene. In this locality, a coarse-grained Bakhtiyari succession several hundred meters thick contains gray marl, limestone, and sandstone with diagnostic marine pelecypod, gastropod, coral, and coralline algae fossils. Foraminiferal and palynological species indicate deposition during early Miocene time. However, the lower Miocene marine interval lies in angular unconformity above ~ 150 m of Bakhtiyari conglomerate that, in turn, unconformably caps an Oligocene marine sequence. These relationships attest to syndepositional deformation and suggest that the oldest Bakhtiyari conglomerate could be Oligocene in age.

The new age information constrains the timing of initial foreland-basin development and proximal Bakhtiyari deposition in the Zagros hinterland. These findings reveal that structural evolution of the High Zagros was underway by early Miocene and probably Oligocene time, earlier than commonly envisioned. The age of the Bakhtiyari Formation in the High Zagros contrasts significantly with the Pliocene–Quaternary Bakhtiyari deposits near the modern deformation front, suggesting a long-term (> 20 Myr) advance of deformation toward the foreland.     

沉积盆地有机质自由基热演化特征及其作为古温标的探索

恢复沉积盆地热历史的方法主要有古地温指标法和动力学模拟法两类.古温标法中最常用的为有机质成熟度和矿物裂变径迹,动力学模拟法中最常用的是盆地的拉张模型和挤压模型.本文针对Ⅰ型有机质自由基浓度作为沉积盆地古温标进行了探索研究.分析了自然演化系列Ⅰ型有机质的自由基的热演化特征;根据自然演化系列的Ⅰ型有机质自由基浓度及时间-温度指数(TTI)数据,对热模拟实验下的Ⅰ型有机质自由基浓度及时间-温度指数(TTI)值进行校正,由校正后的数据初步建立了Ⅰ型有机质的自由基浓度(Ng)与时间-温度指数(TTI)的定量模型.     

伶仃洋沉积动力特点的研究

位于珠江三角洲东侧的伶仃洋，因径流下泄与潮流进退的流向不一，使它各分流口门的出口水道都有主槽和支槽之分，即都有主干水道和分汊水道。陆架高盐海水入侵又使伶仃洋内沉积动力过程在空间分布上发生差异，如沉积物分布有粗－细－稍粗之分；而水体中的密度、速度差异，常常产生锋带，对水下地形的发展有不可忽视的影响。因而全面认识发生在伶仃洋内的沉积动力作用，对深水航道的选线极为重要。     

广西主要港湾余流特征及其对物质输运的影响

广西5个主要海湾各有自己的环流特征：铁山湾存在一个反气旋涡，廉州湾存在一个海湾尺度的气旋环流，钦州湾青菜头南端存在一个气旋式环流，防城湾口有一个气旋环流。珍珠湾余流，在深水航道内自湾内指向湾外。冬季，受北风影响，各湾从湾内流向湾外的余流分量加强。由于反气旋或气旋环流存在，形成和环流相应的泥沙输运。     

渤海湾环流的一次观测和分析

为了对渤海湾的环流形态有一个准确的认识，对渤海湾进行了一次较为全面的海流调查。这次调查以Ｌａｇｒａｎｇｅ观点为指导，采用３种方法：海流计在固定站位连续测流；漂流板轨迹观测和“人工水母”测底层流。本文分别介绍这三种方法及观测过程，并分析得出了渤海湾Ｌａｇｒａｎｇｅ余环流型大致是顺时针的结论