3-D crustal structure in the Agadir region (SW High Atlas,Morocco)

The 1960 Agadir earthquake (Mw 6.0) constitutes the most damaging earthquake event in Morocco. With the expansion of seismic networks during the last decade in Morocco, new seismic data have been collected in this region. The P and S arrivals at 19 stations located in Southern Morocco are used to investigate the lithosphere in the Agadir region. In this study, we use a linearized inversion procedure comprising two steps: (1) finding the minimal 1-D model and simultaneous relocation of hypocentres and (2) determination of local velocity structure using linearized inversion. The model parameterization in this method assumes a continuous velocity field. The resolution tests indicate that the calculated images give near true structure for the studied region from 0- to 45-km depth. The results show that the total crust thickness varies from 30 to 40 km in SW High Atlas and confirm the modest crustal tectonic shortening and thickening in the Atlas Mountains of Morocco. The inferred geological structure reconstructed from the calculated image illustrates the existence of fault-related folding. The evidence for coseismic ruptures in 1960 on the Kasbah anticline combined with the 1960 earthquake hypocentre located in the tomographic image determines the seismic potential of the active fault and related fold. The resulting tomographic image shows a high-velocity anomalies that could be associated with the location of deep active fault (10–30 km) associated with the fold structure. In the South Atlas, theses anomalies could be associated with the South atlas thrust front structure.     

TM/ETM+/LDCM Images for Studying Land Surface Temperature (LST) Interplay with Impervious Surfaces Changes over Time Within the Douala Metropolis,Cameroon

