The vertebral venous system in healthy and scoliotic adolescent spines - a 3D MRI investigation

Introduction. The venous drainage system within vertebral bodies (VBs) has been well documented previously in cadaveric specimens. Advances in 3D imaging and image processing now allow for in vivo quantification of larger venous vessels, such as the basivertebral vein. Differences between healthy and scoliotic VB veins can therefore be investigated. Methods. 20 healthy adolescent controls and 21 AIS patients were recruited (with ethics approval) to undergo 3D MRI, using a 3 Tesla, T1-weighted 3D gradient echo sequence, resulting in 512 slices across the thoraco-lumbar spine, with a voxel size of 0.5x0.5x0.5mm. Using Amira Filament Editor, five transverse slices through the VB were examined simultaneously and the resulting observable vascular network traced. Each VB was assessed, and a vascular network recorded when observable. A local coordinate system was created in the centre of each VB and the vascular networks aligned to this. The length of the vascular network on the left and right sides (with a small central region) of the VB was calculated, and the spatial patterning of the networks assessed level-by-level within each subject. Results. An average of 6 (range 4-10) vascular networks, consistent with descriptions of the basivertebral vein, were identifiable within each subject, most commonly between T10-L1. Differences were seen in the left/right distribution of vessels in the control and AIS subjects. Healthy controls saw a percentage distribution of 29:18:53 across the left:centre:right regions respectively, whereas the AIS subjects had a slightly shifted distribution of 33:25:42. The control group showed consistent spatial patterning of the vascular networks across most levels, but this was not seen in the AIS group. Conclusion. Observation and quantification of the basivertebral vein in vivo is possible using 3D MRI. The AIS group lacked the spatial pattern repetition seen in the control group and minor differences were seen in the left/right distribution of vessels.

Trans-disciplinary research in synthesis of grass pollen aerobiology and its importance for respiratory health in Australasia

Grass pollen is a major trigger for allergic rhinitis and asthma, yet little is known about the timing and levels of human exposure to airborne grass pollen across Australasian urban environments. The relationships between environmental aeroallergen exposure and allergic respiratory disease bridge the fields of ecology, aerobiology, geospatial science and public health. The Australian Aerobiology Working Group comprised of experts in botany, palynology, biogeography, climate change science, plant genetics, biostatistics, ecology, pollen allergy, public and environmental health, and medicine, was established to systematically source, collate and analyse atmospheric pollen concentration data from 11 Australian and six New Zealand sites. Following two week-long workshops, post-workshop evaluations were conducted to reflect upon the utility of this analysis and synthesis approach to address complex multidisciplinary questions. This Working Group described i) a biogeographically dependent variation in airborne pollen diversity, ii) a latitudinal gradient in the timing, duration and number of peaks of the grass pollen season, and iii) the emergence of new methodologies based on trans-disciplinary synthesis of aerobiology and remote sensing data. Challenges included resolving methodological variations between pollen monitoring sites and temporal variations in pollen datasets. Other challenges included “marrying” ecosystem and health sciences and reconciling divergent expert opinion. The Australian Aerobiology Working Group facilitated knowledge transfer between diverse scientific disciplines, mentored students and early career scientists, and provided an uninterrupted collaborative opportunity to focus on a unifying problem globally. The Working Group provided a platform to optimise the value of large existing ecological datasets that have importance for human respiratory health and ecosystems research. Compilation of current knowledge of Australasian pollen aerobiology is a critical first step towards the management of exposure to pollen in patients with allergic disease and provides a basis from which the future impacts of climate change on pollen distribution can be assessed and monitored.

Identification of bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei reveals distinct temperature-driven patterns

