Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Entropy generation for forced convection in a porous saturated circular tube with uniform wall temperature

A numerical study is reported to investigate both the First and the Second Law of Thermodynamics for thermally developing forced convection in a circular tube filled by a saturated porous medium, with uniform wall temperature, and with the effects of viscous dissipation included. A theoretical analysis is also presented to study the problem for the asymptotic region applying the perturbation solution of the Brinkman momentum equation reported by Hooman and Kani [1]. Expressions are reported for the temperature profile, the Nusselt number, the Bejan number, and the dimensionless entropy generation rate in the asymptotic region. Numerical results are found to be in good agreement with theoretical counterparts.

Capillary network models for transport in packed beds: Considerations of pore aspect ratio

The conventional analysis for the estimation of the tortuosity factor for transport in porous media is modified here to account for the effect of pore aspect ratio. Structural models of the porous medium are also constructed for calculating the aspect ratio as a function of porosity. Comparison of the model predictions with the extensive data of Currie (1960) for the effective diffusivity of hydrogen in packed beds shows good agreement with a network model of randomly oriented intersecting pores for porosities upto about 50 percent, which is the region of practical interest. The predictions based on this network model are also found to be in better agreement with the data of Currie than earlier expressions developed for unconsolidated and grainy media.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Serum triglycerides as a risk factor for cardiovascular diseases in the Asia-Pacific Region

Background-The importance of serum triglyceride levels as a risk factor for cardiovascular diseases is uncertain. Methods and Results-We performed an individual participant data meta-analysis of prospective studies conducted in the Asia-Pacific region. Cox models were applied to the combined data from 26 studies to estimate the overall and region-, sex-, and age-specific hazard ratios for major cardiovascular diseases by fifths of triglyceride values. During 796 671 person-years of follow-up among 96 224 individuals, 670 and 667 deaths as a result of coronary heart disease (CHD) and stroke, respectively, were recorded. After adjustment for major cardiovascular risk factors, participants grouped in the highest fifth of triglyceride levels had a 70% (95% CI, 47 to 96) greater risk of CHD death, an 80% (95% CI, 49 to 119) higher risk of fatal or nonfatal CHD, and a 50% (95% CI, 29% to 76%) increased risk of fatal or nonfatal stroke compared with those belonging to the lowest fifth. The association between triglycerides and CHD death was similar across subgroups defined by ethnicity, age, and sex. Conclusions-Serum triglycerides are an important and independent predictor of CHD and stroke risk in the Asia-Pacific region. These results may have clinical implications for cardiovascular risk prediction and the use of lipid-lowering therapy.

Temporal and spatial distribution of anurans in the Pampa Region (Santa Maria, RS).

