Constrains on the structure of Maio Island (Cape Verde) by a three-dimensional gravity model: Imaging partially exhumed magma chambers

We propose a 3-D gravity model for the volcanic structure of the island of Maio (Cape Verde archipelago) with the objective of solving some open questions concerning the geometry and depth of the intrusive Central Igneous Complex. A gravity survey was made covering almost the entire surface of the island. The gravity data was inverted through a non-linear 3-D approach which provided a model constructed in a random growth process. The residual Bouguer gravity field shows a single positive anomaly presenting an elliptic shape with a NWSE trending long axis. This Bouguer gravity anomaly is slightly off-centred with the island but its outline is concordant with the surface exposure of the Central Igneous Complex. The gravimetric modelling shows a high-density volume whose centre of mass is about 4500 m deep. With increasing depth, and despite the restricted gravimetric resolution, the horizontal sections of the model suggest the presence of two distinct bodies, whose relative position accounts for the elongated shape of the high positive Bouguer gravity anomaly. These bodies are interpreted as magma chambers whose coeval volcanic counterparts are no longer preserved. The orientation defined by the two bodies is similar to that of other structures known in the southern group of the Cape Verde islands, thus suggesting a possible structural control constraining the location of the plutonic intrusions.

Gravity Data and Crustal Structure of North-Central Brazil

The Bouguer gravity anomaly of the northwest Ceará state in north-central Brazil was separated into its regional and residual components which were interpreted separately. By assuming that the sources of the regional anomalies are the depth variations of the crust-mantle interface, the mapping of these variations permited identifying crustal thickening zones which may be related to regional structures. The gravity residual sources coincide with occurrences of high-grade rocks (granulites) associated to medium-grade gneisses. Besides, the major strike-slip zones present significant signatures in the gravity data. This geophysical interpretation is compatible with the interpretation that the tectonic framework of the area is related to two crustal blocks conjoined by an A-type suture. The blocks are displaced along an oblique ramp with dextral movement, which played an important role in uplifting high-grade rocks from the lower crust to upper crustal levels. The suture zone corresponds to an imbricated compressive system dipping to the east and complicated by late dextral strike-slip shear zones.

The Southeast Indian Ridge between 88°E and 118°E: Gravity anomalies and crustal accretion at intermediate spreading rates

Although slow spreading ridges characterized by a deep axial valley and fast spreading ridges characterized by an axial bathymetric high have been extensively studied, the transition between these two modes of axial morphology is not well understood. We conducted a geophysical-survey of the intermediate spreading rate Southeast Indian Ridge between 88 degrees E and 118 degrees E, a 2300-km-long section of the ridge located between the Amsterdam hot spot and the Australian-Antarctic Discordance where satellite gravity data suggest that the Southeast Indian Ridge (SEIR) undergoes a change from an axial high in the west to an axial valley in the east. A basic change in axial morphology is found near 103 degrees 30'E in the shipboard data; the axis to the west is marked by an axial high, while a valley is found to the east. Although a well-developed axial high, characteristic of the East Pacific Rise (EPR), is occasionally present, the more common observation is a rifted high that is lower and pervasively faulted, sometimes with significant (> 50 m throw) faults within a kilometer of the axis. A shallow axial valley (< 700 m deep) is observed from 104 degrees E to 114 degrees E with a sudden change to a deep (>1200 m deep) valley across a transform at 114 degrees E. The changes in axial morphology along the SEIR are accompanied by a 500 m increase in near-axis ridge flank depth from 2800 m near 88 degrees E to 3300 m near 114 degrees E and by a 50 mGal increase in the regional level of mantle Bouguer gravity anomalies over the same distance, The regional changes in depth and mantle Bouguer anomaly (MBA) gravity can be both explained by a 1.7-2.4 km change in crustal thickness or by a mantle temperature change of 50 degrees C-90 degrees C. In reality, melt supply (crustal thickness) and mantle temperature are linked, so that changes in both may occur simultaneously and these estimates serve as upper bounds. The along-axis MBA gradient is not uniform. Pronounced steps in the regional level of the MBA gravity occur at 103 degrees 30'E-104 degrees E and at 114 degrees E-116 degrees E and correspond to the changes in the nature of the axial morphology and in the amplitude of abyssal hill morphology suggesting that the different forms of morphology do not grade into each other but rather represent distinctly different forms of axial (s)tructure and tectonics with a sharp transition between them. The change from an axial high to an axial valley requires a threshold effect in which the strength of the lithosphere changes quickly. The presence or absence of a quasi-steady state magma chamber may provide such a mechanism. The different forms of axial morphology are also associated with different intrasegment MBA gravity patterns. Segments with an axial high have an MBA low located at a depth minimum near the center of the segment, At EPR-like segments, the MBA low is about 10 mGal with along-axis gradients of 0.15-0.25 mGal/km, similar to those observed at the EPR, Rifted highs have a shallower low and lower gradients suggesting an attenuated composite magma chamber and a reduced and perhaps episodic melt supply. Segments with a shallow axial valley have very flat along-axis MBA profiles with little correspondence between axial depth and axial MBA gravity.

