Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Is maximizing resilience compatible with established ecological goal functions?

Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] used a simple phytoplanktonzooplankton-nutrient model and a genetic algorithm to determine the parameter values that would maximize the value of certain goal functions. These goal functions were to maximize biomass, maximize flux, maximize flux to biomass ratio, and maximize resilience. It was found that maximizing goal functions maximized resilience. The objective of this study was to investigate whether the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] result was indicative of a general ecosystem principle, or peculiar to the model and parameter ranges used. This study successfully replicated the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] experiment for a number of different model types, however, a different interpretation of the results is made. A new metric, concordance, was devised to describe the agreement between goal functions. It was found that resilience has the highest concordance of all goal functions trialled. for most model types. This implies that resilience offers a compromise between the established ecological goal functions. The parameter value range used is found to affect the parameter versus goal function relationships. Local maxima and minima affected the relationship between parameters and goal functions, and between goal functions. (C) 2003 Elsevier B.V. All rights reserved.

Unsaturated diffusivity functions for concrete derived from NMR images

Deterioration of concrete or reinforcing steel through excessive contaminant concentration is often the result of repeated wetting and drying cycles. At each cycle, the absorption of water carries new contaminants into the unsaturated concrete. Nuclear Magnetic Resonance (NMR) is used with large concrete samples to observe the shape of the wetting profile during a simple one-dimensional wetting process. The absorption of water by dry concrete is modelled by a nonlinear diffusion equation with the unsaturated hydraulic diffusivity being a strongly nonlinear function of the moisture content. Exponential and power functions are used for the hydraulic diffusivity and corresponding solutions of the diffusion equation adequately predict the shape of the experimental wetting profile. The shape parameters, describing the wetting profile, vary little between different blends and are relatively insensitive to subsequent re-wetting experiments allowing universal parameters to be suggested for these concretes.

Usefulness of B-type natriuretic peptide in hypertensive patients with exertional dyspnea and normal left ventricular ejection fraction and correlation with new echocardiographic indexes of systolic and diastolic function

B-type natriuretic peptide (BNP) levels increase in systolic heart failure (HF). However, the value of BNP in hypertensive patients with suspected diastolic HF (symptoms suggestive of HF but normal ejection fraction) and its relation to myocardial function in these patients is unclear. We prospectively studied 72 ambulatory hypertensive subjects (40 women, mean age 58 +/- 8 years) with exertional dyspnea and ejection fraction greater than or equal to50%. Diastolic function was evaluated with transmitral and pulmonary venous Doppler, mitral annular velocities (pulsed-wave tissue Doppler), and flow propagation velocity (color M-mode). Systolic function was assessed with strain and strain rate derived from color tissue Doppler imaging. BNP was related to myocardial function and the presence or absence of global diastolic dysfunction. By conventional Doppler criteria, 34 patients had normal left ventricular diastolic function and 38 had isolated diastolic dysfunction. BNP values were higher in patients with diastolic dysfunction (46 +/- 48 vs 20 +/- 20 pg/ml, p = 0.004) and were related independently to blood pressure, systolic strain rate, left atrial function (p < 0.01 for all), and age (p = 0.015). Patients with diastolic dysfunction and pseudonormal filling had higher BNP levels compared with impaired relaxation (89 +/- 47 vs 35 +/- 42 pg/ml, p = 0.001). However, 79% of patients with diastolic dysfunction had BNP levels within the normal range. We conclude that in ambulatory hypertensive patients with symptoms suggestive of mild HF and normal ejection fraction, BNP is related to atrial and ventricular systolic parameters, blood pressure, and age. Although elevated in the presence of diastolic dysfunction, the BNP level mostly is in the normal range and, therefore, has limited diagnostic value in stable patients with suspected diastolic HF. (C) 2003 by Excerpta Medica, Inc.

Poor correlation of mid-femoral measurements by CT and hip measurements by DXA in the elderly

