Mathematical modeling of PCB bioaccumulation in Perna viridis

In the present work, we built a mathematical model of polychlorinated biphenyl (PCB) bioaccumulation in Perna viridis, namely, a one-compartment model with a time dependent incorporation rate R (μg g⁻¹ lipid per ppb water per day), with positive substrate cooperativity as the underlying physical mechanism. The temporal change of the PCB concentration Q (μg g⁻¹ lipid) in the soft tissues of the mussel depends on the competition of the input rate RW and the output rate kQ, where W is the concentration of PCB in water (ppb water) and k is the elimination rate (per day). From our experimental data, k=0.181±0.017 d⁻¹. The critical concentration in water W_c for positive substrate cooperativity was found to be 2.4 ppb. Below W_c, R is a constant. For a water concentration of 0.5 ppb Aroclor 1254, R=24.0±2.4 μg g⁻¹ lipid ppb⁻¹ d⁻¹. Above W_c, positive substrate cooperativity comes into effect and R becomes a function of time and dependent on the concentration Q in a form R=γQ/(Q+δ). This is the case for a water concentration of 5 ppb Aroclor 1254, where γ=15.1 μg g⁻¹ lipid ppb⁻¹ d¹ and δ≈200 μg g⁻¹ lipid. From this model, the uptake is exponentially increasing when the PCB concentration in the mussel is small compared to 200 μg g⁻¹ lipid, and hyperbolically increasing when the concentration is large compared to 200 μg g⁻¹ lipid, which are consistent with the experimental data. The model is useful for understanding the true processes taking place during the bioaccumulation and for risk assessment with higher confidence. Future experimental data which challenge the present model are anticipated and in fact desirable for improvement and perfection of the model.     

The French Guyana doleritic dykes: geochemical evidence of three populations and new data for the Jurassic Central Atlantic Magmatic Province

A petrographic and geochemical study of 15 Early Jurassic and 7 Proterozoic dolerites of French Guyana, and of one Jurassic dolerite from Ivory-Coast were carried out. The Early Jurassic SSW-NNE trending dykes have doleritic aphyric or gabbroic phyric texture. Their chemical compositions, slightly under-saturated to over-saturated, show moderate to low Mg-ratios (63–36), high TiO₂ contents (1.85–3.56 wt.%), weak rare earth element fractionation [1.8<(La/Yb)_n <4.6], negative Sr-anomalies (0.41<Sr* <0.91) and positive TiO₂-anomalies (1.2<TiO₂*<2.15). They correspond to sub-alkaline high-Ti basalts (HTiB). The La/Nb ratio (1.18±0.23) is close to that of the OIB-type Deccan basalts and incompatible element ratios: Rb/Nb_n (1.07±0.3), Zr/Y (4.6±0.8) are intermediate between those of HTiB in Deccan and those of HTiB in Paraná. These geochemical signatures and recent isotopic data of Bertrand et al. [Bertrand, H., Liegeois, J.P., Deckart, K., Féraud, G., 1999. High-Ti tholeiites in Guinea and Their Connection with the Central Atlantic CFB Province: Elemental and Nd-Sr-Pb Isotopic Evidence for Preferential Zone of Mantle Upwelling in Cause of Rifting. AGU spring meeting (Abst. p 317)] suggest that their magmatic source is different from that of the other basalts of the Central Atlantic Magmatic Province (CAMP). Such signatures are restricted to a central zone coinciding with the Panafrican Rokelides suture. We propose a model of sub-lithospheric preferential channelling of an asthenospheric ascent in this zone. Two other groups of dykes were identified in French Guyana. Compared to the Jurassic ones the Proterozoic dykes have NNW-SSE and E-W trending direction, more important LILE enrichment, low TiO₂ contents (<2 wt%) and Nb–Ta negative anomalies. Their calc-alkaline signature could be the result of a previous subduction and may be related to the 1800 Ma Venturi-Tapajós event, which contaminated the mantle source.     

Modelling of meteorological conditions at an urban scale for the PUMA winter campaign

For the winter 2000 campaign of the Pollution of Urban Midlands Atmosphere project, observation and numerical modelling of meteorological conditions over the West Midlands conurbation, UK, was undertaken. Modelling was performed using the regional atmospheric meteorological system (RAMS). This paper presents a comparison of modelled and observed wind and temperature for 25 and 26 January 2000. The RAMS model uses two nested grids with a mesh size of 2 km for the inner grid which is embedded in the outer grid with a mesh size of 8 km. Statistical evaluation of the model results against the observational data of wind speed, direction and temperature at 10 m was conducted. In general, the modelling results are in a reasonable agreement with observation. The statistical evaluation suggests that model performance is poorer for the inner grid than the outer grid as the model uncertainties (mainly mean bias) transfer from the outer to inner one. The low indices of agreement of temperature and wind are mainly associated with the systematic root-mean-square-difference values. For temperature, the systematic bias may also be affected by representation of cloud amount by the model. For wind, the model tends to have a poor performance for calm conditions, as under a stable anti-cyclonic situation local wind patterns associated with topography may develop, although the topography of the region is relatively flat. The results for the inner grid reveal some subtle spatial patterns at a scale smaller than 10 km near hills and valleys with differences in elevation of a few hundred metres.     

Effective stress analyses of liquefaction-induced deformation in river dikes

