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

Survival analysis statistics applied to threshold data obtained from the ascending forced-choice method of limits

Sensory thresholds are often collected through ascending forced-choice methods. Group thresholds are important for comparing stimuli or populations; yet, the method has two problems. An individual may correctly guess the correct answer at any concentration step and might detect correctly at low concentrations but become adapted or fatigued at higher concentrations. The survival analysis method deals with both issues. Individual sequences of incorrect and correct answers are adjusted, taking into account the group performance at each concentration. The technique reduces the chance probability where there are consecutive correct answers. Adjusted sequences are submitted to survival analysis to determine group thresholds. The technique was applied to an aroma threshold and a taste threshold study. It resulted in group thresholds similar to ASTM or logarithmic regression procedures. Significant differences in taste thresholds between younger and older adults were determined. The approach provides a more robust technique over previous estimation methods.

Folding of graphene slit like pore walls—a simple method of improving CO2 separation from mixtures with CH4 or N2

We report for the first time a detailed procedure for creating a simulation model of energetically stable, folded graphene-like pores and simulation results of CO2/CH4 and CO2/N2 separation using these structures. We show that folding of graphene structures is a very promising method to improve the separation of CO2 from mixtures with CH4 and N2. The separation properties of the analyzed materials are compared with carbon nanotubes having similar diameters or S/V ratio. The presented results have potential importance in the field of CO2 capture and sequestration.

A method of evaluating the accuracy of human body thermoregulation models

Human Body Thermoregulation Models have been widely used in the field of human physiology or thermal comfort studies. However there are few studies on the evaluation method for these models. This paper summarises the existing evaluation methods and critically analyses the flaws. Based on that, a method for the evaluating the accuracy of the Human Body Thermoregulation models is proposed. The new evaluation method contributes to the development of Human Body Thermoregulation models and validates their accuracy both statistically and empirically. The accuracy of different models can be compared by the new method. Furthermore, the new method is not only suitable for the evaluation of Human Body Thermoregulation Models, but also can be theoretically applied to the evaluation of the accuracy of the population-based models in other research fields.

A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

Does fertilizer type and method of application cause significant differences in essential oil yield and composition in rosemary (Rosmarinus officinalis L.)?

Organic fertilizers based on seaweed extract potentially have beneficial effects on many crop plants. Herewe investigate the impact of organic fertilizer on Rosmarinus officinalis measured by both yield and oilquality. Plants grown in a temperature-controlled greenhouse with a natural photoperiod and a controlledirrigation system were treated with seaweed fertilizer and an inorganic fertilizer of matching mineralcomposition but with no organic content. Treatments were either by spraying on to the foliage or wateringdirect to the compost. The essential oil was extracted by hydro-distillation with a Clevenger apparatusand analysed by gas-chromatography mass-spectrometry (GC–MS) and NMR. The chemical composi-tions of the plants were compared, and qualitative differences were found between fertilizer treatmentsand application methods. Thus sprayed seaweed fertilizer showed a significantly higher percentage of�-pinene, �-phellandrene, �-terpinene (monoterpenes) and 3-methylenecycloheptene than other treat-ments. Italicene, �-bisabolol (sesquiterpenes), �-thujene, and E-isocitral (monoterpenes) occurred insignificantly higher percentages for plants watered with the seaweed extract. Each was significantly dif-ferent to the inorganic fertilizer and to controls. The seaweed treatments caused a significant increasein oil amount and leaf area as compared with both inorganic treatments and the control regardless ofapplication method.

A novel method for production of lipid hydroperoxide- or oxysterol-rich low-density lipoprotein

LDL oxidation may be important in atherosclerosis. Extensive oxidation of LDL by copper induces increased uptake by macrophages, but results in decomposition of hydroperoxides, making it more difficult to investigate the effects of hydroperoxides in oxidised LDL on cell function. We describe here a simple method of oxidising LDL by dialysis against copper ions at 4 degrees C, which inhibits the decomposition of hydroperoxides, and allows the production of LDL rich in hydroperoxides (626 +/- 98 nmol/mg LDL protein) but low in oxysterols (3 +/- 1 nmol 7-ketocholesterol/mg LDL protein), whilst allowing sufficient modification (2.6 +/- 0.5 relative electrophoretic mobility) for rapid uptake by macrophages (5.49 +/- 0.75 mu g I-125-labelled hydroperoxide-rich LDL vs. 0.46 +/- 0.04 mu g protein/mg cell protein in 18 h for native LDL). By dialysing under the same conditions, but at 37 degrees C, the hydroperoxides are decomposed extensively and the LDL becomes rich in oxysterols. This novel method of oxidising LDL with high yield to either a hydroperoxide- or oxysterol-rich form by simply altering the temperature of dialysis may provide a useful tool for determining the effects of these different oxidation products on cell function. (C) 2007 Elsevier Ireland Ltd. All rights reserved.

