Unreported fishing, hungry people and political turmoil: the recipe for a food security crisis in Madagascar?

Madagascar, the world's fourth largest island, is one of the world's poorest developing countries, and its people depend heavily on marine resources for subsistence and income. Exports of these resources and foreign fishing access agreements are also important, at least from a large-scale economic perspective. In recent years, concerns have been voiced amongst local fishers and industry groups regarding the growth of the country's fishing effort. Despite these concerns, existing knowledge of the scale, composition and trends of Malagasy fisheries remains poor, and there is negligible information regarding unreported catches and illegal fishing in Madagascar's waters. Small-scale fisheries, which are often substantial in developing countries such as Madagascar, are often unreported or underestimated. Unfortunately, fisheries legislations, management plans and foreign fishing access agreements are often influenced by these incomplete data, leading to serious over-estimations of resource availability. This also appears to be the situation in Madagascar, where the reconstruction of total catches by all Malagasy fisheries sectors conducted here showed that total catches between 1950 and 2008 were twice the volume reported by national fisheries agencies. Most importantly, much of the subsistence sector is missing from official statistics, and signs of decline have already been observed in several stocks, suggesting that current levels of catches are likely to be exceeding sustainable yields. This has profound implications for the economic and ecological sustainability of fisheries, as well as food security in a country where people rely heavily on the ocean for their daily protein needs and livelihoods.     

An integrated zircon geochronological and geochemical investigation into the Miocene plutonic evolution of the Cyclades, Aegean Sea, Greece: Part 1: Geochronology

We use 369 individual U–Pb zircon ages from 14 granitoid samples collected on five islands in the Cyclades in the Aegean Sea, Greece, for constraining the crystallisation history of I- and S-type plutons above the retreating Hellenic subduction zone. Miocene magmatism in the Cyclades extended over a time span from 17 to 11 Ma. The ages for S-type granites are systematically ~2 million years older than those for I-type granites. Considering plutons individually, the zircon data define age spectra ranging from simple and unimodal to complex and multimodal. Seven of the 14 investigated samples yield more than one distinct zircon crystallisation age, with one I-type granodiorite sample from Mykonos Island representing the most complex case with three resolvable age peaks. Two samples from S-type granites on Ikaria appear to have crystallised zircon over 2–3 million years, whereas for the majority of individual samples with multiple zircon age populations the calculated ages deviate by 1–1.5 million years. We interpret our age data to reflect a protracted history involving initial partial melting at deeper lithospheric levels, followed by crystallisation and cooling at shallower crustal levels. Our study corroborates published research arguing that pluton construction is due to incremental emplacement of multiple magma pulses over a few million years. Assuming that multiple age peaks of our 14 samples can indeed serve to quantify time spans for magmatic emplacement, our data suggest that Aegean plutons were constructed over a few million years. Our tectonic interpretation of the U–Pb ages is that the S-type granites resulted from partial melting and migmatisation of the lower crust, possibly starting at ~23 Ma. The I-type granites and associated mafic melts are interpreted to reflect the magmatic arc stage in the Cyclades starting at ~15 Ma.     

TESSIN VISLab—laboratory for scientific visualization

Scientific visualization is an integral part of the modeling workflow, enabling researchers to understand complex or large data sets and simulation results. A high-resolution stereoscopic virtual reality (VR) environment further enhances the possibilities of visualization. Such an environment also allows collaboration in work groups including people of different backgrounds and to present results of research projects to stakeholders or the public. The requirements for the computing equipment driving the VR environment demand specialized software applications which can be run in a parallel fashion on a set of interconnected machines. Another challenge is to devise a useful data workflow from source data sets onto the display system. Therefore, we develop software applications like the OpenGeoSys Data Explorer, custom data conversion tools for established visualization packages such as ParaView and Visualization Toolkit as well as presentation and interaction techniques for 3D applications like Unity. We demonstrate our workflow by presenting visualization results for case studies from a broad range of applications. An outlook on how visualization techniques can be deeply integrated into the simulation process is given and future technical improvements such as a simplified hardware setup are outlined.     

Timing and source constraints on the relationship between mafic and felsic intrusions in the Emeishan large igneous province

