Is China's fifth-largest inland lake to dry-up? Incorporated hydrological and satellite-based methods for forecasting Hulun Lake water levels

Hulun Lake, China's fifth-largest inland lake, experienced severe declines in water level in the period of 2000-2010. This has prompted concerns whether the lake is drying up gradually. A multi-million US dollar engineering project to construct a water channel to transfer part of the river flow from a nearby river to maintain the water level was completed in August 2010. This study aimed to advance the understanding of the key processes controlling the lake water level variation over the last five decades, as well as investigate the impact of the river transfer engineering project on the water level. A water balance model was developed to investigate the lake water level variations over the last five decades, using hydrological and climatic data as well as satellite-based measurements and results from land surface modelling. The investigation reveals that the severe reduction of river discharge (-364±64 mm/yr, ∼70% of the five-decade average) into the lake was the key factor behind the decline of the lake water level between 2000 and 2010. The decline of river discharge was due to the reduction of total runoff from the lake watershed. This was a result of the reduction of soil moisture due to the decrease of precipitation (-49±45 mm/yr) over this period. The water budget calculation suggests that the groundwater component from the surrounding lake area as well as surface run off from the un-gauged area surrounding the lake contributed ∼ net 210 Mm3/yr (equivalent to ∼ 100 mm/yr) water inflows into the lake. The results also show that the water diversion project did prevent a further water level decline of over 0.5 m by the end of 2012. Overall, the monthly water balance model gave an excellent prediction of the lake water level fluctuation over the last five decades and can be a useful tool to manage lake water resources in the future.

Fluxes of organic contaminants from the river catchment into, through and out of the Humber Estuary, UK

This paper summarises the work done on the distribution and reactivity of organic contaminants (simazine, atrazine, lindane, fluoranthene, pyrene, PCB 77, PCB 118) in the Humber Estuary and associated major rivers, as part of the LOIS programme. The preliminary flux calculations show that the most important contributors of selected organic contaminants were the rivers Trent (45% of simazine, 20% of atrazine), Aire (30% of simazine and 33% of atrazine), Don (36 and 37% of fluoranthene and pyrene) and Ouse (18% of fluoranthene and pyrene). For lindane and PCBs, the Aire and Ouse were the key sources. The water flow in all the rivers shows strong seasonal variations, as do the contaminant concentrations. As a result, the mean daily fluxes of these contaminants displayed a strong seasonality. Annual mean concentrations of simazine and atrazine decreased by more than 50% over the period 1994-1995 in most of the rivers, probably as a result of their restricted use in the UK. Mass balance calculations show that the Humber is a sink for atrazine, lindane, PCB 77 and PCB 118, although the degree of removal is generally much lower for atrazine and lindane than for PCB 77 and PCB 118. Mass balance results also show that the Humber can either be a source of fluoranthene and pyrene (in the suspended particulate phase), or a sink (in the dissolved phase), although overall the Humber acts as sink. The budget exercise represents an attempt to quantify the input and output of selected organic contaminants from catchment to ocean. However, due to limited data and assumptions involved in calculations, the estimates should be considered as an order of magnitude approximation. Further improvement both in resolution and accuracy is required.

The white dwarf's carbon fraction as a secondary parameter of Type Ia supernovae

Context. Binary stellar evolution calculations predict thatChandrasekhar-mass carbon/oxygen white dwarfs (WDs) show a radiallyvarying profile for the composition with a carbon depleted core. Manyrecent multi-dimensional simulations of Type Ia supernovae (SNe Ia),however, assume the progenitor WD has a homogeneous chemicalcomposition.
Aims: In this work, we explore the impact ofdifferent initial carbon profiles of the progenitor WD on the explosionphase and on synthetic observables in the Chandrasekhar-mass delayeddetonation model. Spectra and light curves are compared to observationsto judge the validity of the model.
Methods: The explosion phaseis simulated using the finite volume supernova code Leafs, which isextended to treat different compositions of the progenitor WD. Thesynthetic observables are computed with the Monte Carlo radiativetransfer code Artis. Results: Differences in binding energies ofcarbon and oxygen lead to a lower nuclear energy release for carbondepleted material; thus, the burning fronts that develop are weaker andthe total nuclear energy release is smaller. For otherwise identicalconditions, carbon depleted models produce less 56Ni.Comparing different models with similar 56Ni yields showslower kinetic energies in the ejecta for carbon depleted models, butonly small differences in velocity distributions and line velocities inspectra. The light curve width-luminosity relation (WLR) obtained formodels with differing carbon depletion is roughly perpendicular to theobserved WLR, hence the carbon mass fraction is probably only asecondary parameter in the family of SNe Ia.
Tables 3 and 4 are available in electronic form at http://www.aanda.org