Background and aims: Hip fracture is a devastating event in terms of outcome in the elderly, and the best predictor of hip fracture risk is hip bone density, usually measured by dual X-ray absorptiometry (DXA). However, bone density can also be ascertained from computerized tomography (CT) scans, and mid-thigh scans are frequently employed to assess the muscle and fat composition of the lower limb. Therefore, we examined if it was possible to predict hip bone density using mid-femoral bone density. Methods: Subjects were 803 ambulatory white and black women and men, aged 70-79 years, participating in the Health, Aging and Body Composition (Health ABC) Study. Bone mineral content (BMC, g) and volumetric bone mineral density (vBMD, mg/cm(3)) of the mid-femur were obtained by CT, whereas BMC and areal bone mineral density (aBMD, g/cm(2)) of the hip (femoral neck and trochanter) were derived from DXA. Results: In regression analyses stratified by race and sex, the coefficient of determination was low with mid-femoral BMC, explaining 6-27% of the variance in hip BMC, with a standard error of estimate (SEE) ranging from 16 to 22% of the mean. For mid-femur vBMD, the variance explained in hip aBMD was 2-17% with a SEE ranging from 15 to 18%. Adjusting aBMD to approximate volumetric density did not improve the relationships. In addition, the utility of fracture prediction was examined. Forty-eight subjects had one or more fractures (various sites) during a mean follow-up of 4.07 years. In logistic regression analysis, there was no association between mid-femoral vBMD and fracture (all fractures), whereas a 1 SD increase in hip BMD was associated with reduced odds for fracture of similar to60%. Conclusions: These results do not support the use of CT-derived mid-femoral vBMD or BMC to predict DXA-measured hip bone mineral status, irrespective of race or sex in older adults. Further, in contrast to femoral neck and trochanter BMD, mid-femur vBMD was not able to predict fracture (all fractures). (C) 2003, Editrice Kurtis.

Mapping the royal road and other hierarchical functions

In this paper we present a technique for visualising hierarchical and symmetric, multimodal fitness functions that have been investigated in the evolutionary computation literature. The focus of this technique is on landscapes in moderate-dimensional, binary spaces (i.e., fitness functions defined over {0, 1}(n), for n less than or equal to 16). The visualisation approach involves an unfolding of the hyperspace into a two-dimensional graph, whose layout represents the topology of the space using a recursive relationship, and whose shading defines the shape of the cost surface defined on the space. Using this technique we present case-study explorations of three fitness functions: royal road, hierarchical-if-and-only-if (H-IFF), and hierarchically decomposable functions (HDF). The visualisation approach provides an insight into the properties of these functions, particularly with respect to the size and shape of the basins of attraction around each of the local optima.

Correlation between grain nutritional content and pasting properties of pre-gelatinized red rice flour

As with any variety of rice, red rice characteristics are subject to varietal differences, growing conditions, types of processing, and nutritional and rheological properties. This study determined the nutritional characteristics (centesimal composition and minerals) and paste viscosity properties of raw grains of four red rice genotypes (Tradicional MNAPB0405, MNACE0501 and MNACH0501) and the paste viscosity properties of pre-gelatinized flours obtained at different cooking times (20, 30 and 40 min). The main nutritional properties were correlated with the pasting properties of the pre-gelatinized flours. The samples showed differences in nutritional properties and paste viscosity. MNAPB0405 and MNACE0501 showed higher levels of fiber and fat and provided higher caloric energy than Tradicional and MNACH0501, which, in turn, showed higher levels of amylose. MNACH0501 showed higher peak viscosity (2402 cP), higher breakdown viscosity (696 cP) and a greater tendency to retrogradation (1510 cP), while Tradicional, MNAPB0405 and MNACE0501 had pasting profiles with peak viscosities varying between 855 and 1093 cP, breaking viscosity below 85 cP and retrogradation tendency between 376 and 1206 cP. The factors genotype and cooking time influenced the rheological behavior of pre-gelatinized flours, decreasing their pasting properties. The protein and amylose levels are correlated with the pasting properties and can be used as indicators of these properties in different genotypes of red rice, whether raw or processed into pre-gelatinized flours.

Health education: the role and functions of the specialist and the generalist

Education for health is a process in which all public health and medical care personnel are involved. People learn both formally (planned learning experiences) and informally (unplanned learning experiences). Since the patient, the client, the consummer and the community expect public health and medical care personnel to assist them with health and disease issues and problems, the response of the professional "educates" the customer whether the professional intends to educate or not. Therefore, it is incumbent on all public health and medical care professionals to understand their educational functions and their role in health education. It is also important that the role of the specialist in education be clear. The specialist, as to all other specialists, has an in-depth knowledge of his area of expertise, i.e., the teaching/learning process; s/he may function as a consultant to others to enhance the educational potential of their role or s/he may work with a team or with communities or groups of patients. Specific competencies and knowledge are required of the health education specialist; and there is a body of learning and social change theory which provides a frame of reference for planning, implementing and evaluating educational programs. Working with others to enhance their potential to learn and to make informed decisions about health/disease issues is the hallmark of the health education specialist.

