Fitting mixtures of Kent distributions to aid in joint set identification

When examining a rock mass, joint sets and their orientations can play a significant role with regard to how the rock mass will behave. To identify joint sets present in the rock mass, the orientation of individual fracture planer can be measured on exposed rock faces and the resulting data can be examined for heterogeneity. In this article, the expectation-maximization algorithm is used to lit mixtures of Kent component distributions to the fracture data to aid in the identification of joint sets. An additional uniform component is also included in the model to accommodate the noise present in the data.

A systematic approach to error isolation in computerized wastewater simulation models

Activated sludge models are used extensively in the study of wastewater treatment processes. While various commercial implementations of these models are available, there are many people who need to code models themselves using the simulation packages available to them, Quality assurance of such models is difficult. While benchmarking problems have been developed and are available, the comparison of simulation data with that of commercial models leads only to the detection, not the isolation of errors. To identify the errors in the code is time-consuming. In this paper, we address the problem by developing a systematic and largely automated approach to the isolation of coding errors. There are three steps: firstly, possible errors are classified according to their place in the model structure and a feature matrix is established for each class of errors. Secondly, an observer is designed to generate residuals, such that each class of errors imposes a subspace, spanned by its feature matrix, on the residuals. Finally. localising the residuals in a subspace isolates coding errors. The algorithm proved capable of rapidly and reliably isolating a variety of single and simultaneous errors in a case study using the ASM 1 activated sludge model. In this paper a newly coded model was verified against a known implementation. The method is also applicable to simultaneous verification of any two independent implementations, hence is useful in commercial model development.

Surface diffusion and adsorption of hydrocarbons in activated carbon

Adsorption and diffusion in a porous media were studied theoretically and experimentally with a differential transient permeation method. The porous medium is allowed to equilibrate at some specified loading, and then the time trajectory of the permeation process is followed after a small difference between the pressures at the end faces of the porous medium is introduced at time t = 0 +. Such a trajectory us. time would contain adsorption and diffusion characteristics of the system. By studying this for various surface loadings, pore and surface diffusions can be fully characterized. Mathematical modeling of transient permeation is detailed for pure gases or vapors diffusion and adsorption in porous media. Effects of nonlinearity of adsorption isotherm, pressure, temperature and heat effects were considered in the model. Experimental data of diffusion and adsorption of propane, n-butane and n-hexane in activated carbon at different temperatures and loadings show the potential of this method as a useful tool to study adsorption kinetics in porous media. Validity of the model is best tested against the transient data where the kinetics curves exhibit sigmoidal shape, which is a result of the diffusion and adsorption rate during the initial stage of permeation.

A typology of community opportunity and vulnerability in metropolitan Australia

A multivariate model using hierarchical clustering and discriminant analysis is used to identify clusters of community opportunity and community vulnerability across Australia's mega metropolitan regions, Variables used in the model measure aspects of structural economic change, occupational change, human capital, income, unemployment, family/household disadvantage, and housing stress. A nine-cluster solution is used to categorise communities across metropolitan space. Significant between-city variations in the incidence of these clusters of opportunity and vulnerability are apparent, suggesting the emergence of marked differentiation between Australia's mega metropolitan regions in their adjustments to changing economic and social conditions. JEL classification: C49, R11, R12.

Reversal of cardiac and renal fibrosis by pirfenidone and spironolactone in streptozotocin-diabetic rats

