A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back on to the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimize the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.     

The effect of geometry and tidal forcing on hydrodynamics and net sediment transport in semi-enclosed tidal basins

A new depth-averaged exploratory model has been developed to investigate the hydrodynamics and the tidally averaged sediment transport in a semi-enclosed tidal basin. This model comprises the two-dimensional (2DH) dynamics in a tidal basin that consists of a channel of arbitrary length, flanked by tidal flats, in which the water motion is being driven by an asymmetric tidal forcing at the seaward side. The equations are discretized in space by means of the finite element method and solved in the frequency domain. In this study, the lateral variations of the tidal asymmetry and the tidally averaged sediment transport are analyzed, as well as their sensitivity to changes in basin geometry and external overtides. The Coriolis force is taken into account. It is found that the length of the tidal basin and, to a lesser extent, the tidal flat area and the convergence length determine the behaviour of the tidally averaged velocity and the overtides and consequently control the strength and the direction of the tidally averaged sediment transport. Furthermore, the externally prescribed overtides can have a major influence on tidal asymmetry in the basin, depending on their amplitude and phase. Finally, for sufficiently wide tidal basins, the Coriolis force generates significant lateral dynamics.     

The Banhad?o Alkaline Complex, Southeastern Brazil: source and evolution of potassic SiO2-undersaturated high-Ca and low-Ca magmatic series

The Cretaceous Banhad?o alkaline complex in southeastern Brazil presents two potassic SiO₂-undersaturated series. The high-Ca magmatic series consist of initially fractionated olivine (Fo_92-91) + diopside (Wo_48-43En_49-35Ae_0-7), as evidenced by the presence of xenocrysts and xenoliths. In that sequence, diopside (Wo_47-38En_46-37Ae_0-8) + phlogopite + apatite + perovskite (Prv_>92) crystallized to form the phlogopite melteigite and led to the Ca enrichment of the magma. Diopside (Wo_47-41En_32-24 Ae_3-14) continued to crystallize as an early mafic mineral, followed by nepheline (Ne_74.8-70.1Ks_26.3-21.2Qz_7.6-0.9) and leucite (Lc_65-56) and subsequently by melanite and potassic feldspar (Or_85-99Ab_1-7) to form melanite ijolites, wollastonite-melanite urtites and melanite-nepheline syenites. Melanite-pseudoleucite-nepheline syenites are interpreted to be a leucite accumulation. Melanite nephelinite dykes are believed to represent some of the magmatic differentiation steps. The low-Ca magmatic series is representative of a typical fractionation of aegirine-augite (Wo_36-29En_25-4Ae_39-18) + alkali feldspar (Or_57-96Ab_3-43) + nepheline (Ne_76.5-69.0Ks_19.9-14.4Qz_15.1-7.7) + titanite from phonolite magma. The evolution of this series from potassic nepheline syenites to sodic sodalite syenites and sodalitolites is attributed to an extensive fractionation of potassic feldspar, which led to an increase of the NaCl activity in the melt during the final stages forming sodalite-rich rocks. Phonolite dykes followed a similar evolutionary process and also registered some crustal assimilation. The mesocratic nepheline syenites showed interactions with phlogopite melteigites, such as compatible trace element enrichments and the presence of diopside xenocrysts, which were interpreted to be due to a mixing/mingling process of phonolite and nephelinite magmas. The geochemical data show higher TiO₂ and P₂O₅ contents and lower SiO₂ contents for the high-Ca series and different LILE evolution trends and REE chondrite-normalized patterns as compared to the low-Ca series. The ⁸⁷Sr/⁸⁶Sr, ¹⁴³Nd/¹⁴⁴Nd, ²⁰⁶Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb initial ratios for the high-Ca series (0.70407–0.70526, 0.51242–0.51251, 17.782–19.266 and 38.051–39.521, respectively) were slightly different from those of the low-Ca series (0.70542–0.70583, 0.51232–0.51240, 17.758–17.772 and 38.021–38.061, respectively). For both series, a CO₂-rich potassic metasomatized lithospheric mantle enriched the source with rutile-bearing phlogopite clinopyroxenite veins. Kamafugite-like parental magma is attributed to the high-Ca series with major contributions from the melting of the veins. Potassic nephelinite-like parental magma is assigned to the low-Ca series, where the metasomatized wall-rock played a more significant role in the melting process.     

