Landslides induced by the April 2007 Aysén Fjord earthquake, Chilean Patagonia

On the 21st of April 2007, the Aysén Fjord earthquake (Mw 6.2) in southern Chile (45.3° S, 73.0° W) triggered hundreds of landslides in the epicentral area along the fjord coast and surroundings. Some of these landslides induced large tsunami waves within the fjord causing fatalities and damaging several salmon farms, the most important economic activity of the area. The landslides included rock slides and avalanches, rock falls, shallow soil and soil–rock slides, and debris flows. The earthquake was the climax of a seismic swarm that began 3 months earlier. The seismicity is associated with tectonic activity along the Liquiñe–Ofqui fault zone (LOFZ), a major structural feature of the region. The earthquake-induced landslides were mapped and classified from field observations and remote sensing analysis. The landslide areas and epicentral distances are within the expected range for the earthquake magnitude according to worldwide data, while the position of landslides on the slopes strongly suggests topographic amplification effects in triggering the failures. The location of the landslides is also clearly related to some of the main fault branches of the LOFZ. The seismic event has configured a new situation of seismic and landslide hazard in the Aysén region and along the LOFZ, where the presence of towns and economic infrastructure along the coasts of several fjords constitutes a potential risk that was not considered before this seismic event.     

Self-Similar and Self-Affine Properties of Two-Dimensional Fracture Patterns in Rocks

In this paper, we investigate the fractal properties of binary maps of rock fractures at different scales and different geological types, as well as different families of fracture patterns obtained from a two-dimensional Laplacian growth model (LGM). From these analyses we figure out which families of the LGM patterns match the structural properties of the fracture binary maps. The LGM is defined in terms of a nonlinear map that depends on two parameters, λ and \mathfraka\mathfrak{a}, that respectively define the area and shape of the elements of the aggregate that conforms the patterns. The fractal dimension and roughness exponent of the LGM patterns are found to depend on \mathfraka\mathfrak{a}, with 0 < \mathfrak a < 10<\mathfrak {a}<1. From a detailed statistical analysis of these patterns we found that the fractal dimensions of capacity, correlation and information decrease monotonically as \mathfraka\mathfrak{a} increases. We also found that the values of these three fractal dimensions tend to collapse on top of each other as \mathfraka\lessapprox1\mathfrak{a}\lessapprox1. Remarkably, the fractal properties of rock fractures in the scales from millimeters up to a few meters appear to be well represented by the fractal structure of the LGM families of patterns with \mathfraka=0.15\mathfrak{a}=0.15 and 0.30, while the fractal properties of rock fractures in the scale of kilometers seems to be well represented by the LGM family with \mathfraka=0.90\mathfrak{a}=0.90. In addition, the three fractal dimension values of fracture binary maps in the scales from millimeters up to meters were found to be different between them. Nonetheless, for fractures in the scale of kilometers, the values of the three fractal dimensions are very close to each other as an indication of self-similar behavior. Analysis of the corrections to the scaling of the roughness exponent, ζ, suggests that they are negligible for the LGM family of fracture patterns with \mathfraka=0.9\mathfrak{a}=0.9. This finding points to a self-affine structure for this family of patterns. In fact, the calculated roughness exponent results are in the range of values characteristic of rock fractures.     

Orbital structure of self-consistent triaxial stellar systems