Several I- and A-type granite, syenite plutons and spatially associated, giant Fe-Ti-V deposit-bearing mafic-ultramafic layered intrusions occur in the Pan-Xi (Panzhihua-Xichang) area within the inner zone of the Emeishan large igneous province (ELIP). These complexes are interpreted to be related to the Emeishan mantle plume. We present LA-ICP-MS and SIMS zircon U-Pb ages and Hf-Nd isotopic compositions for the gabbros, syenites and granites from these complexes. The dating shows that the age of the felsic intrusive magmatism (256.2 ± 3.0-259.8 ± 1.6 Ma) is indistinguishable from that of the mafic intrusive magmatism (255.4 ± 3.1-259.5 ± 2.7 Ma) and represents the final phase of a continuous magmatic episode that lasted no more than 10 Myr. The upper gabbros in the mafic-ultramafic intrusions are generally more isotopically enriched (lower ε_Nd and ε_Hf) than the middle and lower gabbros, suggesting that the upper gabbros have experienced a higher level of crustal contamination than the lower gabbros. The significantly positive ε_Hf(t) values of the A-type granites and syenites (+4.9 to +10.8) are higher than those of the upper gabbros of the associated mafic intrusion, which shows that they cannot be derived by fractional crystallization of these bodies. They are however identical to those of the mafic enclaves (+7.0 to +11.4) and middle and lower gabbros, implying that they are cogenetic. We suggest that they were generated by fractionation of large-volume, plume-related basaltic magmas that ponded deep in the crust. The deep-seated magma chamber erupted in two stages: the first near a density minimum in the basaltic fractionation trend and the second during the final stage of fractionation when the magma was a low density Fe-poor, Si-rich felsic magma. The basaltic magmas emplaced in the shallow-level magma chambers differentiated to form mafic-ultramafic layered intrusions accompanied by a small amount of crustal assimilation through roof melting. Evolved A-type granites (synenites and syenodiorites) were produced dominantly by crystallization in the deep crustal magma chamber. In contrast, the I-type granites have negative ε_Nd(t) [−6.3 to −7.5] and ε_Hf(t) [−1.3 to −6.7] values, with the Nd model ages () of 1.63−1.67 Ga and Hf model ages () of 1.56−1.58 Ga, suggesting that they were mainly derived from partial melting of Mesoproterozoic crust. In combination with previous studies, this study also shows that plume activity not only gave rise to reworking of ancient crust, but also significant growth of juvenile crust in the center of the ELIP.     

The last interglacial ocean

The final effort of the CLIMAP project was a study of the last interglaciation, a time of minimum ice volume some 122,000 yr ago coincident with the Substage 5e oxygen isotopic minimum. Based on detailed oxygen isotope analyses and biotic census counts in 52 cores across the world ocean, last interglacial sea-surface temperatures (SST) were compared with those today. There are small SST departures in the mid-latitude North Atlantic (warmer) and the Gulf of Mexico (cooler). The eastern boundary currents of the South Atlantic and Pacific oceans are marked by large SST anomalies in individual cores, but their interpretations are precluded by no-analog problems and by discordancies among estimates from different biotic groups. In general, the last interglacial ocean was not significantly different from the modern ocean. The relative sequencing of ice decay versus oceanic warming on the Stage 6/5 oxygen isotopic transition and of ice growth versus oceanic cooling on the Stage 5e/5d transition was also studied. In most of the Southern Hemisphere, the oceanic response marked by the biotic census counts preceded (led) the global ice-volume response marked by the oxygen-isotope signal by several thousand years. The reverse pattern is evident in the North Atlantic Ocean and the Gulf of Mexico, where the oceanic response lagged that of global ice volume by several thousand years. As a result, the very warm temperatures associated with the last interglaciation were regionally diachronous by several thousand years. These regional lead-lag relationships agree with those observed on other transitions and in long-term phase relationships; they cannot be explained simply as artifacts of bioturbational translations of the original signals.     

Comparing mixing-length models of the diabatic wind profile over homogeneous terrain

Models of the diabatic wind profile over homogeneous terrain for the entire atmospheric boundary layer are developed using mixing-length theory and are compared to wind speed observations up to 300 m at the National Test Station for Wind Turbines at Høvsøre, Denmark. The measurements are performed within a wide range of atmospheric stability conditions, which allows a comparison of the models with the average wind profile computed in seven stability classes, showing a better agreement than compared to the traditional surface-layer wind profile. The wind profile is measured by combining cup anemometer and lidar observations, showing good agreement at the overlapping heights. The height of the boundary layer, a parameter required for the wind profile models, is estimated under neutral and stable conditions using surface-layer turbulence measurements, and under unstable conditions based on the aerosol backscatter profile from ceilometer observations.     

