Shaking and splashing—A case study of far-field effects of the Mjølnir asteroid impact on depositional environments in the Barents Sea

The Mjølnir impact crater in the Norwegian Barents Sea features among the 20 largest impact craters listed in the Earth Impact Database. The impact is dated to 142 ± 2.6 Ma, corresponding closely to the Jurassic/Cretaceous boundary in the Boreal stratigraphy. Multidisciplinary studies carried out over the last three decades have suggested that the up to 40 km wide crater was created by a 1–3 km diameter impactor colliding with a shallow epicontinental sea, causing regional havoc and a regional ecological crisis that followed in its wake. Only minor evidence for the consequences of the impact for the surrounding depositional basins has been documented so far. This study describes a large submarine slump penetrated by hydrocarbon exploration well 7121/9-1, located in the southern Hammerfest Basin and approximately 350 km away from the impact site. The slump is dated by a black shale drape, which contains characteristic impact-related biotic assemblages and potential ejecta material. This precise dating enables us to associate the slump with large-scale fault movements and footwall collapse along the basin-bounding Troms-Finnmark Fault Complex, which we conclude were caused by shock waves from the Mjølnir impact and the passage of associated tsunami trains. The draping black shale is interpreted to represent significant reworking of material from the contemporary seabed by tsunamis and currents set up by the impact.     

Scalar advection schemes for ocean modelling on unstructured triangular grids

Several schemes for scalar advection on unstructured triangular grids are assessed for possible use in ocean modelling applications. Finite element, finite volume and finite volume–element approaches are evaluated. A series of tests, including a numerical order of convergence analysis, idealized rotating cone and cylinder experiments, and transport of a tracer through the Stommel Gyre representation of ocean basin-scale circulation, are carried out. Volume element Eulerian–Lagrangian and third-order Runge-Kutta discontinuous Galerkin schemes are recommended for use in tracer studies. Taylor–Galerkin and second-order Runge–Kutta discontinuous Galerkin are found to be robust and accurate second-order schemes. When positivity is required, a fluctuation redistribution scheme was found to be an easily implemented, accurate, and computationally efficient approach. Responsible editor: Phil Dyke     

Velocity variations in carbonate rocks due to dual porosity and wave-induced fluid flow

Stiffness variations in carbonates may be described as resulting from different concentrations of flat compliant pores or cracks, which can have a significant effect on the effective stiffness and acoustic properties (e.g., velocities and attenuations) of dry as well as saturated carbonates, although they carry extremely little porosity. As shown in this paper, the effects of dual porosity and wave-induced fluid flow or pore pressure communication may also play a significant role. On the basis of a previously published T-matrix approach to model the effective viscoelastic properties of cracked porous media, we illustrate the (frequency-dependent) effects of wave-induced fluid flow (mainly squirt flow) or pore pressure communication for a model structure consisting of a mixture of fluid-saturated porous grains and fluid-saturated cavities (vugs, etc.) that are embedded in a solid matrix associated with carbonates. We assume that the pores within the porous grains are decoupled from the pores in the solid matrix (and possibly saturated with different fluids) but that each pore system at the micro and/or mesoscale may or may not be connected. For each of four different connectivity models, we present numerical results for four different cases of microstructure (that emphasize the importance of cracks and flat compliant pores). Our numerical results indicate that the velocity and attenuation spectra of carbonates vary significantly, even when the crack density and all other volume concentrations are constant.     

Glacial and environmental changes over the last 60 000 years in the Polar Ural Mountains,Arctic Russia,inferred from a high‐resolution lake record and other observations from adjacent areas

