Trophic modeling of the Northern Humboldt Current Ecosystem, Part II: Elucidating ecosystem dynamics from 1995 to 2004 with a focus on the impact of ENSO

The Northern Humboldt Current Ecosystem is one of the most productive in the world in terms of fish production. Its location near to the equator permits strong upwelling under relatively low winds, thus creating optimal conditions for the development of plankton communities. These communities ultimately support abundant populations of grazing fish such as the Peruvian anchoveta, Engraulis ringens. The ecosystem is also subject to strong inter-annual environmental variability associated with the El Niño Southern Oscillation (ENSO), which has major effects on nutrient structure, primary production, and higher trophic levels. Here our objective is to model the contributions of several external drivers (i.e. reconstructed phytoplankton changes, fish immigration, and fishing rate) and internal control mechanisms (i.e. predator-prey) to ecosystem dynamics over an ENSO cycle. Steady-state models and time-series data from the Instituto del Mar del Perú (IMARPE) from 1995 to 2004 provide the base data for simulations conducted with the program Ecopath with Ecosim. In simulations all three external drivers contribute to ecosystem dynamics. Changes in phytoplankton quantity and composition (i.e. contribution of diatoms and dino- and silicoflagellates), as affected by upwelling intensity, were important in dynamics of the El Niño of 1997–98 and the subsequent 3 years. The expansion and immigration of mesopelagic fish populations during El Niño was important for dynamics in following years. Fishing rate changes were the most important of the three external drivers tested, helping to explain observed dynamics throughout the modeled period, and particularly during the post-El Niño period. Internal control settings show a mix of predator–prey control settings; however a “wasp-waist” control of the ecosystem by small pelagic fish is not supported.     

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground-penetrating radar

Across the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large-scale storms, while providing important habitats and recreational opportunities. Understanding their long-term evolution is vital in managing their condition in a changing climate. Recently, ground-penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre-scale resolution images of near-surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo-environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground-truth control. Nonetheless, GPR data interpretation can be non-intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi-automated GPR feature-extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo-landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty-organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR-interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.     

An overview of the Gulf of Batabanó (Cuba): Environmental features as revealed by surface sediment characterisation

The main environmental features of the Gulf of Batabanó, Cuba, Caribbean, were investigated through the analyses of surface sediments collected at 23 sites. In order to highlight the potential threats affecting the sedimentary compartment of this area, samples were analysed for: granulometry, mineralogy, heavy metals concentration (As, Cd, Cu, Mn, Ni, Pb, Zn), organic carbon, total nitrogen and radionuclides. Findings were compared with published data and “grey” literature.Results showed: granulometric homogeneity and a widespread carbonatic condition all over the gulf, probably due to stable bathymetry and lack of terrigenous input (except for the La Coloma basin); a rather pristine environment for what concerns heavy metals pollution, except for La Coloma where a large arsenic input was recorded; very low levels of natural and artificial radioactivity; a relevant quantity of sedimentary organic matter, providing biota with useful substrate for feeding and enhancing the food-web development while indirectly supplying lobster fisheries. Combined data highlighted the impact of the Dique Sur in reducing terrigenous input in the coastal area.Future studies should focus on dating of sediment cores for identifying and quantifying the changes acting in the gulf and on investigating the origins of the large arsenic input to La Coloma.     

Taphonomy of guanaco bones in Tierra del Fuego

Guanaco carcasses are deposited in great quantities in Cabo San Pablo, Tierra del Fuego (Argentina), as a result of winter stress. Taphonomic studies indicate that the gnawing action of foxes on guanaco (Lama guanicoe) carcasses produces only very tenuous marks on the bones. Lack of sustained interest in the carcasses by carnivores results in slow disarticulation. The articulated and disarticulated bones are exposed to heavy trampling by guanacos, a process that produces vertical migration of small/dense bones and fracturing of the most weathered bones. An understanding of this ongoing process is important for local archaeology, since modern bones are migrating into archaeological contexts. A regional approach to taphonomy is the most appropriate instrument to solve this and other related problems.     

Polar motion for an elastic Earth model with a homogeneous liquid core using a canonical theory

The rotational motion for an elastic Earth model with a homogeneous liquid core has been obtained using Hamilton's equations. From the canonical equations for the precessional and nutational motions in an inertial frame, the corresponding equations in an Earth fixed frame are deduced. The linearized equations obtained for polar motion and liquid core motion are equivalent to the Sasao-Okubo-Saito equations.     

Kriging of water levels in the Souss aquifer,Morocco

Universal kriging is applied to water table data from the Souss aquifer in central Morocco. The procedure accounts for the spatial variability of the phenomenon to be mapped. With the use of measured elevations of the water table, an experimental variogram is constructed that characterizes the spatial variability of the measured water levels. Spherical and Gaussian variogram models are alternatively used to fit the experimental variogram. The models are used to develop contour maps of water table elevations and corresponding estimation variances. The estimation variances express the reliability of the kriged water table elevation maps. Universal kriging also provides a contour map of the expected elevation of the water table (drift). The differences between the expected and measured water table elevations are called residuals from the drift. Residuals from the drift are compared with residuals obtained by more traditional least-squares analysis.     

Dynamic of the freshwater–saltwater interface in a karstic aquifer under extraordinary recharge action: the Merida Yucatan case study

A hydraulic analysis of the interface between freshwater–saltwater behavior was done in the Merida Yucatan zone, two machines that constantly register the groundwater levels were installed, and three electric conductivity logs were taken from wells. When comparing the measured results with the ones obtained using theoretical equations developed to calculate the freshwater–saltwater interface position, it was proved that in some cases these equations can be applied, and in others not. Two effects that rule the behavior of karst aquifers in extraordinary conditions were found.     

On solution of the inverse problem for confined aquifer flow via maximum likelihood

Joint estimation of transmissivity (T) and storativity (S) in a confined aquifer is done via maximum likelihood (ML). The differential equation of groundwater flow is discretized by the finite-element method, leading to equation _t+x _t=u _t. Elements of matrices and , as well as estimated covariance matrix of noise termu _t, are functions of T and S. By minimizing the negative loglikelihood function corresponding to discretized groundwater flow equation with respect to T and S, ML estimators are obtained. The ML approach is found to yield accurate estimates of T and S (within 9 and 10% of their actual values, respectively) and showed quadratic convergence in Newton's search technique. Prediction of aquifer response, using ML estimators, results in estimated piezometric heads accurate to ±0.5 m from their actual, exact values. Statistical properties of ML estimators are derived and some basic results for statistical inference are given.     

Mode Identification from Line Profiles using the Direct Fitting Technique

The method of moments and the direct fitting method are the onlyspectroscopic methods of mode identification which allow a determination ofall pulsational parameters. The pulsation parameters are required to predictthe light amplitude and phase which can be important discriminants in modeidentification. The direct fitting method has several advantages over themethod of moments. It is not restricted to low spherical harmonic degree or form of the eigenfunction and is not sensitive to the placement of the continuum. In the last few years the method has been applied to several different types of stars. We briefly describe the method and give someexamples of its application.