Characterization of water-soluble organic matter from urban aerosols

Water-soluble organic matter (WSOM) from atmospheric particles comprises a complex array of molecular structures that play an important role on the physic-chemical properties of atmospheric particles and, therefore, are linked to several global-relevant atmospheric processes which impact the climate and public health. Due to the large variety of sources and formation processes, adequate knowledge on WSOM composition and its effects on the properties of atmospheric aerosol are still limited. Therefore, this thesis aims at providing new insights on the molecular composition of WSOM from fine atmospheric aerosols typical of an urban area (Aveiro, Portugal). In a first step, adsorption phenomena of semivolatile organic compounds on quartz fibre filters employed in the collection of atmospheric aerosols were assessed. Afterwards, atmospheric aerosol samples were collected during fifteen months, on a weekly basis. A mass balance of aerosol samples was performed in order to set the relative contribution of elemental carbon, WSOM and water-insoluble organic matter to the aerosol mass collected at the urban area of Aveiro, with a special focus on the assessment of the influence of different meteorological conditions. In order to assess the chemical complexity of the WSOM from urban aerosols, their structural characteristics were studied by means of Fourier transform infrared infrared - Attenuated Total Reflectance (FTIR-ATR) and solid-state cross polarization with magic angle spinning 13C nuclear magnetic resonance (CPMAS 13C NMR) spectroscopies, as well as their elemental composition. The structural characterization of aerosol WSOM samples collected in the urban area highlighted a highly complex mixture of functional groups. It was concluded that aliphatic and aromatic structures, hydroxyl groups and carboxyl groups are characteristic to all samples. The semi-quantitative assessment of the CPMAS 13C NMR data showed different distributions of the various functional groups between the aerosol samples collected at different seasons. Moreover, the presence of signals typical of lignin-derived structures in both CPMAS 13C NMR and FTIR-ATR spectra of the WSOM samples from the colder seasons, highlights the major contribution of biomass burning processes in domestic fireplaces, during low temperature conditions, into the bulk chemical properties of WSOM from urban aerosols. A comprehensive two-dimensional liquid chromatography (LC x LC) method, on-line coupled to a diode array, fluorescence, and evaporative light scattering detectors, was employed for resolving the chemical heterogeneity of the aerosol WSOM samples and, simultaneously, to map the hydrophobicity versus the molecular weight distribution of the samples. The LC x LC method employed a mixed-mode hydrophilic interaction column operating under aqueous reversed phase mode in the first dimension, and a size-exclusion column in the second dimension, which was found to be useful for separating the aerosol WSOM samples into various fractions with distinct molecular weight and hydrophobic features. The estimative of the average molecular weight (Mw) distribution of the urban aerosol WSOM samples ranged from 48 to 942 Da and from 45 to 1241 Da in terms of UV absorption and fluorescence detection, respectively. Findings suggest that smaller Mw group fractions seem to be related to a more hydrophobic nature.

Underway O2/Ar-estiment net community production across the Kuroshio Extension and impact on winter pCO2 drawdown

Senior thesis written for Oceanography 445

Cartografia de zonamento geotécnico no apoio ao planeamento e valorização da exploração do georrecurso: o caso do maciço rochoso da Curviã (Joane)