Last Interglacial synthesis of high-latitude temperature: temperature anomalies and associated errors for 4 time slices

The Last Interglacial (LIG, 129-116 thousand of years BP, ka) represents a test bed for climate model feedbacks in warmer-than-present high latitude regions. However, mainly because aligning different palaeoclimatic archives and from different parts of the world is not trivial, a spatio-temporal picture of LIG temperature changes is difficult to obtain. Here, we have selected 47 polar ice core and sub-polar marine sediment records and developed a strategy to align them onto the recent AICC2012 ice core chronology. We provide the first compilation of high-latitude temperature changes across the LIG associated with a coherent temporal framework built between ice core and marine sediment records. Our new data synthesis highlights non-synchronous maximum temperature changes between the two hemispheres with the Southern Ocean and Antarctica records showing an early warming compared to North Atlantic records. We also observe warmer than present-day conditions that occur for a longer time period in southern high latitudes than in northern high latitudes. Finally, the amplitude of temperature changes at high northern latitudes is larger compared to high southern latitude temperature changes recorded at the onset and the demise of the LIG. We have also compiled four data-based time slices with temperature anomalies (compared to present-day conditions) at 115 ka, 120 ka, 125 ka and 130 ka and quantitatively estimated temperature uncertainties that include relative dating errors. This provides an improved benchmark for performing more robust model-data comparison. The surface temperature simulated by two General Circulation Models (CCSM3 and HadCM3) for 130 ka and 125 ka is compared to the corresponding time slice data synthesis. This comparison shows that the models predict warmer than present conditions earlier than documented in the North Atlantic, while neither model is able to produce the reconstructed early Southern Ocean and Antarctic warming. Our results highlight the importance of producing a sequence of time slices rather than one single time slice averaging the LIG climate conditions.

Seafloor magnetic anomalies and ages observed from shipboard and airborne measurements from New Zealand to the Amundsen Sea (south of 64°S)

Geophysical data acquired using R/V Polarstern constrain the structure and age of the rifted oceanic margin of West Antarctica. West of the Antipodes Fracture Zone, the 145 km wide continent-ocean transition zone (COTZ) of the Marie Byrd Land sector resembles a typical magma-poor margin. New gravity and seismic reflection data indicates initial continental crust of thickness 24 km, that was stretched 90 km. Farther east, the Bellingshausen sector is broad and complex with abundant evidence for volcanism, the COTZ is ~670 km wide, and the nature of crust within the COTZ is uncertain. Margin extension is estimated to be 106-304 km in this sector. Seafloor magnetic anomalies adjacent to Marie Byrd Land near the Pahemo Fracture Zone indicate full-spreading rate during c33-c31 (80-68 Myr) of 60 mm/yr, increasing to 74 mm/yr at c27 (62 Myr), and then dropping to 22 mm/yr by c22 (50 Myr). Spreading rates were lower to the west. Extrapolation towards the continental margin indicates initial oceanic crust formation at around c34y (84 Myr). Subsequent motion of the Bellingshausen plate relative to Antarctica (84-62 Myr) took place east of the Antipodes Fracture Zone at rates <40 mm/yr, typically 5-20 mm/yr. The high extension rate of 30-60 mm/yr during initial margin formation is consistent with steep and symmetrical margin morphology, but subsequent motion of the Bellingshausen plate was slow and complex, and modified rift morphology through migrating deformation and volcanic centers to create a broad and complex COTZ.

An experimental study of discharge partitioning and flow structure at symmetrical bifurcations