Change in temperature is often a major environmental factor in triggering waterborne disease outbreaks. Previous research has revealed temporal and spatial patterns of bacterial population in several aquatic ecosystems. To date, very little information is available on aquaculture environment. Here, we assessed environmental temperature effects on bacterial community composition in freshwater aquaculture system farming of Litopenaeus vannamei (FASFL). Water samples were collected over a one-year period, and aquatic bacteria were characterized by polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE) and 16S rDNA pyrosequencing. Resulting DGGE fingerprints revealed a specific and dynamic bacterial population structure with considerable variation over the seasonal change, suggesting that environmental temperature was a key driver of bacterial population in the FASFL. Pyrosequencing data further demonstrated substantial difference in bacterial community composition between the water at higher (WHT) and at lower (WLT) temperatures in the FASFL. Actinobacteria, Proteobacteria and Bacteroidetes were the highest abundant phyla in the FASFL, however, a large number of unclassified bacteria contributed the most to the observed variation in phylogenetic diversity. The WHT harbored remarkably higher diversity and richness in bacterial composition at genus and species levels when compared to the WLT. Some potential pathogenenic species were identified in both WHT and WLT, providing data in support of aquatic animal health management in the aquaculture industry.

Postglacial Early Permian (late Sakmarian-early Artinskian) shallow-marine carbonate deposition along a 2000 km transect from Timor to west Australia

Late Sakmarian to early Artinskian (Early Permian) carbonate deposition was widespread in the marine intracratonic rift basins that extended into the interior of Eastern Gondwana from Timor in the north to the northern Perth Basin in the south. These basins spanned about 20° of paleolatitude (approximately 35°S to 55°S). This study describes the type section of the Maubisse Limestone in Timor-Leste, and compares this unit with carbonate sections in the Canning Basin (Nura Nura Member of the Poole Sandstone), the Southern Carnarvon Basin (Callytharra Formation) and the northern Perth Basin (Fossil Cliff Member of the Holmwood Shale). The carbonate units have no glacial influence and formed part of a major depositional cycle that, in the southern basins, overlies glacially influenced strata and lies a short distance below mudstone containing marine fossils and scattered dropstones (perhaps indicative of sea ice). In the south marine conditions became more restricted and were replaced by coal measures at the top of the depositional sequence. In the north, the carbonate deposits are possibly bryozoan–crinoidal mounds; whereas in the southern basins they form laterally continuous relatively thin beds, deposited on a very low-gradient seafloor, at the tops of shale–limestone parasequences that thicken upward in parasequence sets. All marine deposition within the sequence took place under very shallow (inner neritic) conditions, and the limestones have similar grain composition. Bryozoan and crinoidal debris dominate the grain assemblages and brachiopod shell fragments, foraminifera and ostracod valves are usually common. Tubiphytes ranged as far south as the Southern Carnarvon Basin, albeit rarely, but is more common to the north. Gastropod and bivalve shell debris, echinoid spines, solitary rugose corals and trilobite carapace elements are rare. The uniformity of the grain assemblage and the lack of tropical elements such as larger fusulinid foraminifera, colonial corals or dasycladacean algae indicate temperate marine conditions with only a small increase in temperature to the north. The depositional cycle containing the studied carbonate deposits represents a warmer phase than the preceding glacially influenced Asselian to early Sakmarian interval and the subsequent cool phase of the “mid” Artinskian that is followed by significant warming during the late Artinskian–early Kungurian. The timing of cooler and warmer intervals in the west Australian basins seems out-of-phase with the eastern Australian succession, but this may be a problem of chronostratigraphic miscorrelation due to endemic faunas and palynofloras.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

A framework for evaluating the impact of obesity prevention strategies on socioeconomic inequalities in weight

We developed a theoretical framework to organize obesity prevention interventions by their likely impact on the socioeconomic gradient of weight. The degree to which an intervention involves individual agency versus structural change influences socioeconomic inequalities in weight. Agentic interventions, such as standalone social marketing, increase socioeconomic inequalities. Structural interventions, such as food procurement policies and restrictions on unhealthy foods in schools, show equal or greater benefit for lower socioeconomic groups. Many obesity prevention interventions belong to the agento–structural types of interventions, and account for the environment in which health behaviors occur, but they require a level of individual agency for behavioral change, including workplace design to encourage exercise and fiscal regulation of unhealthy foods or beverages. Obesity prevention interventions differ in their effectiveness across socioeconomic groups. Limiting further increases in socioeconomic inequalities in obesity requires implementation of structural interventions. Further empirical evaluation, especially of agento–structural type interventions, remains crucial.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

In situ reaction furnace for real-time XRD studies

