TUNPOP, a computer simulation model of the yellowfin tuna population and the surface tuna fishery of the Eastern Pacific Ocean

ENGLISH: Mathematical documentation of TUNP0P, an age-structured computer simulation model of the yellowfin tuna population and surface tuna fishery of the eastern Pacific Ocean, is described. Example runs of the model are presented and discussed, and the sensitivity of the model output to changes in various parameters is examined. SPANISH: Se describe la documentación matemática de TUNP0P, un modelo computador de simulación basado en la edad de la población del atún aleta amarilla y de la pesca atunera epipelágíca del Océano Pacífico oriental. Se presentan y se discuten ejemplos de las pasadas del modelo, y se examina la sensibilidad de los resultados de salida con relación a los cambios de varios parámetros. (PDF contains 47 pages.)

Increase in the density of resting microglia precedes neuritic plaque formation and microglial activation in a transgenic model of Alzheimer's disease

The formation of cerebral senile plaques composed of amyloid beta peptide (A beta) is a fundamental feature of Alzheimer's disease (AD). Glial cells and more specifically microglia become reactive in the presence of A beta. In a triple transgenic model of AD (3 x Tg-AD), we found a significant increase in activated microglia at 12 (by 111%) and 18 (by 88%) months of age when compared with non-transgenic (non-Tg) controls. This microglial activation correlated with A beta plaque formation, and the activation in microglia was closely associated with A beta plaques and smaller A beta deposits. We also found a significant increase in the area density of resting microglia in 3 x Tg-AD animals both at plaque-free stage (at 9 months by 105%) and after the development of A plaques (at 12 months by 54% and at 18 months by 131%). Our results show for the first time that the increase in the density of resting microglia precedes both plaque formation and activation of microglia by extracellular A beta accumulation. We suggest that AD pathology triggers a complex microglial reaction: at the initial stages of the disease the number of resting microglia increases, as if in preparation for the ensuing activation in an attempt to fight the extracellular A beta load that is characteristic of the terminal stages of the disease. Cell Death and Disease (2010) 1, e1; doi:10.1038/cddis.2009.2; published online 14 January 2010

Microscopic Model Of Nano-Scale Particles Removal In High Pressure Co2-Based Solvents

A physical model is presented to describe the kinds of static forces responsible for adhesion of nano-scale copper metal particles to silicon surface with a fluid layer. To demonstrate the extent of particle cleaning, Received in revised form equilibrium separation distance (ESD) and net adhesion force (NAF) of a regulated metal particle with different radii (10-300 nm) on the silicon surface in CO2-based cleaning systems under different pressures were simulated. Generally, increasing the pressure of the cleaning system decreased the net adhesion force between spherical copper particle and silicon surface entrapped with medium. For CO2 + isopropanol cleaning system, the equilibrium separation distance exhibited a maximum at temperature 313.15 K in the Equilibrium separation distance regions of pressure space (1.84-8.02 MPa). When the dimension of copper particle was given, for example, High pressure 50 nm radius particles, the net adhesion force decreased and equilibrium separation distance increased with increased pressure in the CO2 + H2O cleaning system at temperature 348.15 K under 2.50-12.67 MPa pressure range. However, the net adhesion force and equilibrium separation distance both decreased with an increase in surfactant concentration at given pressure (27.6 or 27.5 MPa) and temperature (318 or 298 K) for CO2 + H2O with surfactant PFPE COO-NH4+ or DiF(8)-PO4-Na+. (C) 2008 Elsevier B.V. All rights reserved.

A Trabecular Bone Explant Model of Osteocyte–Osteoblast Co-Culture for Bone Mechanobiology

