Granoblastic olivine aggregates in magnesian chondrules: Planetesimal fragments or thermally annealed solar nebula condensates?

Granoblastic olivine aggregates (GOA) have been discovered in some Type I magnesian chondrules within carbonaceous chondrites by Libourel and Krot [Libourel, G., Krot, A.N., 2007. Evidence for the presence of planetesimal material among the precursors of magnesian chondrules of nebular origin. Earth Planet. Sci. Lett. 254, 1–8], who proposed an origin from pre-existing planetesimals. Amoeboid olivine aggregates (AOA), generally considered as aggregates of solar nebula condensates and found within similar carbonaceous chondrites, display similar equilibrium texture, though on a finer scale. For these reasons, we conducted experiments to determine if annealing of olivine required time scales appropriate to planetesimal or nebular heating. Pressed < 43 µm and < 63 µm San Carlos olivine powder (Fo_88.4) was isothermally heated at temperatures ranging from 1350–1550 °C for 1–100 h. The 100 h runs yield olivine aggregates with well-developed granoblastic texture at all temperatures, manifest as a network of randomly-oriented and sutured olivine grains with 120° triple junctions. Individual olivine grains are 4–6 sided and polygonal by 1450 °C and equigranular texture is developed at high temperature (1500–1550 °C). Melting of olivine commences at 1450–1500 °C and aids in ‘ripening’ and suturing (grain coarsening and grain boundary migration). Textural equilibrium is clearly met at 1550 °C. A planetesimal origin cannot be ruled out; however, the experimental evidence reveals that granoblastic texture can be reproduced in an interval not inconsistent with heating times for nebular objects. GOA may have experienced higher degrees of thermal processing than the finer-grained AOA. If the precursors were the same, grain coarsening would have to be accompanied by modification to bulk and isotopic compositions. However, the precursors could have been olivine condensates formed later than AOA. Annealing may have been a widespread process operating in the primordial solar nebula responsible for thermal processing and formation of GOA prior to their incorporation into chondrules.     

Methylhopane biomarker hydrocarbons in Hamersley Province sediments provide evidence for Neoarchean aerobiosis

Bacteriohopanepolyols with A-ring methyl substituents are bacterial lipid biomarkers that are presently known to originate from just a few specific taxa. When preserved in ancient rocks, the fossilized hydrocarbon skeletons of these molecules have the potential to provide biogeochemical records of those taxa and their respective physiologies. Considering both their occurrences in modern organisms, and their Proterozoic and Phanerozoic sedimentary distributions, hopanes carrying 2-methyl or 3-methyl substituents are proposed to be derived from cyanobacteria and oxygen-respiring methanotrophs, respectively. Here we report the distribution of methylhopanes in 2.72–2.56 billion-year-old, Neoarchean rocks from the Hamersley Province on the Pilbara Craton. The relative abundance of C₃₁-3β-methylhopane, but not that of C₃₁-2α-methylhopane, shows a strong correlation to the carbon isotopic composition of co-occurring kerogen (insoluble particulate organic matter). The unanticipated nature of this relationship provides evidence for a shallow-water locus of carbon cycling through aerobic oxidation of CH₄ and, coincidentally, a means to demonstrate biomarker syngenicity. 2α-Methylhopanes are most abundant in both shale and carbonate from shallow-water sediments, and a positive correlation to carbonate abundance, contrasts to variable, but generally lower, 2α-methylhopane abundances of deeper water facies. These observations are consistent with their origin from cyanobacteria which were likely the cornerstone of microbial communities in shallow-water ecosystems providing molecular oxygen, fixed carbon, and possibly fixed nitrogen.     

Organosilane occurrence in irghizite samples from the Zhamanshin impact crater,Kazakhstan

Abstract— The composition of surface deposits on vesicle walls in irghizites (i.e., impact glasses at site) from the Zhamanshin meteorite crater were studied using time‐of‐flight secondary ion mass spectrometry (TOF‐SIMS). The cavity walls are unique interfaces for condensation of gases from the superheated, high‐silica melt during the impact. Initially, signals from the cavity wall are dominated by hydrocarbon fragments whereas the glass fracture face surrounding the cavity gave only signals corresponding to glass components. After 12 h in ultra high vacuum (UHV), signals from the cavity wall are dominated by peaks corresponding to fragments normally measured from organosilanes and organosiloxanes with the majority of the hydrocarbon signals markedly reduced. Characteristic hydrocarbon fragments are now observed on the glass fracture surface next to the cavity in an annulus around the cavity perimeter. There are also minor signals in this region from organosilanes and organosiloxanes. In contrast, four tektites (Australites) (i.e., glassy distal ejecta) gave no organosilane or organosiloxane signals after the same preparation and vacuum evaporation procedure. These species appear to be formed only at the impact site where higher levels of organic material are likely to be present in soil and are trapped before evaporation. This appears to be the first report of naturally occurring silicon‐organic compounds.     

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.