Major element and mineral composition of gabbroic xenoliths and megacrysts from Jeju Island, South Korea

Gabbroic xenoliths and diverse megacrysts (e.g., clinopyroxenes, amphiboles and plagioclases), which correspond to the lithology ranging from gabbro-norite to gabbro, occur in the Pleisto-Holocene alkali basalts from Jeju Island, South Korea. The gabbroic xenoliths consist primarily of moderate-K2O plagioclase, Ti-Al-rich clinopyroxene and CaO-rich orthopyroxene; additionally, TiO2-rich amphibole (kaersutite) and Ti-Fe oxides might or might not be present. The plagioclase is the most dominant phase (approx. 60-70 vol.%). The xenoliths and megacrysts provide evidence for the modal metasomatism of the lower continental crust by the mafic magmas during the Pleistocene. The coarse grain size (up to 5 mm), moderate Mg# [=100xMg/(Mg+Fe(total)) atomic ratio] of pyroxenes (70-77) and textural features (e.g., poikilitic) indicate that the gabbroic xenoliths are consistent with a cumulus origin. The clinopyroxenes from these xenoliths are enriched in REE with smooth convex-upward MREE patterns, which are expected for cumulus minerals formed from a melt enriched in incompatible trace elements. The strikingly similar major and trace element variations and the patterns of constituent minerals clearly indicate a genetic link between the gabbroic xenoliths (plus megacrysts) and the host basalt, indicating that the xenoliths belong to the Jeju Pleisto-Holocene magma system. On the basis of the textural features, the mineral equilibria and the major and trace element variations, the xenoliths appear to have crystallized from basaltic melts at the reservoir-roof environment within the lower crust (4-7 kbars) above the present Moho estimates beneath Jeju Island, where the xenoliths represent wall rocks. Following the consolidation of the xenolith lithologies, volatile- and incompatible element-enriched melt/fluid, as metasomatic agents, infiltrated through the grain boundaries and/or cracks and reacted with the preexisting anhydrous phases, which produced the metasomatic amphiboles. This volatile-enriched melt/fluid could have evolved from the initially anhydrous compositions to the volatile-saturated compositions by the active fractional crystallization in the Jeju Pleisto-Holocene magma system. This process was significant in that it was a relatively young event and played an important role in the formation of the hydrous minerals and the metasomatization of the lower continental crust, which is a plume-impacted area along the Asian continental margin. The major and trace element analyses of the mineral phases from the xenoliths were performed to define the principal geochemical characteristics of the crustal lithosphere segment represented by the studied xenoliths.

Preliminary ichthyolith biostratigraphy and age profiles at DSDP Sites 91-595 and 91-596

Seventy meters of Cenozoic and Mesozoic pelagic clay cored at DSDP Sites 595 and 596 provide the basis for a preliminary analysis of ichthyolith biostratigraphy in the southwest Pacific. A most likely order of the more reliable ichthyolith events is compared with a synthesis of ichthyolith biostratigraphy in the North Pacific and with dated composite ranges. The resultant preliminary ichthyolith stratigraphy suggests that the Cenozoic is represented by the upper 20 m at Site 596 and 16 to 22 m at Site 595. Mixing of taxa precludes a clear recognition of the Cretaceous/Tertiary boundary at Site 595. The occurrence of 13 newly described subtypes is recorded in Mesozoic sediments at Sites 595 and 596. These new subtypes and previously described Mesozoic forms may be useful for recognizing Mesozoic subdivisions when their occurrences in sequences dated by other microfossils are investigated.

Uppermost Jurassic to lower Cretaceous benthic foraminifera from ODP Hole 123-765C

Benthic foraminifers were studied in 99 samples collected from the lower 200 m of Hole 765C. The studied section ranges from the Tithonian to Aptian, and benthic foraminifers can be subdivided into five assemblages on the basis of faunal diversity and stratigraphic ranges of distinctive species. Compared with deep-water assemblages from Atlantic DSDP sites and Poland, assemblages from the Argo Abyssal Plain display a higher diversity of agglutinated forms, which comprise the autochthonous assemblages. Assemblages at the base of Hole 765C are wholly composed of agglutinated forms, reflecting deposition beneath the carbonate compensation depth (CCD). Most calcareous benthic species are found in turbidite layers, and the presence of an upper Valanginian Praedorothia praehauteriviana Assemblage may indicate deposition at or just below the CCD. The P. praehauteriviana Assemblage from Hole 765C is the temporal equivalent of similar assemblages from DSDP Holes 534A, 416A, 370, 105, and 101 in the Atlantic Ocean and Hole 306 in the Pacific Ocean. Stratigraphic ranges of cosmopolitan agglutinated species at Site 765 generally overlap with their reported ranges in the Atlantic and in the bathyal flysch sequences of the Carpathians; however, several species from Hole 765C have not been previously reported from Uppermost Jurassic to Lower Cretaceous abyssal sediments.

Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Cretaceous-Tertiary ichthyoliths of brown clay intervals