1 Fibrosis leads to chronic impairment of cardiac and renal function and thus reversal of existing fibrosis may improve function and survival. This project has determined whether pirfenidone, a new antifibrotic compound, and spironolactone, an aldosterone antagonist, reverse both deposition of the major extracellular matrix proteins, collagen and fibronectin, and functional changes in the streptozotocin(STZ)-diabetic rat. 2 Streptozotocin (65 mg kg(-1) i.v.)-treated rats given pirfenidone (5-methyl-1-phenyl-2-[1H]-pyridone; approximately-200 mg kg(-1) day(-1) as 0.2-2g l(-1) drinking water) or spironolactone (50 mg kg(-1) day(-1) s.c.) for 4 weeks starting 4 weeks after STZ showed no attenuation of the increased blood glucose concentrations and increased food and water intakes which characterize diabetes in this model. 3 STZ-treatment increased perivascular and interstitial collagen deposition in the left ventricle and kidney, and surrounding the aorta. Cardiac, renal and plasma fibronectin concentrations increased in STZ-diabetic rats. Passive diastolic stiffness increased in isolated hearts from STZ-diabetic rats. Both pirfenidone and spironolactone treatment attenuated these increases without normalizing the decreased + dP/dt(max) of STZ-diabetic hearts. 4 Left ventricular papillary muscles from STZ-treated rats showed decreased maximal positive inotropic responses to noradrenaline, EMD 57033 (calcium sensitizer) and calcium chloride; this was not reversed by pirfenidone or spironolactone treatment. STZ-treatment transiently decreased GFR and urine flow rates in isolated perfused kidneys; pirfenidone but not spironolactone prevented the return to control values. 5 Thus, short-term pirfenidone and spironolactone treatment reversed cardiac and renal fibrosis and attenuated the increased diastolic stiffness without normalizing cardiac contractility or renal function in STZ-diabetic rats.

Proposed modifications to metabolic model for glycogen-accumulating organisms under anaerobic conditions

Filipe et al. (2001) proposed an anaerobic metabolic model for glycogen-accumulating organisms (GAO) in which the succinate-propionate pathway was used to describe the production of propionyl-CoA. However, propionyl-CoA is only an intermediate product in the above pathway. Stopping at propionyl-CoA instead of propionate (the end product of the pathway) results in the consumption of one ATP from succinate to succinyl-CoA, which was not accounted for in the model of Filipe et al. (2001). This resulted in significant errors in the stoichiometric coefficients in the final metabolic model. A modified model is presented in this communication and is shown to fit the experimental data significantly better than the original model. (C) 2002 Wiley Periodicals, Inc.

The Mallory body as an aggresome: In vitro studies

Prior in vivo studies supported the concept that Mallory bodies (MBs) are aggresomes of cytokeratins 8 and 18. However, to test this hypothesis an in vitro model is needed to study the dynamics of MB formation. Such a study is difficult because MBs have never been induced in tissue culture. Therefore, MBs were first induced in vivo in drug-primed mice and then primary cultures of hepatocytes from these mice were studied. Two approaches were utilized: 1. Primary cultures were transfected with plasmids containing the sequence for cytokeratin 18 (CK 18) tagged with green fluorescent protein (GFP). 2. Immunofluorescent staining was used to localize the ubiquitin-proteasome pathway components involved in MB-aggresome complex formation in primary hepatocyte cultures. The cells were double stained with a ubiquitin antibody and one of the following antibodies: CK 8, CK 18, tubulin, mutant ubiquitin (UBB+ 1), transglutaminase, phosphothreonine, and the 20S and 26S proteasome subunits P25 and Tbp7, respectively. In the first approach, fluorescence was observed in keratin filaments and MBs 48 h after the cells were transfected with the CK 18 GFP plasmid. Nascent cytokeratin 18 was preferentially concentrated in MBs. Less fluorescence was observed in the normal keratin filaments. This indicated that MBs continued to form in vitro. The immunofluorescent staining of the hepatocytes showed that CK 8 and 18, ubiquitin, mutant ubiquitin (UBB+ 1), P25, Tbp7, phosphothreonine, tubulin, and transglutaminase were all located at the border or the interior of the MB. These results support the concept that MBs are aggresomes of CK 8 and CK 18 and are a result of inhibition of the ubiquitin-proteasome pathway of protein degradation possibly caused by UBB+ 1. (C) 2002 Elsevier Science.

Spatial structure and habitat variation in a grasshopper hybrid zone

A hybrid zone between the grasshoppers Chorthippus brunneus and C. jacobsi (Orthoptera: Acrididae) in northern Spain has been analyzed for variation in morphology and ecology. These species are readily distinguished by the number of stridulatory pegs on the hind femur. Both sexes are fully winged and inhabit disturbed habitats throughout the study area. We develop a maximum-likelihood approach to fitting a two-dimensional cline to geographical variation in quantitative traits and for estimating associations of population mean with local habitat. This method reveals a cline in peg number approximately 30 km south of the Picos de Europa Mountains that shows substantial deviations in population mean compared with the expectations of simple tension zone models. The inclusion of variation in local vegetation in the model explains a significant proportion of the residual variation in peg number, indicating that habitat-genotype associations contribute to the observed spatial pattern. However, this association is weak, and a number of populations continue to show strong deviations in mean even after habitat is included in the final model. These outliers may be the result of long-distance colonization of sites distant from the cline center or may be due to a patchy pattern of initial contact during postglacial expansion. As well as contrasting with the smooth hybrid zones described for Chorthippus parallelus, this situation also contrasts with the mosaic hybrid zones observed in Gryllus crickets and in parts of the hybrid zone between Bombina toad species, where habitat-genotype associations account for substantial amounts of among-site variation.

