@article {Mel{\'e}ndez2015158,
	title = {{TOMO3D: 3-D joint refraction and reflection traveltime tomography parallel code for active-source seismic data-synthetic test}},
	journal = {Geophysical Journal International},
	volume = {203},
	number = {1},
	year = {2015},
	pages = {158{\textendash}174},
	abstract = {We present a new 3-D traveltime tomography code (TOMO3D) for the modelling of active-source seismic data that uses the arrival times of both refracted and reflected seismic phases to derive the velocity distribution and the geometry of reflecting boundaries in the subsurface. This code is based on its popular 2-D version TOMO2D from which it inherited the methods to solve the forward and inverse problems. The traveltime calculations are done using a hybrid ray-tracing technique combining the graph and bending methods. The LSQR algorithm is used to perform the iterative regularized inversion to improve the initial velocity and depth models. In order to cope with an increased computational demand due to the incorporation of the third dimension, the forward problem solver, which takes most of the run time (\~{}90 per cent in the test presented here), has been parallelized with a combination of multi-processing and message passing interface standards. This parallelization distributes the ray-tracing and traveltime calculations among available computational resources. The code{\textquoteright}s performance is illustrated with a realistic synthetic example, including a checkerboard anomaly and two reflectors, which simulates the geometry of a subduction zone. The code is designed to invert for a single reflector at a time. A data-driven layer-stripping strategy is proposed for cases involving multiple reflectors, and it is tested for the successive inversion of the two reflectors. Layers are bound by consecutive reflectors, and an initial velocity model for each inversion step incorporates the results from previous steps. This strategy poses simpler inversion problems at each step, allowing the recovery of strong velocity discontinuities that would otherwise be smoothened. {\textcopyright} The Authors 2015. Published by Oxford University Press on behalf of The Royal Astronomical Society.},
	doi = {10.1093/gji/ggv292},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-84942105766\&partnerID=40\&md5=95ac8480fd67102b8f5322381f9ff8ef},
	author = {Mel{\'e}ndez, A.a and Korenaga, J.b and Sallar{\`e}s, V.a and Miniussi, A.c and Ranero, C.R.d}
}
@inbook {GomezdelaPenaL.RaneroC.R.GraciaE.BartolomeR.2014,
	title = {{T{\'e}cnicas avanzadas de procesado de datos de s{\'\i}smica de reflexi{\'o}n multicanal aplicadas a mejorar la visualizaci{\'o}n de fallas activas del Mar de Albor{\'a}n}},
	booktitle = {Una aproximaci{\'o}n multidisciplinar al estudio de fallas activas, los terremotos y el riesgo s{\'\i}smico},
	year = {2014},
	pages = {229{\textendash}232},
	isbn = {978-84-617-2049-1},
	url = {http://www.researchgate.net/publication/267574309\_Una\_aproximacin\_multidisciplinar\_al\_estudio\_de\_las\_fallas\_activas\_los\_terremotos\_y\_el\_riesgo\_ssmico},
	author = {G{\'o}mez de la Pe{\~n}a, L. and Ranero, C.R. and Gr{\`a}cia, E. and R., Bartolome and cruise Party, TopoMed}
}
@conference {Bartolome2014,
	title = {{TSUJAL project and marine survey: Crustal characterization of the Rivera Plate-Jalisco Block boundary and its implications for seismic and tsunami hazard assessment5) TSUJAL Working Group: D}},
	booktitle = {2nd Iberian meeting on active faults and paleoseismology {\textendash} IBERFAULT 2014},
	year = {2014},
	month = {oct},
	doi = {10.13140/2.1.3949.1524},
	url = {http://www.researchgate.net/publication/267514007\_TSUJAL\_project\_and\_marine\_survey\_Crustal\_characterization\_of\_the\_Rivera\_Plate-Jalisco\_Block\_boundary\_and\_its\_implications\_for\_seismic\_and\_tsunami\_hazard\_assessment5)\_TSUJAL\_Working\},
	author = {Bartolom{\'e}, R and Cordoba and Nu{\~n}ez-Corn{\'u} and Da{\~n}obeitia and Alonso, J.L. and Cameselle, A.L. and Castellon, A. and Estrada, F. and Perea, H. and Prada, M.}
}
@conference {Institut2013,
	title = {{Tu-02-13 TOMO3D - A New 3-D Joint Refraction and Reflection Travel-time Tomography Code for Active-source Seismic Data}},
	booktitle = {Eage {\textquoteright}13},
	year = {2013},
	month = {jun},
	pages = {10{\textendash}13},
	abstract = {TOMO3D is a code for three-dimensional refraction and reflection travel-time tomography of wide-angle seismic data that inverts for the velocity field and the geometry of a reflector. Ray tracing is performed by a hybrid method combining the graph and bending methods, and the inversion is iteratively solved using an LSQR algorithm. We present a series of benchmark tests with synthetic data for the forward and inverse problems, as well as other more complex tests comparing inversions with refraction travel times only (i.e. first arrivals) to others with refraction and reflection travel times. The combination of refraction and reflection data increases ray coverage and thus velocity resolution, while allowing the determination of major geological interfaces, and reducing the velocity-depth trade-off.},
	issn = {2214-4609},
	doi = {10.3997/2214-4609.20130372},
	url = {http://www.earthdoc.org/publication/publicationdetails/?publication=69247},
	author = {Institut, a Mel{\'e}ndez and Ci{\`e}ncies, De and Ci{\`e}ncies, V Sallar{\`e}s Institut De and Icrea, C R Ranero}
}
@conference {Ranero2012,
	title = {{The tectonic structure of the Tyrrhenian Basin, a complex interaction among faulting and magmatism}},
	booktitle = {Rendiconti Online Societa Geologica Italiana},
	volume = {21},
	number = {PART 1},
	year = {2012},
	pages = {251{\textendash}252},
	keywords = {Crustal structure, Tectonic structure, Tyrrhenian rifting},
	issn = {20358008},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-84872362909\&partnerID=tZOtx3y1},
	author = {Ranero, C. R. and Sallar{\`e}s, V. and Zitellini, N. and Grevemeyer, I. and Guzman, M. and Prada, M. and Moeller, S. and de Franco, R.}
}
@article {Rosas2012,
	title = {{Thrust{\textendash}wrench interference between major active faults in the Gulf of Cadiz (Africa{\textendash}Eurasia plate boundary, offshore SW Iberia): Tectonic implications from coupled analog and numerical modeling}},
	journal = {Tectonophysics},
	volume = {548-549},
	year = {2012},
	month = {jun},
	pages = {1{\textendash}21},
	abstract = {Analog and numerical modeling experiments were carried out to investigate the tectonic interference between intersecting major active strike-slip and thrust faults in the Gulf of Cadiz (Africa-Eurasia plate boundary, offshore SW Iberia). The obtained results show that newly mapped tectonic features located in the fault intersection area (corner zone) consist mostly in oblique (dextral-reverse) faults that accommodate significant strain partitioning. Modeling of this corner-zone faults show that they have endured some degree of rotation, displaying successive evolving geometries and kinematics. Numerical modeling results further show that an interbedded shallow soft layer, accounting for a regional (Late Miocene) gravitational "Chaotic" unit, could explain the mild bathymetric expression of the fault pattern in the corner-zone. Moreover, a recognized depth discrepancy, between the (upper crust) interference fault-pattern and the (lithospheric mantle) seismicity, is interpreted as a manifestation of similar thrust-wrench tectonic interference at different lithospheric depths. Accordingly, an intermediate lower crust-upper mantle aseismic (i.e. softened) depth-domain could be explained by pervasive alteration/serpentinization, prompted by fluid percolation through fault-related fractures associated with the newly revealed corner zone fault-network. Overall obtained results reinforce the relevance of a thrust-wrench multi-rupture seismic scenario as the main cause for the moderate seismicity (Mw < 6.0) in the study area. {\textcopyright} 2012 Elsevier B.V.},
	keywords = {Analog and numerical modeling, Fault interference, Gulf of Cadiz, Lithospheric mantle seismicity, Multi-rupture seismic hazards, Thrust-wrench tectonics},
	issn = {00401951},
	doi = {10.1016/j.tecto.2012.04.013},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-84862115104\&partnerID=tZOtx3y1},
	author = {Rosas, F.M. and Duarte, J.C. and Neves, M.C. and Terrinha, P. and Silva, S. and Matias, L. and Gr{\`a}cia, E. and Bartolome, R.}
}
@booklet {Goldfinger,
	title = {{Turbidite Event History{\textemdash}Methods and Implications for Holocene Paleoseismicity of the Cascadia Subduction Zone}},
	year = {2012},
	publisher = {U. S. Geological Survey},
	abstract = {Appendices. {\textendash} 170 pages, 64 figures, 13 tables.},
	url = {http://digital.csic.es/handle/10261/66386},
	author = {Goldfinger, Chris and Nelson, C. Hans and Morey, Ann E. and Johnson, Joel E. and Patton, Jason and Karabanov, Eugene and Guti{\'e}rrez-Pastor, Julia and Eriksson, Andrew T. and Gr{\'a}cia, Eul{\`a}lia and Dunhill, Gita and Enkin, Randolph J. and Dallimore, Audrey and Vallier, Tracy}
}
@article {Moore2011,
	title = {{Threshold of borehole failure: Breaking in before breaking out, Mississippi fan, Gulf of Mexico}},
	journal = {Geosphere},
	volume = {7},
	number = {3},
	year = {2011},
	pages = {684{\textendash}693},
	abstract = {Resistivity images from Integrated Ocean Drilling Program (IODP) Site U1322 on the Mississippi fan (Gulf of Mexico) show borehole failure as (1) low-resistivity bands interpreted as breakouts and (2) high-resistivity bands. Both features occur as opposing pairs on opposite sides of the borehole, and have similar azimuthal orientations and widths. Failures occur at depths of 90-216 m in sediments very rich in expansive (smectite-illite) clays of 40\%-50\% porosity that are younger than 65 ka. The low-resistivity breakouts resemble similar features in other IODP boreholes from southwest Japan and offshore Oregon. The high-resistivity features are unknown in other boreholes. Estimates of stress magnitudes based on the overburden stress and the extensional tectonic environment in the Gulf of Mexico predict that the borehole was at failure. Experiments were conducted on cores with lithologies equivalent to those of the borehole failure localities from IODP Site U1322 and adjacent Site U1324. These experiments suggest an elastic-plastic deformation with strains of 10\%-15\% before reaching a plastic yielding. In the experiments, strain softening during plastic deformation ranges from 0\% to 20\%. Physically the experimental samples show a combination of lateral bulging and discrete conjugate shears. These experiments suggest that the resistive areas in the borehole are an initial state of bulging, or extrusion, into the borehole. We call these extrusive failures "breakins" to distinguish them from traditional breakouts. Extrusion into borehole decreases the amount of conductive borehole fluid between the bulging sediment and the resistivity tool, increasing the resistivity signal. The high residual strength of the sediment prevents disaggregation and spalling. Where spalling has developed, breakouts occur. This analysis is the first documentation of this incipient stage of borehole failure. {\textcopyright} 2011 Geological Society of America.},
	issn = {1553-040},
	doi = {10.1130/GES00613.1},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-80054764068\&partnerID=tZOtx3y1},
	author = {Moore, J. Casey and Roger Urgeles, R. and Conin, Marianne and Flemings, Peter B. and Iturrino, Gerardo J.}
}
@article {Amblas2011,
	title = {{Transient erosion in the Valencia Trough turbidite systems, NW Mediterranean Basin}},
	journal = {Geomorphology},
	volume = {130},
	number = {3-4},
	year = {2011},
	month = {jul},
	pages = {173{\textendash}184},
	abstract = {Submarine canyons can efficiently drain sediments from continental margins just as river systems do in subaerial catchments. Like in river systems, submarine canyons are often arranged as complex drainage networks that evolve from patterns of erosion and deposition. In the present paper we use a morphometric analysis of submarine canyon-channel long-profiles to study the recent sedimentary history of the Valencia Trough turbidite system (VTTS) in the NW Mediterranean Sea. The VTTS is unique in that it drains sediment from margins with contrasting morphologies through a single "trunk" conduit, the Valencia Channel. The Valencia Channel has been active since the late Miocene, evolving in response to Plio-Quaternary episodes of erosion and deposition. The integrated analysis of long-profiles obtained from high-resolution bathymetric data across the entire turbidite system shows evidence for transient canyon incision in the form of knickpoints and hanging tributaries. Multiple factors appear to have triggered these periods of incision. These include a large debris flow at 11,500. yr. BP that disrupted the upper reaches of the VTTS and glacio-eustatic lowstands that forced shifting of sediment input to the VTTS. Based on these inferences, long-term time-averaged incision rates for the Valencia Channel have been estimated. The evidence we present strongly suggests that Foix Canyon has played a key role in the drainage dynamics of the VTTS in the past.This study builds conceptually on a recent modeling study that provides a morphodynamic explanation for the long-term evolution of submarine canyon thalweg profiles. The procedure and results from this work are of potential application to other submarine sediment drainage systems, past and present, including those containing mid-ocean type valleys like the Valencia Channel. {\textcopyright} 2011 Elsevier B.V.},
	keywords = {Knickpoints, Long-profiles, Submarine canyons, Transient erosion, Turbidity currents},
	issn = {0169555X},
	doi = {10.1016/j.geomorph.2011.03.013},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-79956155543\&partnerID=tZOtx3y1},
	author = {Amblas, David and Gerber, Thomas P. and Canals, Miquel and Pratson, Lincoln F. and Urgeles, Roger and Lastras, Galderic and Calafat, Antoni M.}
}
@article {Garcia2011,
	title = {{Tricyclic Antidepressants-loaded Biodegradable PLGA Nanoparticles: In Vitro Characterization and In Vivo Analgesic and Anti-Allodynic Effect}},
	journal = {Current Nanoscience},
	volume = {7},
	number = {3},
	year = {2011},
	month = {jun},
	pages = {345{\textendash}353},
	abstract = {Tricyclic antidepressants (TCAs) have potent local pain blockade properties that could be of interest in relieving chronic pain states as neuropathic pain. The aim of this work was to reach a persistent control of nociceptive and neuropathic pain by means of an in- jectable controlled release system using lower than usual doses of TCAs. To address this issue, amitriptyline, doxepin and imipramine were encapsulated with poly (lactic-co-glycolic) acid (PLGA) as polymer. Nanoparticles were characterized. The in vitro drug release profile and mechanism was evaluated, and the in vivo analgesic and anti-allodynic activity in front of heat-induced nociceptive pain and sciatic nerve chronic constriction injury, respectively, was tested. The mean{\textpm}SD particle size and drug loadings (\%) of the nanoparticles obtained were 420{\textpm}13, 480{\textpm}73 and 373{\textpm}25nm, and 40.46{\textpm}4.11, 31.09{\textpm}3.02 and 32.20{\textpm}3.20\% for amitriptyline, doxepin and imi- pramine, respectively. According to the Korsmeyer-Peppas model, the release mechanism of doxepin was diffusion controlled, while a combination of Fickian diffusion and polymer relaxation/erosion of the PLGA matrix was involved for amitriptyline and imipramine. Af- ter local infiltration of nanoparticles in rats, the antinociceptive and anti-allodynic activity of the encapsulated drugs were long-lasting and higher than that observed from the solutions. Amitriptyline elicited the lower analgesic effect. Doxepin showed the most outstanding results and its encapsulation led to a 62\% and 229\% increase in antinociceptive and anti-allodynic activity, respectively. So, this drug could be considered as a therapeutical alternative in pain relieving treatments.{\textcopyright} 2011 Bentham Science Publishers Ltd.},
	keywords = {Anti-allodynia, Antinociception, In vitro release, Korsmeyer-peppas, Plga nanoparticles, Tricyclic antidepressants},
	issn = {15734137},
	doi = {10.2174/157341311795542336},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-79956288549\&partnerID=tZOtx3y1},
	author = {Garcia, Xavier and Escribano, Elvira and Colom, Helena and Domenech, Josep and Queralt, Josep}
}
@article {Garcia2011a,
	title = {{A two-phase adaptive finite element method for solid-fluid coupling in complex geometries}},
	journal = {International Journal for Numerical Methods in Fluids},
	volume = {66},
	number = {1},
	year = {2011},
	month = {may},
	pages = {82{\textendash}96},
	abstract = {In this paper we present a method to solve the Navier-Stokes equations in complex geometries, such as porous sands, using a finite-element solver but without the complexity of meshing the porous space. The method is based on treating the solid boundaries as a second fluid and solving a set of equations similar to those used for multi-fluid flow. When combined with anisotropic mesh adaptivity, it is possible to resolve complex geometries starting with an arbitrary coarse mesh. The approach is validated by comparing simulation results with available data in three test cases. In the first we simulate the flow past a cylinder. The second test case compares the pressure drop in flow through random packs of spheres with the Ergun equation. In the last case simulation results are compared with experimental data on the flow past a simplified vehicle model (Ahmed body) at high Reynolds number using large-eddy simulation (LES). Results are in good agreement with all three reference models. {\textcopyright} 2010 John Wiley \& Sons, Ltd..},
	keywords = {Anisotropic mesh adaptivity, Ergun equation, Flow past a cylinder, Flow past sphere packs, Flow past the Ahmed body, Two-fluid approach},
	issn = {02712091},
	doi = {10.1002/fld.2249},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-79953156921\&partnerID=tZOtx3y1},
	author = {Garcia, Xavier and Pavlidis, Dimitrios and Gorman, Gerard J. and Gomes, Jefferson L. M. A. and Piggott, Matthew D. and Aristodemou, Elsa and Mindel, Julian and Latham, John-Paul and Pain, Christopher C. and ApSimon, Helen}
}
@article {Harris2010a,
	title = {{Thermal regime of the Costa Rican convergent margin: 1. Along-strike variations in heat flow from probe measurements and estimated from bottom-simulating reflectors}},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {11},
	number = {12},
	year = {2010},
	month = {dec},
	pages = {n/a{\textendash}n/a},
	abstract = {The thermal structure of convergent margins provides information related to the tectonics, geodynamics, metamorphism, and fluid flow of active plate boundaries. We report 176 heat flow measurements made with a violin bow style probe across the Costa Rican margin at the Middle America Trench. The probe measurements are collocated with seismic reflection lines. These seismic reflection lines show widespread distribution of bottom-simulating reflectors (BSRs). To extend the spatial coverage of heat flow measurements we estimate heat flow from the depth of BSRs. Comparisons between probe measurements and BSR-derived estimates of heat flow are generally within 10\% and improve with distance landward of the deformation front. Together, these determinations provide new information on the thermal regime of this margin. Consistent with previous studies, the margin associated with the northern Nicoya Peninsula is remarkably cool. We define better the southern boundary of the cool region. The northern extent of the cool region remains poorly determined. A regional trend of decreasing heat flow landward of the deformation front is apparent, consistent with the downward advection of heat by the subducting Cocos Plate. High wave number variability at a scale of 5-10 km is significantly greater than the measurement uncertainty and is greater south of the northern Nicoya Peninsula. These heat flow anomalies vary between approximately 20 and 60 mW m-2 and are most likely due to localized fluid flow through mounds and faults on the margin. Simple one-dimensional models show that these anomalies are consistent with flow rates of 7-15 mm yr-1. Across the margin toe variability is significant and likely due to fluid flow through deformation structures associated with the frontal sedimentary prism. Copyright 2010 by the American Geophysical Union.},
	keywords = {fluid flow, heat flow, Middle America Trench, subduction zones},
	issn = {15252027},
	doi = {10.1029/2010GC003272},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-78650544794\&partnerID=tZOtx3y1},
	author = {Harris, Robert N. and Grevemeyer, Ingo and Ranero, Cesar R. and Villinger, Heinrich and Barckhausen, Udo and Henke, Thomas and Mueller, Christian and Neben, Soenke}
}
@article {Harris2010,
	title = {{Thermal regime of the Costa Rican convergent margin: 2. Thermal models of the shallow Middle America subduction zone offshore Costa Rica}},
	journal = {Geochemistry, Geophysics, Geosystems},
	volume = {11},
	number = {12},
	year = {2010},
	month = {dec},
	pages = {n/a{\textendash}n/a},
	abstract = {At the Costa Rica margin along the Middle America Trench along-strike variations in heat flow are well mapped. These variations can be understood in terms of either ventilated fluid flow, where exposed basement allows fluids to freely advect heat between the crustal aquifer and ocean, or insulated fluid flow where continuous sediment cover restricts heat advection to within the crustal aquifer. We model fluid flow within the subducting aquifer using Nusselt number approximations coupled with finite element models of subduction and explore its effect on temperatures along the subduction thrust. The sensitivity of these models to the initial thermal state of the plate and styles of fluid flow, either ventilated or insulated, is explored. Heat flow measurements on cool crust accreted at the East Pacific Rise are consistent with ventilated hydrothermal cooling that continues with subduction. These models yield much cooler temperatures than predicted from simulations initialized with conductive predictions and without hydrothermal circulation. Heat flow transects on warm crust accreted at the Cocos-Nazca spreading center are consistent with models of insulated hydrothermal circulation that advects heat updip within the subducting crustal aquifer. Near the trench these models are warmer than conductive predictions and cooler than conductive predictions downdip of the trench. Comparisons between microseismicity and modeled isotherms suggest that the updip limit of microseismicity occurs at temperatures warmer than 100{\textdegree}C and that the downdip extent of microseismicity is bounded by the intersection of the subduction thrust with the base of the overriding crust. Copyright 2010 by the American Geophysical Union.},
	keywords = {fluid flow, subduction zones, thermal model},
	issn = {15252027},
	doi = {10.1029/2010GC003273},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-78650528793\&partnerID=tZOtx3y1},
	author = {Harris, Robert N. and Spinelli, Glenn and Ranero, Cesar R. and Grevemeyer, Ingo and Villinger, Heinrich and Barckhausen, Udo}
}
@article {Simancas2009,
	title = {{Tectonic relationships of Southwest Iberia with the allochthons of Northwest Iberia and the Moroccan Variscides}},
	journal = {Comptes Rendus Geoscience},
	volume = {341},
	number = {2-3},
	year = {2009},
	month = {feb},
	pages = {103{\textendash}113},
	abstract = {The Iberian Massif poses a problem of relationships between its northwestern and southern parts. Suture terranes (ophiolites and high-pressure rocks) crop out in NW Iberia but only as allochthonous units, unconnected from their root zone. Sutures cropping out in SW Iberia are discussed in order to relate them to the unknown root of the NW Iberia allochthons. On the other hand, the Moroccan Variscides are very briefly presented with a view to propose their correlation with the Iberian zones. Particularly important is the transition from the Variscides to the Paleoproterozoic basement in Morocco, which is a key argument for palaeogeographic reconstructions. {\textcopyright} 2008 Acad{\'e}mie des sciences.},
	keywords = {Correlations, Iberia, Morocco, Sutures, Variscan orogen},
	issn = {16310713},
	doi = {10.1016/j.crte.2008.11.003},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-63049089339\&partnerID=tZOtx3y1},
	author = {Simancas, Jos{\'e} Fernando and Azor, Antonio and Mart{\'\i}nez-Poyatos, David and Tahiri, Abdelfatah and El Hadi, Hassan and Gonz{\'a}lez-Lodeiro, Francisco and P{\'e}rez-Esta{\'u}n, Andr{\'e}s and Carbonell, Ram{\'o}n}
}
@article {Gutscher2009,
	title = {{Tectonic shortening and gravitational spreading in the Gulf of Cadiz accretionary wedge: Observations from multi-beam bathymetry and seismic profiling}},
	journal = {Marine and Petroleum Geology},
	volume = {26},
	number = {5},
	year = {2009},
	month = {may},
	pages = {647{\textendash}659},
	abstract = {The Gulf of Cadiz lies astride the complex plate boundary between Africa and Eurasia west of the Betic-Rif mountain belt. We report on the results of recent bathymetric swathmapping and multi-channel seismic surveys carried out here. The seafloor is marked by contrasting morphological provinces, spanning the SW Iberian and NW Moroccan continental margins, abyssal plains and an elongate, arcuate, accretionary wedge. A wide variety of tectonic and gravitational processes appear to have shaped these structures. Active compressional deformation of the wedge is suggested by folding and thrusting of the frontal sedimentary layers as well as basal duplexing in deeper internal units. There is evidence for simultaneous gravitational spreading occurring upslope. The very shallow mean surface and basal slopes of the accretionary wedge (1{\textdegree} each) indicate a very weak decollement layer, geometrically similar to the Mediterranean Ridge accretionary complex. Locally steep slopes (up to 10{\textdegree}) indicate strongly focused, active deformation and potential gravitational instabilities. The unusual surface morphology of the upper accretionary wedge includes "raft-tectonics" type fissures and abundant sub-circular depressions. Dissolution and/or diapiric processes are proposed to be involved in the formation of these depressions. {\textcopyright} 2008 Elsevier Ltd. All rights reserved.},
	keywords = {Accretionary wedge, Decollement, Gravitational spreading, Iberia, Morocco},
	issn = {02648172},
	doi = {10.1016/j.marpetgeo.2007.11.008},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-67349128699\&partnerID=tZOtx3y1},
	author = {Gutscher, Marc-Andr{\'e} and Dominguez, Stephane and Westbrook, Graham K. and Gente, Pascal and Babonneau, Nathalie and Mulder, Thierry and Gonthier, Eliane and Bartolom{\'e}, Rafael and Luis, Joaquim and Rosas, Filipe and Terrinha, Pedro}
}
@article {Latham2008,
	title = {{Three-dimensional particle shape acquisition and use of shape library for DEM and FEM/DEM simulation}},
	journal = {Minerals Engineering},
	volume = {21},
	number = {11},
	year = {2008},
	month = {oct},
	pages = {797{\textendash}805},
	abstract = {Numerical simulation that will capture the complex behaviour of rock fragment systems, e.g., in mining and civil engineering, requires both the computational mechanics capability to model particle interactions between complex shapes and an associated means to represent the kind of arbitrary or angular geometry relevant to problems involving rock fragments. This paper is concerned less with the modelling and more with the representation. Here, we focus on representation geared to {\textquoteright}soft contact{\textquoteright} modelling using either combined finite-discrete element (FEM/DEM) methods, or non-spherical DEM methods such as multi-sphere approximations of irregular geometry. 3D laser ranging (LADAR) is used to capture astonishingly realistic rock aggregate geometries. We report on the work flow procedures to generate computationally meshed virtual particles for modelling. The design of a shape library and a suggested procedure for selecting virtual particles for input to FEM/DEM or DEM models is discussed together with the use of inertia moments for shape descriptors. Use of the shape library for shape descriptor analysis is also illustrated. DEM simulations of packing using irregular particles from the shape library are presented. {\textcopyright} 2008 Elsevier Ltd. All rights reserved.},
	keywords = {Discrete element modelling, Gravity concentration, Mineral processing, Simulation, Sorting methods},
	issn = {08926875},
	doi = {10.1016/j.mineng.2008.05.015},
	url = {http://www.scopus.com/inward/record.url?eid=2-s2.0-50149086948\&partnerID=tZOtx3y1},
	author = {Latham, John-Paul and Munjiza, Antonio and Garcia, Xavier and Xiang, Jiansheng and Guises, Romain}
}
