Hierarchical location-allocation models for congested systems

In this paper we address the issue of locating hierarchical facilities in the presence of congestion. Two hierarchical models are presented, where lower level servers attend requests first, and then, some of the served customers are referred to higher level servers. In the first model, the objective is to find the minimum number of servers and theirlocations that will cover a given region with a distance or time standard. The second model is cast as a Maximal Covering Location formulation. A heuristic procedure is then presented together with computational experience. Finally, some extensions of these models that address other types of spatial configurations are offered.

Deep crustal electromagnetic structure of Bhuj earthquake region (India) and its implications

The existence of fluids and partial melt in the lower crust of the seismically active Kutch rift basin (on the western continental margin of India) owing to underplating has been proposed in previous geological and geophysical studies. This hypothesis is examined using magnetotelluric (MT) data acquired at 23 stations along two profiles across Kutch Mainland Uplift and Wagad Uplift. A detailed upper crustal structure is also presented using twodimensional inversion of MT data in the Bhuj earthquake (2001) area. The prominent boundaries of reflection in the upper crust at 5, 10 and 20 km obtained in previous seismic reflection profiles correlate with conductive structures in our models. The MT study reveals 1-2 km thick Mesozoic sediments under the Deccan trap cover. The Deccan trap thickness in this region varies from a few meters to 1.5 km. The basement is shallow on the northern side compared to the south and is in good agreement with geological models as well as drilling information. The models for these profiles indicate that the thickness of sediments would further increase southwards into the Gulf of Kutch. Significant findings of the present study indicate 1) the hypocentre region of the earthquake is devoid of fluids, 2) absence of melt (that is emplaced during rifting as suggested from the passive seismological studies) in the lower crust and 3) a low resistive zone in the depth range of 5-20 km. The present MT study rules out fluidsand melt (magma) as the causative factors that triggered the Bhuj earthquake. The estimated porosity value of 0.02% will explain 100-500 ohm·m resistivity values observed in the lower crust. Based on the seismic velocities and geochemical studies, presence of garnet is inferred. The lower crust consists of basalts - probably generated by partial melting of metasomatised garnet peridotite at deeper depths in the lithosphere - and their composition might be modified by reaction with the spinel peridotites.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

Empirical models for predicting the production of wild mushrooms in Scots pine (Pinus sylvestris L.) forests in the Central Pyrenees

Mushroom picking has become a widespread autumn recreational activity in the Central Pyrenees and other regions of Spain. Predictive models that relate mushroom production or fungal species richness with forest stand and site characteristics are not available. This study used mushroom production data from 24 Scots pine plots over 3 years to develop a predictive model that could facilitate forest management decisions when comparing silvicultural options in terms of mushroom production. Mixed modelling was used to model the dependence of mushroom production on stand and site factors. The results showed that productions were greatest when stand basal area was approximately 20 m2 ha-1. Increasing elevation and northern aspect increased total mushroom production as well as the production of edible and marketed mushrooms. Increasing slope decreased productions. Marketed Lactarius spp., the most important group collected in the region, showed similar relationships. The annual variation in mushroom production correlated with autumn rainfall. Mushroom species richness was highest when the total production was highest.

Particle shape and orientation in laser diffraction and static image analysis size distribution analysis of micrometer sized rectangular particles

Laser diffraction (LD) and static image analysis (SIA) of rectangular particles [United States Pharmacopeia, USP30-NF25, General Chapter <776>, Optical Miroscopy.] have been systematically studied. To rule out sample dispersion and particle orientation as the root cause of differences in size distribution profiles, we immobilize powder samples on a glass plate by means of a dry disperser. For a defined region of the glass plate, we measure the diffraction pattern as induced by the dispersed particles, and the 2D dimensions of the individual particles using LD and optical microscopy, respectively. We demonstrate a correlation between LD and SIA, with the scattering intensity of the individual particles as the dominant factor. In theory, the scattering intensity is related to the square of the projected area of both spherical and rectangular particles. In traditional LD the size distribution profile is dominated by the maximum projected area of the particles (A). The diffraction diameters of a rectangular particle with length L and breadth B as measured by the LD instrument approximately correspond to spheres of diameter ØL and ØB respectively. Differences in the scattering intensity between spherical and rectangular particles suggest that the contribution made to the overall LD volume probability distribution by each rectangular particle is proportional to A2/L and A2/B. Accordingly, for rectangular particles the scattering intensity weighted diffraction diameter (SIWDD) explains an overestimation of their shortest dimension and an underestimation of their longest dimension. This study analyzes various samples of particles whose length ranges from approximately 10 to 1000 μm. The correlation we demonstrate between LD and SIA can be used to improve validation of LD methods based on SIA data for a variety of pharmaceutical powders all with a different rectangular particle size and shape.

Trust region versus line search for computing the optical flow

We consider the numerical treatment of the optical flow problem by evaluating the performance of the trust region method versus the line search method. To the best of our knowledge, the trust region method is studied here for the first time for variational optical flow computation. Four different optical flow models are used to test the performance of the proposed algorithm combining linear and nonlinear data terms with quadratic and TV regularization. We show that trust region often performs better than line search; especially in the presence of non-linearity and non-convexity in the model.

The summer North Atlantic Oscillation in CMIP3 models and related uncertainties in projected summer drying in Europe

This paper discusses uncertainties in model projections of summer drying in the Euro-Mediterranean region related to errors and uncertainties in the simulation of the summer NAO (SNAO). The SNAO is the leading mode of summer SLP variability in the North Atlantic/European sector and modulates precipitation not only in the vicinity of the SLP dipole (northwest Europe) but also in the Mediterranean region. An analysis of CMIP3 models is conducted to determine the extent to which models reproduce the signature of the SNAO and its impact on precipitation and to assess the role of the SNAO in the projected precipitation reductions. Most models correctly simulate the spatial pattern of the SNAO and the dry anomalies in northwest Europe that accompany the positive phase. The models also capture the concurrent wet conditions in the Mediterranean, but the amplitude of this signal is too weak, especially in the east. This error is related to the poor simulation of the upper-level circulation response to a positive SNAO, namely the observed trough over the Balkans that creates potential instability and favors precipitation. The SNAO is generally projected to trend upwards in CMIP3 models, leading to a consistent signal of precipitation reduction in NW Europe, but the intensity of the trend varies greatly across models, resulting in large uncertainties in the magnitude of the projected drying. In the Mediterranean, because the simulated influence of the SNAO is too weak, no precipitation increase occurs even in the presence of a strong SNAO trend, reducing confidence in these projections.