The Gibraltar Arc seismogenic zone (part 2): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by tsunami modeling and seismic intensity

The Great Lisbon earthquake has the largest documented felt area of any shallow earthquake and an estimated magnitude of 8.5–9.0. The associated tsunami ravaged the coast of SW Portugal and the Gulf of Cadiz, with run-up heights reported to have reached 5–15 m. While several source regions offshore SW Portugal have been proposed (e.g.— Gorringe Bank, Marquis de Pombal fault), no single source appears to be able to account for the great seismic moment as well as all the historical tsunami amplitude and travel time observations. A shallow east dipping fault plane beneath the Gulf of Cadiz associated with active subduction beneath Gibraltar, represents a candidate source for the Lisbon earthquake of 1755.Here we consider the fault parameters implied by this hypothesis, with respect to total slip, seismic moment, and recurrence interval to test the viability of this source. The geometry of the seismogenic zone is obtained from deep crustal studies and can be represented by an east dipping fault plane with mean dimensions of 180 km (N–S) × 210 km (E–W). For 10 m of co-seismic slip an Mw 8.64 event results and for 20 m of slip an Mw 8.8 earthquake is generated. Thus, for convergence rates of about 1 cm/yr, an event of this magnitude could occur every 1000–2000 years. Available kinematic and sedimentological data are in general agreement with such a recurrence interval. Tsunami wave form modeling indicates a subduction source in the Gulf of Cadiz can partly satisfy the historical observations with respect to wave amplitudes and arrival times, though discrepancies remain for some stations. A macroseismic analysis is performed using site effect functions calculated from isoseismals observed during instrumentally recorded strong earthquakes in the region (M7.9 1969 and M6.8 1964). The resulting synthetic isoseismals for the 1755 event suggest a subduction source, possibly in combination with an additional source at the NW corner of the Gulf of Cadiz can satisfactorily explain the historically observed seismic intensities. Further studies are needed to sample the turbidites in the adjacent abyssal plains to better document the source region and more precisely calibrate the chronology of great earthquakes in this region.     

Late Pleistocene highstand and recession of a small, high-altitude pluvial lake, Jakes Valley, central Great Basin, USA

Models of factors controlling late Pleistocene pluvial lake-level fluctuations in the Great Basin are evaluated by dating lake levels in Jakes Valley. “Jakes Lake” rose to a highstand at 13,870 ± 50 ¹⁴C Yr B.P., receded to a stillstand at 12,440 ± 50 ¹⁴C yr B.P., and receded steadily to desiccation thereafter. The Jakes Lake highstand is roughly coincident with highstands of lakes Bonneville, Lahontan and Russell. The rise to highstand and recession of Jakes Lake were most likely controlled by a storm track steered by the polar jet stream. The final stillstand of Jakes Lake helps constrain timing of northward retreat of the polar jet stream during the Pleistocene-Holocene transition.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Are Natural Climate Forcings Able to Counteract the Projected Anthropogenic Global Warming?

A two-dimensional global climate model is used to assessthe climatic changes associated with the new IPCC SRES emissions scenarios and to determine which kind of changes in total solar irradiance and volcanic perturbations could mask the projected anthropogenic global warming associated to the SRES scenarios. Our results suggest that only extremely unlikely changes in total solar irradiance and/or volcanic eruptions would be able to overcome the simulated anthropogenic global warming over the century. Nevertheless, in the critical interval of the next two decades the externally-driven natural climate variability might possibly confuse the debate about temperature trends and impede detection of the anthropogenic climate change signal.     

Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Simultaneous use of geochemical and geophysical methods to characterise abandoned landfills

. Groundwater sampling and geophysical methods determined a serious contamination problem associated with refilled exploitation sites at the Cal Dimoni area, Llobregat delta, Barcelona, Spain. To characterise this process, hydrogeochemical analyses were performed and showed the following modifications to groundwater chemical composition: increasing pH values, changing redox conditions, significant increases in total organic carbon (TOC) and certain trace elements, and high groundwater conductivity values. Major ion content accumulations were found under the refilled area. In contrast, elements involved in the oxidation–reduction processes, such as iron, manganese and nitrates, clearly diminished. Electromagnetic prospecting methods were also performed and delineated the contamination plume extent. These methods also showed separate sources of contamination, one clearly related to the groundwater–refilled zone leachate interaction, another as a consequence of the manure–accumulation surface site. Geochemical and geophysical methods have shown similar results for locating groundwater contamination sources, and for determining leachate generation mechanisms and flow paths.     