Update of a Footprint-Based Approach for the Characterisation of Complex Measurement Sites

Horizontal heterogeneity can significantly affect the flux data quality at monitoring sites in complex terrain. In heterogeneous conditions, the adoption of the eddy-covariance technique is contraindicated by the lack of horizontal homogeneity and presence of advective conditions. In addition, uncertainty concerning the sources or sinks influencing a measurement compromises the data interpretation. The consideration of the spatial context of a measurement, defined by a footprint analysis, can therefore provide an important tool for data quality assessment. This study presents an update of an existing footprint-based quality evaluation concept for flux measurement sites in complex terrain. The most significant modifications in the present version are the use of a forward Lagrangian stochastic trajectory model for the determination of the spatial context of the measurements, and the determination of effective roughness lengths with a flux aggregation model in a pre-processing step. Detailed terrain data gathered by remote sensing methods are included. This approach determines spatial structures in the quality of flux data for varying meteorological conditions. The results help to identify terrain influences affecting the quality of flux data, such as dominating obstacles upwind of the site, or slopes biasing the wind field, so that the most suitable footprint regions for the collection of high-quality datasets can be identified. Additionally, the approach can be used to evaluate the performance of a coordinate rotation procedure, and to check to what extent the measured fluxes are representative for a target land-use type.     

Evolution of a post-batholith dike swarm in central coastal Queensland, Australia: arc-front to backarc?

A swarm of felsic and mafic dikes cuts a Late Carboniferous–Permian batholith called the Urannah Suite in central coastal Queensland. Late Permian–Triassic westward thrusting (Hunter–Bowen Orogeny) exposed this mid-crustal Suite and the crosscutting dikes, thus enabling examination of dikes that range in age from syn- to post-batholithic. Although both mafic and felsic dikes have the same dominant northerly strike, field, geochronologic and geochemical examination reveal that the swarm is composite; felsic dikes are older (285 Ma) and geochemically and isotopically distinct from mafic dikes (273–229 Ma). Dike compositions are compared to those of the host plutonic rocks, and to volcanic rocks the same age as the dikes. Whereas the felsic (older) dikes are compositionally similar to their host granites (initial ⁸⁷Sr/⁸⁶Sr>0.7045), the mafic (younger) dikes are isotopically (Sr, Nd, Pb) less radiogenic. Moreover, several different types of mafic dikes are evident based on geochemistry; most of these have mixed characteristics in terms of tectonic classification. All but two dikes of basalt and basaltic andesite composition classify as ‘with-in plate' on Ti–Zr–Y tectonic classification plots. Many of the basalts have high TiO₂ contents (1.5–3.0 wt.%). Most of these have REE and spider diagram patterns like calc-alkaline tholeiites, the exceptions being a few alkali basalts recognized by their alkali content, and high Ti, Ce, Nb and Zr contents. When put into the context of all plutonic rocks in the area (late Paleozoic and Mesozoic), these dikes record a transition at 280 Ma, after which time, all magmatism in the region is less isotopically evolved (initial ⁸⁷Sr/⁸⁶Sr=0.7033–0.7044). A model of slab retreat and hinge movement to the east in the latest Permian explains the change of geochemical signature from arc-front to backarc at about 280 Ma.     

Paleokarst system development in the San Andres Formation, Permian Basin, revealed by seismic characterization

Paleokarst systems are one of the major factors resulting in carbonate reservoir heterogeneity and compartmentalization. Nevertheless, few effective workflows have been proposed to map the 3D distribution of such systems. We describe a detailed seismic characterization approach integrating core, well log and rock physics analysis, to reveal a complex subsurface paleokarst system in the San Andres Formation, Permian basin, West Texas. In the area of high volume production, the collapsed paleokarst system is characterized by irregularly developed crackle and fracture breccias, mosaic breccias and cave fillings in the Upper San Andres Formation, which are delineated using seismic acoustic impedance. Along the transition from platform to basin, the paleokarst system is marked by a linear collapse including sags and small vertical faults that are recognizable in seismic imaging. Production data indicates that tight paleokarst zones cause reservoir compartmentalization and influence fluid communication between wells. The complex paleokarst system development is explained using a carbonate platform hydrological model, an outcrop analogue similar to modern marine hydrological environments within carbonate islands. Our method of model development for complex subsurface paleokarst systems may be applicable to other paleoenvironments.