Este trabalho centra-se no estudo do aproveitamento expectável do maciço rochoso da pedreira da Curviã N.o 2 (Joane, Vila Nova de Famalicão, no N Portugal), através da obtenção de um bloco unitário tipo que forneça indicações para a exploração do recurso geológico para fins industriais e/ou ornamentais. Desta forma, investiga-se se num dado limite de zona geotécnica do maciço rochoso e propicio a obtenção de blocos com dimensão, avaliados apos o processo de transformação, nomeadamente, para enrocamento em obras marítimas ou balastro em obras ferroviárias. Foram seleccionados diversos afloramentos, tendo-se recorrido a técnica de amostragem linear as superfícies expostas do maciço. Esta técnica e uma das formas mais expeditas de coligir dados geológico-geotécnicos relativos as descontinuidades. Procedeu-se, ainda, a um tratamento estatístico das descontinuidades, bem como dos parâmetros geológico-geotécnicos e geomecânicos a estas associadas, propostos pela Sociedade Internacional de Mecânica das Rochas (ISRM). Todos os dados foram representados cartograficamente numa base apoiada pelos Sistemas de Informação Geográfica (SIG) e utilizadas as ferramentas de geologia estrutural, analise morfotectónica, modelação digital de terreno e cartografia de zonamento geotécnico. O zonamento geotécnico do maciço granítico foi realizado sempre em estreita ligação com o conhecimento das características do maciço ”in situ”. Pretende-se que esta metodologia contribua para um melhor conhecimento da compartimentação dos maciços rochosos em geral e, em particular, do modelo geotécnico comportamental do maciço rochoso da Curviã N.o2.

Projeto de uma unidade de produção de biodiesel num reator ultrassónico pela via supercrítica

A concretização deste estágio na INCBIO teve como principal objetivo o projeto e dimensionamento de uma unidade de produção de biodiesel por transesterificação nãocatalítica num reator ultrassónico com condições supercríticas de pressão e temperatura. Com vista à concretização do trabalho, iniciou-se a realização do estudo do estado da arte relativo à produção de biodiesel, com particular enfoque na produção de biodiesel por via catalítica e por via supercrítica e na produção de biodiesel com uso de tecnologia ultrassónica. Conclui-se que nenhum estudo contempla a combinação simultânea da produção de biodiesel por via supercrítica através da tecnologia ultrassónica. Este estudo do estado da arte permitiu ainda definir as condições de temperatura, pressão e rácio mássico (250 °C, 95 bar e 1:1 respetivamente) a considerar no projeto da unidade de produção de biodiesel deste trabalho. Com base no estudo do estado da arte efetuado e com base nas características da matéria procedeu-se à definição do processo de produção de biodiesel. Para a definição do processo começou-se por elaborar o diagrama de blocos do processo (BFD) e o diagrama de fluxo do processo (PFD). Com base nos diagramas e na composição da matéria-prima, procedeu-se à quantificação dos reagentes (metanol) com base na estequiometria das reações envolvidas e ao cálculo do balanço de massa. O balanço de massa foi calculado com base na estequiometria das reações envolvidas e foi também calculado através do software de simulação ASPEN PLUS. Após o cálculo do balanço de massa elaborou-se o diagrama de tubulação e instrumentação (P&ID), que contém todos os equipamentos, válvulas, instrumentação e tubagens existentes na unidade. Após a definição do processo e cálculo do balanço de massa procedeu-se ao dimensionamento mecânico e cálculo hidráulico dos tanques, tubagem, bombas, permutador de calor, reator ultrassónico, válvulas de controlo e instrumentação de acordo com as normas ASME. Nesta fase do trabalho foram consultados diversos fornecedores possíveis para a compra de todo o material necessário. O dimensionamento mecânico e cálculo hidráulico efetuados permitiram, entre outras informações relevantes, obter as dimensões necessárias à construção do layout e à elaboração do desenho 3D. Com os resultados obtidos e desenhos elaborados, é possível avançar com a construção da unidade, pelo que pode-se inferir que o objetivo de projetar uma unidade de produção de biodiesel num reator ultrassónico com condições supercríticas foi alcançado. Por fim, efetuou-se uma análise económica detalhada que possibilita a comparação de uma unidade de produção de biodiesel por via catalítica (unidade de produção da INCBIO) com a unidade de produção de biodiesel num reator ultrassónico com condições supercríticas (unidade projetada neste trabalho). Por questões de confidencialidade, não foram revelados maior parte dos resultados da análise económica efetuada. No entanto, conclui-se que os custos de construção da unidade de produção de biodiesel num reator ultrassónico com condições supercríticas são mais baixos cerca de 35 a 40%, quando comparados com os custos de construção da unidade de produção de biodiesel por via catalítica, evidenciando assim que a combinação em simultâneo das condições supercríticas com a tecnologia ultrassónica possibilita a diminuição dos custos de produção.