The pattern of ichthyolith distribution established in sequences with stratigraphies based on calcareous or siliceous microfossils is used to provide age correlations for three deep-sea pelagic clay intervals that lack the better known microfossils. At Site 637, approximately 25 m of brown clay in Cores 103-637A-21R through 103-637A-23R underlies upper Miocene sediments and is of Paleocene to early Eocene age. At Site 639, 1.7 m of brown clay in Core 103-639C-2R is Eocene to Oligocene. At Site 640, 3.5 m of clay in Cores 103-640A-1R and 103-640A-2R contains a Cretaceous to Paleocene sequence, with the Cretaceous/Tertiary boundary between 84 and 103 cm in Section 103-640A-2R-1.

Model experiment on optimal size of a fish stock with environmental uncertainties

We analyze the effect of environmental uncertainties on optimal fishery management in a bio-economic fishery model. Unlike most of the literature on resource economics, but in line with ecological models, we allow the different biological processes of survival and recruitment to be affected differently by environmental uncertainties. We show that the overall effect of uncertainty on the optimal size of a fish stock is ambiguous, depending on the prudence of the value function. For the case of a risk-neutral fishery manager, the overall effect depends on the relative magnitude of two opposing effects, the 'convex-cost effect' and the 'gambling effect'. We apply the analysis to the Baltic cod and the North Sea herring fisheries, concluding that for risk neutral agents the net effect of environmental uncertainties on the optimal size of these fish stocks is negative, albeit small in absolute value. Under risk aversion, the effect on optimal stock size is positive for sufficiently high coefficients of constant relative risk aversion.

Analyses of ferromanganese micronodules from the central Arctic Ocean

In order to determine geochemical compositions of Late Cenozoic Arctic seawater, oxide fractions were chemically separated from 15 samples of hand-picked ferromanganese micronodules (50-300 mu m). The success of the chemical separation is indicated by the fact that >97% of the Sr in the oxide fraction is seawater-derived. Rare-earth element (REE) abundances of the Arctic micronodule oxide fractions are much lower than those of bulk Fe-Mn nodules from other ocean basins of the world (e.g., 33 vs. 145 ppm Nd), but the Arctic oxides are enriched in Ce relative to Nd (Ce-N/Nd-N=2.2+/-0.5) and have convex-upward, shale-normalized REE patterns (Nd-N/Gd-N=0.61+/-0.06, Gd-N/Yb-N = 1.5+/-0.2, Nd-N/Yb-N = 0.9+/-0.2), typical of other hydrogenous and diagenetic marine Fe-Mn-oxides. Bulk sediment samples from the central Arctic Ocean have REE abundances and patterns that are characteristic of those of post-Archean shale. Non-detrital fractions (calcite + oxide coatings) of Recent Arctic foraminifera have REE abundances and patterns similar to those of Recent foraminifera from the Atlantic Ocean. Electron microprobe analyses (n=178) of transition elements in 29 Arctic Fe-Mn micronodules from five different stratigraphic intervals of Late Cenozoic sediment indicate that oxide accretion occurred as a result of hydrogenetic and diagenetic processes close to the sediment-seawater interface. Transition element ratios suggest that no oxide accretion occurred during transitions from oxic to suboxic diagenetic conditions. Only K is correlated with Si and Al, and ratios of these elements suggest that they are associated with illite or phillipsite. Ca and Mg are correlated with Mn, which indicates variable substitution of these elements from seawater into the manganate phase. The geochemical characteristics of Arctic Fe-Mn micronodules indicate that the REEs of the oxide fractions were ultimately derived from seawater. However, because of minute contributions of Sr from siliciclastic detritus during diagenesis or during the chemical leaching procedure, Sr isotope compositions of the oxide fractions cannot be used to trace temporal changes in the Sr-87/Sr-86 ratio of Arctic seawater or to improve the chronostratigraphy.

Middle Miocene to Quaternary diatoms in the Nankai Trough and Japan Trench

During Leg 87 of the Deep Sea Drilling Project, eleven holes were drilled at Sites 582 and 583 in the Nankai Trough, off Shikoku, southwestern Honshu, and three holes at Site 584 in the Japan Trench, off northeastern Honshu, Japan. In the former area, a low-latitude diatom zone called the Pseudoeunotia doliolus Zone is recognized in thick upper Quaternary sediments, which yield rare but characteristic admixtures of marine planktonic, marine tychopelagic-tobenthic, and nonmarine diatoms. In the latter area, all the sediments recovered contain abundant to common diatoms, allowing recognition of 12 continuous diatom zones from upper Quaternary through lower middle Miocene. Three hiatuses occur in this area around the Pleistocene/Pliocene boundary and in the upper and middle Miocene. In addition, 19 modified diatom zones for a lower Miocene through upper Quaternary interval are proposed. These middle-to-highlatitude zones are numerically coded (NPD1-NPD12) and represent the entire North Pacific. The establishment of these zones is based primarily on Leg 87 data and other DSDP materials and partially on several Japanese subaerial sequences. Correlation of the new zonal framework with previously established frameworks is attempted by the evaluation of operational usefulness of previously used datums. Resting spores of Chaetoceros and its related forms are recorded with specific intent for the first time, and possible ramifications of its frequency variation are presented. Nine new species are proposed: Delphineis sheshukovae Akiba n. sp., Denticulopsis praelauta Akiba and Koizumi n. sp., Kisseleviella ezoensis Akiba n. sp., Nitzschia umaoiensis Akiba n. sp., Thalassiosira jouseae Akiba n. sp., T. praenidulus Akiba n. sp., T. sancettae Akiba n. sp., T. umaoiensis Akiba n. sp., and T. urahoroensis Akiba n. sp. Transfers of systematic positions of the following four taxa are also proposed: Delphineis simonsenii (Mertz) Akiba n. comb., Ikebea tenuis (Brun) Akiba n. comb., Thalassiosira delicata (Barron) Akiba n. stat., and Thalassiothrixrobusta (Schrader) Akiba n. comb.