Equilibrium self-assembly of colloids with distinct interaction sites: Thermodynamics, percolation, and cluster distribution functions

We calculate the equilibrium thermodynamic properties, percolation threshold, and cluster distribution functions for a model of associating colloids, which consists of hard spherical particles having on their surfaces three short-ranged attractive sites (sticky spots) of two different types, A and B. The thermodynamic properties are calculated using Wertheim's perturbation theory of associating fluids. This also allows us to find the onset of self-assembly, which can be quantified by the maxima of the specific heat at constant volume. The percolation threshold is derived, under the no-loop assumption, for the correlated bond model: In all cases it is two percolated phases that become identical at a critical point, when one exists. Finally, the cluster size distributions are calculated by mapping the model onto an effective model, characterized by a-state-dependent-functionality (f) over bar and unique bonding probability (p) over bar. The mapping is based on the asymptotic limit of the cluster distributions functions of the generic model and the effective parameters are defined through the requirement that the equilibrium cluster distributions of the true and effective models have the same number-averaged and weight-averaged sizes at all densities and temperatures. We also study the model numerically in the case where BB interactions are missing. In this limit, AB bonds either provide branching between A-chains (Y-junctions) if epsilon(AB)/epsilon(AA) is small, or drive the formation of a hyperbranched polymer if epsilon(AB)/epsilon(AA) is large. We find that the theoretical predictions describe quite accurately the numerical data, especially in the region where Y-junctions are present. There is fairly good agreement between theoretical and numerical results both for the thermodynamic (number of bonds and phase coexistence) and the connectivity properties of the model (cluster size distributions and percolation locus).

On differentiability and analyticity of positive definite functions

We derive a set of differential inequalities for positive definite functions based on previous results derived for positive definite kernels by purely algebraic methods. Our main results show that the global behavior of a smooth positive definite function is, to a large extent, determined solely by the sequence of even-order derivatives at the origin: if a single one of these vanishes then the function is constant; if they are all non-zero and satisfy a natural growth condition, the function is real-analytic and consequently extends holomorphically to a maximal horizontal strip of the complex plane.

Stratification of the earth beneath the Azores from P and S receiver functions

Seismic recordings of IRIS/IDA/GSN station CMLA and of several temporary stations in the Azores archipelago are processed with P and S receiver function (PRF and SRF) techniques. Contrary to regional seismic tomography these methods provide estimates of the absolute velocities and of the Vp/Vs ratio up to a depth of similar to 300 km. Joint inversion of PRFs and SRFs for a few data sets consistently reveals a division of the subsurface medium into four zones with a distinctly different Vp/Vs ratio: the crust similar to 20 km thick with a ratio of similar to 1.9 in the lower crust, the high-Vs mantle lid with a strongly reduced VpNs velocity ratio relative to the standard 1.8, the low-velocity zone (LVZ) with a velocity ratio of similar to 2.0, and the underlying upper-mantle layer with a standard velocity ratio. Our estimates of crustal thickness greatly exceed previous estimates (similar to 10 km). The base of the high-Vs lid (the Gutenberg discontinuity) is at a depth of-SO km. The LVZ with a reduction of S velocity of similar to 15% relative to the standard (IASP91) model is terminated at a depth of similar to 200 km. The average thickness of the mantle transition zone (TZ) is evaluated from the time difference between the S410p and SKS660p, seismic phases that are robustly detected in the S and SKS receiver functions. This thickness is practically similar to the standard IASP91 value of 250 km. and is characteristic of a large region of the North Atlantic outside the Azores plateau. Our data are indicative of a reduction of the S-wave velocity of several percent relative to the standard velocity in a depth interval from 460 to 500 km. This reduction is found in the nearest vicinities of the Azores, in the region sampled by the PRFs, but, as evidenced by SRFs, it is missing at a distance of a few hundred kilometers from the islands. We speculate that this anomaly may correspond to the source of a plume which generated the Azores hotspot. Previously, a low S velocity in this depth range was found with SRF techniques beneath a few other hotspots.

From efficiency to optimality in proportional reinsurance under group correlation

Based on our recent discovery of closed form formulae of efficient Mean Variance retentions in variable quota-share proportional reinsurance under group correlation, we analyzed the influence of different combination of correlation and safety loading levels on the efficient frontier, both in a single period stylized problem and in a multiperiod one.