Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.

Modeling harvest rates and numbers from age and sex ratios: A demonstration for elephant populations

Illegal harvest rates of wildlife populations are often unknown or difficult to estimate from field data due to under-reporting or incomplete detection of carcasses. This is especially true for elephants that are killed for ivory or in conflicts with people. We describe a method to infer harvest rates from coarse field data of three population parameters, namely, adult female to male ratio, male old-adult to young-adult ratio, and proportion of adult males in the population using Jensen's (2000) 2-sex, density-dependent Leslie matrix model. The specific combination of male and female harvest rates and numbers can be determined from the history of harvest and estimate of population size. We applied this technique to two populations of elephants for which data on age structure and records of mortality were available-a forest-dwelling population of the Asian elephant (at Nagarahole, India) and an African savannah elephant population (at Samburu, Kenya) that had experienced male-biased harvest regimes over 2-3 decades. For the Nagarahole population, the recorded numbers of male and female elephants killed illegally during 1981-2000 were 64% and 88% of the values predicted by the model, respectively, implying some non-detection or incomplete reporting while for the Samburu population the recorded and modeled numbers of harvest during 1990-1999 closely matched. This technique, applicable to any animal population following logistic growth model, can be especially useful for inferring illegal harvest numbers of forest elephants in Africa and Asia.

Integrating biodiversity and conservation with modern agricultural landscapes

To achieve food security and meet the demands of the ever-growing human populations, farming systems have assumed unsustainable practices to produce more from a finite land area. This has been cause for concern mainly due to the often-irreversible damage done to the otherwise productive agricultural landscapes. Agro-ecology is proclaimed to be deteriorating due to eroding integrity of connected ecological mosaics and vulnerability to climate change. This has contributed to declining species diversity, loss of buffer vegetation, fragmentation of habitats, and loss of natural pollinators or predators, which eventually leads to decline in ecosystem services. Currently, a hierarchy of conservation initiatives is being considered to restore ecological integrity of agricultural landscapes. However, the challenge of identifying a suitable conservation strategy is a daunting task in view of socio-ecological factors that may constrain the choice of available strategies. One way to mitigate this situation and integrate biodiversity with agricultural landscapes is to implement offset mechanisms, which are compensatory and balancing approaches to restore the ecological health and function of an ecosystem. This needs to be tailored to the history of location specific agricultural practices, and the social, ecological and environmental conditions. The offset mechanisms can complement other initiatives through which farmers are insured against landscape-level risks such as droughts, fire and floods. For countries in the developing world with significant biodiversity and extensive agriculture, we should promote a comprehensive model of sustainable agricultural landscapes and ecosystem services, replicable at landscape to regional scales. Arguably, the model can be a potential option to sustain the integrity of biodiversity mosaic in agricultural landscapes.

Identifying Priority Areas for Conservation and Management in Diverse Tropical Forests

The high concentration of the world's species in tropical forests endows these systems with particular importance for retaining global biodiversity, yet it also presents significant challenges for ecology and conservation science. The vast number of rare and yet to be discovered species restricts the applicability of species-level modelling for tropical forests, while the capacity of community classification approaches to identify priorities for conservation and management is also limited. Here we assessed the degree to which macroecological modelling can overcome shortfalls in our knowledge of biodiversity in tropical forests and help identify priority areas for their conservation and management. We used 527 plant community survey plots in the Australian Wet Tropics to generate models and predictions of species richness, compositional dissimilarity, and community composition for all the 4,313 vascular plant species recorded across the region (>1.3 million communities (grid cells)). We then applied these predictions to identify areas of tropical forest likely to contain the greatest concentration of species, rare species, endemic species and primitive angiosperm families. Synthesising these alternative attributes of diversity into a single index of conservation value, we identified two areas within the Australian wet tropics that should be a high priority for future conservation actions: the Atherton Tablelands and Daintree rainforest. Our findings demonstrate the value of macroecological modelling in identifying priority areas for conservation and management actions within highly diverse systems, such as tropical forests.

Forest area estimation and reporting: implications for conservation, management and REDD

Periodic estimation, monitoring and reporting on area under forest and plantation types and afforestation rates are critical to forest and biodiversity conservation, sustainable forest management and for meeting international commitments. This article is aimed at assessing the adequacy of the current monitoring and reporting approach adopted in India in the context of new challenges of conservation and reporting to international conventions and agencies. The analysis shows that the current mode of monitoring and reporting of forest area is inadequate to meet the national and international requirements. India could be potentially over-reporting the area under forests by including many non-forest tree categories such as commercial plantations of coconut, cashew, coffee and rubber, and fruit orchards. India may also be under-reporting deforestation by reporting only gross forest area at the state and national levels. There is a need for monitoring and reporting of forest cover, deforestation and afforestation rates according to categories such as (i) natural/primary forest, (ii) secondary/degraded forests, (iii) forest plantations, (iv) commercial plantations, (v) fruit orchards and (vi) scattered trees.

Hydrogen bonding, halogen bonding and lithium bonding: an atoms in molecules and natural bond orbital perspective towards conservation of total bond order, inter- and intra-molecular bonding