Our knowledge about the glaciation history in the Russian Arctic has to a large extent been based on geomorphological mapping supplemented by studies of short stratigraphical sequences found in exposed sections. Here we present new geochronological data from the Polar Ural Mountains along with a high‐resolution sediment record from Bolshoye Shchuchye, the largest and deepest lake in the mountain range. Seismic profiles show that the lake contains a 160‐m‐thick sequence of unconsolidated lacustrine sediments. A well‐dated 24‐m‐long core from the southern end of the lake spans the last 24 cal. ka. From downward extrapolation of sedimentation rates we estimate that sedimentation started about 50–60 ka ago, most likely just after a large glacier had eroded older sediments from the basin. Terrestrial cosmogenic nuclide (TCN) exposure dating (¹⁰Be) of boulders and Optically Stimulated Luminescence (OSL) dating of sediments indicate that this part of the Ural Mountains was last covered by a coherent ice‐field complex during Marine Isotope Stage (MIS) 4. A regrowth of the glaciers took place during a late stage of MIS 3, but the central valleys remained ice free until the present. The presence of small‐ and medium‐sized glaciers during MIS 2 is reflected by a sequence of glacial varves and a high sedimentation rate in the lake basin and likewise from ¹⁰Be dating of glacial boulders. The maximum extent of the mountain glaciers during MIS 2 was attained prior to 24 cal. ka BP. Some small present‐day glaciers, which are now disappearing completely due to climate warming, were only slightly larger during the Last Glacial Maximum (LGM) as compared to AD 1953. A marked decrease in sedimentation rate around 18–17 cal. ka BP indicates that the glaciers then became smaller and probably disappeared altogether around 15–14 cal. ka BP.     

Geodynamic significance of ophiolites within the Calabrian Arc

Abstract Slices of oceanic lithosphere belonging to the neo‐Tethys realm crop out discontinuously in the northern Calabrian Arc, Southern Apennines. They consist of high‐pressure–low‐temperature metamorphic ophiolitic sequences formed from metaultramafics, metabasites and alternating metapelites, metarenites, marbles and calcschist. Ophiolites occupy an intermediate position in the northern Calabrian Arc nappe pile, situated between overlying Hercynian continental crust and the underlying Apenninic limestone units. In the literature, these ophiolitic sequences are subdivided into several tectonometamorphic units. Geochemical characteristics indicate that metabasites were derived from subalkaline basalts with tholeiitic affinity (transitional mid‐oceanic ridge basalt type), and a harzburgitic‐lherzolitic protolith is suggested for the serpentinites. The pressure–temperature‐deformation paths of the metabasites from different outcrops display similar features: (i) the prograde segment follows a typical Alpine geothermal gradient up to a metamorphic climax at 350°C and 0.9 GPa and crystallization of the high‐pressure mineral assemblage occurs along a pervasive foliation developed during a compressive tectonic event; and (ii) the retrogression path can be subdivided in two segments, the first is characterized by nearly isothermal decompression to approximately 400°C and 0.3 GPa and the second follows a cooling trajectory. During low‐pressure conditions, a second deformation event produces millimetric to decametric scale asymmetric folds that describe west‐verging major structures. The third deformation event is characterized by brittle extensional structures. The tectonometamorphic evolution of the ophiolitic sequences from the different outcrops is similar. Both thermobarometric modeling and tectonic history indicate that the studied rocks underwent Alpine subduction and exhumation processes as tectonic slices inside a west‐verging accretionary wedge. The subduction of oceanic lithosphere was towards the present east; therefore, the Hercynian continental crust, overthrusted on the ophiolitic accretionary wedge after the neo‐Tethys closure, was part of the African paleomargin or a continental microplate between Africa and Europe.     

Eddy Diffusion Coefficients of Suspended Particulate Matter: Effects of Wind Energy Transfer and Stratification