Shear stress and sediment transport calculations for sheet flow under waves

A simple method is provided for calculating transport rates of not too fine (d(50) greater than or equal to 0.20 mm) sand under sheet flow conditions. The method consists of a Meyer-Peter-type transport formula operating on a time-varying Shields parameter, which accounts for both acceleration-asymmetry and boundary layer streaming. While velocity moment formulae, e.g.., = Constant x calibrated against U-tube measurements, fail spectacularly under some real waves (Ribberink, J.S., Dohmen-Janssen, C.M., Hanes, D.M., McLean, S.R., Vincent, C., 2000. Near-bed sand transport mechanisms under waves. Proc. 27th Int. Conf. Coastal Engineering, Sydney, ASCE, New York, pp. 3263-3276, Fig. 12), the new method predicts the real wave observations equally well. The reason that the velocity moment formulae fail under these waves is partly the presence of boundary layer streaming and partly the saw-tooth asymmetry, i.e., the front of the waves being steeper than the back. Waves with saw-tooth asymmetry may generate a net landward sediment transport even if = 0, because of the more abrupt acceleration under the steep front. More abrupt accelerations are associated with thinner boundary layers and greater pressure gradients for a given velocity magnitude. The two real wave effects are incorporated in a model of the form Q(s)(t) = Q(s)[theta(t)] rather than Q(S)(t) = Q(S)[u(infinity)(t)], i.e., by expressing the transport rate in terms of an instantaneous Shields parameter rather than in terms of the free stream velocity, and accounting for both streaming and accelerations in the 0(t) calculations. The instantaneous friction velocities u(*)(t) and subsequently theta(t) are calculated as follows. Firstly, a linear filter incorporating the grain roughness friction factor f(2.5) and a phase angle phi(tau) is applied to u(infinity)(t). This delivers u(*)(t) which is used to calculate an instantaneous grain roughness Shields parameter theta(2.5)(t). Secondly, a constant bed shear stress is added which corresponds to the streaming related bed shear stress -rho ($) over bar((u) over tilde(w) over tilde)(infinity) . The method can be applied to any u(infinity)(t) time series, but further experimental validation is recommended before application to conditions that differ strongly from the ones considered below. The method is not recommended for rippled beds or for sheet flow with typical prototype wave periods and d(50) < 0.20 turn. In such scenarios, time lags related to vertical sediment movement become important, and these are not considered by the present model. (C) 2002 Elsevier Science B.V. All rights reserved.

Age does not influence the bone response to treadmill exercise in female rats

Purpose: Because it is believed that bone may respond to exercise differently at different ages, we compared bone responses in immature and mature rats after 12 wk of treadmill running. Methods: Twenty-two immature (5-wk-old) and 21 mature (17-wk-old) female Sprague Dawley rats were randomized into a running (trained, N = 10 immature, 9 mature) or a control group (controls, N 12 immature, 12 mature) before sacrifice 12 wk later. Rats ran on a treadmill five times per week for 60-70 min at speeds up to 26 m.min(-1). Both at baseline and after intervention, we measured total body, lumbar spine, and proximal femoral bone mineral, as well as total body soft tissue composition using dual-energy x-ray absorptiometry (DXA) in vivo. After sacrificing the animals, we measured dynamic and static histomorphometry and three-point bending strength of the tibia. Results: Running training was associated with greater differences in tibial subperiosteal area, cortical cross-sectional area, peak load, stiffness, and moment of inertia in immature and mature rats (P < 0.05). The trained rats had greater periosteal bone formation rates (P < 0.01) than controls, but there was no difference in tibial trabecular bone histomorphometry. Similar running-related gains were seen in DXA lumbar spine area (P = 0.04) and bone mineral content (BMC; P = 0.03) at both ages. For total body bone area and BMC, the immature trained group increased significantly compared with controls (P < 0.05), whereas the mature trained group gained less than did controls (P < 0.01). Conclusion: In this in vivo model, where a similar physical training program was performed by immature and mature female rats, we demonstrated that both age groups were sensitive to loading and that bone strength gains appeared to result more from changes in bone geometry than from improved material properties.