One hundred complexes have been investigated exhibiting D-X center dot center dot center dot A interactions, where X = H, Cl or Li and DX is the `X bond' donor and A is the acceptor. The optimized structures of all these complexes have been used to propose a generalized `Legon-Millen rule' for the angular geometry in all these interactions. A detailed Atoms in Molecules (AIM) theoretical analysis confirms an important conclusion, known in the literature: there is a strong correlation between the electron density at the X center dot center dot center dot A bond critical point (BCP) and the interaction energy for all these interactions. In addition, we show that extrapolation of the fitted line leads to the ionic bond for Li-bonding (electrostatic) while for hydrogen and chlorine bonding, it leads to the covalent bond. Further, we observe a strong correlation between the change in electron density at the D-X BCP and that at the X center dot center dot center dot A BCP, suggesting conservation of the bond order. The correlation found between penetration and electron density at BCP can be very useful for crystal structure analysis, which relies on arbitrary van der Waals radii for estimating penetration. Various criteria proposed for shared-and closed-shell interactions based on electron density topology have been tested for H/Cl/Li bonded complexes. Finally, using the natural bond orbital (NBO) analysis it is shown that the D-X bond weakens upon X bond formation, whether it is ionic (DLi) or covalent (DH/DCl) and the respective indices such as ionicity or covalent bond order decrease. Clearly, one can think of conservation of bond order that includes ionic and covalent contributions to both D-X and X center dot center dot center dot A bonds, for not only X = H/Cl/Li investigated here but also any atom involved in intermolecular bonding.

Distribution, Demography, and Conservation of Lion-tailed Macaques (Macaca silenus) in the Anamalai Hills Landscape, Western Ghats, India

The status of the endemic and endangered lion-tailed macaque (Macaca silenus) has not been properly assessed in several regions of the Western Ghats of southern India. We conducted a study in Parambikulam Forest Reserve in the state of Kerala to determine the distribution, demography, and status of lion-tailed macaques. We laid 5km(2) grid cells on the map of the study area (644km(2)) and made four replicated walks in each grid cell using GPS. We gathered data on lion-tailed macaque group locations, demography, and site covariates including trail length, duration of walk, proportion of evergreen forest, height of tallest trees, and human disturbance index. We also performed occupancy modeling using PRESENCE ver. 3.0. We estimated a minimum of 17 groups of macaques in these hills. Low detection and occupancy probabilities indicated a low density of lion-tailed macaques in the study area. Height of the tallest trees correlated positively whereas human disturbance and proportion of evergreen forest correlated negatively with occupancy in grid cells. We also used data from earlier studies carried out in the surrounding Anamalai Tiger Reserve and Nelliyampathy Hills to discuss the conservation status in the large Anamalai Hills Landscape. This landscape harbors an estimated population of 1108 individuals of lion-tailed macaques, which is about one third of the entire estimated wild population of this species. A conservation plan for this landscape could be used as a model for conservation in other regions of the Western Ghats.

Identifying regions for conservation of sloth bears through occupancy modelling in north-eastern Karnataka, India

In the absence of information on species in decline with contracting ranges, management should emphasize remaining populations and protection of their habitats. Threatened by anthropogenic pressure including habitat degradation and loss, sloth bears (Melursus ursinus) in India have become limited in range, habitat, and population size. We identified ecological and anthropogenic determinants of occurrence within an occupancy framework to evaluate habitat suitability of non-protected regions (with sloth bears) in northeastern Karnataka, India. We employed a systematic sampling methodology to yield presence absence data to examine a priori hypotheses of determinants that affected occupancy. These covariates were broadly classified as habitat or anthropogenic factors. Mean number of termite mounds and trees positively influenced sloth bear occupancy, and grazing pressure expounded by mean number of livestock dung affected it negatively. Also, mean percentage of shrub coverage had no impact on bear inhabitance. The best fitting model further predicted habitats in Bukkasagara, Agoli, and Benakal reserved forests to have 38%, 75%, and 88%, respectively, of their sampled grid cells with high occupancies (>0.70) albeit little or no legal protection. We recommend a conservation strategy that includes protection of vegetation stand-structure, maintenance of soil moisture, and enrichment of habitat for the long-term welfare of this species.

Conservation Properties of the Trapezoidal Rule in Linear Time Domain Analysis of Acoustics and Structures

The trapezoidal rule, which is a special case of the Newmark family of algorithms, is one of the most widely used methods for transient hyperbolic problems. In this work, we show that this rule conserves linear and angular momenta and energy in the case of undamped linear elastodynamics problems, and an ``energy-like measure'' in the case of undamped acoustic problems. These conservation properties, thus, provide a rational basis for using this algorithm. In linear elastodynamics problems, variants of the trapezoidal rule that incorporate ``high-frequency'' dissipation are often used, since the higher frequencies, which are not approximated properly by the standard displacement-based approach, often result in unphysical behavior. Instead of modifying the trapezoidal algorithm, we propose using a hybrid finite element framework for constructing the stiffness matrix. Hybrid finite elements, which are based on a two-field variational formulation involving displacement and stresses, are known to approximate the eigenvalues much more accurately than the standard displacement-based approach, thereby either bypassing or reducing the need for high-frequency dissipation. We show this by means of several examples, where we compare the numerical solutions obtained using the displacement-based and hybrid approaches against analytical solutions.