The new furnace at the Materials Characterization by X-ray Diffraction beamline at Elettra has been designed for powder diffraction measurements at high temperature (up to 1373 K at the present state). Around the measurement region the geometry of the radiative heating element assures a negligible temperature gradient along the capillary and can accommodate either powder samples in capillary or small flat samples. A double capillary holder allows flow-through of gas in the inner sample capillary while the outer one serves as the reaction chamber. The furnace is coupled to a translating curved imaging-plate detector, allowing the collection of diffraction patterns up to 2[theta] [asymptotically equal to] 130°.

In situ investigation of explosive crystallization in a-Ge: Experimental determination of the interface response function using dynamic transmission electron microscopy

The crystallization of amorphous semiconductors is a strongly exothermic process. Once initiated the release of latent heat can be sufficient to drive a self-sustaining crystallization front through the material in a manner that has been described as explosive. Here, we perform a quantitative in situ study of explosive crystallization in amorphous germanium using dynamic transmission electron microscopy. Direct observations of the speed of the explosive crystallization front as it evolves along a laser-imprinted temperature gradient are used to experimentally determine the complete interface response function (i.e., the temperature-dependent front propagation speed) for this process, which reaches a peak of 16 m/s. Fitting to the Frenkel-Wilson kinetic law demonstrates that the diffusivity of the material locally/immediately in advance of the explosive crystallization front is inconsistent with those of a liquid phase. This result suggests a modification to the liquid-mediated mechanism commonly used to describe this process that replaces the phase change at the leading amorphous-liquid interface with a change in bonding character (from covalent to metallic) occurring in the hot amorphous material.

Regolith mass balance inferred from combined mineralogical, geochemical and geophysical studies: Mule Hole gneissic watershed, South India

The aim of this study is to propose a method to assess the long-term chemical weathering mass balance for a regolith developed on a heterogeneous silicate substratum at the small experimental watershed scale by adopting a combined approach of geophysics, geochemistry and mineralogy. We initiated in 2003 a study of the steep climatic gradient and associated geomorphologic features of the edge of the rifted continental passive margin of the Karnataka Plateau, Peninsular India. In the transition sub-humid zone of this climatic gradient we have studied the pristine forested small watershed of Mule Hole (4.3 km(2)) mainly developed on gneissic substratum. Mineralogical, geochemical and geophysical investigations were carried out (i) in characteristic red soil profiles and (ii) in boreholes up to 60 m deep in order to take into account the effect of the weathering mantle roots. In addition, 12 Electrical Resistivity Tomography profiles (ERT), with an investigation depth of 30 m, were generated at the watershed scale to spatially characterize the information gathered in boreholes and soil profiles. The location of the ERT profiles is based on a previous electromagnetic survey, with an investigation depth of about 6 m. The soil cover thickness was inferred from the electromagnetic survey combined with a geological/pedological survey. Taking into account the parent rock heterogeneity, the degree of weathering of each of the regolith samples has been defined using both the mineralogical composition and the geochemical indices (Loss on Ignition, Weathering Index of Parker, Chemical Index of Alteration). Comparing these indices with electrical resistivity logs, it has been found that a value of 400 Ohm m delineates clearly the parent rocks and the weathered materials, Then the 12 inverted ERT profiles were constrained with this value after verifying the uncertainty due to the inversion procedure. Synthetic models based on the field data were used for this purpose. The estimated average regolith thickness at the watershed scale is 17.2 m, including 15.2 m of saprolite and 2 m of soil cover. Finally, using these estimations of the thicknesses, the long-term mass balance is calculated for the average gneiss-derived saprolite and red soil. In the saprolite, the open-system mass-transport function T indicates that all the major elements except Ca are depleted. The chlorite and biotite crystals, the chief sources for Mg (95%), Fe (84%), Mn (86%) and K (57%, biotite only), are the first to undergo weathering and the oligoclase crystals are relatively intact within the saprolite with a loss of only 18%. The Ca accumulation can be attributed to the precipitation of CaCO3 from the percolating solution due to the current and/or the paleoclimatic conditions. Overall, the most important losses occur for Si, Mg and Na with -286 x 10(6) mol/ha (62% of the total mass loss), -67 x 10(6) mol/ha (15% of the total mass loss) and -39 x 10(6) mol/ha (9% of the total mass loss), respectively. Al, Fe and K account for 7%, 4% and 3% of the total mass loss, respectively. In the red soil profiles, the open-system mass-transport functions point out that all major elements except Mn are depleted. Most of the oligoclase crystals have broken down with a loss of 90%. The most important losses occur for Si, Na and Mg with -55 x 10(6) mol/ha (47% of the total mass loss), -22 x 10(6) mol/ha (19% of the total mass loss) and -16 x 10(6) mol/ha (14% of the total mass loss), respectively. Ca, Al, K and Fe account for 8%, 6%, 4% and 2% of the total mass loss, respectively. Overall these findings confirm the immaturity of the saprolite at the watershed scale. The soil profiles are more evolved than saprolite but still contain primary minerals that can further undergo weathering and hence consume atmospheric CO2.