Rare element concentrations measured on filtered seawater of water bottle samples during POLARSTRERN cruise ANT-XXIV/3 (GIPY5)

The concentrations of rare earth elements (REEs) in seawater display systematic variations related to weathering inputs, particle scavenging and water mass histories. Here we investigate the REE concentrations of water column profiles in the Atlantic sector of the Southern Ocean, a key region of the global circulation and primary production. The data reveal a pronounced contrast between the vertical profiles in the Antarctic Circumpolar Current (ACC) and those to the south of the ACC in the Weddell Gyre (WG). The ACC profiles exhibit the typical increase of REE concentrations with water depth and a change in the shape of the profiles from near linear for the light REEs to more convex for the heavy REEs. In contrast, the WG profiles exhibit high REE concentrations throughout the water column with only the near surface samples showing slightly reduced concentrations indicative of particle scavenging. Seawater normalised REE patterns reveal the strong remineralisation signal in the ACC with the light REEs preferentially removed in surface waters and the mirror image pattern of their preferential release in deep waters. In the WG the patterns are relatively homogenous reflecting the prevalence of well-mixed Lower Circumpolar Deep Water (LCDW) that follows shoaling isopycnals in the region. In the WG particle scavenging of REEs is comparatively small and limited to the summer months by light limitation and winter sea ice cover. Considering the surface water depletion compared to LCDW and that the surface waters of the WG are replaced every few years, the removal rate is estimated to be on the order of 1 nmol/m3/yr for La and Nd. The negative cerium anomalies observed in deep waters are some of the strongest found globally with only the deepest waters in parts of the Pacific having stronger anomalies. These deep waters have been isolated from fresh continental REE inputs during their long journey through the abyssal Indo-Pacific ocean and suggests that the high REE concentrations found in the ACC and WG reflect contributions from old deep waters.

Integrating tracking in concentrating photovoltaics using laterally moving optics

In this work the concept of integrating tracking in concentrating photovoltaics is briefly summarized and possible fields of application are classified. A previously proposed system setup relies on the use of two rotational symmetric laterally moving plano-convex lenses to achieve 500× concentration over an angular range of ±24 ◦ . However, the circular lens apertures are less suitable for application in lens array structures. A new design algorithm based on the Simultaneous Multiple Surface algorithm in three dimensions (SMS3D) demonstrates the ability to address this problem. Performance simulations show that the resulting non-rotational symmetric design outperforms its conventional rotational symmetric counterpart

Distributed collaborative processing in wireless sensor networks with application to target localization and beamforming

