An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

Regional magma plumbing and emplacement mechanisms of the Faroe‐Shetland Sill Complex: implications for magma transport and petroleum systems within sedimentary basins

The movement of magma through the shallow crust and the impact of subsurface sill complexes on the hydrocarbon systems of prospective sedimentary basins has long been an area of interest and debate. Based on 3D seismic reflection and well data, we present a regional analysis of the emplacement and magmatic plumbing system of the Palaeogene Faroe‐Shetland Sill Complex (FSSC), which is intruded into the Mesozoic and Cenozoic sequences of the Faroe‐Shetland Basin (FSB). Identification of magma flow directions through detailed seismic interpretation of approximately 100 sills indicates that the main magma input zones into the FSB were controlled primarily by the NE–SW basin structure that compartmentalise the FSB into its constituent sub‐basins. An analysis of well data shows that potentially up to 88% of sills in the FSSC are <40 m in thickness, and thus below the vertical resolution limit of seismic data at depths at which most sills occur. This resolution limitation suggests that caution needs to be exercised when interpreting magmatic systems from seismic data alone, as a large amount of intrusive material could potentially be missed. The interaction of the FSSC with the petroleum systems of the FSB is not well understood. Given the close association between the FSSC and potential petroleum migration routes into some of the oil/gas fields (e.g. Tormore), the role the intrusions may have played in compartmentalisation of basin fill needs to be taken fully into account to further unlock the future petroleum potential of the FSB.     

Modelling landscape evolution on geological time scales: a new method based on irregular spatial discretization

We present simulations of large-scale landscape evolution on tectonic time scales obtained from a new numerical model which allows for arbitrary spatial discretization. The new method makes use of efficient algorithms from the field of computational geometry to compute the set of natural neighbours of any irregular distribution of points in a plane. The natural neighbours are used to solve geomorphic equations that include erosion/deposition by channelled flow and diffusion. The algorithm has great geometrical flexibility, which makes it possible to solve problems involving complex boundaries, radially symmetrical uplift functions and horizontal tectonic transport across strike-slip faults. The algorithm is also ideally suited for problems which require large variations in spatial discretization and/or self-adaptive meshing. We present a number of examples to illustrate the power of the new approach and its advantages over more 'classical' models based on regular (rectangular) discretization. We also demonstrate that the synthetic river networks and landscapes generated by the model obey the laws of network composition and have scaling properties similar to those of natural landscapes. Finally we explain how orographically controlled precipitation and flexural isostasy may be easily incorporated in the model without sacrificing efficiency.
     

Volcanism in the Sumisu Rift, I. Major element, volatile, and stable isotope geochemistry

A bimodal volcanic suite with KAr ages of 0.05–1.40 Ma was collected from the Sumisu Rift using alvin. These rocks are contemporaneous with island arc tholeiite lavas of the Izu-Ogasawara arc 20 km to the east, and provide a present day example of volcanism associated with arc rifting and back-arc basin initiation. Major element geochemistry of the basalts is most similar to that of basalts found in other, more mature back-arc basins, which indicates that back-arc basins need not begin their magmatic evolution with lavas bearing strong arc signatures.Volatile concentrations distinguish Sumisu Rift basalts from island arc basalts and MORB. H₂O contents, which are at least four times greater than in MORB, suppress plagioclase crystallization. This suppression results in a more mafic fractionating assemblage, which prevents Al₂O₃ depletion and delays the initiation of Fe₂O_3^(tot) and TiO₂ enrichment. However, unlike arc basalts,Fe^3+/ΣFe ratios are only slightly higher than in MORB and are insufficient to cause magnetite saturation early enough to suppress Fe₂O_3^(tot) and TiO₂ enrichment. Thus, major element trends are more similar to those of MORB than arcs.H₂O, CO₂ and S are undersaturated relative to pure phase solubility curves, indicating exsolution of an H₂O-rich mixed gas phase. HighH₂O/S, highδD, and low (MORB-like)δ³⁴S ratios are considered primary and distinctive of the back-arc basin setting.     

Climatic and halokinetic controls on alluvial–lacustrine sedimentation during compressional deformation,Andean forearc,northern Chile

The Salar de Atacama forms one of a series of forearc basins developed along the western flank of the Central Andes. Exposed along the northwest margin of the basin, a salt‐cored range, the Cordillera de la Sal, records the Mid‐Miocene to recent sedimentological and structural development of this basin. Sediments of the Mid‐Miocene Vilama Formation record the complex interaction between regional/local climate change, halokinesis and compressional deformation. This study reveals how these factors have controlled the facies development and distribution within the Salar de Atacama. Detailed sedimentary logging, cross‐sections and present day geomorphology through the northern Cordillera de la Sal have been used to establish a lithostratigraphy, chronostratigraphy and the regional distribution of the Vilama Formation. The Vilama Formation documents an increase in aridity with a hiatus in sedimentation from Mid‐Miocene to 9 Ma with initial uplift of the Cordillera de la Sal. From 9 Ma to 8.5 Ma deposition of a meandering fluvial system is recorded followed by a rapid decrease in sedimentation till 6 Ma. From 6 to 2 Ma, the deposition of extensive palustrine carbonates and distal alluvial–mudflat–lacustrine demonstrates the existence of an extensive lake within the Salar de Atacama. Post 2 Ma, the lake decreased in size and braided alluvial gravels associated with alluvial fans were widespread through the region suggesting a final shift to hyperarid conditions. By comparing the Vilama Formation with similar age facies throughout northern Chile and southern Peru, several shifts in climate are recognized. Climate signatures within northern Chile appear to be largely diachronous with the last regional event in the Mid‐Miocene. Since that time, humid events have been restricted to either Precordillerian basins or the Central Atacama. Within the Central Atacama, the final switch to hyperarid conditions was not till the earliest Pleistocene, much later than previously estimated within the region.     

A nested numerical tidal model of the southern New England Bight

A nested, barotropic numerical tidal model has been applied to the southern New England Bight to include Long Island Sound (LIS), Block Island Sound (BIS), Rhode Island Sound (RIS), Buzzards Bay (BB) and a portion of the southern New England shelf. An explicit finite difference representation of the two-dimensional vertically averaged equations of motion and continuity, employing forward time and centered spatial differences with the bottom friction term being evaluated at both time levels, is applied to the study area to predict tidal circulation dynamics. Horizontal friction and advective terms have been neglected. Tidal height data along the open boundaries is specified from available data. Model predictions of corange and high and low cophase lines agree with available data and show the strongly co-oscillatory behavior of LIS. Tidal current ellipses agree reasonably well with available data showing rectilinear motion in LIS, BIS and BB and rotary motion in RIS.