The osteocyte network is recognized as the major mechanical sensor in the bone remodeling process, and osteocyte-osteoblast communication acts as an important mediator in the coordination of bone formation and turnover. In this study, we developed a novel 3D trabecular bone explant co-culture model that allows live osteocytes situated in their native extracellular matrix environment to be interconnected with seeded osteoblasts on the bone surface. Using a low-level medium perfusion system, the viability of in situ osteocytes in bone explants was maintained for up to 4 weeks, and functional gap junction intercellular communication (GJIC) was successfully established between osteocytes and seeded primary osteoblasts. Using this novel co-culture model, the effects of dynamic deformational loading, GJIC, and prostaglandin E-2 (PGE(2)) release on functional bone adaptation were further investigated. The results showed that dynamical deformational loading can significantly increase the PGE(2) release by bone cells, bone formation, and the apparent elastic modulus of bone explants. However, the inhibition of gap junctions or the PGE(2) pathway dramatically attenuated the effects of mechanical loading. This 3D trabecular bone explant co-culture model has great potential to fill in the critical gap in knowledge regarding the role of osteocytes as a mechano-sensor and how osteocytes transmit signals to regulate osteoblasts function and skeletal integrity as reflected in its mechanical properties.

Absolute measurement of roughness and lateral-correlation length of random surfaces by use of the simplified model of image-speckle contrast

We present a method of image-speckle contrast for the nonprecalibration measurement of the root-mean-square roughness and the lateral-correlation length of random surfaces with Gaussian correlation. We use the simplified model of the speckle fields produced by the weak scattering object in the theoretical analysis. The explicit mathematical relation shows that the saturation value of the image-speckle contrast at a large aperture radius determines the roughness, while the variation of the contrast with the aperture radius determines the lateral-correlation length. In the experimental performance, we specially fabricate the random surface samples with Gaussian correlation. The square of the image-speckle contrast is measured versus the radius of the aperture in the 4f system, and the roughness and the lateral-correlation length are extracted by fitting the theoretical result to the experimental data. Comparison of the measurement with that by an atomic force microscope shows our method has a satisfying accuracy. (C) 2002 Optical Society of America.

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

A model of the olfactory bulb and beyond

The olfactory bulb of mammals aids in the discrimination of odors. A mathematical model based on the bulbar anatomy and electrophysiology is described. Simulations of the highly non-linear model produce a 35-60 Hz modulated activity, which is coherent across the bulb. The decision states (for the odor information) in this system can be thought of as stable cycles, rather than as point stable states typical of simpler neuro-computing models. Analysis shows that a group of coupled non-linear oscillators are responsible for the oscillatory activities. The output oscillation pattern of the bulb is determined by the odor input. The model provides a framework in which to understand the transformation between odor input and bulbar output to the olfactory cortex. This model can also be extended to other brain areas such as the hippocampus, thalamus, and neocortex, which show oscillatory neural activities. There is significant correspondence between the model behavior and observed electrophysiology.

It has also been suggested that the olfactory bulb, the first processing center after the sensory cells in the olfactory pathway, plays a role in olfactory adaptation, odor sensitivity enhancement by motivation, and other olfactory psychophysical phenomena. The input from the higher olfactory centers to the inhibitory cells in the bulb are shown to be able to modulate the response, and thus the sensitivity, of the bulb to odor input. It follows that the bulb can decrease its sensitivity to a pre-existing and detected odor (adaptation) while remaining sensitive to new odors, or can increase its sensitivity to discover interesting new odors. Other olfactory psychophysical phenomena such as cross-adaptation are also discussed.

A GIS-based model of soil erosion and transport

Soil erosion is a natural process that occurs when the force of wind, raindrops or running water on the soil surface exceeds the cohesive forces that bind the soil together. In general, vegetation cover protects the soil from the effects of these erosive forces. However, land management activities such as ploughing, burning or heavy grazing may disturb this protective layer, exposing the underlying soil. The decision making process in rural catchment management is often supported by the predictive modelling of soil erosion and sediment transport processes within the catchment, using established techniques such as the Universal Soil Loss Equation [USLE] and the Agricultural Nonpoint Source pollution model [AGNPS]. In this article, the authors examine the range of erosion models currently available and describe the application of one of these to the Burrishoole catchment on the north-west coast of Ireland, which has suffered heavy erosion of blanket peat in recent years.