Finite Element Solution Of Surface-Tension Driven Flows In Laser Surface-Melting

Lasers are very efficient in heating localized regions and hence they find a wide application in surface treatment processes. The surface of a material can be selectively modified to give superior wear and corrosion resistance. In laser surface-melting and welding problems, the high temperature gradient prevailing in the free surface induces a surface-tension gradient which is the dominant driving force for convection (known as thermo-capillary or Marangoni convection). It has been reported that the surface-tension driven convection plays a dominant role in determining the melt pool shape. In most of the earlier works on laser-melting and related problems, the finite difference method (FDM) has been used to solve the Navier Stokes equations [1]. Since the Reynolds number is quite high in these cases, upwinding has been used. Though upwinding gives physically realistic solutions even on a coarse grid, the results are inaccurate. McLay and Carey have solved the thermo-capillary flow in welding problems by an implicit finite element method [2]. They used the conventional Galerkin finite element method (FEM) which requires that the pressure be interpolated by one order lower than velocity (mixed interpolation). This restricts the choice of elements to certain higher order elements which need numerical integration for evaluation of element matrices. The implicit algorithm yields a system of nonlinear, unsymmetric equations which are not positive definite. Computations would be possible only with large mainframe computers.Sluzalec [3] has modeled the pulsed laser-melting problem by an explicit method (FEM). He has used the six-node triangular element with mixed interpolation. Since he has considered the buoyancy induced flow only, the velocity values are small. In the present work, an equal order explicit FEM is used to compute the thermo-capillary flow in the laser surface-melting problem. As this method permits equal order interpolation, there is no restriction in the choice of elements. Even linear elements such as the three-node triangular elements can be used. As the governing equations are solved in a sequential manner, the computer memory requirement is less. The finite element formulation is discussed in this paper along with typical numerical results.

Model Predictive Static Programming: A Computationally Efficient Technique For Suboptimal Control Design

Combining the philosophies of nonlinear model predictive control and approximate dynamic programming, a new suboptimal control design technique is presented in this paper, named as model predictive static programming (MPSP), which is applicable for finite-horizon nonlinear problems with terminal constraints. This technique is computationally efficient, and hence, can possibly be implemented online. The effectiveness of the proposed method is demonstrated by designing an ascent phase guidance scheme for a ballistic missile propelled by solid motors. A comparison study with a conventional gradient method shows that the MPSP solution is quite close to the optimal solution.

A probabilistic constrained nonlinear optimization framework to optimize RED parameters

The random early detection (RED) technique has seen a lot of research over the years. However, the functional relationship between RED performance and its parameters viz,, queue weight (omega(q)), marking probability (max(p)), minimum threshold (min(th)) and maximum threshold (max(th)) is not analytically availa ble. In this paper, we formulate a probabilistic constrained optimization problem by assuming a nonlinear relationship between the RED average queue length and its parameters. This problem involves all the RED parameters as the variables of the optimization problem. We use the barrier and the penalty function approaches for its Solution. However (as above), the exact functional relationship between the barrier and penalty objective functions and the optimization variable is not known, but noisy samples of these are available for different parameter values. Thus, for obtaining the gradient and Hessian of the objective, we use certain recently developed simultaneous perturbation stochastic approximation (SPSA) based estimates of these. We propose two four-timescale stochastic approximation algorithms based oil certain modified second-order SPSA updates for finding the optimum RED parameters. We present the results of detailed simulation experiments conducted over different network topologies and network/traffic conditions/settings, comparing the performance of Our algorithms with variants of RED and a few other well known adaptive queue management (AQM) techniques discussed in the literature.