Abstract The proliferation of wireless sensor networks and the variety of envisioned applications associated with them has motivated the development of distributed algorithms for collaborative processing over networked systems. One of the applications that has attracted the attention of the researchers is that of target localization where the nodes of the network try to estimate the position of an unknown target that lies within its coverage area. Particularly challenging is the problem of estimating the target’s position when we use received signal strength indicator (RSSI) due to the nonlinear relationship between the measured signal and the true position of the target. Many of the existing approaches suffer either from high computational complexity (e.g., particle filters) or lack of accuracy. Further, many of the proposed solutions are centralized which make their application to a sensor network questionable. Depending on the application at hand and, from a practical perspective it could be convenient to find a balance between localization accuracy and complexity. Into this direction we approach the maximum likelihood location estimation problem by solving a suboptimal (and more tractable) problem. One of the main advantages of the proposed scheme is that it allows for a decentralized implementation using distributed processing tools (e.g., consensus and convex optimization) and therefore, it is very suitable to be implemented in real sensor networks. If further accuracy is needed an additional refinement step could be performed around the found solution. Under the assumption of independent noise among the nodes such local search can be done in a fully distributed way using a distributed version of the Gauss-Newton method based on consensus. Regardless of the underlying application or function of the sensor network it is al¬ways necessary to have a mechanism for data reporting. While some approaches use a special kind of nodes (called sink nodes) for data harvesting and forwarding to the outside world, there are however some scenarios where such an approach is impractical or even impossible to deploy. Further, such sink nodes become a bottleneck in terms of traffic flow and power consumption. To overcome these issues instead of using sink nodes for data reporting one could use collaborative beamforming techniques to forward directly the generated data to a base station or gateway to the outside world. In a dis-tributed environment like a sensor network nodes cooperate in order to form a virtual antenna array that can exploit the benefits of multi-antenna communications. In col-laborative beamforming nodes synchronize their phases in order to add constructively at the receiver. Some of the inconveniences associated with collaborative beamforming techniques is that there is no control over the radiation pattern since it is treated as a random quantity. This may cause interference to other coexisting systems and fast bat-tery depletion at the nodes. Since energy-efficiency is a major design issue we consider the development of a distributed collaborative beamforming scheme that maximizes the network lifetime while meeting some quality of service (QoS) requirement at the re¬ceiver side. Using local information about battery status and channel conditions we find distributed algorithms that converge to the optimal centralized beamformer. While in the first part we consider only battery depletion due to communications beamforming, we extend the model to account for more realistic scenarios by the introduction of an additional random energy consumption. It is shown how the new problem generalizes the original one and under which conditions it is easily solvable. By formulating the problem under the energy-efficiency perspective the network’s lifetime is significantly improved. Resumen La proliferación de las redes inalámbricas de sensores junto con la gran variedad de posi¬bles aplicaciones relacionadas, han motivado el desarrollo de herramientas y algoritmos necesarios para el procesado cooperativo en sistemas distribuidos. Una de las aplicaciones que suscitado mayor interés entre la comunidad científica es la de localization, donde el conjunto de nodos de la red intenta estimar la posición de un blanco localizado dentro de su área de cobertura. El problema de la localization es especialmente desafiante cuando se usan niveles de energía de la seal recibida (RSSI por sus siglas en inglés) como medida para la localization. El principal inconveniente reside en el hecho que el nivel de señal recibida no sigue una relación lineal con la posición del blanco. Muchas de las soluciones actuales al problema de localization usando RSSI se basan en complejos esquemas centralizados como filtros de partículas, mientas que en otras se basan en esquemas mucho más simples pero con menor precisión. Además, en muchos casos las estrategias son centralizadas lo que resulta poco prácticos para su implementación en redes de sensores. Desde un punto de vista práctico y de implementation, es conveniente, para ciertos escenarios y aplicaciones, el desarrollo de alternativas que ofrezcan un compromiso entre complejidad y precisión. En esta línea, en lugar de abordar directamente el problema de la estimación de la posición del blanco bajo el criterio de máxima verosimilitud, proponemos usar una formulación subóptima del problema más manejable analíticamente y que ofrece la ventaja de permitir en¬contrar la solución al problema de localization de una forma totalmente distribuida, convirtiéndola así en una solución atractiva dentro del contexto de redes inalámbricas de sensores. Para ello, se usan herramientas de procesado distribuido como los algorit¬mos de consenso y de optimización convexa en sistemas distribuidos. Para aplicaciones donde se requiera de un mayor grado de precisión se propone una estrategia que con¬siste en la optimización local de la función de verosimilitud entorno a la estimación inicialmente obtenida. Esta optimización se puede realizar de forma descentralizada usando una versión basada en consenso del método de Gauss-Newton siempre y cuando asumamos independencia de los ruidos de medida en los diferentes nodos. Independientemente de la aplicación subyacente de la red de sensores, es necesario tener un mecanismo que permita recopilar los datos provenientes de la red de sensores. Una forma de hacerlo es mediante el uso de uno o varios nodos especiales, llamados nodos “sumidero”, (sink en inglés) que actúen como centros recolectores de información y que estarán equipados con hardware adicional que les permita la interacción con el exterior de la red. La principal desventaja de esta estrategia es que dichos nodos se convierten en cuellos de botella en cuanto a tráfico y capacidad de cálculo. Como alter¬nativa se pueden usar técnicas cooperativas de conformación de haz (beamforming en inglés) de manera que el conjunto de la red puede verse como un único sistema virtual de múltiples antenas y, por tanto, que exploten los beneficios que ofrecen las comu¬nicaciones con múltiples antenas. Para ello, los distintos nodos de la red sincronizan sus transmisiones de manera que se produce una interferencia constructiva en el recep¬tor. No obstante, las actuales técnicas se basan en resultados promedios y asintóticos, cuando el número de nodos es muy grande. Para una configuración específica se pierde el control sobre el diagrama de radiación causando posibles interferencias sobre sis¬temas coexistentes o gastando más potencia de la requerida. La eficiencia energética es una cuestión capital en las redes inalámbricas de sensores ya que los nodos están equipados con baterías. Es por tanto muy importante preservar la batería evitando cambios innecesarios y el consecuente aumento de costes. Bajo estas consideraciones, se propone un esquema de conformación de haz que maximice el tiempo de vida útil de la red, entendiendo como tal el máximo tiempo que la red puede estar operativa garantizando unos requisitos de calidad de servicio (QoS por sus siglas en inglés) que permitan una decodificación fiable de la señal recibida en la estación base. Se proponen además algoritmos distribuidos que convergen a la solución centralizada. Inicialmente se considera que la única causa de consumo energético se debe a las comunicaciones con la estación base. Este modelo de consumo energético es modificado para tener en cuenta otras formas de consumo de energía derivadas de procesos inherentes al funcionamiento de la red como la adquisición y procesado de datos, las comunicaciones locales entre nodos, etc. Dicho consumo adicional de energía se modela como una variable aleatoria en cada nodo. Se cambia por tanto, a un escenario probabilístico que generaliza el caso determinista y se proporcionan condiciones bajo las cuales el problema se puede resolver de forma eficiente. Se demuestra que el tiempo de vida de la red mejora de forma significativa usando el criterio propuesto de eficiencia energética.

