Diurnal warm-layer events in the western Mediterranean and European shelf seas

We characterize near-surface ocean diurnal warm-layer events, using satellite observations and fields from numerical weather forecasting. The study covers April to September, 2006, over the area 11°W to 17°E and 35°N to 57°N, with 0.1° cells. We use hourly satellite SSTs from which peak amplitudes of diurnal cycles in SST (dSSTs) can be estimated with error ∼0.3 K. The diurnal excursions of SST observed are spatially and temporally coherent. The largest dSSTs exceed 6 K, affect 0.01% of the surface, and are seen in the Mediterranean, North and Irish Seas. There is an anti-correlation between the magnitude and the horizontal length scale of dSST events. Events wherein dSST exceeds 4 K have length scales of ≤40 km. From the frequency distribution of different measures of wind-speed minima, we infer that extreme dSST maxima arise where conditions of low wind speed are sustained from early morning to mid afternoon.

Horizontal and vertical structure of the Eyjafjallajökull ash cloud over the UK: a comparison of airborne lidar observations and simulations

During April and May 2010 the ash cloud from the eruption of the Icelandic volcano Eyjafjallajökull caused widespread disruption to aviation over northern Europe. The location and impact of the eruption led to a wealth of observations of the ash cloud were being obtained which can be used to assess modelling of the long range transport of ash in the troposphere. The UK FAAM (Facility for Airborne Atmospheric Measurements) BAe-146-301 research aircraft overflew the ash cloud on a number of days during May. The aircraft carries a downward looking lidar which detected the ash layer through the backscatter of the laser light. In this study ash concentrations derived from the lidar are compared with simulations of the ash cloud made with NAME (Numerical Atmospheric-dispersion Modelling Environment), a general purpose atmospheric transport and dispersion model. The simulated ash clouds are compared to the lidar data to determine how well NAME simulates the horizontal and vertical structure of the ash clouds. Comparison between the ash concentrations derived from the lidar and those from NAME is used to define the fraction of ash emitted in the eruption that is transported over long distances compared to the total emission of tephra. In making these comparisons possible position errors in the simulated ash clouds are identified and accounted for. The ash layers seen by the lidar considered in this study were thin, with typical depths of 550–750 m. The vertical structure of the ash cloud simulated by NAME was generally consistent with the observed ash layers, although the layers in the simulated ash clouds that are identified with observed ash layers are about twice the depth of the observed layers. The structure of the simulated ash clouds were sensitive to the profile of ash emissions that was assumed. In terms of horizontal and vertical structure the best results were obtained by assuming that the emission occurred at the top of the eruption plume, consistent with the observed structure of eruption plumes. However, early in the period when the intensity of the eruption was low, assuming that the emission of ash was uniform with height gives better guidance on the horizontal and vertical structure of the ash cloud. Comparison of the lidar concentrations with those from NAME show that 2–5% of the total mass erupted by the volcano remained in the ash cloud over the United Kingdom.

Evaluation of boundary-layer type in a weather forecast model utilising long-term Doppler lidar observations

Many studies evaluating model boundary-layer schemes focus either on near-surface parameters or on short-term observational campaigns. This reflects the observational datasets that are widely available for use in model evaluation. In this paper we show how surface and long-term Doppler lidar observations, combined in a way to match model representation of the boundary layer as closely as possible, can be used to evaluate the skill of boundary-layer forecasts. We use a 2-year observational dataset from a rural site in the UK to evaluate a climatology of boundary layer type forecast by the UK Met Office Unified Model. In addition, we demonstrate the use of a binary skill score (Symmetric Extremal Dependence Index) to investigate the dependence of forecast skill on season, horizontal resolution and forecast leadtime. A clear diurnal and seasonal cycle can be seen in the climatology of both the model and observations, with the main discrepancies being the model overpredicting cumulus capped and decoupled stratocumulus capped boundary-layers and underpredicting well mixed boundary-layers. Using the SEDI skill score the model is most skillful at predicting the surface stability. The skill of the model in predicting cumulus capped and stratocumulus capped stable boundary layer forecasts is low but greater than a 24 hr persistence forecast. In contrast, the prediction of decoupled boundary-layers and boundary-layers with multiple cloud layers is lower than persistence. This process based evaluation approach has the potential to be applied to other boundary-layer parameterisation schemes with similar decision structures.