Vertical structure and horizontal gradients of aerosol extinction coefficients over coastal India inferred from airborne lidar measurements during the Integrated Campaign for Aerosol, Gases and Radiation Budget (ICARB) field campaign

Quantitative estimates of the vertical structure and the spatial gradients of aerosol extinction coefficients have been made from airborne lidar measurements across the coastline into offshore oceanic regions along the east and west coasts of India. The vertical structure revealed the presence of strong, elevated aerosol layers in the altitude region of similar to 2-4 km, well above the atmospheric boundary layer (ABL). Horizontal gradients also showed a vertical structure, being sharp with the e(-1) scaling distance (D-0H) as small as similar to 150 km in the well-mixed regions mostly under the influence of local source effects. Above the ABL, where local effects are subdued, the gradients were much shallower (similar to 600-800 km); nevertheless, they were steep compared to the value of similar to 1500-2500 km reported for columnar AOD during winter. The gradients of these elevated layers were steeper over the east coast of India than over the west coast. Near-simultaneous radio sonde (Vaisala, Inc., Finland) ascents made over the northern Bay of Bengal showed the presence of convectively unstable regions, first from surface to similar to 750-1000 m and the other extending from 1750 to 3000 m separated by a stable region in between. These can act as a conduit for the advection of aerosols and favor the transport of continental aerosols in the higher levels (> 2 km) into the oceans without entering the marine boundary layer below. Large spatial gradient in aerosol optical and hence radiative impacts between the coastal landmass and the adjacent oceans within a short distance of < 300 km (even at an altitude of 3 km) during summer and the premonsoon is of significance to the regional climate.

Axially homogeneous, zero mean flow buoyancy-driven turbulence in a vertical pipe

We report an experimental study of a new type of turbulent flow that is driven purely by buoyancy. The flow is due to an unstable density difference, created using brine and water, across the ends of a long (length/diameter = 9) vertical pipe. The Schmidt number Sc is 670, and the Rayleigh number (Ra) based on the density gradient and diameter is about 10(8). Under these conditions the convection is turbulent, and the time-averaged velocity at any point is `zero'. The Reynolds number based on the Taylor microscale, Re-lambda, is about 65. The pipe is long enough for there to be an axially homogeneous region, with a linear density gradient, about 6-7 diameters long in the midlength of the pipe. In the absence of a mean flow and, therefore, mean shear, turbulence is sustained just by buoyancy. The flow can be thus considered to be an axially homogeneous turbulent natural convection driven by a constant (unstable) density gradient. We characterize the flow using flow visualization and particle image velocimetry (PIV). Measurements show that the mean velocities and the Reynolds shear stresses are zero across the cross-section; the root mean squared (r.m.s.) of the vertical velocity is larger than those of the lateral velocities (by about one and half times at the pipe axis). We identify some features of the turbulent flow using velocity correlation maps and the probability density functions of velocities and velocity differences. The flow away from the wall, affected mainly by buoyancy, consists of vertically moving fluid masses continually colliding and interacting, while the flow near the wall appears similar to that in wall-bound shear-free turbulence. The turbulence is anisotropic, with the anisotropy increasing to large values as the wall is approached. A mixing length model with the diameter of the pipe as the length scale predicts well the scalings for velocity fluctuations and the flux. This model implies that the Nusselt number would scale as (RaSc1/2)-Sc-1/2, and the Reynolds number would scale as (RaSc-1/2)-Sc-1/2. The velocity and the flux measurements appear to be consistent with the Ra-1/2 scaling, although it must be pointed out that the Rayleigh number range was less than 10. The Schmidt number was not varied to check the Sc scaling. The fluxes and the Reynolds numbers obtained in the present configuration are Much higher compared to what would be obtained in Rayleigh-Benard (R-B) convection for similar density differences.