A laboratory incubation method for determining the rate of microbiological degradation of skeletal muscle tissue in soil

A controlled laboratory experiment is described, in principle and practice, which can be used for the of determination the rate of tissue decomposition in soil. By way of example, an experiment was conducted to determine the effect of temperature (12°C, 22°C) on the aerobic decomposition of skeletal muscle tissue (Organic Texel × Suffolk lamb (Ovis aries)) in a sandy loam soil. Measurements of decomposition processes included muscle tissue mass loss, microbial CO2 respiration, and muscle tissue carbon (C) and nitrogen (N). Muscle tissue mass loss at 22°C always was greater than at 12°C (p < 0.001). Microbial respiration was greater in samples incubated at 22°C for the initial 21 days of burial (p < 0.01). All buried muscle tissue samples demonstrated changes in C and N content at the end of the experiment. A significant correlation (p < 0.001) was demonstrated between the loss of muscle tissue-derived C (C1) and microbially-respired C (Cm) demonstrating CO2 respiration may be used to predict mass loss and hence biodegradation. In this experiment Q10 (12°C - 22°C) = 2.0. This method is recommended as a useful tool in determining the effect of environmental variables on the rate of decomposition of various tissues and associated materials.

Adsorption, orientation and thermal decomposition of copper(II) hexafluoroacetylacetonate on rutile TiO2(110)

We have investigated the adsorption and thermal decomposition of copper hexafluoroacetylacetonate (Cu-11(hfaC)(2)) on single crystal rutile TiO2(110). Low energy electron diffraction shows that room temperature saturation coverage of the Cu-II(hfac)(2) adsorbate forms an ordered (2 x 1) over-layer. X-ray and ultra-violet photoemission spectroscopy of the saturated surface were recorded as the sample was annealed in a sequential manner to reveal decomposition pathways. The results show that the molecule dissociatively adsorbs by detachment of one of the two ligands to form hfac and Cu-1(hfac) which chemisorb to the substrate at 298 K. These ligands only begin to decompose once the surface temperature exceeds 473 K where Cu core level shifts indicate metallisation. This reduction from Cu(I) to Cu(0) takes place in the absence of an external reducing agent and without disproportionation and is accompanied by the onset of decomposition of the hfac ligands. Finally, C K-edge near edge X-ray absorption fine structure experiments indicate that both the ligands adsorb aligned in the < 001 > direction and we propose a model in which the hfac ligands adsorb on the 5-fold coordinated Ti atoms and the Cu-1(hfac) moiety attaches to the bridging O atoms in a square planar geometry. The calculated tilt angle for these combined geometries is approximately 10 degrees to the surface normal.

Calibration of uncertain flood inundation models using remotely sensed water levels

A traditional method of validating the performance of a flood model when remotely sensed data of the flood extent are available is to compare the predicted flood extent to that observed. The performance measure employed often uses areal pattern-matching to assess the degree to which the two extents overlap. Recently, remote sensing of flood extents using synthetic aperture radar (SAR) and airborne scanning laser altimetry (LIDAR) has made more straightforward the synoptic measurement of water surface elevations along flood waterlines, and this has emphasised the possibility of using alternative performance measures based on height. This paper considers the advantages that can accrue from using a performance measure based on waterline elevations rather than one based on areal patterns of wet and dry pixels. The two measures were compared for their ability to estimate flood inundation uncertainty maps from a set of model runs carried out to span the acceptable model parameter range in a GLUE-based analysis. A 1 in 5-year flood on the Thames in 1992 was used as a test event. As is typical for UK floods, only a single SAR image of observed flood extent was available for model calibration and validation. A simple implementation of a two-dimensional flood model (LISFLOOD-FP) was used to generate model flood extents for comparison with that observed. The performance measure based on height differences of corresponding points along the observed and modelled waterlines was found to be significantly more sensitive to the channel friction parameter than the measure based on areal patterns of flood extent. The former was able to restrict the parameter range of acceptable model runs and hence reduce the number of runs necessary to generate an inundation uncertainty map. A result of this was that there was less uncertainty in the final flood risk map. The uncertainty analysis included the effects of uncertainties in the observed flood extent as well as in model parameters. The height-based measure was found to be more sensitive when increased heighting accuracy was achieved by requiring that observed waterline heights varied slowly along the reach. The technique allows for the decomposition of the reach into sections, with different effective channel friction parameters used in different sections, which in this case resulted in lower r.m.s. height differences between observed and modelled waterlines than those achieved by runs using a single friction parameter for the whole reach. However, a validation of the modelled inundation uncertainty using the calibration event showed a significant difference between the uncertainty map and the observed flood extent. While this was true for both measures, the difference was especially significant for the height-based one. This is likely to be due to the conceptually simple flood inundation model and the coarse application resolution employed in this case. The increased sensitivity of the height-based measure may lead to an increased onus being placed on the model developer in the production of a valid model