Parametrizing the horizontal inhomogeneity of ice water content using CloudSat data products

In order to calculate unbiased microphysical and radiative quantities in the presence of a cloud, it is necessary to know not only the mean water content but also the distribution of this water content. This article describes a study of the in-cloud horizontal inhomogeneity of ice water content, based on CloudSat data. In particular, by focusing on the relations with variables that are already available in general circulation models (GCMs), a parametrization of inhomogeneity that is suitable for inclusion in GCM simulations is developed. Inhomogeneity is defined in terms of the fractional standard deviation (FSD), which is given by the standard deviation divided by the mean. The FSD of ice water content is found to increase with the horizontal scale over which it is calculated and also with the thickness of the layer. The connection to cloud fraction is more complicated; for small cloud fractions FSD increases as cloud fraction increases while FSD decreases sharply for overcast scenes. The relations to horizontal scale, layer thickness and cloud fraction are parametrized in a relatively simple equation. The performance of this parametrization is tested on an independent set of CloudSat data. The parametrization is shown to be a significant improvement on the assumption of a single-valued global FSD

Impact of Laptev Sea flaw polynyas on the atmospheric boundary layer and ice production using idealized mesoscale simulations

The interaction between polynyas and the atmospheric boundary layer is examined in the Laptev Sea using the regional, non-hydrostatic Consortium for Small-scale Modelling (COSMO) atmosphere model. A thermodynamic sea-ice model is used to consider the response of sea-ice surface temperature to idealized atmospheric forcing. The idealized regimes represent atmospheric conditions that are typical for the Laptev Sea region. Cold wintertime conditions are investigated with sea-ice–ocean temperature differences of up to 40 K. The Laptev Sea flaw polynyas strongly modify the atmospheric boundary layer. Convectively mixed layers reach heights of up to 1200 m above the polynyas with temperature anomalies of more than 5 K. Horizontal transport of heat expands to areas more than 500 km downstream of the polynyas. Strong wind regimes lead to a more shallow mixed layer with strong near-surface modifications, while weaker wind regimes show a deeper, well-mixed convective boundary layer. Shallow mesoscale circulations occur in the vicinity of ice-free and thin-ice covered polynyas. They are forced by large turbulent and radiative heat fluxes from the surface of up to 789 W m−2, strong low-level thermally induced convergence and cold air flow from the orographic structure of the Taimyr Peninsula in the western Laptev Sea region. Based on the surface energy balance we derive potential sea-ice production rates between 8 and 25 cm d−1. These production rates are mainly determined by whether the polynyas are ice-free or covered by thin ice and by the wind strength.

Relationships between dendritic morphology, spatial distribution and firing patterns in rat layer 1 neurons

The cortical layer 1 contains mainly small interneurons, which have traditionally been classified according to their axonal morphology. The dendritic morphology of these cells, however, has received little attention and remains ill defined. Very little is known about how the dendritic morphology and spatial distribution of these cells may relate to functional neuronal properties. We used biocytin labeling and whole cell patch clamp recordings, associated with digital reconstruction and quantitative morphological analysis, to assess correlations between dendritic morphology, spatial distribution and membrane properties of rat layer 1 neurons. A total of 106 cells were recorded, labeled and subjected to morphological analysis. Based on the quantitative patterns of their dendritic arbor, cells were divided into four major morphotypes: horizontal, radial, ascendant, and descendant cells. Descendant cells exhibited a highly distinct spatial distribution in relation to other morphotypes, suggesting that they may have a distinct function in these cortical circuits. A significant difference was also found in the distribution of firing patterns between each morphotype and between the neuronal populations of each sublayer. Passive membrane properties were, however, statistically homogeneous among all subgroups. We speculate that the differences observed in active membrane properties might be related to differences in the synaptic input of specific types of afferent fibers and to differences in the computational roles of each morphotype in layer 1 circuits. Our findings provide new insights into dendritic morphology and neuronal spatial distribution in layer 1 circuits, indicating that variations in these properties may be correlated with distinct physiological functions.

