Dynamics of settling charged particles in turbulence : theory and experiments

It has been proposed that inertial clustering may lead to an increased collision rate of water droplets in clouds. Atmospheric clouds and electrosprays contain electrically charged particles embedded in turbulent flows, often under the influence of an externally imposed, approximately uniform gravitational or electric force. In this thesis, we present the investigation of charged inertial particles embedded in turbulence. We have developed a theoretical description for the dynamics of such systems of charged, sedimenting particles in turbulence, allowing radial distribution functions to be predicted for both monodisperse and bidisperse particle size distributions. The governing parameters are the particle Stokes number (particle inertial time scale relative to turbulence dissipation time scale), the Coulomb-turbulence parameter (ratio of Coulomb ’terminalar speed to turbulence dissipation velocity scale), and the settling parameter (the ratio of the gravitational terminal speed to turbulence dissipation velocity scale). For the monodispersion particles, The peak in the radial distribution function is well predicted by the balance between the particle terminal velocity under Coulomb repulsion and a time-averaged ’drift’ velocity obtained from the nonuniform sampling of fluid strain and rotation due to finite particle inertia. The theory is compared to measured radial distribution functions for water particles in homogeneous, isotropic air turbulence. The radial distribution functions are obtained from particle positions measured in three dimensions using digital holography. The measurements support the general theoretical expression, consisting of a power law increase in particle clustering due to particle response to dissipative turbulent eddies, modulated by an exponential electrostatic interaction term. Both terms are modified as a result of the gravitational diffusion-like term, and the role of ’gravity’ is explored by imposing a macroscopic uniform electric field to create an enhanced, effective gravity. The relation between the radial distribution functions and inward mean radial relative velocity is established for charged particles.

Dynamics of the last four glacial terminations recorded in Lake Van, Turkey

A well-dated suite of Lake Van climate-proxy data covering the last 360 ka documents environmental changes over 4 glacial/interglacial cycles in Eastern Anatolia, Turkey. The picture of cold and dry glacials and warm and wet interglacials emerging from pollen, organic carbon, authigenic carbonate content, elemental profiling by XRF and lithological analyses is inconsistent with classical interpretation of ox- ygen isotopic composition of carbonates pointing to a more complex pattern in Lake Van region. Detailed analysis of glacial terminations allows for the constraining of a depositional model explaining different patterns observed in all the proxies. We hypothesize that variations in relative contribution of rainfall, snowmelt and glacier meltwater recharging the basin have a very important role for all sedimentary processes in Lake Van. Lake level of glacial Lake Van, predominantly fed by snowmelt, was low, the water column was oxic, and carbonates precipitating in the epilimnion recorded the light isotopic signature of inflow. During terminations, increasing rainfall and significant supply of mountain glaciers' meltwater contributed to lake level rise. Increased rainfall enhanced density gradients in the water column, and hindered mixing leading to development of bottom-water anoxia. Carbonates precipitating during terminations show large fluctuations in their isotopic composition. Full interglacial conditions in Lake Van are characterized by high or slowly falling lake level. Rainfall and snowmelt feed the lake but due to re-established mixing, the isotopic composition of authigenic carbonates is heavier and closer to that of evaporation-influenced lake water than that of runoff representing snowmelt and atmospheric precipitation.

DYPTOP: a cost-efficient TOPMODEL implementation to simulate sub-grid spatio-temporal dynamics of global wetlands and peatlands

Simulating the spatio-temporal dynamics of inundation is key to understanding the role of wetlands under past and future climate change. Earlier modelling studies have mostly relied on fixed prescribed peatland maps and inundation time series of limited temporal coverage. Here, we describe and assess the the Dynamical Peatland Model Based on TOPMODEL (DYPTOP), which predicts the extent of inundation based on a computationally efficient TOPMODEL implementation. This approach rests on an empirical, grid-cell-specific relationship between the mean soil water balance and the flooded area. DYPTOP combines the simulated inundation extent and its temporal persistency with criteria for the ecosystem water balance and the modelled peatland-specific soil carbon balance to predict the global distribution of peatlands. We apply DYPTOP in combination with the LPX-Bern DGVM and benchmark the global-scale distribution, extent, and seasonality of inundation against satellite data. DYPTOP successfully predicts the spatial distribution and extent of wetlands and major boreal and tropical peatland complexes and reveals the governing limitations to peatland occurrence across the globe. Peatlands covering large boreal lowlands are reproduced only when accounting for a positive feedback induced by the enhanced mean soil water holding capacity in peatland-dominated regions. DYPTOP is designed to minimize input data requirements, optimizes computational efficiency and allows for a modular adoption in Earth system models.