We used a multipolar code to create, through the dissipationless collapses of systems of 1,000,000 particles, three self-consistent triaxial stellar systems with axial ratios corresponding to those of E4, E5 and E6 galaxies. The E5 and E6 models have small, but significant, rotational velocities although their total angular momenta are zero, that is, they exhibit figure rotation; the rotational velocity decreases with decreasing flattening of the models and for the E4 model it is essentially zero. Except for minor changes, probably caused by unavoidable relaxation effects, the systems are highly stable. The potential of each system was subsequently approximated with interpolating formulae yielding smooth potentials, stationary for the non-rotating model and stationary in the rotating frame for the rotating ones. The Lyapunov exponents could then be computed for randomly selected samples of the bodies that make up the different systems, allowing the recognition of regular and partially and fully chaotic orbits. Finally, the regular orbits were Fourier analyzed and classified using their locations on the frequency map. As it could be expected, the percentages of chaotic orbits increase with the flattening of the system. As one goes from E6 through E4, the fraction of partially chaotic orbits relative to that of fully chaotic ones increases, with the former surpassing the latter in model E4; the likely cause of this behavior is that triaxiality diminishes from E6 through E4, the latter system being almost axially symmetric. We especulate that some of the partially chaotic orbits may obey a global integral akin to the long axis component of angular momentum. Our results show that is perfectly possible to have highly stable triaxial models with large fractions of chaotic orbits, but such systems cannot have constant axial ratios from center to border: a slightly flattened reservoir of highly chaotic orbits seems to be mandatory for those systems.     

Geochemical processes and mobilization of toxic metals and metalloids in an As-rich base metal waste pile in Zimapán,Central Mexico

The geochemistry and mineralogy of samples collected along depth profiles from an As-rich tailing deposit with abundant calcite was studied to determine the processes that influence the mobility of Fe, Zn, Cu, Ni, Cd, As, Sb, Cr and Tl. In spite of their near neutral pH, almost all of them are acid potential generators. Total concentrations decreased as: Fe > As > Zn > Pb > Cu > Sb > Cd > Cr > Ni > Tl. Soluble contents were lower and followed a slightly different order. Mobility decreased as: Tl > Cd, Zn, Cu, Sb, Ni, As > Fe, Pb > Cr. Higher soluble concentrations of Fe, Cu, Zn, As, Pb, and Ni were found in low-pH samples and of Sb and Tl in near-neutral samples. Sulfide oxidation processes are developing in the tailing’s dam. These processes do not have a trend with depth but occur mainly in acid layers. Near neutral layers formed by primary sulfides and calcite probably correspond to wastes produced from the processing of ore coming mainly from pods within the skarn, and acid layers with abundant secondary minerals from material mined from chimneys and mantos. The presence of calcite influences speciation, neutralizes acid mine drainage (AMD), and decreases the mobility of most toxic metals and metalloids (TMMs). However, a hard-pan layer was not observed in the studied profiles. Retention of TMM within tailings probably occurs through the formation of low solubility metal carbonates and from elevation of pH that promotes Fe hydroxides precipitation that may retain As, Sb and metals. Calcite occurrence promotes As, Cd, Cu, Fe, Zn, Pb, Cd and Cr retention, does not play a role on Tl and Ni mobilization, and increases Sb release.     

Determining Hydraulic Conductivity Using Pumping Data from Low-Flow Sampling

Hydraulic conductivity values computed using the steady-state discharge and drawdown attained while low-flow sampling were evaluated to determine if they were equivalent to those determined from slug testing. Based on testing 12 wells, it was found that the results were statistically equivalent. Conductivity values computed using low-flow sampling parameters were also evaluated as to their reproducibility in actual practice by analyzing consultant data for three wells sampled over three quarterly monitoring periods by four field technicians. The results were found to be reproducible within about a factor of 2 or better. Since the method is based on only one pair of parameters, diligence is required in attaining steady state and in accurately measuring the flow rate and drawdown. Conductivity values computed using this approach can enhance the use of low-flow data gathered in water quality sampling, avoid the need for slug testing in a subsequent phase of investigation, and help reduce the cost of characterizing sites when multilevel samplers are used. Given the practical range of discharge in low-flow sampling, the method was found to be applicable at conductivity values somewhat greater than 10⁻⁶ cm/s. Given the typical accuracy of water level meters and pressure transducers and a maximum discharge of 1 L/min, as mandated by regulatory guidance, the method has a calculated upper conductivity limit in the range of 10⁻³ to 10⁻² cm/s.     

The continent‐ocean transition on the northwestern South China Sea