Approche cartésienne pour le calcul du vent en terrain complexe avec application à la propagation des feux de forêt

La méthode de projection et l'approche variationnelle de Sasaki sont deux techniques permettant d'obtenir un champ vectoriel à divergence nulle à partir d'un champ initial quelconque. Pour une vitesse d'un vent en haute altitude, un champ de vitesse sur une grille décalée est généré au-dessus d'une topographie donnée par une fonction analytique. L'approche cartésienne nommée Embedded Boundary Method est utilisée pour résoudre une équation de Poisson découlant de la projection sur un domaine irrégulier avec des conditions aux limites mixtes. La solution obtenue permet de corriger le champ initial afin d'obtenir un champ respectant la loi de conservation de la masse et prenant également en compte les effets dûs à la géométrie du terrain. Le champ de vitesse ainsi généré permettra de propager un feu de forêt sur la topographie à l'aide de la méthode iso-niveaux. L'algorithme est décrit pour le cas en deux et trois dimensions et des tests de convergence sont effectués.

Études des fuites excitoniques dans des familles de boîtes quantiques d'InAs/InP par PLRT par addition de fréquences

Ce mémoire porte sur les mécanismes de relaxation et de fuite des excitons dans des systèmes de boîtes quantiques(BQs) d’InAs/InP. Les systèmes sont composés d’un sub- strat volumique d’InP, appelé matrice (M), d’un puits quantique d’InAs, nommé couche de mouillage (CM), et des familles de BQs d’InAs. La distinction entre les familles est faite par le nombre de monocouche d’épaisseur des boîtes qui sont beaucoup plus larges que hautes. Une revue de littérature retrace les principaux mécanismes de relaxation et de fuite des excitons dans les systèmes. Ensuite, différents modèles portant sur la fuite thermique des excitons des BQs sont comparés. Les types de caractérisations déjà produites et les spécifications des croissances des échantillons sont présentés. L’approche adoptée pour ce mémoire a été de caractériser temporellement la dynamique des BQs avec des mesures d’absorbtion transitoire et de photoluminescence résolue en temps (PLRT) par addition de fréquences. L’expérience d’absorption transitoire n’a pas fait ressortir de résultats très probants, mais elle est expliquée en détails. Les mesures de PLRT ont permis de suivre en température le temps de vie effectif des excitons dans des familles de BQs. Ensuite, avec un modèle de bilan détaillé, qui a été bien explicité, il a été possible d’identifier le rôle de la M et de la CM dans la relaxation et la fuite des excitons dans les BQs. Les ajustements montrent plus précisément que la fuite de porteurs dans les BQs se fait sous la forme de paires d’électrons-trous corrélées.

Développement d’une méthode d’extraction des contaminants émergents dans les solides particulaires par LDTD-APCI-MS/MS