Fundamentación del conjunto tractor-apero para la labor de surcado y su influencia en el consumo energético en cultivos de raíces y tubérculos

El cumplimiento de los objetivos agrotécnicos de los cultivos de raíces y tubérculos exige que la labor de surcado configure un caballón que permita el trabajo de los posteriores aperos que participan en las demás labores de cultivo. En la agricultura cubana se emplean diversos tipos de surcadores que difieren en las características de los elementos aporcadores. Este hecho afecta a la capacidad de trabajo del conjunto tractor-apero y por extensión en el consumo de combustible y el rendimiento del cultivo. El objetivo de este trabajo ha sido fundamentar el empleo del apero más adecuado para la labor de surcado en cultivos de raíces y tubérculos (patata, boniato, yuca) en Fluvisoles, teniendo en cuenta el comportamiento de las dimensiones y forma del caballón, cualidades tecnológicas explotativas del conjunto tractor-apero, consumo energético del proceso productivo y el rendimiento agrícola. Los tratamientos considerados estuvieron determinados por la utilización del surcador fertilizador IIMA-BALDAN 4.500 (T1), y los surcadores aporcadores SA-5 (T2) y SA-3 (T3). En los tres tratamientos se alcanzó un perfil de la sección transversal del caballón de configuración cóncavo-convexo. El surcador BALDAN logró el mejor comportamiento de los indicadores tecnológico-explotativos lo que se tradujo en la mayor capacidad de trabajo de 1,36 ha h-1 y el menor consumo de combustible de 5,38 L ha-1. Este mismo tratamiento proporcionó el mayor rendimiento medio de los cultivos superando en un 43,2% al del surcador SA-5 y en un 57,9% al del surcador SA-3. Desde el punto de vista energético T1, obtuvo una mejor eficiencia energética superior en un 43,5% a la del tratamiento T2, y en un 57,1% a la del tratamiento T3. ABSTRACT The achievement of agrotechnical objectives in the cultivation of roots and tubers demands the plow through activity creates a ridge to allow the work of subsequent farming tools participating in other cultivation works. In Cuban Agriculture various kinds of furrowers are used which differ in the characteristics of their earth up elements. This fact affects work capacity of the tractor-farming tool combination and by extension to fuel consumption and cultivation yield. The aim of this work has been to support the use of the most accurate farming tool in the plow through activity for the cultivation of roots and tubers (potato, sweet potato, cassava) in Fluvisoils, taking into account the behaviour of the ridge´s dimensions and shape, technologic-explotative qualities of the tractor-farming tool combination, energy consumption of the productive process and agricultural yield. The treatments considered were determined by the use of fertlizing furrower IIMA-BALDAN 4.500 (T1), and the earth up furrowers SA-5 (T2) and SA-3 (T3). During the three treatments a transverse section profile of the ridge under concave-convex configuration was obtained. The furrower BALDAN achieved the best behaviour of the technologic-explotative indicators, which was translated into the largest work capacity of 1,36 ha h-1 and the smallest fuel consumption of 5,38 L ha-1. This treatment also provided the biggest average yield of the cultivation, overcoming in a 43,2% the furrower SA-5 and furrower SA-3 in a 57,9%. From the energetic point of view T1 gained better energetic efficiency than treatment T2 in a 43,5% and treatment T3 in a 57,1%.

Aplicaciones de interés forestal de las imágenes obtenidas con el sensor Multi-angle Imaging SpectroRadiometer (MISR)