Does the plankton community follow the horizontal water quality heterogeneity in a tropical urban reservoir (Guarapiranga reservoir, Sao Paulo, Brazil)?

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The importance of vertical and horizontal dimensions of the sediment matrix in structuring nematodes across spatial scales

Intensive surveys have been conducted to unravel spatial patterns of benthic infauna communities. Although it has been recognized that benthic organisms are spatially structured along the horizontal and vertical dimensions of the sediment, little is known on how these two dimensions interact with each other. In this study we investigated the interdependence between the vertical and horizontal dimensions in structuring marine nematodes assemblages. We tested whether the similarity in nematode species composition along the horizontal dimension was dependent on the vertical layer of the sediment. To test this hypothesis, three-cm interval sediment samples (15 cm depth) were taken independently from two bedforms in three estuaries. Results indicated that assemblages living in the top layers are more abundant, species rich and less variable, in terms of species presence/absence and relative abundances, than assemblages living in the deeper layers. Results showed that redox potential explained the greatest amount (12%) of variability in species composition, more than depth or particle size. The fauna inhabiting the more oxygenated layers were more homogeneous across the horizontal scales than those from the reduced layers. In contrast to previous studies, which suggested that reduced layers are characterized by a specific set of tolerant species, the present study showed that species assemblages in the deeper layers are more causal (characterized mainly by vagrant species). The proposed mechanism is that at the superficial oxygenated layers, species have higher chances of being resuspended and displaced over longer distances by passive transport, while at the deeper anoxic layers they are restricted to active dispersal from the above and nearby sediments. Such restriction in the dispersal potential together with the unfavorable environmental conditions leads to randomness in the presence of species resulting in the high variability between assemblages along the horizontal dimension.

On the seasonal variations of salinity of the tropical Atlantic mixed layer

The physical processes controlling the mixed layer salinity (MLS) seasonal budget in the tropical Atlantic Ocean are investigated using a regional configuration of an ocean general circulation model. The analysis reveals that the MLS cycle is generally weak in comparison of individual physical processes entering in the budget because of strong compensation. In evaporative regions, around the surface salinity maxima, the ocean acts to freshen the mixed layer against the action of evaporation. Poleward of the southern SSS maxima, the freshening is ensured by geostrophic advection, the vertical salinity diffusion and, during winter, a dominant contribution of the convective entrainment. On the equatorward flanks of the SSS maxima, Ekman transport mainly contributes to supply freshwater from ITCZ regions while vertical salinity diffusion adds on the effect of evaporation. All these terms are phase locked through the effect of the wind. Under the seasonal march of the ITCZ and in coastal areas affected by river (7°S:15°N), the upper ocean freshening by precipitations and/or runoff is attenuated by vertical salinity diffusion. In the eastern equatorial regions, seasonal cycle of wind forced surface currents advect freshwaters, which are mixed with subsurface saline water because of the strong vertical turbulent diffusion. In all these regions, the vertical diffusion presents an important contribution to the MLS budget by providing, in general, an upwelling flux of salinity. It is generally due to vertical salinity gradient and mixing due to winds. Furthermore, in the equator where the vertical shear, associated to surface horizontal currents, is developed, the diffusion depends also on the sheared flow stability.

Biomass of seston in the benthic layer of the Southwestern Indian Ocean

