Tsunami vulnerability assessment of Casablanca-Morocco using numerical modelling and GIS tools

Earthquakes and tsunamis along Morocco's coasts have been reported since historical times. The threat posed by tsunamis must be included in coastal risk studies. This study focuses on the tsunami impact and vulnerability assessment of the Casablanca harbour and surrounding area using a combination of tsunami inundation numerical modelling, field survey data and geographic information system. The tsunami scenario used here is compatible with the 1755 Lisbon event that we considered to be the worst case tsunami scenario. Hydrodynamic modelling was performed with an adapted version of the Cornell Multigrid Coupled Tsunami Model from Cornell University. The simulation covers the eastern domain of the Azores-Gibraltar fracture zone corresponding to the largest tsunamigenic area in the North Atlantic. The proposed vulnerability model attempts to provide an insight into the tsunami vulnerability of building stock. Results in the form of a vulnerability map will be useful for decision makers and local authorities in preventing the community resiliency for tsunami hazards.

Toplological Entropy in the Syncronization of Piecewise Linear and Monotone Maps. Coupled Duffing Oscillators

In this paper is presented a relationship between the synchronization and the topological entropy. We obtain the values for the coupling parameter, in terms of the topological entropy, to achieve synchronization of two unidirectional and bidirectional coupled piecewise linear maps. In addition, we prove a result that relates the synchronizability of two m-modal maps with the synchronizability of two conjugated piecewise linear maps. An application to the unidirectional and bidirectional coupled identical chaotic Duffing equations is given. We discuss the complete synchronization of two identical double-well Duffing oscillators, from the point of view of symbolic dynamics. Working with Poincare cross-sections and the return maps associated, the synchronization of the two oscillators, in terms of the coupling strength, is characterized.

Paleomagnetic study of the Great Foum Zguid dyke (southern Morocco): A positive contact test related to metasomatic processes

When a paleomagnetic pole is sought for in an igneous body, the host rocks should be subjected to a contact test to assure that the determined paleopole has the age of the intrusion. If the contact test is positive, it precludes the possibility that the measured magnetization is a later effect. Therefore, we investigated the variations of the remanent magnetization along cross-sections of rocks hosting the Foum Zguid dyke (southern Morocco) and the dyke itself. A positive contact test was obtained, but it is mainly related with Chemical/Crystalline Remanent Magnetization due to metasomatic processes in the host-rocks during magma intrusion and cooling, and not only with Thermo-Remanent Magnetization as ordinarily assumed in standard studies. Paleomagnetic data obtained within the dyke then reflect the Earth magnetic field during emplacement of this well-dated (196.9 +/- 1.8 Ma) intrusion.

Preliminary results of a study of magnetic properties in the Foum-Zguid dyke (Morocco)

This work focuses on the study of flow and propagation of magma using rock magnetic analyses along sections across the thick Jurassic dyke of Foum-Zguid (Southern Morocco). Thermomagnetic data show that Ti-poor titanomagnetite is the main magnetic carrier. Petrographic analysis shows that the main Ti phase (ilmenite) occurs either as lamellae within spinel (center of the dyke) or as isolated grains (dyke margin). Bulk magnetic properties display distinct behavior according to the distance to the dyke margin; grain size of the main magnetic carrier decreases towards the center of the dyke, while the natural remanent magnetization and the bulk magnetic susceptibility increase. Only the magnetic susceptibility ellipsoid close to the dyke margin corresponds to that usually found in thin dykes, with the magnetic foliation sub parallel to dyke margins. Maximum principal axis is in most cases either parallel or perpendicular to the intersection between the planes of magnetic foliation and dyke wall. Moreover, when this axis is perpendicular to the intersection it is associated with a more oblate magnetic susceptibility ellipsoid shape, indicating the presence of complex magnetic fabrics. The studied magnetic properties show that, in this 100 m wide thick dyke, flow structures related with dyke propagation are only preserved close to the quickly cooled dyke margins.

Thick dyke emplacement and internal flow: A structural and magnetic fabric study of the deep-seated dolerite dyke of Foum Zguid (southern Morocco)