Application of a computational model for complex fluvial ecosystems: The population dynamics of zebra mussel Dreissena polymorpha as a case study

The potential and adaptive flexibility of population dynamic P-systems (PDP) to study population dynamics suggests that they may be suitable for modelling complex fluvial ecosystems, characterized by a composition of dynamic habitats with many variables that interact simultaneously. Using as a model a reservoir occupied by the zebra mussel Dreissena polymorpha, we designed a computational model based on P systems to study the population dynamics of larvae, in order to evaluate management actions to control or eradicate this invasive species. The population dynamics of this species was simulated under different scenarios ranging from the absence of water flow change to a weekly variation with different flow rates, to the actual hydrodynamic situation of an intermediate flow rate. Our results show that PDP models can be very useful tools to model complex, partially desynchronized, processes that work in parallel. This allows the study of complex hydroecological processes such as the one presented, where reproductive cycles, temperature and water dynamics are involved in the desynchronization of the population dynamics both, within areas and among them. The results obtained may be useful in the management of other reservoirs with similar hydrodynamic situations in which the presence of this invasive species has been documented.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Dynamics of a local Alpine flooding event in October 2011: moisture source and large-scale circulation

Alpine heavy precipitation events often affect small catchments, although the circulation pattern leading to the event extends over the entire North Atlantic. The various scale interactions involved are particularly challenging for the numerical weather prediction of such events. Unlike previous studies focusing on the southern Alps, here a comprehensive study of a heavy precipitation event in the northern Alps in October 2011 is presented with particular focus on the role of the large-scale circulation in the North Atlantic/European region. During the event exceptionally high amounts of total precipitable water occurred in and north of the Alps. This moisture was initially transported along the flanks of a blocking ridge over the North Atlantic. Subsequently, strong and persistent northerly flow established at the upstream flank of a trough over Europe and steered the moisture towards the northern Alps. Lagrangian diagnostics reveal that a large fraction of the moisture emerged from the West African coast where a subtropical upper-level cut-off low served as an important moisture collector. Wave activity flux diagnostics show that the ridge was initiated as part of a low-frequency, large-scale Rossby wave train while convergence of fast transients helped to amplify it locally in the North Atlantic. A novel diagnostic for advective potential vorticity tendencies sheds more light on this amplification and further emphasizes the role of the ridge in amplifying the trough over Europe. Operational forecasts misrepresented the amplitude and orientation of this trough. For the first time, this study documents an important pathway for northern Alpine flooding, in which the interaction of synoptic-scale to large-scale weather systems and of long-range moisture transport from the Tropics are dominant. Moreover, the trapping of moisture in a subtropical cut-off near the West African coast is found to be a crucial precursor to the observed European high-impact weather.

Comparative carbon cycle dynamics of the present and last interglacial

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO₂ remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO₂ fertilization, land use, wild fire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO₂ dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO₂ and δ¹³ CO₂ changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO₂ dynamics from 8ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO₂ changes after 122 ka BP. This failure to simulate late-Eemian CO₂ dynamics could be a result of the imposed forcings such as prescribed CaCO₃ accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO₂ dynamics eshallow water CaCO₃ accumulation, peat and permafrost carbon dynamics are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO₂ dynamics.

Molecular Machanics of Water Monomer-Maple Code

This is the Maple code to support the molecular dynamics of a water monomer molecule, allowing investigation of the classical vibrations of this molecule.

Molecular Mechanics of Water Monomer