Application of the maximum entropy method to the (2 + 1)D four-fermion model

We investigate spectral functions extracted using the maximum entropy method from correlators measured in lattice simulations of the (2+1)-dimensional four-fermion model. This model is particularly interesting because it has both a chirally broken phase with a rich spectrum of mesonic bound states and a symmetric phase where there are only resonances. In the broken phase we study the elementary fermion, pion, sigma, and massive pseudoscalar meson; our results confirm the Goldstone nature of the π and permit an estimate of the meson binding energy. We have, however, seen no signal of σ→ππ decay as the chiral limit is approached. In the symmetric phase we observe a resonance of nonzero width in qualitative agreement with analytic expectations; in addition the ultraviolet behavior of the spectral functions is consistent with the large nonperturbative anomalous dimension for fermion composite operators expected in this model.

Electrolyte transport in the mouse trachea: No evidence for a contribution of luminal K+ conductance

Recent studies on frog skin acini have challenged the question whether Cl- secretion or Na+ absorption in the airways is driven by luminal K+ channels in series to a basolateral K+ conductance. We examined the possible role of luminal K+ channels in electrolyte transport in mouse trachea in Ussing-chamber experiments. Tracheas of both normal and CFTR (-/-) mice showed a dominant amiloride-sensitive Na+ absorption under both, control conditions and after cAMP-dependent stimulation. The lumen-negative transepithelial voltage was enhanced after application of IBMX and forskolin and Cl- secretion was activated. Electrolyte secretion induced by IBMX and forskolin was inhibited by luminal glibenclamide and the blocker of basolateral Na(+)2Cl(-)K(+) cotransporter azosemide. Similarly, the compound 29313, a blocker of basolateral KCNQ1/KCNE3 K+ channels effectively blocked Cl- secretion when applied to either the luminal or basolateral side of the epithelium. RT-PCR analysis suggested expression of additional K+ channels in tracheal epithelial cells such as Slo1 and Kir6.2. However, we did not detect any functional evidence for expression of luminal K+ channels in mouse airways, using luminal 29313, clotrimazole and Ba2+ or different K+ channel toxins such as charybdotoxin, apamin and alpha-dendrotoxin. Thus, the present study demonstrates Cl- secretion in mouse airways, which depends on basolateral Na(+)2Cl(-)K(+) cotransport and luminal CFTR and non-CFTR Cl- channels. Cl- secretion is maintained by the activity of basolateral K+ channels, while no clear evidence was found for the presence of a luminal K+ conductance.

In vivo and in vitro models demonstrate a role for caveolin-1 in the pathogenesis of ischaemic acute renal failure

Caveolae and their proteins, the caveolins, transport macromolecules; compartmentalize signalling molecules; and are involved in various repair processes. There is little information regarding their role in the pathogenesis of significant renal syndromes such as acute renal failure (ARF). In this study, an in vivo rat model of 30 min bilateral renal ischaemia followed by reperfusion times from 4 h to 1 week was used to map the temporal and spatial association between caveolin-1 and tubular epithelial damage (desquamation, apoptosis, necrosis). An in vitro model of ischaemic ARF was also studied, where cultured renal tubular epithelial cells or arterial endothelial cells were subjected to injury initiators modelled on ischaemia-reperfusion (hypoxia, serum deprivation, free radical damage or hypoxia-hyperoxia). Expression of caveolin proteins was investigated using immunohistochemistry, immunoelectron microscopy, and immunoblots of whole cell, membrane or cytosol protein extracts. In vivo, healthy kidney had abundant caveolin-1 in vascular endothelial cells and also some expression in membrane surfaces of distal tubular epithelium. In the kidneys of ARF animals, punctate cytoplasmic localization of caveolin-1 was identified, with high intensity expression in injured proximal tubules that were losing basement membrane adhesion or were apoptotic, 24 h to 4 days after ischaemia-reperfusion. Western immunoblots indicated a marked increase in caveolin-1 expression in the cortex where some proximal tubular injury was located. In vitro, the main treatment-induced change in both cell types was translocation of caveolin-1 from the original plasma membrane site into membrane-associated sites in the cytoplasm. Overall, expression levels did not alter for whole cell extracts and the protein remained membrane-bound, as indicated by cell fractionation analyses. Caveolin-1 was also found to localize intensely within apoptotic cells. The results are indicative of a role for caveolin-1 in ARF-induced renal injury. Whether it functions for cell repair or death remains to be elucidated. Copyright (C) 2003 John Wiley Sons, Ltd.