Rifted margins are created as a result of stretching and breakup of continental lithosphere that eventually leads to oceanic spreading and formation of a new oceanic basin. A cornerstone for understanding what processes control the final transition to seafloor spreading is the nature of the continent‐ocean transition (COT). We reprocessed multichannel seismic profiles and use available gravity data to study the structure and variability of the COT along the Northwest subbasin (NWSB) of the South China Sea. We have interpreted the seismic images to discern continental from oceanic domains. The continental‐crust domain is characterized by tilted fault blocks generally overlain by thick syn‐rift sedimentary units, and underlain by fairly continuous Moho reflections typically at 8–10 s twtt. The thickness of the continental crust changes greatly across the basin, from ~20 to 25 km under the shelf and uppermost slope, to ~9–6 km under the lower slope. The oceanic‐crust domain is characterized by a highly reflective top of basement, little faulting, no syntectonic strata and fairly constant thickness (over tens to hundreds of km) of typically 6 km, but ranging from 4 to 8 km. The COT is imaged as a ~5–10 km wide zone where oceanic‐type features directly abut or lap on continental‐type structures. The South China margin continental crust is cut by abundant normal faults. Seismic profiles show an along‐strike variation in the tectonic structure of the continental margin. The NE‐most lines display ~20–40 km wide segments of intense faulting under the slope and associated continental‐crust thinning, giving way to a narrow COT and oceanic crust. Towards the SW, faulting and thinning of the continental crust occurs across a ~100–110 km wide segment with a narrow COT and abutting oceanic crust. We interpret this 3D structural variability and the narrow COT as a consequence of the abrupt termination of continental rifting tectonics by the NE to SW propagation of a spreading centre. We suggest that breakup occurred abruptly by spreading centre propagation rather than by thinning during continental rifting. We propose a kinematic evolution for the oceanic domain of the NWSB consisting of a southward spreading centre propagation followed by a first narrow ridge jump to the north, and then a younger larger jump to the SE, to abandon the NWSB and create the East subbasin of the South China Sea.     

Large‐scale margin collapse during Messinian early sea‐level drawdown: the SW Valencia trough,NW Mediterranean

In this study, detailed mapping of the ‘Messinian markers’ and examination of their geometrical relationships in the SW Valencia trough (Western Mediterranean) have revealed the style and depositional processes associated with emersion of continental margins during the Messinian Salinity Crisis (MSC). Based on multichannel seismic profiles and well data, this article evidences the existence of two Messinian depositional units in intermediate basins (Complex Unit and Upper Unit) and four main Messinian erosional surfaces (Margin Erosion Surface, Bottom Surface, Top (Erosion) Surface and Intermediate Surface). Results show that (1) initial rapid sea‐level drawdown and exposure of the shelf and upper slope of the Valencia margin induced large‐scale destabilization of the continental slope and deposition of large detrital bodies at the base‐of‐slope in the form of major mass‐transport deposits (MTD); (2) as sea level continued to drop, the development of the Margin Erosion Surface attained full development on the margins and eroded the clastic units (MTDs) deposited during initial drawdown. At the same time, a submarine drainage network formed in the deepwater Valencia trough; (3) persistent lowstand and restrictive conditions in the area resulted in deposition of the evaporites that form the Upper Unit in the SW Valencia trough.     

Fish remains from the Formative Period (1000 BC–AD 400) of Lake Titicaca, Bolivia: Zooarchaeology and taphonomy

This paper presents a detailed zooarchaeological study of fish remains from the site of Kala Uyuni, Bolivia. The fish remains of 31 flotation samples from different archaeological contexts and dated to the Formative Period (1000 BC–AD 400) are described in terms of frequency, weight, taxa distribution, sizes, skeletal representations, and exposure to heat. The results confirm the importance of fish in the diet of the inhabitants of the site and provide insights about the particularities of their exploitation, consumption, and differential discard across the site throughout time. A strong decline in the procurement and consumption of fish is observed for the Late Formative which seems to be a consequence of both environmental and social processes.