Dynamics of convectively driven banded jets in the laboratory

The banded organization of clouds and zonal winds in the atmospheres of the outer planets has long fascinated observers. Several recent studies in the theory and idealized modeling of geostrophic turbulence have suggested possible explanations for the emergence of such organized patterns, typically involving highly anisotropic exchanges of kinetic energy and vorticity within the dissipationless inertial ranges of turbulent flows dominated (at least at large scales) by ensembles of propagating Rossby waves. The results from an attempt to reproduce such conditions in the laboratory are presented here. Achievement of a distinct inertial range turns out to require an experiment on the largest feasible scale. Deep, rotating convection on small horizontal scales was induced by gently and continuously spraying dense, salty water onto the free surface of the 13-m-diameter cylindrical tank on the Coriolis platform in Grenoble, France. A “planetary vorticity gradient” or “β effect” was obtained by use of a conically sloping bottom and the whole tank rotated at angular speeds up to 0.15 rad s−1. Over a period of several hours, a highly barotropic, zonally banded large-scale flow pattern was seen to emerge with up to 5–6 narrow, alternating, zonally aligned jets across the tank, indicating the development of an anisotropic field of geostrophic turbulence. Using particle image velocimetry (PIV) techniques, zonal jets are shown to have arisen from nonlinear interactions between barotropic eddies on a scale comparable to either a Rhines or “frictional” wavelength, which scales roughly as (β/Urms)−1/2. This resulted in an anisotropic kinetic energy spectrum with a significantly steeper slope with wavenumber k for the zonal flow than for the nonzonal eddies, which largely follows the classical Kolmogorov k−5/3 inertial range. Potential vorticity fields show evidence of Rossby wave breaking and the presence of a “hyperstaircase” with radius, indicating instantaneous flows that are supercritical with respect to the Rayleigh–Kuo instability criterion and in a state of “barotropic adjustment.” The implications of these results are discussed in light of zonal jets observed in planetary atmospheres and, most recently, in the terrestrial oceans.

Global observed changes in daily climate extremes of temperature and precipitation

A suite of climate change indices derived from daily temperature and precipitation data, with a primary focus on extreme events, were computed and analyzed. By setting an exact formula for each index and using specially designed software, analyses done in different countries have been combined seamlessly. This has enabled the presentation of the most up-to-date and comprehensive global picture of trends in extreme temperature and precipitation indices using results from a number of workshops held in data-sparse regions and high-quality station data supplied by numerous scientists world wide. Seasonal and annual indices for the period 1951-2003 were gridded. Trends in the gridded fields were computed and tested for statistical significance. Results showed widespread significant changes in temperature extremes associated with warming, especially for those indices derived from daily minimum temperature. Over 70% of the global land area sampled showed a significant decrease in the annual occurrence of cold nights and a significant increase in the annual occurrence of warm nights. Some regions experienced a more than doubling of these indices. This implies a positive shift in the distribution of daily minimum temperature throughout the globe. Daily maximum temperature indices showed similar changes but with smaller magnitudes. Precipitation changes showed a widespread and significant increase, but the changes are much less spatially coherent compared with temperature change. Probability distributions of indices derived from approximately 200 temperature and 600 precipitation stations, with near-complete data for 1901-2003 and covering a very large region of the Northern Hemisphere midlatitudes (and parts of Australia for precipitation) were analyzed for the periods 1901-1950, 1951-1978 and 1979-2003. Results indicate a significant warming throughout the 20th century. Differences in temperature indices distributions are particularly pronounced between the most recent two periods and for those indices related to minimum temperature. An analysis of those indices for which seasonal time series are available shows that these changes occur for all seasons although they are generally least pronounced for September to November. Precipitation indices show a tendency toward wetter conditions throughout the 20th century.

Synthesis of galacto-oligosaccharide from lactose using beta-galactosidase from Kluyveromyces lactis: studies on batch and continuous UF membrane-fitted bioreactors

A study of galacto-oligosaccharides (GOS) synthesis from lactose with beta-galactosidase from Kluyveromyces lactis (Maxilact(R) L2000) was carried out. The synthesis was performed using various initial lactose concentrations ranging from 220 to 400 mg/mL and enzyme concentrations ranging from 3 to 9 U/mL, and was investigated at 40degreesC and pH 7, in a stirred-tank reactor. In the experimental range examined, the results showed the amount of GOS formed depended on lactose concentration but not on enzyme concentration. Galactose was a competitive inhibitor, while glucose was a non-competitive inhibitor. In a further study, a laboratory-scale reactor system, fitted with a 10-kDa NMWCO composite regenerated cellulose membrane, was used in a continuous process. The reactor was operated in cross-flow mode. The effect of operating pressures on flux and productivity was investigated by applying different transmembrane pressures to the system. The continuous process showed better production performance compared to the batch synthesis with the same lactose and enzyme concentrations at 40degreesC, pH 7. Comparison of product structures from batch and continuous processes, analyzed by HPAEPAD and methylation analysis, showed similarities but differed from the structures found in a commercial GOS product (Vivinal(R)GOS). (C) 2004 Wiley Periodicals, Inc.