Douala, the most important metropolis of Cameroon, is a sub-Saharan wet coastal environment of which the anarchic urbanization is a socio-economic and environmental problem, significantly influencing the local climate. In this study, three Landsat images from 1986 (TM), 2007 (ETM+) and 2016 (LDCM), were utilized to investigate the effect of this urbanization on the increasing land surface temperature (LST) between these dates. Thus, the urban indices (UI), determined from the Landsat Visible and NIR channels were used to identify impervious areas (Urban Fabric and bare soil) of urban area. It has been shown from the UI images that, impervious areas have been increased from 1986 to 2016. The LST images derived have a continual expansion of zones and points of heat throughout these dates. The correlation analysis of LST and UI, at the pixel-scale, indicated the positive relationship between these parameters, which could show a real impact of urbanization on the increasing temperature in the area. These correlations are fairly low in 1986 (maximum R-square value is about 0.35) and in 2007 (maximum R-square value is about 0.44. In 2016, a high positive correlation (maximum R-square value is about 0.77) confirm that, the impervious areas strengthen the temperature and the Urban Heat Island effect in Douala urban zone. Overall, the earth observation images and the geographic information system techniques were effective approaches for aiming at environment monitoring and analyzing urban growth patterns and evaluating their impacts on urban climates.     

Long-term seismicity in regions of present day low seismic activity: the example of western Europe

In western Europe, the knowledge of long-term seismicity is based on reliable historical seismicity and covers a time period of less than 700 years. Despite the fact that the seismic activity is considered as low in the region extending from the Lower Rhine Embayment to England, historical information collected recently suggests the occurrence of three earthquakes with magnitude around 6.0 or greater. These events are a source of information for the engineer or the scientist involved in mitigation against large earthquakes. We provide information relevant to this aspect for the Belgian earthquake of September 18, 1692. The severity of the damage described in original sources indicates that its epicentral intensity could be IX (EMS-98 scale) and that the area with intensity VII and greater than VII has at least a mean radius of 45 km. Following relationships between average macroseismic radii and magnitude for earthquakes in stable continental regions, its magnitude M_s is estimated as between 6.0 and 6.5. To extend in time our knowledge of the seismic activity, we conducted paleoseismic investigations in the Roer Graben to address the question of the possible occurrence of large earthquakes with coseismic surface ruptures. Our study along the Feldbiss fault (the western border of the graben) demonstrates its recent activity and provides numerous lines of evidence of Holocene and Late Pleistocene large earthquakes. It suggests that along the 10 km long Bree fault scarp, the return period for earthquakes with magnitude from 6.2 to 6.7 ranges from 10,000 to 20,000 years during the last 50,000 years. Considering as possible the occurrence of similar earthquakes along all the Quaternary faults in the Lower Rhine Embayment, a large earthquake could occur there each 500–1000 years. These results are important in two ways. (i) The evidence that large earthquakes occur in western Europe in the very recent past which is not only attested by historical sources, but also suggested by paleoseismic investigations in the Roer Graben. (ii) The existence of a scientific basis to better evaluate the long-term seismicity in this part of Europe (maximal magnitude and return period) in the framework of seismic hazard assessment.     

Wave climate simulation for southern region of the South China Sea

This study investigates long-term variability and wave characteristic trends in the southern region of the South China Sea (SCS). We implemented the state-of-the art WAVEWATCH III spectral wave model to simulate a 31-year wave hindcast. The simulation results were used to assess the inter-annual variability and long-term changes in the SCS wave climate for the period 1979 to 2009. The model was forced with Climate Forecast System Reanalysis winds and validated against altimeter data and limited available measurements from an Acoustic Wave and Current recorder located offshore of Terengganu, Malaysia. The mean annual significant wave height and peak wave period indicate the occurrence of higher wave heights and wave periods in the central SCS and lower in the Sunda shelf region. Consistent with wind patterns, the wave direction also shows southeasterly (northwesterly) waves during the summer (winter) monsoon. This detailed hindcast demonstrates strong inter-annual variability of wave heights, especially during the winter months in the SCS. Significant wave height correlated negatively with Niño3.4 index during winter, spring and autumn seasons but became positive in the summer monsoon. Such correlations correspond well with surface wind anomalies over the SCS during El Nino events. During El Niño Modoki, the summer time positive correlation extends northeastwards to cover the entire domain. Although significant positive trends were found at 95 % confidence levels during May, July and September, there is significant negative trend in December covering the Sunda shelf region. However, the trend appears to be largely influenced by large El Niño signals.     

Nonextensive dust acoustic waves in a charge varying dusty plasma

Our recent analysis on nonlinear nonextensive dust-acoustic waves (DA) [Amour and Tribeche in Phys. Plasmas 17:063702, 2010] is extended to include self-consistent nonadiabatic grain charge fluctuation. The appropriate nonextensive electron charging current is rederived based on the orbit-limited motion theory. Our results reveal that the amplitude, strength and nature of the nonlinear DA waves (solitons and shocks) are extremely sensitive to the degree of ion nonextensivity. Stronger is the electron correlation, more important is the charge variation induced nonlinear wave damping. The anomalous dissipation effects may prevail over that dispersion as the electrons evolve far away from their Maxwellian equilibrium. Our investigation may be of wide relevance to astronomers and space scientists working on interstellar dusty plasmas where nonthermal distributions are turning out to be a very common and characteristic feature.     

RETRACTED ARTICLE: Numerical analysis of multiphase flow in porous material

Recent developments in the application of x-ray micro-tomography in laboratory geomechanics have allowed all the individual grains of sand in a test sample to be seen and identified uniquely in 3D. Combining such imaging capabilities with experiments carried out “in situ” within an imaging set-up has led to the possibility of directly observing the mechanisms of deformation as they happen. The challenge has thus become extracting pertinent, quantified information from these rich time-lapse 3D images to elucidate the mechanics at play. This paper presents a new approach (ID-Track) for the quantification of individual grain kinematics (displacements and rotations) of large quantities of sand grains (tens of thousands) in a test sample undergoing loading. With ID-Track, grains are tracked between images based on some geometrical feature(s) that allow their unique identification and matching between images. This differs from Digital Image Correlation (DIC), which makes measurements by recognising patterns between images. Since ID-Track does not use the image of a grain for tracking, it is significantly faster than DIC. The technique is detailed in the paper, and is shown to be fast and simple, giving good measurements of displacements, but suffering in the measurement of rotations when compared with Discrete DIC. Subsequently, results are presented from successful applications of ID-track to triaxial tests on two quite different sands: the angular Hostun sand and the rounded Caicos Ooids. This reveals details on the performance of the technique for different grain shapes and insight into the differences in the grain-scale mechanisms occurring in these two sands as they exhibit strain localisation under triaxial loading.     

Dimensional Reduction of Pattern-Based Simulation Using Wavelet Analysis

A pattern-based simulation technique using wavelet analysis is proposed for the simulation (wavesim) of categorical and continuous variables. Patterns are extracted by scanning a training image with a template and then storing them in a pattern database. The dimension reduction of patterns in the pattern database is performed by wavelet decomposition at certain scale and the approximate sub-band is used for pattern database classification. The pattern database classification is performed by the k-means clustering algorithm and classes are represented by a class prototype. For the simulation of categorical variables, the conditional cumulative density function (ccdf) for each class is generated based on the frequency of the individual categories at the central node of the template. During the simulation process, the similarity of the conditioning data event with the class prototypes is measured using the L ₂-norm. When simulating categorical variables, the ccdf of the best matched class is used to draw a pattern from a class. When continuous variables are simulated, a random pattern is drawn from the best matched class. Several examples of conditional and unconditional simulation with two- and three- dimensional data sets show that the spatial continuity of geometric features and shapes is well reproduced. A comparative study with the filtersim algorithm shows that the wavesim performs better than filtersim in all examples. A full-field case study at the Olympic Dam base metals deposit, South Australia, simulates the lithological rock-type units as categorical variables. Results show that the proportions of various rock-type units in the hard data are well reproduced when similar to those in the training image; when rock-type proportions between the training image and hard data differ, the results show a compromise between the two.     

The LATEA metacraton (Central Hoggar, Tuareg shield, Algeria): behaviour of an old passive margin during the Pan-African orogeny

Historically, the Tuareg shield is divided into three parts bordered by mega-shear zones with the centre, the Central Polycyclic Hoggar, characterized by Archaean and Palaeoproterozoic lithologies. Nearly 10 years ago, the Tuareg shield was shown to be composed of 23 displaced terranes [Geology 22 (1994) 641] whose relationships were deciphered in Aïr to the SE [Precambr. Res. 67 (1994) 59]. The Polycyclic Central Hoggar terranes were characterized by the presence of well preserved Archaean/Palaeoproterozoic and Neoproterozoic lithologies.We show here that the terranes from Central Hoggar (Laouni, Azrou-n-Fad, Tefedest, Egéré-Aleksod) belonged to a single old passive margin, to which we gave the acronym name LATEA, which behaved as a craton during the Mesoproterozoic and the Early-Middle Neoproterozoic but was partly destabilized and dissected during the Late Neoproterozoic as a consequence of its involvement as a passive margin in the Pan-African orogen.An early Pan-African phase consisted of thrust sheets including garnet-bearing lithologies (eclogite, amphibolite, gneiss) that can be mapped and correlated in three LATEA terranes. In the Tin Begane area, P–T–t paths have been established from>15 kbar––790 °C (eclogite) to 4 kbar––500 °C (greenschist retrogression) through 12 kbar––830 °C (garnet amphibolite) and 8 kbar––700 °C (garnet gneiss), corresponding to the retrograde path of a Franciscan-type loop. Sm–Nd geochronology on minerals and laser ablation ICP-MS on garnet show the mobility of REE, particularly LREE, during the retrograde greenschist facies that affects, although slightly, some of these rocks. The amphibolite-facies metamorphism has been dated at 685 ± 19 Ma and the greenschist facies at 522 ± 27 Ma. During the thrust phase, the Archaean–Palaeoproterozoic basement was only locally affected by the Pan-African tectonics. LATEA behaved as a craton. Other juvenile terranes were also thrust early onto LATEA: the Iskel island arc at ≈850 Ma to the west of LATEA, the Serouenout terrane in the 700–620 Ma age range to the east. No subduction-related magmas have intruded LATEA during this epoch, which behaved as a passive margin.During the main Pan-African phase (625–580 Ma), LATEA was dissected by mega-shear zones that induced several hundreds km of relative displacement and allowed the emplacement of high-K calc-alkaline batholiths. Smaller movements continued till 525 Ma, accompanied by the emplacement of subcircular plutons with alkaline affinity. Here is dated the Ounane granodiorite (624 ± 15 Ma; ⁸⁷Sr/⁸⁶Sr_i=0.70839 ± 0.00016; 6WR, MSWD=0.87) and the Tisselliline granite (552 ± 15 Ma; ⁸⁷Sr/⁸⁶Sr_i=0.7074 ± 0.0001; 5WR, MSWD=1.4). Nd isotopes indicate a preponderant Palaeoproterozoic crustal source for these two plutons: _Nd=−14 to −21 at 624 Ma and T_DM=1650–2320 Ma for Ounane and _Nd=−13 to −15 at 555 Ma and T_DM=1550–1720 Ma for Tisselliline. Our model links these intrusions to a linear lithospheric delamination along mega-shear zones, allowing the hot asthenosphere to rise, melt by adiabatic pressure release and inducing the melting of the Palaeoproterozoic and Archaean lower crust.The LATEA cratonic microcontinent remained however sufficiently rigid to preserve Archaean and Palaeoproterozoic lithologies as well as Middle Neoproterozoic oceanic thrust sheets. This corresponds to the notion of metacraton [J. African Earth Sci. 34 (2002) 119], i.e. a craton that has been remobilized during an orogenic event but is still recognizable dominantly through its rheological, geochronological, isotopic and sometimes petrological characteristics.     

Rapid analysis of risk assessment using developed simulation of chemical industrial accidents software package

The environmental consequences are defined as consequences of accidental release of hazardous substances to the natural environment. This release can lead to many hazards depending on the material stored. The consequences of these hazards to the environment are widespread and have significant importance to human communities living in the surroundings. The mathematical models are extremely useful tools to predict the impacts of chemical process accidents. The objective of this paper is to develop a software package for accident simulation and damage potential estimation. The software is coded in visual basic and is compatible with windows working environments. The software is called Simulation of chemical industrial accident. This application is a comprehensive software package which can be integrated with geographical information system to predict and display the consequence of chemical hazards. The software is a user-friendly and effective tool for evaluating the consequences of major chemical accidents, process decision making for land-use planning, namely locating suitable hazardous installations, hazardous waste disposal areas and emergency response plan.     

Numerical modeling of the flow and transport of radionuclides in heterogeneous porous media

This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and stable, as confirmed by numerical simulations.