A series of effective stress analyses is carried out on the seismic performance of river dikes based on the case histories during the 1993 Hokkaido-Nansei-oki and 1995 Hyogoken-Nambu earthquakes in Japan. Seven case histories selected for the analyses involve a crest settlement ranging from none to 2.7 m in the dikes 3–6 m high with evidence of liquefaction at foundation soil. The effective stress model used is based on a multiple shear mechanism and was developed by one of the authors. The soil parameters are evaluated based on the site investigation and laboratory test results. The results of the analyses are basically consistent with the observed performance of the river dikes. In particular, the effective stress model shows a reasonable capability to reproduce the varying degree of settlements depending on the geotechnical conditions of foundation soils beneath the dikes. The analyses also indicate that the effect of a cohesive soil layer mixed with the liquefiable sand layers beneath the dikes can be a primary factor for reducing the liquefaction-induced deformation of dikes.     

Optimal placement of dampers and actuators based on stochastic approach

A general method is developed for optimal application of dampers and actuators by installing them at optimal location on seismic-resistant structures. The study includes development of a statistical criterion, formulation of a general optimization problem and establishment of a solution procedure. Numerical analysis of the seismic response in time-history of controlled structures is used to verify the proposed method for optimal device application and to demonstrate the effectiveness of seismic response control with optimal device location. This study shows that the proposed method for the optimal device application is simple and general, and that the optimally applied dampers and actuators are very efficient for seismic response reduction.     

Models of passive margin inversion: implications for the Rhenohercynian fold-and-thrust belt, Belgium and Germany

Positive tectonic inversion is related to the transmission of compressional stresses along a décollement into the foreland of an orogenic zone. This stress and strain concentration in regions remote from the main orogenic front is commonly related to the presence of pre-existing rheological heterogeneities such as normal syn-depositional faults. During inversion, these pre-existing normal faults are reactivated as reverse faults. Tectonic inversion in the Rhenohercynian fold-and-thrust belt during the Variscan Orogeny shows that inversion is likely synchronous with the onset of collision in the hinterland. Here, we present the results of a simplified thermo-mechanical model (STM) which allows one to study strain partitioning between two orogenic zones. We show that, if the two orogenic zones have the same mechanical properties, the viscosity of the décollement, which links them, controls the initial strain partitioning. During subsequent finite shortening, erosional processes determine the partitioning of strain rate. The presence of a weak structure in the inverted zone and of a low-viscosity décollement leads to initial strain concentration in the inverted track rather than in the collision zone and a progressive decrease in strain partitioning between the two orogenic zones. The STM results are in good agreement with results of a 2D finite-element model. We conclude that, in the western part of the Rhenohercynian Massif, simultaneous uplift and deformation within the Mid-German Crystalline Rise (the main collision zone) and the Ardenne Anticlinorium (the inverted zone) lead to interpreting this orogenic event as a case of vice tectonic rather than the propagation of a ‘wave of folding’ towards the Variscan front, as suggested by previous authors.     

The null distribution of sample serial correlation coefficient

 The null distribution of the lag-k sample serial correlation coefficient (r _k , k=1,2,3) was investigated by Monte Carlo simulation. For a time series with normal, exponential, Pearson 3, EV1 (Gumbel), or generalized Pareto (GP) distribution type, the null distribution of its r _k can be approximated by the normal distribution with mean −1/(n−k) and variance 1/(n−1). But for a time series with the lognormal, EV2 or EV3 (Weibull) distribution type, the null distribution of r _k is skewed distributed. In such cases, a simulation technique is suggested to construct percentile confidence intervals at a given significance level.     

Model-independent Travel-time Attributes for 2-D, Finite-offset Multicoverage Reflections

— In this paper, we provide a 5-parameter stacking formula to transform 2-D prestack data into a particular common-offset section. This requires the knowledge of the near-surface velocity only and it is expected that ray theory holds to describe primary reflections. The earth model can be arbitrarily inhomogeneous. The new stacking approach can be viewed as a generalization of the 3-parameter common-reflection-surface (CRS) stack, by which 2-D multicoverage data are stacked into a simulated zero-offset section. The new 5-parameter formula can handle P-P, P-S and S-S reflections.     

Microseismicity Induced by Heavy Rainfall Around Flooded Vertical Ore Veins

—?Microseismicity (M?M?M?M??6 were also monitored at a crustal movement monitoring station located several hundred meters from the veins. It was found that the opening of the vertical ore veins primary led to significant strain and tilt, but not to seismicity, because the delay and the longer duration of the seismicity were significant. Most seismic events involve thrusting mechanisms that are consistent with the present stress state of E-W-oriented tectonic compression, but are not consistent with the opening of the deepest ore vein. Interstingly, all the events within a few months of the heavy rainfall occurred near the faults that offset the deepest ore veins, wheareas all those events located away from the deepest ore veins occurred many months after the heavy rainfall. Consequently, the delayed diffusion of water appears to have played a dominant role in reducing rock strength, which led to seismicity in the Ikuno mine.     

Power of the Mann–Whitney test for detecting a shift in median or mean of hydro-meteorological data

 The non-parametric Mann–Whitney (MW) statistic test has been popularly used to assess the significance of a shift in median or mean of hydro-meteorological time series. It has been considered that the test is more suitable for non-normally distributed data and it may be not sensitive to the distribution type of sample data. However, no evidence has been provided to demonstrate these. This study investigates the power of the test in various circumstances by means of Monte Carlo simulation. Simulation results demonstrate that the power of the test is very sensitive to various properties of sample data. The power depends on the pre-assigned significance level, magnitude of a shift, sample size, and its occurrence position within a time series; and it is also strongly affected by the variation, skewness, and distribution type of a time series. The bigger the magnitude of a shift, the more powerful the test is; the larger the sample size, the more powerful the test is; and the bigger the variation within a time series, the less power the test has. The test has the highest power if a shift occurs at the midpoint of a time series. For the samples with different distribution types, the power of the test is dramatically different. The test has the highest power for time series with the extreme value type III (EV3) distribution while it indicates the lowest power for time series with the lognormal distribution.