Recent research has examined the factors controlling the geometrical configuration of bifurcations, determined the range of stability conditions for a number of bifurcation types and assessed the impact of perturbations on bifurcation evolution. However, the flow division process and the parameters that influence flow and sediment partitioning are still poorly characterized. To identify and isolate these parameters, three-dimensional velocities were measured at 11 cross-sections in a fixed-walled experimental bifurcation. Water surface gradients were controlled, and systematically varied, using a weir in each distributary. As may be expected, the steepest distributary conveyed the most discharge ( was dominant) while the mildest distributary conveyed the least discharge ( was subordinate). A zone of water surface super-elevation was co-located with the bifurcation in symmetric cases or displaced into the subordinate branch in asymmetric cases. Downstream of a relatively acute-angled bifurcation, primary velocity cores were near to the water surface and against the inner banks, with near-bed zones of lower primary velocity at the outer banks. Downstream of an obtuse-angled bifurcation, velocity cores were initially at the outer banks, with near-bed zones of lower velocities at the inner banks, but patterns soon reverted to match the acute-angled case. A single secondary flow cell was generated in each distributary, with water flowing inwards at the water surface and outwards at the bed. Circulation was relatively enhanced within the subordinate branch, which may help explain why subordinate distributaries remain open, may play a role in determining the size of commonly-observed topographic features, and may thus exert some control on the stability of asymmetric bifurcations. Further, because larger values of circulation result from larger gradient disadvantages, the length of confluence-diffluence units in braided rivers or between diffluences within delta distributary networks may vary depending upon flow structures inherited from upstream and whether, and how, they are fed by dominant or subordinate distributaries. Copyright (C) 2011 John Wiley & Sons, Ltd.

The relationship between stem-form, stand-closure and site -conditions : the influence of environmental conditions on tree allometry and forest structure in west-central Alberta

Changes in the configuration of a tree stern result insignificant differences in its total volume and in the proportion of that volume that is merchantable timber. Tree allometry, as represented by stem-fo~, is the result of the vertical force of gravity and the horizontal force of wind. The effect of wind force is demonstrated in the relationship between stem-form, standclosure and site-conditions. An increase in wind force on the individual tree due to a decrease in stand density should produce a more tapered tree. The density of the stand is determined by the conditions that the trees are growing under. The ability of the tree to respond to increased wind force may also be a function of these conditions . This stem-form/stand-closure/site-conditions relationship was examined using a pre-existing database from westcentral Alberta. This database consisted of environmental, vegetation, soils and timber data covering a wide range of sites. There were 653 sample trees with 82 variables that formed the basis of the analysis. There were eight tree species consisting of Pinus contorta, Picea mariana, Picea engelmannii x glauca, Abies lasiocarpa, Larix laricina, Populus tremuloides, Betula papyrifera and Populus balsamifera plus a comprehensive all-species data set. As the actual conformation of the stern is very individual, stem-fo~was represented by the diameter at breast height to total height r~tio. The four stand-closure variables, crown closure, total basal area, total volume and total number of stems were reduced to total basal area and total number of stems utilizing a bivariate correlation matrix by species. Site-conditions were subdivided into macro, meso and micro variables and reduced in number 3 using cross-tabulations, bivariate correlation and principal components analysis as screening tools. The stem-fo~/stand-closure relationship was examined using bivariate correlation coefficients for stem-fo~ with total number of stems and stem-fo~ with total basal area. The stem-fo~/site-conditions and the stand-closure/site- conditions relationships were examined using multiple correlation coefficients. The stem-form/stand-closure/site-conditions relationship was examined using multiple correlation coefficients in separate analyses for both total number of stems and total basal area. An increase in stand-closure produced a decrease in stem-form for both total number of stems and total basal area for most species. There was a significant relationship between stem-form and site-conditions and between stand-closure and site-conditions for both total number of stems and total basal area for most species. There was a significant relationship between the stemform and site-conditions, including the stand-closure, for most species; total number of stems was involved independently of the site-conditions in the prediction of stem-form and total basal area was not. Larix laricina and Betula papyrifera were the exceptions to the trends observed with most species. The influence of both stand-closure (total number of stems in particular) and site-conditions (elevation in particular) suggest that forest management practices should include these- ecological parameters in determining appropriate restocking levels.

Tracking and mean structure of easterly waves over the Intra-Americas Sea

Easterly waves (EWs) are prominent features of the intertropical convergence zone (ITCZ), found in both the Atlantic and Pacific during the Northern Hemisphere summer and fall, where they commonly serve as precursors to hurricanes over both basins.Alarge proportion of Atlantic EWs are known to form over Africa, but the origin of EWs over the Caribbean and east Pacific in particular has not been established in detail. In this study reanalyses are used to examine the coherence of the large-scale wave signatures and to obtain track statistics and energy conversion terms for EWs across this region. Regression analysis demonstrates that some EW kinematic structures readily propagate between the Atlantic and east Pacific, with the highest correlations observed across Costa Rica and Panama. Track statistics are consistent with this analysis and suggest that some individual waves are maintained as they pass from the Atlantic into the east Pacific, whereas others are generated locally in the Caribbean and east Pacific. Vortex anomalies associated with the waves are observed on the leeward side of the Sierra Madre, propagating northwestward along the coast, consistent with previous modeling studies of the interactions between zonal flow and EWs with model topography similar to the Sierra Madre. An energetics analysis additionally indicates that the Caribbean low-level jet and its extension into the east Pacific—known as the Papagayo jet—are a source of energy for EWs in the region. Two case studies support these statistics, as well as demonstrate the modulation of EW track and storm development location by the MJO.

