Basalt–trachybasalt samples in Gale Crater,Mars

The ChemCam instrument on the Mars Science Laboratory (MSL) rover, Curiosity, observed numerous igneous float rocks and conglomerate clasts, reported previously. A new statistical analysis of single‐laser‐shot spectra of igneous targets observed by ChemCam shows a strong peak at ~55 wt% SiO₂ and 6 wt% total alkalis, with a minor secondary maximum at 47–51 wt% SiO₂ and lower alkali content. The centers of these distributions, together with the rock textures, indicate that many of the ChemCam igneous targets are trachybasalts, Mg# = 27 but with a secondary concentration of basaltic material, with a focus of compositions around Mg# = 54. We suggest that all of these igneous rocks resulted from low‐pressure, olivine‐dominated fractionation of Adirondack (MER) class‐type basalt compositions. This magmatism has subalkaline, tholeiitic affinities. The similarity of the basalt endmember to much of the Gale sediment compositions in the first 1000 sols of the MSL mission suggests that this type of Fe‐rich, relatively low‐Mg#, olivine tholeiite is the dominant constituent of the Gale catchment that is the source material for the fine‐grained sediments in Gale. The similarity to many Gusev igneous compositions suggests that it is a major constituent of ancient Martian magmas, and distinct from the shergottite parental melts thought to be associated with Tharsis and the Northern Lowlands. The Gale Crater catchment sampled a mixture of this tholeiitic basalt along with alkaline igneous material, together giving some analogies to terrestrial intraplate magmatic provinces.     

Mapping an aggregate extraction site off the Eastern English Channel: A methodology in support of monitoring and management

Each year approximately 23–28 million tonnes of sand and gravel are removed from offshore sediments around England and Wales. This study was located in a licensed marine aggregate extraction site off Shoreham in the Eastern English Channel (EEC thereafter). Results from the multibeam survey showed the presence of dredged pits created by suction hopper dredging and elongated furrows created by trailer suction hopper dredging in the area where sand and gravel had been excavated. Electronic Monitoring System (EMS) contained the dredging intensity recorded annually at the site; this information was combined with particle size data providing interpreted maps, which informed the status of the sediments at the site.     

Anisotropic P-wave attenuation measured from a multi-azimuth surface seismic reflection survey

A system of aligned vertical fractures produces azimuthal variations in stacking velocity and amplitude variation with offset, characteristics often reported in seismic reflection data for hydrocarbon exploration. Studies of associated attenuation anisotropy have been mostly theoretical, laboratory or vertical seismic profiling based. We used an 11 common‐midpoint‐long portion of each of four marine surface‐seismic reflection profiles, intersecting each other at 45° within circa 100 m of a common location, to measure the azimuthal variation of effective attenuation, Q⁻¹_eff and stacking velocity, in a shallow interval, about 100 m thick, in which consistently orientated vertical fracturing was expected due to an underlying salt diapirism. We found qualitative and quantitative consistency between the azimuthal variation in the attenuation and stacking velocity, and published amplitude variation with offset results. The 135° azimuth line showed the least apparent attenuation (1000 Q⁻¹_eff= 16 ± 7) and the fastest stacking velocity, hence we infer it to be closest to the fracture trend: the orthogonal 45° line showed the most apparent attenuation (1000Q⁻¹_eff= 52 ± 15) and slowest stacking velocity. The variation of Q⁻¹_eff with azimuth φ is well fitted by 1000Q⁻¹_eff = 34 − 18cos[2(φ+40°)] giving a fracture direction of 140 ± 23° (±1SD, derived from ‘bootstrapping’ fits to all 11⁴ combinations of individual common‐midpoint/azimuth measurements), compared to 134 ± 47° from published amplitude variation with offset data. The effects of short‐window spectral estimation and choices of spectral ratio bandwidth and offset ranges used in attenuation analysis, individually give uncertainties of up to ±13° in fracture direction. This magnitude of azimuthal variation can be produced by credible crack geometries (e.g., dry cracks, radius 6.5 m, aspect ratio 3 × 10⁻⁵, crack density 0.2) but we do not claim these to be the actual properties of the interval studied, because of the lack of well control (and its consequences for the choice of theoretical model and host rock physical properties) and the small number of azimuths available here.     

The 15 August 2007 Peru Earthquake and Tsunami: Influence of the Source Characteristics on the Tsunami Heights

The tsunami caused by the 2007 Peru earthquake (M_w 8.0) provoked less damage than by the seismic shaking itself (numerous casualties due to the earthquake in the vicinity of Pisco). However, it propagated across the Pacific Ocean and small waves were observed on one tide gauge in Taiohae Bay (Nuku Hiva, Marquesas, French Polynesia). We invert seismological data to recover the rupture pattern in two steps. The first step uses surface waves to find a solution for the moment tensor, and the second step uses body waves to compute the slip distribution in the source area. We find the slip distribution to consist of two main slip patches in the source area. The inversion of surface waves yields a scalar moment of 8.9 10²⁰ Nm, and body-wave inversion gives 1.4 10²¹ Nm. The inversion of tsunami data recorded on a single deep ocean sensor also can be used to compute a fault slip pattern (yielding a scalar moment of 1.1 10²¹ Nm). We then use these different sources to model the tsunami propagation across the Pacific Ocean, especially towards Nuku Hiva. While the source model taken from the body-wave inversion yields computed tsunami waves systematically too low with respect to observations (on the central Pacific Ocean DART buoy as on the Polynesian tide gauge), the source model established from the surface-wave inversion is more efficient to fit the observations, confirming that the tsunami is sensitive to the low frequency component of the source. Finally we also discuss the modeling of the late tsunami arrivals in Taiohae Bay using several friction coefficients for the sea bottom.     