Simulation of the classical molecular dynamics of a water molecule can be useful in explaining normal modes of motion, Fourier Transforms, and fundamental frequencies of vibration, as illustrated herein.

Water, bottom deposits, and macrobenthos in the south part of the East Siberian Sea in September 2003

Macrobenthos biomass and bottom biocoenoses were studied in the sublittoral zone of the southern East Siberian Sea. The macrobenthos is characterized by relatively high abundance (from 30 to 2680 #/m**2), biomass (from 0.25 to 578.8 g/m**2), and diversity (83 species in total). Lateral distribution of macrobenthos biomass correlates with a substrate type and salinity and is substantially higher in areas washed by the Arctic water mass than in estuaries with mixed fresh and Arctic waters and shows a tendency to decreasing in the convergence zone of different water masses. The highest macrobenthos biomass is observed in cores of water masses in the Long Strait area and in the eastern part of the sea.

Composition and fluxes of sinking particulate material and bottom sediments in the northeast Black Sea in 1998-1999

For the first time deep-sea mooring stations with sediment traps were deployed in the northeast Black Sea. One sediment trap for long-term studies was located at Station 1 (44°15'N, 37°43'E, deployment depth 1800 m, depth 1900 m). The trap collected sinking sedimentary material from January to May 1998. Material collectors were changed every 15 days. Other stations with sediment traps for short-term studies (September-October 1999) were located on the shelf: Station 2 (44°16'N, 38°37'E, deployment depth 45 m, depth 50 m) and on the bottom of the canyon: Station 3 (44°16'N, 38°22'E, deployment depth 1145 m, depth 1150 m), Station 4 (44°11'N, 38°21'E, deployment depths 200, 1550, 1650 m, depth 1670 m). Collected material indicates that vertical particle fluxes are controlled by seasonal changes of in situ production and by dynamics of terrigenous matter input. Higher vertical particle flux of carbonate and biogenic silica was in spring due to bloom of plankton organisms. Maximum of coccolith bloom is in April-May. Bloom of diatoms begins in March. In winter and autumn lithogenic material dominates in total flux. Its amount strongly depends on storms and river run-off. Suspended particle material differs from surface shelf sediments by finer particles (mainly clay fraction) and high content of clay minerals and biogenic silica. This material may form lateral fluxes with higher concentration of particles transported along the bottom of deep-sea canyons from the shelf to the deep basin within the nepheloid layer. In winter such transportation of sedimentary material is more intensive due to active vertical circulation of water masses.

The hydrographic conditions at depth below 1000 m water depth in the region of the Romanche Fracture Zones

During the "Atlantic Expedition" in1965 (IQSY) a comprehensive bathymetric survey and a few hydrographic stations were made by R.V. "Meteor" in the equatorial region of the Mid-Atlantic Ridge. The survey results are shown in a bythymetric chart covering the western parts of the Romanche- and Chain Fracture Zones. West of the original Romanche Trench another deep trench with a medium depth of 6000 m was discovered. The maximum sounding obtained was 7028 m. Both trenches apparently belong to the same fracture zone, but are distinctly separated from each other. The estern boundary of the trench against the Brasil Basin is formed by a sill rising to a depth of about 4400 m. The serial hydrographic observations give some indications of the flow of the cold Westatlantic deep water in the fracture zone area and its influence on the hydrographic conditions in the East-Atlantic Basin. The upper limit of the nearly homogenious Westatlantic bottom water with an Antarctic components lies about 4400 m. The water mass entering the system of trenches of the Romanche Fracture Zone over the western sill originates from the lower part of the discontinuity layer lying above the bottom water. Potential temperatures of 0.6°C were the lowest observed by "Meteor" in the western trench. There seems to be a remarkable tongue of relatively high salinity and a minimum of oxygen in the deep water of this trench. At present we can only speculate upon the origin of this highly saline deep water tongue underneath the eastward moving relatively thin layer of less saline Westatlantic deep water. In the range of the sill separating both trenches a lee wave is indicated by the distribution of salinity and oxygen, which implies a vertical transport of water masses. Caused by this transport it is assumed that relatively cold water may be lifted temporarily to a depth, where it can pass the northbounding ridge, thus getting directly into the Sierra Leone Basin. In the original Romanche Trench the cold Westatlantic deep water seems to fill the whole trough, but its extension remains limited to the trench itself. The water masses found east of the sill separating the trench from the East-Atlantic Basin originate from the lower part of the discontinuity layer. With potential temperatures of about 1.3°C they are much warmer than those observed in the Romanche Trench bottom water.

Size variations of Discoaster multiradiatus across the Paleocene-Eocene thermal maximum

Size measurements of the calcareous nannofossil taxon Discoaster multiradiatus were carried out across the Paleocene-Eocene Thermal Maximum (PETM) in Ocean Drilling Program Holes 690B (Maud Rise, Weddell Sea) and 1209B (Shatsky Rise, Pacific Ocean). Morphometric investigations show that D. multiradiatus specimens are generally larger at ODP Site 1209 than at ODP Site 690. A limited increase in size of D. multiradiatus is recorded at ODP Site 1209, whereas significant enlargements characterize ODP Site 690. Preservation is comparable at both sites: nannofossils are moderately preserved with some evidence of etching/overgrowth in the PETM interval. Yet, D. multiradiatus variations do not correlate with preservation state and morphometric data most likely represent primary signals rather than diagenetic artifacts. There is a direct relationship between D. multiradiatus size and paleotemperatures: largest specimens are coeval with global warming associated with the PETM, inferred to result from excess atmospheric CO2 due to (partial) oxidation of massive quantities of methane. Size increases and largest specimens of D. multiradiatus occur at different stratigraphic levels within PETM at ODP Sites 690 and 1209. A marked shift in diameter size was observed at the onset and peak of the Carbon Isotopic Excursion (CIE) at ODP Site 690, but only at the end of CIE and initial recovery interval at ODP Site 1209. This diachroneity is puzzling, but indeed correlates well with reconstructed changes in surface and thermocline water masses temperature and salinity in the PETM interval at low and high latitudes. The presumed high concentrations of carbon dioxide seem to have not influenced the morphometry of D. multiradiatus. The major size increase of D. multiradiatus in the CIE of ODP Site 690 could represent the migration of larger-sized allochtonus specimens that moved from peri-equatorial/subtropical areas to higher latitudes during the warmest interval of the PETM, although no direct evidence of distinct populations/subpopulations has been obtained from the frequency diagrams. As a result, we infer that D. multiradiatus is a proxy of water masses stratification and might be used for deriving temperature-salinity-nutrient conditions in the mixed layer and thermocline and their dynamics.

Spatial and seasonal distribution of Thaumarchaeota in the water column and sediment of the southern North Sea

We have examined the spatial and seasonal distribution of Thaumarchaeota in the water column and sediment of the southern North Sea using the specific intact polar lipid (IPL) hexose, phosphohexose (HPH) crenarchaeol, as well as thaumarchaeotal 16S rRNA gene abundances and expression. In the water column, a higher abundance of Thaumarchaeota was observed in the winter season than in the summer, which is in agreement with previous studies, but this was not the case in the sediment where Thaumarchaeota were most abundant in spring and summer. This observation corresponds well with the idea that ammonia availability is a key factor in thaumarchaeotal niche determination. In the surface waters of the southern North Sea, we observed a spatial variability in HPH crenarchaeol, thaumarchaeotal 16S rRNA gene abundance and transcriptional activity that corresponded well with the different water masses present. In bottom waters, a clear differentiation based on water masses was not observed; instead, we suggest that observed differences in thaumarchaeotal abundance with depth may be related to resuspension from the sediment. This could be due to suspension of benthic Thaumarchaeota to the water column or due to delivery of e.g. resuspended sediment or ammonium to the water column, which could be utilized by pelagic Thaumarchaeota. This study has shown that the seasonality of Thaumarchaeota in water and sediment is different and highlights the importance of water masses, currents and sedimentary processes in determining the spatial abundance of Thaumarchaeota in the southern North Sea.