Las aplicaciones de la teledetección al seguimiento de lo que ocurre en la superficie terrestre se han ido multiplicando y afinando con el lanzamiento de nuevos sensores por parte de las diferentes agencias espaciales. La necesidad de tener información actualizada cada poco tiempo y espacialmente homogénea, ha provocado el desarrollo de nuevos programas como el Earth Observing System (EOS) de la National Aeronautics and Space Administration (NASA). Uno de los sensores que incorpora el buque insignia de ese programa, el satélite TERRA, es el Multi-angle Imaging SpectroRadiometer (MISR), diseñado para capturar información multiangular de la superficie terrestre. Ya desde los años 1970, se conocía que la reflectancia de las diversas ocupaciones y usos del suelo variaba en función del ángulo de observación y de iluminación, es decir, que eran anisotrópicas. Tal variación estaba además relacionada con la estructura tridimensional de tales ocupaciones, por lo que se podía aprovechar tal relación para obtener información de esa estructura, más allá de la que pudiera proporcionar la información meramente espectral. El sensor MISR incorpora 9 cámaras a diferentes ángulos para capturar 9 imágenes casi simultáneas del mismo punto, lo que permite estimar con relativa fiabilidad la respuesta anisotrópica de la superficie terrestre. Varios trabajos han demostrado que se pueden estimar variables relacionadas con la estructura de la vegetación con la información que proporciona MISR. En esta Tesis se ha realizado una primera aplicación a la Península Ibérica, para comprobar su utilidad a la hora de estimar variables de interés forestal. En un primer paso se ha analizado la variabilidad temporal que se produce en los datos, debido a los cambios en la geometría de captación, es decir, debido a la posición relativa de sensores y fuente de iluminación, que en este caso es el Sol. Se ha comprobado cómo la anisotropía es mayor desde finales de otoño hasta principios de primavera debido a que la posición del Sol es más cercana al plano de los sensores. También se ha comprobado que los valores máximo y mínimo se van desplazando temporalmente entre el centro y el extremo angular. En la caracterización multiangular de ocupaciones del suelo de CORINE Land Cover que se ha realizado, se puede observar cómo la forma predominante en las imágenes con el Sol más alto es convexa con un máximo en la cámara más cercana a la fuente de iluminación. Sin embargo, cuando el Sol se encuentra mucho más bajo, ese máximo es muy externo. Por otra parte, los datos obtenidos en verano son mucho más variables para cada ocupación que los de noviembre, posiblemente debido al aumento proporcional de las zonas en sombra. Para comprobar si la información multiangular tiene algún efecto en la obtención de imágenes clasificadas según ocupación y usos del suelo, se han realizado una serie de clasificaciones variando la información utilizada, desde sólo multiespectral, a multiangular y multiespectral. Los resultados muestran que, mientras para las clasificaciones más genéricas la información multiangular proporciona los peores resultados, a medida que se amplían el número de clases a obtener tal información mejora a lo obtenido únicamente con información multiespectral. Por otra parte, se ha realizado una estimación de variables cuantitativas como la fracción de cabida cubierta (Fcc) y la altura de la vegetación a partir de información proporcionada por MISR a diferentes resoluciones. En el valle de Alcudia (Ciudad Real) se ha estimado la fracción de cabida cubierta del arbolado para un píxel de 275 m utilizando redes neuronales. Los resultados muestran que utilizar información multiespectral y multiangular puede mejorar casi un 20% las estimaciones realizadas sólo con datos multiespectrales. Además, las relaciones obtenidas llegan al 0,7 de R con errores inferiores a un 10% en Fcc, siendo éstos mucho mejores que los obtenidos con el producto elaborado a partir de datos multiespectrales del sensor Moderate Resolution Imaging Spectroradiometer (MODIS), también a bordo de Terra, para la misma variable. Por último, se ha estimado la fracción de cabida cubierta y la altura efectiva de la vegetación para 700.000 ha de la provincia de Murcia, con una resolución de 1.100 m. Los resultados muestran la relación existente entre los datos espectrales y los multiangulares, obteniéndose coeficientes de Spearman del orden de 0,8 en el caso de la fracción de cabida cubierta de la vegetación, y de 0,4 en el caso de la altura efectiva. Las estimaciones de ambas variables con redes neuronales y diversas combinaciones de datos, arrojan resultados con R superiores a 0,85 para el caso del grado de cubierta vegetal, y 0,6 para la altura efectiva. Los parámetros multiangulares proporcionados en los productos elaborados con MISR a 1.100 m, no obtienen buenos resultados por sí mismos pero producen cierta mejora al incorporarlos a la información espectral. Los errores cuadráticos medios obtenidos son inferiores a 0,016 para la Fcc de la vegetación en tanto por uno, y 0,7 m para la altura efectiva de la misma. Regresiones geográficamente ponderadas muestran además que localmente se pueden obtener mejores resultados aún mejores, especialmente cuando hay una mayor variabilidad espacial de las variables estimadas. En resumen, la utilización de los datos proporcionados por MISR ofrece una prometedora vía de mejora de resultados en la media-baja resolución, tanto para la clasificación de imágenes como para la obtención de variables cuantitativas de la estructura de la vegetación. ABSTRACT Applications of remote sensing for monitoring what is happening on the land surface have been multiplied and refined with the launch of new sensors by different Space Agencies. The need of having up to date and spatially homogeneous data, has led to the development of new programs such as the Earth Observing System (EOS) of the National Aeronautics and Space Administration (NASA). One of the sensors incorporating the flagship of that program, the TERRA satellite, is Multi-angle Imaging Spectroradiometer (MISR), designed to capture the multi-angle information of the Earth's surface. Since the 1970s, it was known that the reflectance of various land covers and land uses varied depending on the viewing and ilumination angles, so they are anisotropic. Such variation was also related to the three dimensional structure of such covers, so that one could take advantage of such a relationship to obtain information from that structure, beyond which spectral information could provide. The MISR sensor incorporates 9 cameras at different angles to capture 9 almost simultaneous images of the same point, allowing relatively reliable estimates of the anisotropic response of the Earth's surface. Several studies have shown that we can estimate variables related to the vegetation structure with the information provided by this sensor, so this thesis has made an initial application to the Iberian Peninsula, to check their usefulness in estimating forest variables of interest. In a first step we analyzed the temporal variability that occurs in the data, due to the changes in the acquisition geometry, i.e. the relative position of sensor and light source, which in this case is the Sun. It has been found that the anisotropy is greater from late fall through early spring due to the Sun's position closer to the plane of the sensors. It was also found that the maximum and minimum values are displaced temporarily between the center and the ends. In characterizing CORINE Land Covers that has been done, one could see how the predominant form in the images with the highest sun is convex with a maximum in the camera closer to the light source. However, when the sun is much lower, the maximum is external. Moreover, the data obtained for each land cover are much more variable in summer that in November, possibly due to the proportional increase in shadow areas. To check whether the information has any effect on multi-angle imaging classification of land cover and land use, a series of classifications have been produced changing the data used, from only multispectrally, to multi-angle and multispectral. The results show that while for the most generic classifications multi-angle information is the worst, as there are extended the number of classes to obtain such information it improves the results. On the other hand, an estimate was made of quantitative variables such as canopy cover and vegetation height using information provided by MISR at different resolutions. In the valley of Alcudia (Ciudad Real), we estimated the canopy cover of trees for a pixel of 275 m by using neural networks. The results showed that using multispectral and multiangle information can improve by almost 20% the estimates that only used multispectral data. Furthermore, the relationships obtained reached an R coefficient of 0.7 with errors below 10% in canopy cover, which is much better result than the one obtained using data from the Moderate Resolution Imaging Spectroradiometer (MODIS), also onboard Terra, for the same variable. Finally we estimated the canopy cover and the effective height of the vegetation for 700,000 hectares in the province of Murcia, with a spatial resolution of 1,100 m. The results show a relationship between the spectral and the multi-angle data, and provide estimates of the canopy cover with a Spearman’s coefficient of 0.8 in the case of the vegetation canopy cover, and 0.4 in the case of the effective height. The estimates of both variables using neural networks and various combinations of data, yield results with an R coefficient greater than 0.85 for the case of the canopy cover, and 0.6 for the effective height. Multi-angle parameters provided in the products made from MISR at 1,100 m pixel size, did not produce good results from themselves but improved the results when included to the spectral information. The mean square errors were less than 0.016 for the canopy cover, and 0.7 m for the effective height. Geographically weighted regressions also showed that locally we can have even better results, especially when there is high spatial variability of estimated variables. In summary, the use of the data provided by MISR offers a promising way of improving remote sensing performance in the low-medium spatial resolution, both for image classification and for the estimation of quantitative variables of the vegetation structure.