A stable one-step-ahead predictive control of non-linear systems

In this paper stability of one-step ahead predictive controllers based on non-linear models is established. It is shown that, under conditions which can be fulfilled by most industrial plants, the closed-loop system is robustly stable in the presence of plant uncertainties and input–output constraints. There is no requirement that the plant should be open-loop stable and the analysis is valid for general forms of non-linear system representation including the case out when the problem is constraint-free. The effectiveness of controllers designed according to the algorithm analyzed in this paper is demonstrated on a recognized benchmark problem and on a simulation of a continuous-stirred tank reactor (CSTR). In both examples a radial basis function neural network is employed as the non-linear system model.

Combined bias and outlier identification in dynamic data reconciliation

Measured process data normally contain inaccuracies because the measurements are obtained using imperfect instruments. As well as random errors one can expect systematic bias caused by miscalibrated instruments or outliers caused by process peaks such as sudden power fluctuations. Data reconciliation is the adjustment of a set of process data based on a model of the process so that the derived estimates conform to natural laws. In this paper, techniques for the detection and identification of both systematic bias and outliers in dynamic process data are presented. A novel technique for the detection and identification of systematic bias is formulated and presented. The problem of detection, identification and elimination of outliers is also treated using a modified version of a previously available clustering technique. These techniques are also combined to provide a global dynamic data reconciliation (DDR) strategy. The algorithms presented are tested in isolation and in combination using dynamic simulations of two continuous stirred tank reactors (CSTR).

Advancements in the estimation of ice particle fall speeds using laboratory and field measurements

Accurate estimates for the fall speed of natural hydrometeors are vital if their evolution in clouds is to be understood quantitatively. In this study, laboratory measurements of the terminal velocity vt for a variety of ice particle models settling in viscous fluids, along with wind-tunnel and field measurements of ice particles settling in air, have been analyzed and compared to common methods of computing vt from the literature. It is observed that while these methods work well for a number of particle types, they fail for particles with open geometries, specifically those particles for which the area ratio Ar is small (Ar is defined as the area of the particle projected normal to the flow divided by the area of a circumscribing disc). In particular, the fall speeds of stellar and dendritic crystals, needles, open bullet rosettes, and low-density aggregates are all overestimated. These particle types are important in many cloud types: aggregates in particular often dominate snow precipitation at the ground and vertically pointing Doppler radar measurements. Based on the laboratory data, a simple modification to previous computational methods is proposed, based on the area ratio. This new method collapses the available drag data onto an approximately universal curve, and the resulting errors in the computed fall speeds relative to the tank data are less than 25% in all cases. Comparison with the (much more scattered) measurements of ice particles falling in air show strong support for this new method, with the area ratio bias apparently eliminated.

A comparative study of multilayered and single layered neural network based predictive controllers

The authors compare the performance of two types of controllers one based on the multilayered network and the other based on the single layered CMAC network (cerebellar model articulator controller). The neurons (information processing units) in the multi-layered network use Gaussian activation functions. The control scheme which is considered is a predictive control algorithm, along the lines used by Willis et al. (1991), Kambhampati and Warwick (1991). The process selected as a test bed is a continuous stirred tank reactor. The reaction taking place is an irreversible exothermic reaction in a constant volume reactor cooled by a single coolant stream. This reactor is a simplified version of the first tank in the two tank system given by Henson and Seborg (1989).

Hydrochemical processes in lowland rivers: insights from in situ, high-resolution monitoring

This paper introduces new insights into the hydrochemical functioning of lowland river systems using field-based spectrophotometric and electrode technologies. The streamwater concentrations of nitrogen species and phosphorus fractions were measured at hourly intervals on a continuous basis at two contrasting sites on tributaries of the River Thames – one draining a rural catchment, the River Enborne, and one draining a more urban system, The Cut. The measurements complement those from an existing network of multi-parameter water quality sondes maintained across the Thames catchment and weekly monitoring based on grab samples. The results of the sub-daily monitoring show that streamwater phosphorus concentrations display highly complex dynamics under storm conditions dependent on the antecedent catchment wetness, and that diurnal phosphorus and nitrogen cycles occur under low flow conditions. The diurnal patterns highlight the dominance of sewage inputs in controlling the streamwater phosphorus and nitrogen concentrations at low flows, even at a distance of 7 km from the nearest sewage treatment works in the rural River Enborne. The time of sample collection is important when judging water quality against ecological thresholds or standards. An exhaustion of the supply of phosphorus from diffuse and multiple septic tank sources during storm events was evident and load estimation was not improved by sub-daily monitoring beyond that achieved by daily sampling because of the eventual reduction in the phosphorus mass entering the stream during events. The results highlight the utility of sub-daily water quality measurements and the discussion considers the practicalities and challenges of in situ, sub-daily monitoring.