Higher oxides of chlorine: absorption cross-sections of Cl2O6 and Cl2O4, the decomposition of Cl2O6, and the reactions of OClO with O and O-3

The absorption cross-sections of Cl2O6 and Cl2O4 have been obtained using a fast flow reactor with a diode array spectrometer (DAS) detection system. The absorption cross-sections at the wavelengths of maximum absorption (lambda(max)) determined in this study are those of Cl2O6: (1.47 +/- 0.15) x 10(-17) cm(2) molecule(-1), at lambda(max) = 276 nm and T = 298 K; and Cl2O4: (9.0 +/- 2.0) x 10(-19) cm(2) molecule(-1), at lambda(max) = 234 nm and T = 298 K. Errors quoted are two standard deviations together with estimates of the systematic error. The shapes of the absorption spectra were obtained over the wavelength range 200-450 nm for Cl2O6 and 200-350 nm for Cl2O4, and were normalized to the absolute cross-sections obtained at lambda(max) for each oxide, and are presented at 1 nm intervals. These data are discussed in relation to previous measurements. The reaction of O with OCIO has been investigated with the objective of observing transient spectroscopic absorptions. A transient absorption was seen, and the possibility is explored of identifying the species with the elusive sym-ClO3 or ClO4, both of which have been characterized in matrices, but not in the gas-phase. The photolysis of OCIO was also re-examined, with emphasis being placed on the products of reaction. UV absorptions attributable to one of the isomers of the ClO dimer, chloryl chloride (ClClO2) were observed; some Cl2O4 was also found at long photolysis times, when much of the ClClO2 had itself been photolysed. We suggest that reports of Cl2O6 formation in previous studies could be a consequence of a mistaken identification. At low temperatures, the photolysis of OCIO leads to the formation of Cl2O3 as a result of the addition of the ClO primary product to OCIO. ClClO2 also appears to be one product of the reaction between O-3 and OCIO, especially when the reaction occurs under explosive conditions. We studied the kinetics of the non-explosive process using a stopped-flow technique, and suggest a value for the room-temperature rate coefficient of (4.6 +/- 0.9) x 10(-19) cm(3) molecule(-1) s(-1) (limit quoted is 2sigma random errors). The photochemical and thermal decomposition of Cl2O6 is described in this paper. For photolysis at k = 254 nm, the removal of Cl2O6 is not accompanied by the build up of any other strong absorber. The implications of the results are either that the photolysis of Cl2O6 produces Cl-2 directly, or that the initial photofragments are converted rapidly to Cl-2. In the thermal decomposition of Cl2O6, Cl2O4 was shown to be a product of reaction, although not necessarily the major one. The kinetics of decomposition were investigated using the stopped-flow technique. At relatively high [OCIO] present in the system, the decay kinetics obeyed a first-order law, with a limiting first-order rate coefficient of 0.002 s(-1). (C) 2004 Elsevier B.V. All rights reserved.

Approximate Riemann solutions of the two-dimensional shallow-water equations

A finite-difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow-water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gasdynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. An extension to the two-dimensional equations with source terms, is included. The scheme is applied to a dam-break problem with cylindrical symmetry.

Solutions of a two-dimensional dam-break problem

Solutions of a two-dimensional dam break problem are presented for two tailwater/reservoir height ratios. The numerical scheme used is an extension of one previously given by the author [J. Hyd. Res. 26(3), 293–306 (1988)], and is based on numerical characteristic decomposition. Thus approximate solutions are obtained via linearised problems, and the method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations.

Prediction of supercritical flow in open channels

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow water equations in open channels. A linearised problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.