Tidal currents,energy flux and bottom boundary layer thickness in the Clyde Sea and North Channel of the Irish Sea

A high-resolution three-dimensional model of the Clyde Sea and the adjacent North Channel of the Irish Sea is used to compute the major diurnal and semidiurnal tides in the region, the associated energy fluxes and thickness of the bottom boundary layer. Initially, the accuracy of the model is assessed by performing a detailed comparison of computed tidal elevations and currents in the region, against an extensive database that exists for the M₂, S₂, N₂, K₁ and O₁ tides. Subsequently, the model is used to compute the tidal energy flux vectors in the region. These show that the major energy flux is confined to the North Channel region, with little energy flux into the Clyde Sea. Comparison with the observed energy flux in the North Channel shows that its across-channel distribution and its magnitude are particularly sensitive to the phase difference between elevation and current. Consequently, small changes in the computed values of these parameters due to slight changes of the order of the uncertainty in the open-boundary values to the model, can significantly influence the computed energy flux. The thickness of the bottom boundary layer in the region is computed using a number of formulations. Depending upon the definition adopted, the empirical coefficient C used to determine its thickness varies over the range 0.1 to 0.3, in good agreement with values found in the literature. In the North Channel, the boundary layer thickness occupies the whole water depth, and hence tidal turbulence produced at the sea bed keeps the region well mixed. In the Clyde Sea, the boundary layer thickness is a small fraction of the depth, and hence the region stratifies.Responsible Editor: Phil Dyke     

Fast lagrangian particle tracking in unstructured ocean model grids

Ocean Dynamics - Lagrangian particle tracking, based on currents derived from hydrodynamic models, is an important tool in quantifying bio-physical transports in the ocean. Particle tracking in the...     

Temperature evolution between 50 K and 320 K of the thermal expansion tensor of gypsum derived from neutron powder diffraction data

The unit cell parameters, extracted from Rietveld analysis of neutron powder diffraction data collected between 4.2 K and 320 K, have been used to calculate the temperature evolution of the thermal expansion tensor for gypsum for 50 ≤ T ≤ 320 K. At 300 K the magnitudes of the principal axes are α ₁₁  = 1.2(6) × 10⁻⁶ K⁻¹, α ₂₂  = 36.82(1) × 10⁻⁶ K⁻¹ and α ₃₃  = 25.1(5) × 10⁻⁶ K⁻¹. The maximum axis, α ₂₂ , is parallel to b, and using Institution of Radio Engineers (IRE) convention for the tensor orthonormal basis, the axes α ₁₁ and α ₃₃ have directions equal to (−0.979, 0, 0.201) and (0.201, 0, 0.979) respectively. The orientation and temperature dependent behaviour of the thermal expansion tensor is related to the crystal structure in the I2/a setting. Received 12 February 1998 / Revised, accepted 19 October 1998     

Geological evidence for neotectonic activity during deglaciation of the southern Sperrin Mountains,Northern Ireland

In the southern Sperrin Mountains, Northern Ireland, stacked glacigenic sequences that accumulated during deglaciation (ca. 17000–13000 yr BP) overlie a basement of jointed and mascerated bedrock. The glacigenic sequences comprise interbedded glaciotectonic shears, diamictic breccias and rock rafts. At one site a normal fault with a metre‐scale vertical displacement cuts through part of the sequence and is overlain by a glacial diamict. Sediments at an adjacent site show that faulting and associated hydrothermal activity was related to neotectonic reactivation of pre‐existing Caledonian lineaments caused by ice unloading. From stratigraphical and directional evidence, fault reactivation occurred early in the deglaciation after north central Ireland ice had retreated southwards into lowland areas, but before Sperrin Mountain ice readvanced from the north. This relationship provides evidence for the relative timing of neotectonic activity in Northern Ireland, and demonstrates the effects of glacio‐isostatic unloading near ice‐sheet centres. Copyright © 1999 John Wiley & Sons, Ltd.     

Identification of the principal patterns of summer moisture transport in South America and their representation by WCRP/CMIP3 global climate models

The goal of this study is to assess the ability of a set of global climate models (GCMs) to represent the main regional spatial patterns of austral summer moisture transport in South America in order to evaluate if one of the possible causes of GCMs misestimating summer precipitation in this region could be associated with an erroneous representation of these patterns. For this purpose, NCEP/NCAR reanalysis and 20 GCMs from the WCRP/CMIP3 multi-model dataset for the period 1960–1999 were considered. Extreme cases of moisture transport patterns were selected to assess their association with rainfall anomalies. Results obtained indicate that only some aspects of water vapor transport and convergence in South America as well as the associated precipitation anomalies can be reproduced adequately by GCMs. Finally, a case study is presented showing that one of the moisture transport patterns identified was observed during December 2008–February 2009.     

NMR Logging to Estimate Hydraulic Conductivity in Unconsolidated Aquifers