The present study aimed to determine the richness, occurrence constancy, reproductive modes. standard of abundance distribution, season of vocalization and to test correlation among climatic variables and activity of vocalization of anurans in a region of the Pampa Biome, Santa Maria, Rio Grande do Sul State. During the period of Novernber/2001 to October/2002 monthly collections were carried out utilizing the `survey at breeding site` method and examination of specimens kept in the Colecao Herpetologica do Setor de Zoologia da Universidade Federal de Santa Maria (ZUFSM). Tire Occurrence of 25 species of anurans was recorded. The anurofauna recorded represents 30% of the species known to Occur in Rio Grande do Sul, and comprises species generally associated with grasslands in this state and neighboring countries. Four reproductive modes were recorded: mode 1 (14 species: 58.3%) mode 11 and 30 (9 species` 37.5%) and mode 24 (1 species; 4.2%). The low diversification of reproductive modes is likely related to the homogeneity of the grassland habitat. Most species were constant or accessory in the Study area and the species abundance distribution patterns fit in the Broken Stick and Log-normal models. characterized by homogeneity of species abundance distribution. Most species showed great plasticity in habitat. but few were plastic in vocalization sites use. There was a weak positive correlation between species richness and precipitation. There was also a weak positive correlation between the abundance of species calling activity and maximum average temperatures. These correlations indicated that, in the study area. the abundance of calling males is more affected by the temperature, and species richness is more affected by precipitation, despite the fact that significantly higher species richness occurs during the hottest period of the year. These results showed that the climatological variables examined were not enough to explain the seasonal occurrence of species, thus the influence of other environmental variables merit to be tested in future studies.

(Review) Body mass index and risk of diabetes mellitus in the Asia-Pacific region

Few prospective data from the Asia Pacific region are available relating body mass index to the risk of diabetes. Our objective was to provide reliable age, sex and region specific estimates of the associations between body mass index and diabetes. Twenty-seven cohort studies from Asia, New Zealand and Australia, including 154,989 participants, contributed 1,244,793 person-years of follow-up. Outcome data included a combination of incidence of diabetes (based on blood glucose measurements) and fatal diabetes events. Hazard ratios were calculated from Cox models, stratified by sex and cohort, and adjusted for age at risk and smoking. During follow-up (mean = 8 years), 75 fatal diabetes events and 242 new cases of diabetes were documented. There were continuous positive associations between baseline body mass index and risk of diabetes with each 2 kg/m(2) lower body mass index associated with a 27% (23-30%) lower risk of diabetes. The associations were stronger in younger age groups, and regional comparisons demonstrated slightly stronger associations in Asian than in Australasian cohorts (P = 0.04). This overview provides evidence of a strong continuous association between body mass index and diabetes in the Asia Pacific region. The results indicate considerable potential for reduction in incidence of diabetes with population-wide lowering of body mass index in this region.

Conformational analysis of LYS(11-36), a peptide derived from the beta-sheet region of T4 lysozyme, in TFE and SDS

The solution conformation of a peptide LYS(11-36), which corresponds to the beta-sheet region in T4 lysozyme, has been examined in aqueous solution, TFE, and SDS micelles by CD and H-1 NMR spectroscopy. Secondary structure predictions suggest some beta-sheet and turn character in aqueous solution but predict a helical conformation in a more hydrophobic environment. The predictions were supported by the CD and NMR studies which showed the peptide to be relatively unstructured in aqueous solution, although there was some evidence of a beta-turn conformer which was maintained in 200 mM SDS and, to a lesser extent, in 50% TFE. The peptide was significantly helical in the presence of either 50% TFE or 200 mM SDS. TFE and SDS titrations showed that the peptide could form helical, sheet, or extended structure depending on the TFE or SDS concentration. The studies indicate that peptide environment is the determining factor in secondary structure adopted by LYS(11-36).

Regeneration of the Corneal Epithelium after Debridement of its Central Region: An Autoradiographic Study on Rabbits

Purpose: To investigate the proliferative behavior of the corneal and limbal epithelia after debridement on the central region of the rabbit cornea. Methods: After scraping a circular epithelial area, 5 mm in diameter, in the center of the cornea, (3)H-thymidine ((3)H-TdR) was injected intravitreally, and the rabbits killed from 1 to 49 days afterward. The cornea, together with the adjacent conjunctiva, was processed for autoradiography. Results: The regenerating epithelium at the center of the cornea exhibited high frequencies of labeled nuclei when compared to controls. The mitotic indexes for the limbus were comparable in experimental and control eyes. The unique basal stratum of the limbal epithelium exhibited quick proliferation and vertical migration in all eyes. Cells that remained labeled for four weeks or more were observed throughout the corneal epithelium, including its basal stratum, and this did not depend on epithelial damage. Conclusion: Corneal epithelium wounds are healed by sliding and proliferation of cells surrounding the epithelial gap without any evidence for the participation of the limbal epithelium. Daughter cells labeled with (3)H-TdR were visualized in all layers of the corneal epithelium up to 7 weeks after the DNA precursor injection. However, at this long interval, the only labeled cells in the limbus were in the suprabasal layers.

The plausibility of long-wavelength stress correlation or stress magnitude as a mechanism for precursory seismicity: Results from two simple elastic models

Observations of accelerating seismic activity prior to large earthquakes in natural fault systems have raised hopes for intermediate-term eartquake forecasting. If this phenomena does exist, then what causes it to occur? Recent theoretical work suggests that the accelerating seismic release sequence is a symptom of increasing long-wavelength stress correlation in the fault region. A more traditional explanation, based on Reid's elastic rebound theory, argues that an accelerating sequence of seismic energy release could be a consequence of increasing stress in a fault system whose stress moment release is dominated by large events. Both of these theories are examined using two discrete models of seismicity: a Burridge-Knopoff block-slider model and an elastic continuum based model. Both models display an accelerating release of seismic energy prior to large simulated earthquakes. In both models there is a correlation between the rate of seismic energy release with the total root-mean-squared stress and the level of long-wavelength stress correlation. Furthermore, both models exhibit a systematic increase in the number of large events at high stress and high long-wavelength stress correlation levels. These results suggest that either explanation is plausible for the accelerating moment release in the models examined. A statistical model based on the Burridge-Knopoff block-slider is constructed which indicates that stress alone is sufficient to produce accelerating release of seismic energy with time prior to a large earthquake.

Load-unload response ratio and accelerating moment/energy release critical region scaling and earthquake prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

Critical comparison of Julius Kruttschnitt Mineral Research Centre (JKMRC) blast fragmentation models

Blast fragmentation can have a significant impact on the profitability of a mine. An optimum run of mine (ROM) size distribution is required to maximise the performance of downstream processes. If this fragmentation size distribution can be modelled and controlled, the operation will have made a significant advancement towards improving its performance. Blast fragmentation modelling is an important step in Mine to Mill™ optimisation. It allows the estimation of blast fragmentation distributions for a number of different rock mass, blast geometry, and explosive parameters. These distributions can then be modelled in downstream mining and milling processes to determine the optimum blast design. When a blast hole is detonated rock breakage occurs in two different stress regions - compressive and tensile. In the-first region, compressive stress waves form a 'crushed zone' directly adjacent to the blast hole. The second region, termed the 'cracked zone', occurs outside the crush one. The widely used Kuz-Ram model does not recognise these two blast regions. In the Kuz-Ram model the mean fragment size from the blast is approximated and is then used to estimate the remaining size distribution. Experience has shown that this model predicts the coarse end reasonably accurately, but it can significantly underestimate the amount of fines generated. As part of the Australian Mineral Industries Research Association (AMIRA) P483A Mine to Mill™ project, the Two-Component Model (TCM) and Crush Zone Model (CZM), developed by the Julius Kruttschnitt Mineral Research Centre (JKMRC), were compared and evaluated to measured ROM fragmentation distributions. An important criteria for this comparison was the variation of model results from measured ROM in the-fine to intermediate section (1-100 mm) of the fragmentation curve. This region of the distribution is important for Mine to Mill™ optimisation. The comparison of modelled and Split ROM fragmentation distributions has been conducted in harder ores (UCS greater than 80 MPa). Further work involves modelling softer ores. The comparisons will be continued with future site surveys to increase confidence in the comparison of the CZM and TCM to Split results. Stochastic fragmentation modelling will then be conducted to take into account variation of input parameters. A window of possible fragmentation distributions can be compared to those obtained by Split . Following this work, an improved fragmentation model will be developed in response to these findings.

Linearization of a naturally occurring circular protein maintains structure but eliminates hemolytic activity

Cyclotides are a recently discovered family of disulfide rich proteins from plants that contain a circular protein backbone. They are exceptionally stable, as exemplified by their use in native medicine of the prototypic cyclotide kalata B1. The peptide retains uterotonic activity after the plant from which it is derived is boiled to make a medicinal tea. The circular backbone is thought to be in part responsible for the stability of the cyclotides, and to investigate its role in determining structure and biological activity, an acyclic derivative, des-(24-28)-kalata B1, was chemically synthesized and purified. This derivative has five residues removed from the 29-amino acid circular backbone of kalata B1 in a loop region corresponding to a processing site in the biosynthetic precursor protein. Two-dimensional NMR spectra of the peptide were recorded, assigned, and used to identify a series of distance, angle, and hydrogen bonding restraints. These were in turn used to determine a representative family of solution structures. Of particular interest was a determination of the structural similarities and differences between des-(2428)-kalata B1 and native kalata B1. Although the overall three-dimensional fold remains very similar to that of the native circular protein, removal of residues 24-28 of kalata B1 causes disruption of some structural features that are important to the overall stability. Furthermore, loss of hemolytic activity is associated with backbone truncation and linearization.