Crystallization conditions of Los Angeles, a basaltic Martian meteorite

The texture of Los Angeles (stone 1) is dominated by relatively large (0.5−2.0 mm) anhedral to subhedral grains of pyroxene, and generally subhedral to euhedral shocked plagioclase feldspar (maskelynite). Minor phases include subhedral titanomagnetite and ilmenite, Fe-rich olivine, olivine+augite-dominated symplectites [some of which include a Si-rich phase and some which do not], pyrrhotite, phosphate(s), and an impact shock-related alkali- and silica-rich glass closely associated with anhedral to euhedral silica grains. Observations and model calculations indicate that the initial crystallization of Mg-rich pigeonitic pyroxenes at ≤1150 °C, probably concomitantly with plagioclase, was followed by pigeonitic and augitic compositions between 1100 and 1050 °C whereas between 1050 and 920 to 905 °C pyroxene of single composition crystallized. Below 920 to 905 °C, single composition Fe-rich clinopyroxene exsolved to augite and pigeonite. Initial appearance of titanomagnetite probably occurred near 990 °C and FMQ-1.5 whereas at and below 990 °C and ≥FMQ-1.5 titanomagnetite and single composition Fe-rich clinopyroxene may have started to react, producing ilmenite and olivine. However, judging from the most common titanomagnetite compositions, we infer that most of this reaction likely occurred between 950 and 900 °C at FMQ-1.0±0.2 and nearly simultaneously with pyroxene exsolution, thus producing assemblages of pigeonite, titanomagnetite, olivine, ilmenite, and augite. We deem this reaction as the most plausible explanation for the formation of the olivine+augite-dominated symplectites in Los Angeles. But we cannot preclude possible contributions to the symplectites from the shock-related alkali- and silica-rich glass or shocked plagioclase, and the breakdown of Fe-rich pigeonite compositions to olivine+augite+silica below 900 °C. Reactions between Fe-Ti oxides and silicate minerals in Los Angeles and other similar basaltic Martian meteorites can control the T-fO₂ equilibration path during cooling, which may better explain the relative differences in fO₂ among the basaltic Martian meteorites.     

Spatial and temporal variation of the nekton and hyperbenthos from a temperate European estuary with regulated freshwater inflow

The aquatic macrofauna of the Guadalquivir estuary were sampled (1 mm mesh persiana net) at 5 sampling sites located along the entire (except the tidal freshwater region) estuarine gradient of salinity (outer 50 km). A total of 134 fish and macroinvertebrate species was collected but only 62 were considered common or regularly present in the estuary. Univariate measures of the community structure showed statistically significant differences among sampling sites: species richness, abundance, and biomass decreased in the upstream direction, being positively correlated with the salinity. Temporal differences of these three variables were also statistically significant. While a clear seasonal pattern (minimum densities in winter and maximum in spring-summer) was observed for abundance and biomass, no such pattern existed for the number of species. Mysids was the most dominant group throughout the estuary (96% to 99% of abundance; 49% to 85% of biomass), although fish biomass was also important at the outer estuary (36% to 38%). Multivariate analyses indicated highly significant spatial variation in the macrofaunal communities observed along the salinity gradient. These analyses suggest that the underlying structure was a continuum with more or less overlapping distributions of the species dependent on their ability to tolerate different physicochemical conditions. There were also significant temporal (intermonthly + interannual) variation of the estuarine community; the relative multivariate dispersion indicated that monthly variation was more considerable (relative multivariate dispersion >1) at the outer part of the estuary during the wet year (last 20 km) and was higher in the inner stations during the dry year (32 to 50 km from the river mouth). Since a clear negative exponential relationship was observed between the freshwater input (from a dam located 110 km upstream) and water salinity at all sampling stations, it is concluded that the human freshwater management is probably affecting the studied estuarine communities. While the higher seasonal (long-term) stability of the salinity gradient, due to the human control of the freshwater input, may facilitate the recruitment of marine species juveniles during the meteorologically unstable early-spring, the additional (short-term) salinity fluctuations during the warm period may negatively affect species that complete their lifecycle within the estuary.