El problema del rozamiento en el análisis de estructuras de fábrica mediante modelos de sólidos rígidos

La evaluación de la seguridad de estructuras antiguas de fábrica es un problema abierto.El material es heterogéneo y anisótropo, el estado previo de tensiones difícil de conocer y las condiciones de contorno inciertas. A comienzos de los años 50 se demostró que el análisis límite era aplicable a este tipo de estructuras, considerándose desde entonces como una herramienta adecuada. En los casos en los que no se produce deslizamiento la aplicación de los teoremas del análisis límite estándar constituye una herramienta formidable por su simplicidad y robustez. No es necesario conocer el estado real de tensiones. Basta con encontrar cualquier solución de equilibrio, y que satisfaga las condiciones de límite del material, en la seguridad de que su carga será igual o inferior a la carga real de inicio de colapso. Además esta carga de inicio de colapso es única (teorema de la unicidad) y se puede obtener como el óptimo de uno cualquiera entre un par de programas matemáticos convexos duales. Sin embargo, cuando puedan existir mecanismos de inicio de colapso que impliquen deslizamientos, cualquier solución debe satisfacer tanto las restricciones estáticas como las cinemáticas, así como un tipo especial de restricciones disyuntivas que ligan las anteriores y que pueden plantearse como de complementariedad. En este último caso no está asegurada la existencia de una solución única, por lo que es necesaria la búsqueda de otros métodos para tratar la incertidumbre asociada a su multiplicidad. En los últimos años, la investigación se ha centrado en la búsqueda de un mínimo absoluto por debajo del cual el colapso sea imposible. Este método es fácil de plantear desde el punto de vista matemático, pero intratable computacionalmente, debido a las restricciones de complementariedad 0 y z 0 que no son ni convexas ni suaves. El problema de decisión resultante es de complejidad computacional No determinista Polinomial (NP)- completo y el problema de optimización global NP-difícil. A pesar de ello, obtener una solución (sin garantía de exito) es un problema asequible. La presente tesis propone resolver el problema mediante Programación Lineal Secuencial, aprovechando las especiales características de las restricciones de complementariedad, que escritas en forma bilineal son del tipo y z = 0; y 0; z 0 , y aprovechando que el error de complementariedad (en forma bilineal) es una función de penalización exacta. Pero cuando se trata de encontrar la peor solución, el problema de optimización global equivalente es intratable (NP-difícil). Además, en tanto no se demuestre la existencia de un principio de máximo o mínimo, existe la duda de que el esfuerzo empleado en aproximar este mínimo esté justificado. En el capítulo 5, se propone hallar la distribución de frecuencias del factor de carga, para todas las soluciones de inicio de colapso posibles, sobre un sencillo ejemplo. Para ello, se realiza un muestreo de soluciones mediante el método de Monte Carlo, utilizando como contraste un método exacto de computación de politopos. El objetivo final es plantear hasta que punto está justificada la busqueda del mínimo absoluto y proponer un método alternativo de evaluación de la seguridad basado en probabilidades. Las distribuciones de frecuencias, de los factores de carga correspondientes a las soluciones de inicio de colapso obtenidas para el caso estudiado, muestran que tanto el valor máximo como el mínimo de los factores de carga son muy infrecuentes, y tanto más, cuanto más perfecto y contínuo es el contacto. Los resultados obtenidos confirman el interés de desarrollar nuevos métodos probabilistas. En el capítulo 6, se propone un método de este tipo basado en la obtención de múltiples soluciones, desde puntos de partida aleatorios y calificando los resultados mediante la Estadística de Orden. El propósito es determinar la probabilidad de inicio de colapso para cada solución.El método se aplica (de acuerdo a la reducción de expectativas propuesta por la Optimización Ordinal) para obtener una solución que se encuentre en un porcentaje determinado de las peores. Finalmente, en el capítulo 7, se proponen métodos híbridos, incorporando metaheurísticas, para los casos en que la búsqueda del mínimo global esté justificada. Abstract Safety assessment of the historic masonry structures is an open problem. The material is heterogeneous and anisotropic, the previous state of stress is hard to know and the boundary conditions are uncertain. In the early 50's it was proven that limit analysis was applicable to this kind of structures, being considered a suitable tool since then. In cases where no slip occurs, the application of the standard limit analysis theorems constitutes an excellent tool due to its simplicity and robustness. It is enough find any equilibrium solution which satisfy the limit constraints of the material. As we are certain that this load will be equal to or less than the actual load of the onset of collapse, it is not necessary to know the actual stresses state. Furthermore this load for the onset of collapse is unique (uniqueness theorem), and it can be obtained as the optimal from any of two mathematical convex duals programs However, if the mechanisms of the onset of collapse involve sliding, any solution must satisfy both static and kinematic constraints, and also a special kind of disjunctive constraints linking the previous ones, which can be formulated as complementarity constraints. In the latter case, it is not guaranted the existence of a single solution, so it is necessary to look for other ways to treat the uncertainty associated with its multiplicity. In recent years, research has been focused on finding an absolute minimum below which collapse is impossible. This method is easy to set from a mathematical point of view, but computationally intractable. This is due to the complementarity constraints 0 y z 0 , which are neither convex nor smooth. The computational complexity of the resulting decision problem is "Not-deterministic Polynomialcomplete" (NP-complete), and the corresponding global optimization problem is NP-hard. However, obtaining a solution (success is not guaranteed) is an affordable problem. This thesis proposes solve that problem through Successive Linear Programming: taking advantage of the special characteristics of complementarity constraints, which written in bilinear form are y z = 0; y 0; z 0 ; and taking advantage of the fact that the complementarity error (bilinear form) is an exact penalty function. But when it comes to finding the worst solution, the (equivalent) global optimization problem is intractable (NP-hard). Furthermore, until a minimum or maximum principle is not demonstrated, it is questionable that the effort expended in approximating this minimum is justified. XIV In chapter 5, it is proposed find the frequency distribution of the load factor, for all possible solutions of the onset of collapse, on a simple example. For this purpose, a Monte Carlo sampling of solutions is performed using a contrast method "exact computation of polytopes". The ultimate goal is to determine to which extent the search of the global minimum is justified, and to propose an alternative approach to safety assessment based on probabilities. The frequency distributions for the case study show that both the maximum and the minimum load factors are very infrequent, especially when the contact gets more perfect and more continuous. The results indicates the interest of developing new probabilistic methods. In Chapter 6, is proposed a method based on multiple solutions obtained from random starting points, and qualifying the results through Order Statistics. The purpose is to determine the probability for each solution of the onset of collapse. The method is applied (according to expectations reduction given by the Ordinal Optimization) to obtain a solution that is in a certain percentage of the worst. Finally, in Chapter 7, hybrid methods incorporating metaheuristics are proposed for cases in which the search for the global minimum is justified.

A practical approach for outdoors distributed target localization in wireless sensor networks

Wireless sensor networks are posed as the new communication paradigm where the use of small, low-complexity, and low-power devices is preferred over costly centralized systems. The spectra of potential applications of sensor networks is very wide, ranging from monitoring, surveillance, and localization, among others. Localization is a key application in sensor networks and the use of simple, efficient, and distributed algorithms is of paramount practical importance. Combining convex optimization tools with consensus algorithms we propose a distributed localization algorithm for scenarios where received signal strength indicator readings are used. We approach the localization problem by formulating an alternative problem that uses distance estimates locally computed at each node. The formulated problem is solved by a relaxed version using semidefinite relaxation technique. Conditions under which the relaxed problem yields to the same solution as the original problem are given and a distributed consensusbased implementation of the algorithm is proposed based on an augmented Lagrangian approach and primaldual decomposition methods. Although suboptimal, the proposed approach is very suitable for its implementation in real sensor networks, i.e., it is scalable, robust against node failures and requires only local communication among neighboring nodes. Simulation results show that running an additional local search around the found solution can yield performance close to the maximum likelihood estimate.