Optical characterization of the growth and spatial structure of a substorm onset arc

We present a detailed case study of the characteristics of auroral forms that constitute the first ionospheric signatures of substorm expansion phase onset. Analysis of the optical frequency and along-arc (azimuthal) wave number spectra provides the strongest constraint to date on the potential mechanisms and instabilities in the near-Earth magnetosphere that accompany auroral onset and which precede poleward arc expansion and auroral breakup. We evaluate the frequency and growth rates of the auroral forms as a function of azimuthal wave number to determine whether these wave characteristics are consistent with current models of the substorm onset mechanism. We find that the frequency, spatial scales, and growth rates of the auroral forms are most consistent with the cross-field current instability or a ballooning instability, most likely triggered close to the inner edge of the ion plasma sheet. This result is supportive of a near-Earth plasma sheet initiation of the substorm expansion phase. We also present evidence that the frequency and phase characteristics of the auroral undulations may be generated via resonant processes operating along the geomagnetic field. Our observations provide the most powerful constraint to date on the ionospheric manifestation of the physical processes operating during the first few minutes around auroral substorm onset.

Time-lapse geophysical investigations over known archaeological features using electrical resistivity imaging and earth resistance

Electrical methods of geophysical survey are known to produce results that are hard to predict at different times of the year, and under differing weather conditions. This is a problem which can lead to misinterpretation of archaeological features under investigation. The dynamic relationship between a ‘natural’ soil matrix and an archaeological feature is a complex one, which greatly affects the success of the feature’s detection when using active electrical methods of geophysical survey. This study has monitored the gradual variation of measured resistivity over a selection of study areas. By targeting difficult to find, and often ‘missing’ electrical anomalies of known archaeological features, this study has increased the understanding of both the detection and interpretation capabilities of such geophysical surveys. A 16 month time-lapse study over 4 archaeological features has taken place to investigate the aforementioned detection problem across different soils and environments. In addition to the commonly used Twin-Probe earth resistance survey, electrical resistivity imaging (ERI) and quadrature electro-magnetic induction (EMI) were also utilised to explore the problem. Statistical analyses have provided a novel interpretation, which has yielded new insights into how the detection of archaeological features is influenced by the relationship between the target feature and the surrounding ‘natural’ soils. The study has highlighted both the complexity and previous misconceptions around the predictability of the electrical methods. The analysis has confirmed that each site provides an individual and nuanced situation, the variation clearly relating to the composition of the soils (particularly pore size) and the local weather history. The wide range of reasons behind survey success at each specific study site has been revealed. The outcomes have shown that a simplistic model of seasonality is not universally applicable to the electrical detection of archaeological features. This has led to the development of a method for quantifying survey success, enabling a deeper understanding of the unique way in which each site is affected by the interaction of local environmental and geological conditions.

Pelagic molybdenum concentration anomalies and the impact of sediment resuspension on the molybdenum budget in two tidal systems of the North Sea

The seasonal dynamics of molybdenum (Mo) were studied in the water column of two tidal basins of the German Wadden Sea (Sylt-Rømø and Spiekeroog) between 2007 and 2011. In contrast to its conservative behaviour in the open ocean, both, losses of more than 50% of the usual concentration level of Mo in seawater and enrichments up to 20% were observed repeatedly in the water column of the study areas. During early summer, Mo removal by adsorption on algae-derived organic matter (e.g. after Phaeocystis blooms) is postulated to be a possible mechanism. Mo bound to organic aggregates is likely transferred to the surface sediment where microbial decomposition enriches Mo in the pore water. First δ98/95Mo data of the study area disclose residual Mo in the open water column being isotopically heavier than MOMo (Mean Ocean Molybdenum) during a negative Mo concentration anomaly, whereas suspended particulate matter shows distinctly lighter values. Based on field observations a Mo isotope enrichment factor of ε = −0.3‰ has been determined which was used to argue against sorption on metal oxide surfaces. It is suggested here that isotope fractionation is caused by biological activity and association to organic matter. Pelagic Mo concentration anomalies exceeding the theoretical salinity-based concentration level, on the other hand, cannot be explained by replenishment via North Sea waters alone and require a supply of excess Mo. Laboratory experiments with natural anoxic tidal flat sediments and modelled sediment displacement during storm events suggest fast and effective Mo release during the resuspension of anoxic sediments in oxic seawater as an important process for a recycling of sedimentary sulphide bound Mo into the water column.

Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.