Experimental investigation of water impact on axisymmetric bodies

The results of an elaborate experimental investigation on bottom slamming of axisymmetric objects are presented. Drop tests have been performed on a hemisphere and two conical shapes with different deadrise angles. The test setup is designed so as to prevent small rotations of the test objects which cause scatter in the measurement data. The pressure distribution and evolution as well as the body motion parameters are measured during impact. By means of a high speed camera the water uprise is visualized and the wetting factor is determined for the cones. The results are compared with a three-dimensional asymptotic theory for axisymmetric rigid bodies with constant entry velocity. The ratio between the registered peak pressures and the asymptotic theory are in accordance with comparable experiments in the literature. The asymptotic theory, however, is found to be quite conservative, since the measured peak pressure levels appear to be approximately 50% to 75% of the theoretical levels.     

通辽CK－9井水位动态及影响因素的分析

本文在对通辽ＣＫ－９井、水文地质条件和水动态特征初步分析的基础上，采用ＰＣ－１５００计算机，对该井的观测数据与气压、蒸发、降雨、开采、河水位、固体潮、地震波和潜水位等相互关系，进行了定量的相互分析．研究所得结果，对该井水位动态变化和影响因素的关系，提出了新的结论．     

A water retention model accounting for void ratio changes in double porosity clays

This paper presents a constitutive model that predicts the water retention behaviour of compacted clays with evolving bimodal pore size distributions. In line with previous research, the model differentiates between the water present inside the saturated pores of the clay aggregates (the microstructure) and the water present inside the pores between clay aggregates (the macrostructure). A new formulation is then introduced to account for the effect of the macrostructural porosity changes on the retention behaviour of the soil, which results in a consistent evolution of the air-entry value of suction with volumetric deformations. Data from wetting tests on three different active clays (i.e. MX-80 bentonite, FEBEX bentonite, and Boom clay), subjected to distinct mechanical restraints, were used to formulate, calibrate, and validate the proposed model. Results from free swelling tests were also modelled by using both the proposed double porosity model and a published single porosity model, which confirmed the improvement in the predictions of degree of saturation by the present approach. The proposed retention model might be applied, for example, to the simulation of the hydromechanical behaviour of engineered bentonite barriers in underground nuclear waste repositories, where compacted active clays are subjected to changes of both suction and porosity structure under restrained volume conditions.

     

江苏金坛地质遗迹特征及旅游开发建议

金坛市拥有丰富的地质遗迹资源，包括地层剖面、古生物化石产地与古文化遗址、地质构造形迹、地质地貌景观、水体景观、岩溶洞穴、矿产产地与观赏石景点8大类。对其应遵循原态保护的原则，以茅山风景名胜区为重点，因地制宜地开展有特色的大众地学科普旅游，为地方经济和社会发展作出贡献。     

Stress-induced viscosity changes and the existence of dominant wavelengths in folds

The theory of folding in stratified media presented by Biot and Ramberg has been extended by considering a more general type of material response. The model consists of a viscous layer embedded in a less viscous medium, compressed parallel to the layering. A transition from Newtonian to non-Newtonian behavior is considered and an approximate solution to the governing equation is discussed. The results give the effect of local, stress-induced changes in the viscosity on the profile of the fold. The results indicate that as the transition to non-Newtonian behavior proceeds, the wavelength selection process described by Biot and Ramberg breaks down; the wavelength of the fold which develops probably will not be the “dominant” wavelength defined by Biot.