Magnetic topology and current channels in plasmas with toroidal current density inversions

The equilibrium magnetic field inside axisymmetric plasmas with inversions on the toroidal current density is studied. Structurally stable non-nested magnetic surfaces are considered. For any inversion in the internal current density the magnetic families define several positive current channels about a central negative one. A general expression relating the positive and negative currents is derived in terms of a topological anisotropy parameter. Next, an analytical local solution for the poloidal magnetic flux is derived and shown compatible with current hollow magnetic pitch measurements shown in the literature. Finally, the analytical solution exhibits non-nested magnetic families with positive anisotropy, indicating that the current inside the positive channels have at least twice the magnitude of the central one.

SIS: a program to generate draft genome sequence scaffolds for prokaryotes

Background: Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Results: We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of SIS, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that SIS has overall better performance. Conclusions: SIS is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of SIS in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

Methane airborne measurements and comparison to global models during BARCA

Tropical regions, especially the Amazon region, account for large emissions of methane (CH4). Here, we present CH4 observations from two airborne campaigns conducted within the BARCA (Balanco Atmosferico Regional de Carbono na Amazonia) project in the Amazon basin in November 2008 (end of the dry season) and May 2009 (end of the wet season). We performed continuous measurements of CH4 onboard an aircraft for the first time in the Amazon region, covering the whole Amazon basin with over 150 vertical profiles between altitudes of 500 m and 4000 m. The observations support the finding of previous ground-based, airborne, and satellite measurements that the Amazon basin is a large source of atmospheric CH4. Isotope analysis verified that the majority of emissions can be attributed to CH4 emissions from wetlands, while urban CH4 emissions could be also traced back to biogenic origin. A comparison of five TM5 based global CH4 inversions with the observations clearly indicates that the inversions using SCIAMACHY observations represent the BARCA observations best. The calculated CH4 flux estimate obtained from the mismatch between observations and TM5-modeled CH4 fields ranges from 36 to 43 mg m(-2) d(-1) for the Amazon lowland region.