Gross sedimentation rates (GSR) have been measured using sediment traps placed at nine different levels above the bed (0·3, 0·5, 0·8, 1·0, 2·0, 4·0, 6·0, 8·0 and 10·0 m). The sediment traps were deployed for 1·25 years and recovered 28 times during the study period. Low average GSR values of 5·5 g m^-2 day^-1 were obtained at 10·0 m, and high average GSR values of 114·8 g m^-2 day^-1 were obtained at 0·3 m. An expression for the eddy diffusion coefficient of suspended particulate matter (K_s), based on the measured GSR is given. The expression has been used for modelling of K_s at the different trap levels above the bed. High values (≈42 cm² s^-1) of K_s were obtained at the upper traps, whereas low values (≈2 cm² s^-1) were obtained near the bed. Comparison between level of turbulent energy in terms of shear stress at the boundaries of the water column, i.e. from the wind and the bed flow, showed that wind energy exceeded that of the bed flow by a factor 16. At 5·0 m K_s was positively correlated (r=0·66) to the eddy diffusion coefficient of momentum (K_m) derived from the wind energy transfer to the water, giving an average β of 0·5 for K_s =βK_m. The density difference between surface and bottom waters has been designated a parameter of stratification, and is discussed in relation to variations of K_s and K_m .     

Assessment of Erythemal UV Level in Nepal Based on Solar UV Estimates from Total Ozone Mapping Spectrometer

Nepal lies on the southern slope of Himalaya in Asia. In a width ranging between 150 and 250 km, the altitude varies greatly from about 100 m at its southern border to a maximum of 8848 min the northern part. Like the variation in altitude, climatic condition varies quite a lot. Long-term monthly mean erythemal UV daily dose values for Nepal are evaluated using Total Ozone Mapping Spectrometer (TOMS) estimation from the time of its overpass between 1996 and 2003. The results are presented as summer and winter maps of mean UV levels in each satellite grid. The mean winter erythemal UV daily dose ranges between 2.1 and 3.6 kJ m-2 whereas summer values are found to lie between 4.6 and 9.7 kJ m-2. The altitude variation increases the UV levels by about 0.2 kJ km-1 in winter months, and 0.9 kJ km-1 in summer. A multiyear monthly average erythemal daily dose in most of the areas shows that the summer value is about three times higher than that in winter. Although year-to-year variation is not pronounced in high- and mid-elevation regions, UV levels seemed to decrease from 1997 to 2002 in the southern part of the country in the low elevation region by about 5.35%. Due to the combined effects of the altitude, low ozone concentration in the troposphere, and thin air, surface UV radiation at higher altitudes is found to be higher than in the surrounding regions.     

The Zn abundance and isotopic composition of diatom frustules,a proxy for Zn availability in ocean surface seawater

We have developed cleaning methods for extracting diatomopal from bulk marine sediment samples, for measurement of both zinc (Zn) abundance and isotope composition. This cleaning technique was then applied to a set of Holocene core-top samples from the Southern Ocean. The measured δ⁶⁶Zn (reported relative to the JMC_Lyon standard) and Zn/Si ratios from the Southern Ocean diatom_opal samples range from 0.7 to 1.5‰, and from 14 to 0.9 μmol/mol, respectively. The Zn abundance and isotope composition data show a clear correlation with opal burial rates and other oceanographic parameters. In common with previous work, we interpret the systematic changes in the Zn/Si ratio to be linked to the variability in the concentrations of bioavailable Zn in the ambient surface seawater where the diatom opal is formed. This variability is likely to be primarily controlled by the degree to which Zn is taken up into phytoplankton biomass. The observed systematic pattern in the δ⁶⁶Zn compositions of the diatom_opal core-top samples is, similarly, likely to reflect changes in the δ⁶⁶Zn composition of the ambient Zn in the surface waters above the core-top sites, which is progressively driven towards isotopically heavier values by preferential incorporation of the lighter isotopes into phytoplankton organic material. Thus, the systematic relationship between Zn isotopes and abundance observed in the core-top diatom_opal samples suggests a potential tool for investigating the biogeochemical cycling of Zn in the past surface ocean for down-core diatom_opal material. In this respect, it may be possible to test hypotheses that attribute variations in atmospheric CO₂ on glacial–interglacial timescales to the degree to which trace metals limited primary productivity in HNLC zones.