A unified model of flames as gasdynamic discontinuities

Viewed on a hydrodynamic scale, flames in experiments are often thin so that they may be described as gasdynamic discontinuities separating the dense cold fresh mixture from the light hot burned products. The original model of a flame as a gasdynamic discontinuity was due to Darrieus and to Landau. In addition to the fluid dynamical equations, the model consists of a flame speed relation describing the evolution of the discontinuity surface, and jump conditions across the surface which relate the fluid variables on the two sides of the surface. The Darrieus-Landau model predicts, in contrast to observations, that a uniformly propagating planar flame is absolutely unstable and that the strength of the instability grows with increasing perturbation wavenumber so that there is no high-wavenumber cutoff of the instability. The model was modified by Markstein to exhibit a high-wavenumber cutoff if a phenomenological constant in the model has an appropriate sign. Both models are postulated, rather than derived from first principles, and both ignore the flame structure, which depends on chemical kinetics and transport processes within the flame. At present, there are two models which have been derived, rather than postulated, and which are valid in two non-overlapping regions of parameter space. Sivashinsky derived a generalization of the Darrieus-Landau model which is valid for Lewis numbers (ratio of thermal diffusivity to mass diffusivity of the deficient reaction component) bounded away from unity. Matalon & Matkowsky derived a model valid for Lewis numbers close to unity. Each model has its own advantages and disadvantages. Under appropriate conditions the Matalon-Matkowsky model exhibits a high-wavenumber cutoff of the Darrieus-Landau instability. However, since the Lewis numbers considered lie too close to unity, the Matalon-Matkowsky model does not capture the pulsating instability. The Sivashinsky model does capture the pulsating instability, but does not exhibit its high-wavenumber cutoff. In this paper, we derive a model consisting of a new flame speed relation and new jump conditions, which is valid for arbitrary Lewis numbers. It captures the pulsating instability and exhibits the high-wavenumber cutoff of all instabilities. The flame speed relation includes the effect of short wavelengths, not previously considered, which leads to stabilizing transverse surface diffusion terms.

Constitutively active signal transducer and activator of transcription 5 can replace the requirement for growth hormone in adipogenesis of 3T3-F442A preadipocytes

Although it is the best characterized in vitro model of GH action, the mechanisms used by GH to induce differentiation of murine 3T3-F442A preadipocytes remain unclear. Here we have examined the role of three transcriptional regulators in adipogenesis. These regulators are either rapidly induced in response to GH [Stra13, signal transducer and activator of transcription (Stat) 3] or of central importance to GH signaling (Stat5). Retroviral transfection of 3T3-F442A preadipocytes was used to increase expression of Stra13, Stat3, and Stat5a. Only Stat5a transfection increased the expression of adipogenic markers peroxisome proliferator-activated receptor gamma, CCAAT enhancer binding protein (C/EBP)alpha, and adipose protein 2/fatty acid-binding protein in response to GH, as determined by quantitative RT-PCR. Transfection with constitutively active Stat3 and Stat5a revealed that constitutively active Stat5a but not Stat3 was able to replace the GH requirement for adipogenesis. Constitutively active Stat5a but not Stat3 was able to increase the formation of lipid droplets and expression of alpha-glycerol phosphate dehydrogenase toward levels seen in mature adipocytes. Constitutively active Stat5a was also able to increase the expression of transcripts for C/EBPalpha to similar levels as GH, and of C/EBPbeta, peroxisome proliferator-activated receptor gamma, and adipose protein 2/fatty acid-binding protein transcripts to a lesser extent. An in vivo role for GH in murine adipogenesis is supported by significantly decreased epididymal fat depot size in young GH receptor-deleted mice, before manifestation of the lipolytic actions of GH. We conclude that Stat5 is a critical factor in GH-induced, and potentially prolactin-induced, murine adipogenesis.