Douze contaminants émergents (composés pharmaceutiques, pesticides et hormones) ont été quantifiés et extraits de l'eau de rivières et d’échantillons d'eaux usées municipales. La séparation des solides en suspension est effectuée par filtration des échantillons d'eau. L'effet de filtration sur les concentrations de contaminants dissous a été évaluée afin de minimiser les pertes de composés cibles. Les échantillons ont été lyophilisés et ont été extraits en deux cycles par ultrasons combinés avec une étape de nettoyage sur cartouche d’extraction de type C18. La quantification a été réalisée en utilisant la spectrométrie de masse. Les recouvrements de la méthode pour tous les composés ont varié de 68 à 112% dans toutes les matrices étudiées, sauf pour le sulfaméthoxazole et le diclofénac où des recouvrements plus modestes ont été obtenus (38 à 85%). Les limites de détection pour les 12 analytes dans les sédiments et particules en suspension (SPM) de la rivière variaient de 0,7 à 9,4 ng g-1 et de 21 à 92 ng g-1, pour les échantillons SPM de station d'épuration. Tous les contaminants émergents cibles ont été détectés à des concentrations variant de 3 à 5440 ng g-1 dans les matrices étudiées, les concentrations les plus élevées ont été observées dans les échantillons SPM de stations d'épuration. Une partie importante de certains de ces contaminants est clairement associée aux sédiments de rivière ou aux particules en suspension. L’optimisation des processus de traitement de l'eau et le devenir environnemental doit absolument tenir compte de la fraction de contaminants qui liée à des particules si on espère avoir un bilan de masse raisonnable.

Numerical Simulation of Two-Part Underwater Towing System

The motion instability is an important issue that occurs during the operation of towed underwater vehicles (TUV), which considerably affects the accuracy of high precision acoustic instrumentations housed inside the same. Out of the various parameters responsible for this, the disturbances from the tow-ship are the most significant one. The present study focus on the motion dynamics of an underwater towing system with ship induced disturbances as the input. The study focus on an innovative system called two-part towing. The methodology involves numerical modeling of the tow system, which consists of modeling of the tow-cables and vehicles formulation. Previous study in this direction used a segmental approach for the modeling of the cable. Even though, the model was successful in predicting the heave response of the tow-body, instabilities were observed in the numerical solution. The present study devises a simple approach called lumped mass spring model (LMSM) for the cable formulation. In this work, the traditional LMSM has been modified in two ways. First, by implementing advanced time integration procedures and secondly, use of a modified beam model which uses only translational degrees of freedoms for solving beam equation. A number of time integration procedures, such as Euler, Houbolt, Newmark and HHT-α were implemented in the traditional LMSM and the strength and weakness of each scheme were numerically estimated. In most of the previous studies, hydrodynamic forces acting on the tow-system such as drag and lift etc. are approximated as analytical expression of velocities. This approach restricts these models to use simple cylindrical shaped towed bodies and may not be applicable modern tow systems which are diversed in shape and complexity. Hence, this particular study, hydrodynamic parameters such as drag and lift of the tow-system are estimated using CFD techniques. To achieve this, a RANS based CFD code has been developed. Further, a new convection interpolation scheme for CFD simulation, called BNCUS, which is blend of cell based and node based formulation, was proposed in the study and numerically tested. To account for the fact that simulation takes considerable time in solving fluid dynamic equations, a dedicated parallel computing setup has been developed. Two types of computational parallelisms are explored in the current study, viz; the model for shared memory processors and distributed memory processors. In the present study, shared memory model was used for structural dynamic analysis of towing system, distributed memory one was devised in solving fluid dynamic equations.

Jmodels: Java models of earth system regulation

The JModel suite consists of a number of models of aspects of the Earth System. They can all be run from the JModels website. They are written in the Java language for maximum portability, and are capable of running on most computing platforms including Windows, MacOS and Unix/Linux. The models are controlled via graphical user interfaces (GUI), so no knowledge of computer programming is required to run them. The models currently available from the JModels website are: Ocean phosphorus cycle Ocean nitrogen and phosphorus cycles Ocean silicon and phosphorus cycles Ocean and atmosphere carbon cycle Energy radiation balance model (under development) The main purpose of the models is to investigate how material and energy cycles of the Earth system are regulated and controlled by different feedbacks. While the central focus is on these feedbacks and Earth System stabilisation, the models can also be used in other ways. These resources have been developed by: National Oceanography Centre, Southampton project led by Toby Tyrrell and Andrew Yool, focus on how the Earth system works.