Nitrogen speciation and phosphorus fractionation dynamics in a lowland chalk catchment

A detailed analysis of temporal and spatial trends in nitrogen (N) speciation and phosphorus (P) fractionation in the Wylye, a lowland Chalk sub-catchment of the Hampshire Avon, UK is presented, identifying the sources contributing to nutrient enrichment, and temporal variability in the fractionation of nutrients in transit from headwaters to lower reaches of the river. Samples were collected weekly from ten monitoring stations with daily sampling at three further sites over one year, and monthly inorganic N and total reactive P (TRP) concentrations at three of the ten weekly monitoring stations over a ten year period are also presented. The data indicate significant daily and seasonal variation in nutrient fractionation in the water column, resulting from plant uptake of dissolved organic and inorganic nutrient fractions in the summer months, increased delivery of both N and P from diffuse sources in the autumn to winter period and during high flow events, and lack of dilution of point source discharges to the Wylye from septic tank, small package Sewage Treatment Works (STW) and urban Waste Water Treatment Works (WwTW) during the summer low flow period. Weekly data show that contributing source areas vary along the river with headwater N and P strongly influenced by diffuse inorganic N and particulate P fluxes, and SRP and organic-rich point source contributions from STW and WwTW having a greater influence in the lower reaches. Long-term data show a decrease in TRP concentrations at all three monitoring stations, with the most pronounced decrease occurring downstream from Warminster WwTW, following the introduction of P stripping at the works in 2001. Inorganic N demonstrates no statistically significant change over the ten year period of record in the rural headwaters, but an increase in the lower reaches downstream from the WwTW which may be due to urban expansion in the lower catchment.

The consolidation and bond strength of rafted sea ice

The consolidation and bond strength of rafted sea ice were investigated through a series of experiments undertaken in the Ice Physics Laboratory at the UCL. To simulate a section of rafted sea ice, blocks of laboratory grown saline ice were stacked in an insulated tank with spacers between adjacent blocks to allow saline water to flood in. The rate of consolidation was then monitored using a combination of temperature readings recorded in the ice and liquid layer, salinity measurements of the liquid layer, and cores taken at specific times of interest. Two states of consolidation were observed: thermodynamic consolidation where the ice blocks were physically bonded but the bond strength was weak, and mechanical consolidation where the bond had reached full strength. Results showed that the rafted ice had physically bonded in less than a day, however it took many more days (6 to 30 depending on the environmental conditions) for the bond to reach maximum strength. Increasing the thickness of the ice, the salinity of the water and the inter-block gap size all increased the consolidation time. Once consolidated, ice cores were taken and sheared using the asymmetric four-point bending method to measure the strength of the bond between the ice blocks. These were then compared to the shear strength of solid ice blocks simulating level sea ice. Our results show that the shear strength of the bond between the rafted ice blocks is about 30% weaker than that of level ice.

A rate and state friction law for saline ice

Sea ice friction models are necessary to predict the nature of interactions between sea ice floes. These interactions are of interest on a range of scales, for example, to predict loads on engineering structures in icy waters or to understand the basin-scale motion of sea ice. Many models use Amonton's friction law due to its simplicity. More advanced models allow for hydrodynamic lubrication and refreezing of asperities; however, modeling these processes leads to greatly increased complexity. In this paper we propose, by analogy with rock physics, that a rate- and state-dependent friction law allows us to incorporate memory (and thus the effects of lubrication and bonding) into ice friction models without a great increase in complexity. We support this proposal with experimental data on both the laboratory (∼0.1 m) and ice tank (∼1 m) scale. These experiments show that the effects of static contact under normal load can be incorporated into a friction model. We find the parameters for a first-order rate and state model to be A = 0.310, B = 0.382, and μ0 = 0.872. Such a model then allows us to make predictions about the nature of memory effects in moving ice-ice contacts.

Method of determining a geophysical-scale sea ice rheology from laboratory measurements

We present a methodology that allows a sea ice rheology, suitable for use in a General Circulation Model (GCM), to be determined from laboratory and tank experiments on sea ice when combined with a kinematic model of deformation. The laboratory experiments determine a material rheology for sea ice, and would investigate a nonlinear friction law of the form τ ∝ σ n⅔, instead of the more familiar Amonton's law, τ = μσn (τ is the shear stress, μ is the coefficient of friction and σ n is the normal stress). The modelling approach considers a representative region R containing ice floes (or floe aggregates), separated by flaws. The deformation of R is imposed and the motion of the floes determined using a kinematic model, which will be motivated from SAR observations. Deformation of the flaws is inferred from the floe motion and stress determined from the material rheology. The stress over R is then determined from the area-weighted contribution from flaws and floes