Numerical Study of Tsunami Generated by Multiple Submarine Slope Failures in Resurrection Bay, Alaska, during the M W 9.2 1964 Earthquake

We use a viscous slide model of Jiang and LeBlond (1994) coupled with nonlinear shallow water equations to study tsunami waves in Resurrection Bay, in south-central Alaska. The town of Seward, located at the head of Resurrection Bay, was hit hard by both tectonic and local landslide-generated tsunami waves during the M _W 9.2 1964 earthquake with an epicenter located about 150 km northeast of Seward. Recent studies have estimated the total volume of underwater slide material that moved in Resurrection Bay during the earthquake to be about 211 million m³. Resurrection Bay is a glacial fjord with large tidal ranges and sediments accumulating on steep underwater slopes at a high rate. Also, it is located in a seismically active region above the Aleutian megathrust. All these factors make the town vulnerable to locally generated waves produced by underwater slope failures. Therefore it is crucial to assess the tsunami hazard related to local landslide-generated tsunamis in Resurrection Bay in order to conduct comprehensive tsunami inundation mapping at Seward. We use numerical modeling to recreate the landslides and tsunami waves of the 1964 earthquake to test the hypothesis that the local tsunami in Resurrection Bay has been produced by a number of different slope failures. We find that numerical results are in good agreement with the observational data, and the model could be employed to evaluate landslide tsunami hazard in Alaska fjords for the purposes of tsunami hazard mitigation.     

A meshless method to simulate solute transport in heterogeneous porous media

We derive a meshless numerical method based on smoothed particle hydrodynamics (SPH) for the simulation of conservative solute transport in heterogeneous geological formations. We demonstrate that the new proposed scheme is stable, accurate, and conserves global mass. We evaluate the performance of the proposed method versus other popular numerical methods for the simulation of one- and two-dimensional dispersion and two-dimensional advective–dispersive solute transport in heterogeneous porous media under different Pèclet numbers. The results of those benchmarks demonstrate that the proposed scheme has important advantages over other standard methods because of its natural ability to control numerical dispersion and other numerical artifacts. More importantly, while the numerical dispersion affecting traditional numerical methods creates artificial mixing and dilution, the new scheme provides numerical solutions that are “physically correct”, greatly reducing these artifacts.     

First lacustrine varve chronologies from Mexico: impact of droughts,ENSO and human activity since AD 1840 as recorded in maar sediments from Valle de Santiago

We present varve chronologies for sediments from two maar lakes in the Valle de Santiago region (Central Mexico): Hoya La Alberca (AD 1852–1973) and Hoya Rincón de Parangueo (AD 1839–1943). These are the first varve chronologies for Mexican lakes. The varved sections were anchored with tephras from Colima (1913) and Paricutín (1943/1944) and ²¹⁰Pb ages. We compare the sequences using the thickness of seasonal laminae and element counts (Al, Si, S, Cl, K, Ti, Mn, Fe, and Sr) determined by micro X-ray fluorescence spectrometry. The formation of the varve sublaminae is attributed to the strongly seasonal climate regime. Limited rainfall and high evaporation rates in winter and spring induce precipitation of carbonates (high Ca, Sr) enriched in ¹³C and ¹⁸O, whereas rainfall in summer increases organic and clastic input (plagioclase, quartz) with high counts of lithogenic elements (K, Al, Ti, and Si). Eolian input of Ti occurs also in the dry season. Moving correlations (5-yr windows) of the Ca and Ti counts show similar development in both sequences until the 1930s. Positive correlations indicate mixing of allochthonous Ti and autochthonous Ca, while negative correlations indicate their separation in sublaminae. Negative excursions in the correlations correspond with historic and reconstructed droughts, El Niño events, and positive SST anomalies. Based on our data, droughts (3–7 year duration) were severe and centred around the following years: the early 1850s, 1865, 1880, 1895, 1905, 1915 and the late 1920s with continuation into the 1930s. The latter dry period brought both lake systems into a critical state making them susceptible to further drying. Groundwater overexploitation due to the expansion of irrigation agriculture in the region after 1940 induced the transition from calcite to aragonite precipitation in Alberca and halite infiltration in Rincón. The proxy data indicate a faster response to increased evaporation for Rincón, the lake with the larger maar dimensions, solar radiation receipt and higher conductivity, whereas the smaller, steeper Alberca maar responded rapidly to increased precipitation.     

First incoherent scatter radar observations of ionospheric heating on the second electron gyro-harmonic

We report first results from a unique experiment performed at the HIPAS ionospheric modification facility in conjunction with the Poker Flat incoherent scatter radar in Alaska. High-power radio waves at 2.85 MHz, which corresponds to the second electron gyro-harmonic at ~245 km altitude, were transmitted into the nighttime ionosphere. Clear evidence of F-region ionospheric electron temperature enhancements were found, for the first time at this pump frequency, maximizing when the pump frequency is close to the second gyro-harmonic and double resonance. This is consistent with previous pump-enhanced artificial optical observations. We estimate the plasma heating efficiency to be approximately double that for higher pump frequencies.