A glacier inventory for the Buordakh Massif, Cherskiy Range, northeast Siberia, and evidence for recent glacier recession

The Buordakh Massif, in the Cherskiy Range of northeast Siberia, contains mountains over 3000 in and, despite its and climate, numerous glaciers. This paper presents a glacier inventory for the region and documents some 80 glaciers, which range in size from 0.1 to 10.4 km(2) (total glacierized area is ca. 70 km(2)). The inventory is based on mapping derived from Landsat 7 ETM+ satellite imagery from August 2001, augmented with data from field investigations obtained at that time. The glaciers in this region are of the 'firn-less,' cold, continental type, and their mass balance relies heavily on the formation of superimposed ice. The most recent glacier maximum extents have also been delineated, and these are believed to date from the Little Ice Age (ca. A.D. 1550-1850). Glacier areal extent has reduced by some 14.8 km(2) (ca. 17%) since this most. recent maximum. Of the 80 glaciers catalogued, 49 have undergone a measurable retreat from their most recent maximum extent.

Laurentide ice streaming on the Canadian Shield: A conflict with the soft-bedded ice stream paradigm?

Elucidating the controls on the location and vigor of ice streams is crucial to understanding the processes that lead to fast disintegration of ice flows and ice sheets. In the former North American Laurentide ice sheet, ice stream occurrence appears to have been governed by topographic troughs or areas of soft-sediment geology. This paper reports robust evidence of a major paleo-ice stream over the northwestern Canadian Shield, an area previously assumed to be incompatible with fast ice flow because of the low relief and relatively hard bedrock. A coherent pattern of subglacial bedforms (drumlins and megascalle glacial lineations) demarcates the ice stream flow set, which exhibits a convergent onset zone, a narrow main trunk with abrupt lateral margins, and a lobate terminus. Variations in bedform elongation ratio within the flow set match theoretical expectations of ice velocity. In the center of the ice stream, extremely parallel megascalle glacial lineations tens of kilometers long with elongation ratios in excess of 40:1 attest to a single episode of rapid ice flow. We conclude that while bed properties are likely to be influential in determining the occurrence and vigor of ice streams, contrary to established views, widespread soft-bed geology is not an essential requirement for those ice streams without topographic control. We speculate that the ice stream acted as a release valve on ice-sheet mass balance and was initiated by the presence of a proglacial lake that destabilized the ice-sheet margin and propagated fast ice flow through a series of thermomechanical feedbacks involving ice flow and temperature.

Late-20th-century changes in glacier extent in the Caucasus Mountains, Russia/Georgia

Glaciers occupy an area of similar to 1600 km(2) in the Caucasus Mountains. There is widespread evidence of retreat since the Little Ice Age, but an up-to-date regional assessment of glacier change is lacking. In this paper, satellite imagery (Landsat Thematic Mapper and Enhanced Thematic Mapper Plus) is used to obtain the terminus position of 113 glaciers in the central Caucasus in 1985 and 2000, using a manual delineation process based on a false-colour composite (bands 5, 4, 3). Measurements reveal that 94% of the glaciers have retreated, 4% exhibited no overall change and 2% advanced. The mean retreat rate equates to similar to 8 m a(-1), and maximum retreat rates approach similar to 38 m a(-1). The largest (>10 km(2)) glaciers retreated twice as much (similar to 12 m a(-1)) as the smallest (<1 km(2)) glaciers (similar to 6 m a(-1)), and glaciers at lower elevations generally retreated greater distances. Supraglacial debris cover has increased in association with glacier retreat, and the surface area of bare ice has reduced by similar to 10% between 1985 and 2000. Results are compared to declassified Corona imagery from the 1960s and 1970s and detailed field measurements and mass-balance data for Djankuat glacier, central Caucasus. It is concluded that the decrease in glacier area appears to be primarily driven by increasing temperatures since the 1970s and especially since the mid-1990s. Continued retreat could lead to considerable changes in glacier runoff, with implications for regional water resources.