Knowledge on forced magma injection and magma flow in dykes is crucial for the understanding of how magmas migrate through the crust to the Earth's surface. Because many questions still persist, we used the long, thick, and deep-seated Foum Zguid dyke (Morocco) to investigate dyke emplacement and internal flow by means of magnetic methods, structural analysis, petrography, and scanning electron microscopy. We also investigated how the host rocks accommodated the intrusion. Regarding internal flow: 1. Important variations of the rock magnetic properties and magnetic fabric occur with distance from dyke wall; 2. anisotropy of anhysteretic remanent magnetization reveals that anisotropy of magnetic susceptibility (AMS) results mainly from the superposition of subfabrics with distinct coercivities and that the imbrication between magnetic foliation and dyke plane is more reliable to deduce flow than the orientation of the AMS maximum principal axis; and 3. a dominant upward flow near the margins can be inferred. The magnetic fabric closest to the dyke wall likely records magma flow best due to fast cooling, whereas in the core the magnetic properties have been affected by high-temperature exsolution and metasomatic effects due to slow cooling. Regarding dyke emplacement, this study shows that the thick forceful intrusion induced deformation by homogeneous flattening and/or folding of the host sedimentary strata. Dewatering related to heat, as recorded by thick quartz veins bordering the dyke in some localities, may have also helped accommodating dyke intrusion. The spatial arrangement of quartz veins and their geometrical relationship with the dyke indicate a preintrusive to synintrusive sinistral component of strike slip.

Variation of magnetic properties in sedimentary rocks hosting the Foum Zguid dyke (southern Morocco): Combined effects of re-crystallization and Fe-metasomatism

The effects of dyke intrusion on the magnetic properties of host sedimentary rocks are still poorly understood. Therefore, we have evaluated bulk magnetic parameters of standard palaeomagnetic samples collected along several sections across the sediments hosting the Foum Zguid dyke in southern Morocco. The study has been completed with the evaluation of the magnetic fabric after laboratory application of sequential heating experiments. The present study shows that: (1) close to Fourn Zguid dykes, the variations of the bulk magnetic parameters and of the magnetic fabric is strongly related with re-crystallization and Fe-metasomatism intensity. (2) The thermal experiments on AMS of samples collected farther from the dyke and, thus, less affected by heating during dyke emplacement, indicate that 300-400 degrees C is the minimum experimental temperature necessary to trigger appreciable transformations of the pre-existing magnetic fabrics. For temperatures higher than ca. 580 degrees C, the magnetic fabric transformations are fully realized, with complete transposition of the initial fabric to a fabric similar to that of samples collected close to the dyke. Therefore, measured variations of the magnetic fabric can be used to evaluate re-crystallization temperatures experienced by the host sedimentary rock during dyke emplacement. The distinct magnetic behaviour observed along the cross-sections strongly suggests that samples collected farther from the dyke margins did not experience thermal episodes with temperatures higher than 300 degrees C after dyke emplacement. (3) AMS data shows a gradual variation of the magnetic fabric with distance from the dyke margin, from sub-horizontal K-3 away from the dyke to vertical K3 close to the dyke. Experimental heating shows that heat alone can be responsible for this strong variation. Therefore, such orientation changes should not be unequivocally interpreted as the result of a stress field (resulting from the emplacement of the dyke, for instance). (4) Magnetic studies prove to be a very sensitive tool to assess rock magnetic transformations, thermally and chemically induced by dyke intrusion in hosting sediments.

Genus and Braid Index Associated to Sequences of Renormalizable Lorenz Maps

We describe the Lorenz links generated by renormalizable Lorenz maps with reducible kneading invariant (K(f)(-), = K(f)(+)) = (X, Y) * (S, W) in terms of the links corresponding to each factor. This gives one new kind of operation that permits us to generate new knots and links from the ones corresponding to the factors of the *-product. Using this result we obtain explicit formulas for the genus and the braid index of this renormalizable Lorenz knots and links. Then we obtain explicit formulas for sequences of these invariants, associated to sequences of renormalizable Lorenz maps with kneading invariant (X, Y) * (S,W)*(n), concluding that both grow exponentially. This is specially relevant, since it is known that topological entropy is constant on the archipelagoes of renormalization.

Performance of coastal sea-defense infrastructure at El Jadida (Morocco) against tsunami threat: lessons learned from the Japanese 11 March 2011 tsunami

This paper seeks to investigate the effectiveness of sea-defense structures in preventing/reducing the tsunami overtopping as well as evaluating the resulting tsunami impact at El Jadida, Morocco. Different tsunami wave conditions are generated by considering various earthquake scenarios of magnitudes ranging from M-w = 8.0 to M-w = 8.6. These scenarios represent the main active earthquake faults in the SW Iberia margin and are consistent with two past events that generated tsunamis along the Atlantic coast of Morocco. The behavior of incident tsunami waves when interacting with coastal infrastructures is analyzed on the basis of numerical simulations of near-shore tsunami waves' propagation. Tsunami impact at the affected site is assessed through computing inundation and current velocity using a high-resolution digital terrain model that incorporates bathymetric, topographic and coastal structures data. Results, in terms of near-shore tsunami propagation snapshots, waves' interaction with coastal barriers, and spatial distributions of flow depths and speeds, are presented and discussed in light of what was observed during the 2011 Tohoku-oki tsunami. Predicted results show different levels of impact that different tsunami wave conditions could generate in the region. Existing coastal barriers around the El Jadida harbour succeeded in reflecting relatively small waves generated by some scenarios, but failed in preventing the overtopping caused by waves from others. Considering the scenario highly impacting the El Jadida coast, significant inundations are computed at the sandy beach and unprotected areas. The modeled dramatic tsunami impact in the region shows the need for additional tsunami standards not only for sea-defense structures but also for the coastal dwellings and houses to provide potential in-place evacuation.