Nuclear magnetic resonance (NMR) logging provides a new means of estimating the hydraulic conductivity (K) of unconsolidated aquifers. The estimation of K from the measured NMR parameters can be performed using the Schlumberger‐Doll Research (SDR) equation, which is based on the Kozeny–Carman equation and initially developed for obtaining permeability from NMR logging in petroleum reservoirs. The SDR equation includes empirically determined constants. Decades of research for petroleum applications have resulted in standard values for these constants that can provide accurate estimates of permeability in consolidated formations. The question we asked: Can standard values for the constants be defined for hydrogeologic applications that would yield accurate estimates of K in unconsolidated aquifers? Working at 10 locations at three field sites in Kansas and Washington, USA, we acquired NMR and K data using direct‐push methods over a 10‐ to 20‐m depth interval in the shallow subsurface. Analysis of pairs of NMR and K data revealed that we could dramatically improve K estimates by replacing the standard petroleum constants with new constants, optimal for estimating K in the unconsolidated materials at the field sites. Most significant was the finding that there was little change in the SDR constants between sites. This suggests that we can define a new set of constants that can be used to obtain high resolution, cost‐effective estimates of K from NMR logging in unconsolidated aquifers. This significant result has the potential to change dramatically the approach to determining K for hydrogeologic applications.     

The genus Pseudo-nitzschia (Bacillariophyceae) in New Zealand: analysis of the last decade's monitoring data

The diatom genus Pseudo-nitzschia blooms throughout New Zealand's coastal waters. More than a decade of monitoring data were analysed in this study and seasonality was a significant factor in the appearance of Pseudo-nitzschia blooms. Austral spring and summer were the main bloom periods, although there were regional differences. Between January 2000 and August 2011, 8.4% of 29,000 seawater samples analysed for Pseudo-nitzschia exceeded the voluntary trigger level for biotoxin testing (10×10⁴ cells l^?1), and these were analysed by fluorescent in situ hybridisation (FISH) assays. Pseudo-nitzschia fraudulenta and P. pseudodelicatissima (low domoic acid [DA] producers) were the dominant bloom formers throughout New Zealand (each contributing 25%) followed by the highly toxic P. australis (10%). Shellfish flesh testing for DA was triggered on 8477 occasions between 2000 and 2011, but no samples exceeded the regulatory limit (20 mg/kg); the maximum concentration was 13 mg/kg (Marlborough Sounds, mid-winter 2010).     

Rethinking communication in risk interpretation and action

We have studied the attenuation characteristics of eastern Himalaya and southern Tibet by using local earthquake data set that consists of 123 well-located events, recorded by the Himalayan Nepal Tibet Seismic Experiment operated during 2001–2003. We have used single backscattering model to calculate frequency-dependent values of coda Q (\(Q_\mathrm{c}\)). The estimation of \(Q_\mathrm{c}\) is made at central frequencies 2, 4, 8 and 12 Hz through five lapse time windows from 10 to 50 s starting at double the travel time of the S-wave. The observed \(Q_\mathrm{c}\) is found to be strongly frequency-dependent and follows a similar trend as observed in other tectonically active parts of the Himalaya. The trend of variation of \(Q_\mathrm{c}\) with lapse time and the corresponding apparent depths is also studied. Increase in \(Q_\mathrm{c}\) values with the lapse time suggests that the deeper part of the study region is less heterogeneous than the shallower part. The observed values of \(Q_0\) (\(Q_\mathrm{c}\) at 1 Hz) and frequency parameter n indicate that the medium beneath the study area is highly heterogeneous and tectonically very active. A regionalization of the estimated \(Q_0\) is carried out, and a contour map is prepared for the whole region. Some segments of Lesser Himalaya and Sub-Himalaya exhibit very low \(Q_0\) , while the whole Tethyan Himalaya and some parts of Greater Himalaya are characterized by low \(Q_0\) values. Our results are comparable with those obtained from tectonically active regions in the world.     

The coalescence and organization of lahars at Semeru volcano, Indonesia

We present multi-parameter geophysical measurements of rainfall-induced lahars at Semeru Volcano, East Java, using two observation sites 510 m apart, 11.5 km from the summit. Our study site in the Curah Lengkong channel is composed of a 30-m wide box-valley, with a base of gravel and lava bedrock, representing an ideal geometry for high density measurements of active lahars. Instrumentation included pore-pressure sensors (stage), a broad-band seismograph (arrival times, vibrational energy, and turbulence), video footage, and direct bucket sampling. A total of 8 rainfall-induced lahars were recorded, with durations of 1–3 h, heights 0.5–2 m, and peak velocities 3–6 m/s. Flow types ranged from dilute to dense hyperconcentrated flows. These recorded flows were commonly composed of partly coalesced, discrete and unsteady gravity current packets, represented by multiple peaks within each lahar. These packets most likely originate from multiple lahar sources, and can be traced between instrument sites. Those with the highest concentrations and greatest wetted areas were often located mid-lahar at our measured reach, accelerating towards the flow front. As these lahars travel downstream, the individual packets thus coalesce and the flow develops a more organised structure. Observations of different degrees of coalescence between these discrete flow packets illustrate that a single mature debris flow may have formed from multiple dynamically independent lahars, each with different origins.