At stations to 1530 m depth in the Mozambique Channel and on the Saya-de-Malha and Walters banks seston biomass 2 m above the bottom was lower than at 30 m. Above the Walters shoal this difference was 13.2 mg/m**3 and was not equal to zero for P < 0.001. These results contradict previous ideas of biomass increase in benthic layers. The most likely cause of the observed impoverishment of plankton may be predominant consumption of living zooplankton component of seston by bottom and near-bottom predators. In the area of the Walters shoal this consumption is estimated as being about 300 mg/m**2 per day. Animals inhabiting this area live mainly on plankton brought in by horizontal advection, so that existence of faunal assemblages even on shallow-water submarine elevations is supported not mainly by local photosynthesis, but by primary production of surrounding waters.

Instabilities in two-liquid layers subject to a horizontal temperature gradient

We study theoretically the stability of two superposed fluid layers heated laterally. The fluids are supposed to be immiscible, the interface undeformable and of infinite horizontal extension. Combined thermocapillary and buoyancy forces give rise to a basic flow when a temperature difference is applied. The calculations are performed for a melt of GaAs under a layer of molten B2 O3 , a configuration of considerable technological importance. Four dif- ferent flow patterns and five temperature configurations are found for the basic state in this system. A linear stability analysis shows that the basic state may be destabilized by oscilla- tory motions leading to the so-called hydrothermal waves. Depending on the relative height of the two layers these hydrothermal waves propagate parallel or perpendicular to the temperature gradient. This analysis reveals that these perturbations can alter significantly the liquid flow in the liquid-encapsulated crystal growth techniques.

(Table 2) Abundance of the main groups of planktonic animals in the 250-900 m water layer at the Broken Spur and Lost City sites according to visual counts from the Mir manned submersibles in 2002 and 2005

Three dives of the Mir manned submersibles with plankton counts and two vertical plankton hauls with a BR net were carried out above the Lost City (Atlantis underwater massif) and the Broken Spur hydrothermal fields during cruise 50 of R/V Akademik Mstislav Keldysh. Above the Atlantis seamount no significant increase in plankton concentration was found. Above the Lost City field horizontal heterogeneity of plankton distribution in the near-bottom layer and in overlying water layers was shown. Near-bottom aggregations of euphausiids and amphipods previously reported by other scientists seem to be related to attraction of these animals by the submersible's headlights rather than represent a natural phenomenon.

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.

Complementary methods to investigate the development of clogging within a horizontal sub-surface flow tertiary treatment wetland

A combination of experimental methods was applied at a clogged, horizontal subsurface flow (HSSF) municipal wastewater tertiary treatment wetland (TW) in the UK, to quantify the extent of surface and subsurface clogging which had resulted in undesirable surface flow. The three dimensional hydraulic conductivity profile was determined, using a purpose made device which recreates the constant head permeameter test in-situ. The hydrodynamic pathways were investigated by performing dye tracing tests with Rhodamine WT and a novel multi-channel, data-logging, flow through Fluorimeter which allows synchronous measurements to be taken from a matrix of sampling points. Hydraulic conductivity varied in all planes, with the lowest measurement of 0.1 md1 corresponding to the surface layer at the inlet, and the maximum measurement of 1550 md1 located at a 0.4m depth at the outlet. According to dye tracing results, the region where the overland flow ceased received five times the average flow, which then vertically short-circuited below the rhizosphere. The tracer break-through curve obtained from the outlet showed that this preferential flow-path accounted for approximately 80% of the flow overall and arrived 8 h before a distinctly separate secondary flow-path. The overall volumetric efficiencyof the clogged system was 71% and the hydrology was simulated using a dual-path, dead-zone storage model. It is concluded that uneven inlet distribution, continuous surface loading and high rhizosphere resistance is responsible for the clog formation observed in this system. The average inlet hydraulic conductivity was 2 md1, suggesting that current European design guidelines, which predict that the system will reach an equilibrium hydraulic conductivity of 86 md1, do not adequately describe the hydrology of mature systems.