Arrow plot and CA maps on microarray preprocessing methods

Microarray allow to monitoring simultaneously thousands of genes, where the abundance of the transcripts under a same experimental condition at the same time can be quantified. Among various available array technologies, double channel cDNA microarray experiments have arisen in numerous technical protocols associated to genomic studies, which is the focus of this work. Microarray experiments involve many steps and each one can affect the quality of raw data. Background correction and normalization are preprocessing techniques to clean and correct the raw data when undesirable fluctuations arise from technical factors. Several recent studies showed that there is no preprocessing strategy that outperforms others in all circumstances and thus it seems difficult to provide general recommendations. In this work, it is proposed to use exploratory techniques to visualize the effects of preprocessing methods on statistical analysis of cancer two-channel microarray data sets, where the cancer types (classes) are known. For selecting differential expressed genes the arrow plot was used and the graph of profiles resultant from the correspondence analysis for visualizing the results. It was used 6 background methods and 6 normalization methods, performing 36 pre-processing methods and it was analyzed in a published cDNA microarray database (Liver) available at http://genome-www5.stanford.edu/ which microarrays were already classified by cancer type. All statistical analyses were performed using the R statistical software.

Synchronization and Basins of Synchronized in Two-Dimensional Piecewise Maps Via Coupling Three Pieces of One-Dimensional Maps

This paper is devoted to the synchronization of a dynamical system defined by two different coupling versions of two identical piecewise linear bimodal maps. We consider both local and global studies, using different tools as natural transversal Lyapunov exponent, Lyapunov functions, eigenvalues and eigenvectors and numerical simulations. We obtain theoretical results for the existence of synchronization on coupling parameter range. We characterize the synchronization manifold as an attractor and measure the synchronization speed. In one coupling version, we give a necessary and sufficient condition for the synchronization. We study the basins of synchronization and show that, depending upon the type of coupling, they can have very different shapes and are not necessarily constituted by the whole phase space; in some cases, they can be riddled.

Scaling law in Saddle-node bifurcations for one-dimensional maps: A complex variable approach

The study of transient dynamical phenomena near bifurcation thresholds has attracted the interest of many researchers due to the relevance of bifurcations in different physical or biological systems. In the context of saddle-node bifurcations, where two or more fixed points collide annihilating each other, it is known that the dynamics can suffer the so-called delayed transition. This phenomenon emerges when the system spends a lot of time before reaching the remaining stable equilibrium, found after the bifurcation, because of the presence of a saddle-remnant in phase space. Some works have analytically tackled this phenomenon, especially in time-continuous dynamical systems, showing that the time delay, tau, scales according to an inverse square-root power law, tau similar to (mu-mu (c) )(-1/2), as the bifurcation parameter mu, is driven further away from its critical value, mu (c) . In this work, we first characterize analytically this scaling law using complex variable techniques for a family of one-dimensional maps, called the normal form for the saddle-node bifurcation. We then apply our general analytic results to a single-species ecological model with harvesting given by a unimodal map, characterizing the delayed transition and the scaling law arising due to the constant of harvesting. For both analyzed systems, we show that the numerical results are in perfect agreement with the analytical solutions we are providing. The procedure presented in this work can be used to characterize the scaling laws of one-dimensional discrete dynamical systems with saddle-node bifurcations.

Enhanced PCA-based localization using depth maps with missing data

In this paper a new method for self-localization of mobile robots, based on a PCA positioning sensor to operate in unstructured environments, is proposed and experimentally validated. The proposed PCA extension is able to perform the eigenvectors computation from a set of signals corrupted by missing data. The sensor package considered in this work contains a 2D depth sensor pointed upwards to the ceiling, providing depth images with missing data. The positioning sensor obtained is then integrated in a Linear Parameter Varying mobile robot model to obtain a self-localization system, based on linear Kalman filters, with globally stable position error estimates. A study consisting in adding synthetic random corrupted data to the captured depth images revealed that this extended PCA technique is able to reconstruct the signals, with improved accuracy. The self-localization system obtained is assessed in unstructured environments and the methodologies are validated even in the case of varying illumination conditions.