Wet and dry model granulates under mechanical load : a confocal microscopy study

Granular matter, also known as bulk solids, consists of discrete particles with sizes between micrometers and meters. They are present in many industrial applications as well as daily life, like in food processing, pharmaceutics or in the oil and mining industry. When handling granular matter the bulk solids are stored, mixed, conveyed or filtered. These techniques are based on observations in macroscopic experiments, i.e. rheological examinations of the bulk properties. Despite the amply investigations of bulk mechanics, the relation between single particle motion and macroscopic behavior is still not well understood. For exploring the microscopic properties on a single particle level, 3D imaging techniques are required.rnThe objective of this work was the investigation of single particle motions in a bulk system in 3D under an external mechanical load, i.e. compression and shear. During the mechanical load the structural and dynamical properties of these systems were examined with confocal microscopy. Therefor new granular model systems in the wet and dry state were designed and prepared. As the particles are solid bodies, their motion is described by six degrees of freedom. To explore their entire motion with all degrees of freedom, a technique to visualize the rotation of spherical micrometer sized particles in 3D was developed. rnOne of the foci during this dissertation was a model system for dry cohesive granular matter. In such systems the particle motion during a compression of the granular matter was investigated. In general the rotation of single particles was the more sensitive parameter compared to the translation. In regions with large structural changes the rotation had an earlier onset than the translation. In granular systems under shear, shear dilatation and shear zone formation were observed. Globally the granular sediments showed a shear behavior, which was known already from classical shear experiments, for example with Jenike cells. Locally the shear zone formation was enhanced, when near the applied load a pre-diluted region existed. In regions with constant volume fraction a mixing between the different particle layers occurred. In particular an exchange of particles between the current flowing region and the non-flowing region was observed. rnThe second focus was on model systems for wet granular matter, where an additional binding liquid is added to the particle suspension. To examine the 3D structure of the binding liquid on the micrometer scale independently from the particles, a second illumination and detection beam path was implemented. In shear and compression experiments of wet clusters and bulk systems completely different dynamics compared to dry cohesive models systems occured. In a Pickering emulsion-like system large structural changes predominantly occurred in the local environment of binding liquid droplets. These large local structural changes were due to an energy interplay between the energy stored in the binding droplet during its deformation and the binding energy of particles at the droplet interface. rnConfocal microscopy in combination with nanoindentation gave new insights into the single particle motions and dynamics of granular systems under a mechanical load. These novel experimental results can help to improve the understanding of the relationship between bulk properties of granular matter, such as volume fraction or yield stress and the dynamics on a single particle level.rnrn

Volatile organic compounds at the air-sea interface : gas exchange rates, oceanic emissions and the effect of ocean acidification

Oceans are key sources and sinks in the global budgets of significant atmospheric trace gases, termed Volatile Organic Compounds (VOCs). Despite their low concentrations, these species have an important role in the atmosphere, influencing ozone photochemistry and aerosol physics. Surprisingly, little work has been done on assessing their emissions or transport mechanisms and rates between ocean and atmosphere, all of which are important when modelling the atmosphere accurately.rnA new Needle Trap Device (NTD) - GC-MS method was developed for the effective sampling and analysis of VOCs in seawater. Good repeatability (RSDs <16 %), linearity (R2 = 0.96 - 0.99) and limits of detection in the range of pM were obtained for DMS, isoprene, benzene, toluene, p-xylene, (+)-α-pinene and (-)-α-pinene. Laboratory evaluation and subsequent field application indicated that the proposed method can be used successfully in place of the more usually applied extraction techniques (P&T, SPME) to extend the suite of species typically measured in the ocean and improve detection limits. rnDuring a mesocosm CO2 enrichment study, DMS, isoprene and α-pinene were identified and quantified in seawater samples, using the above mentioned method. Based on correlations with available biological datasets, the effects of ocean acidification as well as possible ocean biological sources were investigated for all examined compounds. Future ocean's acidity was shown to decrease oceanic DMS production, possibly impact isoprene emissions but not affect the production of α-pinene. rnIn a separate activity, ocean - atmosphere interactions were simulated in a large scale wind-wave canal facility, in order to investigate the gas exchange process and its controlling mechanisms. Air-water exchange rates of 14 chemical species (of which 11 VOCs) spanning a wide range of solubility (dimensionless solubility, α = 0:4 to 5470) and diffusivity (Schmidt number in water, Scw = 594 to 1194) were obtained under various turbulent (wind speed at ten meters height, u10 = 0:8 to 15ms-1) and surfactant modulated (two different sized Triton X-100 layers) surface conditions. Reliable and reproducible total gas transfer velocities were obtained and the derived values and trends were comparable to previous investigations. Through this study, a much better and more comprehensive understanding of the gas exchange process was accomplished. The role of friction velocity, uw* and mean square slope, σs2 in defining phenomena such as waves and wave breaking, near surface turbulence, bubbles and surface films was recognized as very significant. uw* was determined as the ideal turbulent parameter while σs2 described best the related surface conditions. A combination of both uw* and σs2 variables, was found to reproduce faithfully the air-water gas exchange process. rnA Total Transfer Velocity (TTV) model provided by a compilation of 14 tracers and a combination of both uw* and σs2 parameters, is proposed for the first time. Through the proposed TTV parameterization, a new physical perspective is presented which provides an accurate TTV for any tracer within the examined solubility range. rnThe development of such a comprehensive air-sea gas exchange parameterization represents a highly useful tool for regional and global models, providing accurate total transfer velocity estimations for any tracer and any sea-surface status, simplifying the calculation process and eliminating inevitable calculation uncertainty connected with the selection or combination of different parameterizations.rnrn

Megistobenthic faunal diversity of the Antalya Gulf: Crustacea

The spatio-temporal distribution of megistobenthic crustacean assemblages from the Antalya Gulf, located in the Levantine Sea is described. In order to provide a comprehensive overview of the spatio-temporal patterns of the crustacean community, 3 transect including depth of 10, 25, 75, 125 and 200 m, were studied between 2014 and 2015 to investigate their association with a set of environmental parameters in representative months of each season (spring, summer, autumn and winter). For its economic importance in Levantine waters, a focus analysis of deep-water rose shrimp Parapenaeus longirostris (Lucas, 1846) was done, to investigate the length frequency composition of the population of the Antalya Gulf. A total of 58 crustacean species were encountered in the study area, of these species identified, 18 species were recognized as alien species in the Mediterranean Sea. Throughout the year the most frequent species of the study were the hermit crab Pagurus prideaux (Leach, 1815) and Parapenaeus longirostris (Lucas, 1846) followed by the Indo-Pacific swimming crab Charybdis longicollis (Leene, 1938) and by the invasive shrimp Marsupenaeus japonicus (Spence Bate, 1888). Few species contributing to a high amount to the total biomass were found throughout the year. These species were Charybdis longicollis and Parapenaeus longirostris. Stations of the study area showed similar values of diversity indices of benthic crustacean community among the three transect. The highest values of faunistic indices were detected in autumn and winter (October and February), and also varied along the depth gradient, with the highest values found between 25 and 75 meters. The multivariate analyses conducted on the abundance data point out major differences between depths and between seasons. Therefore, according to cluster analysis and ordination over abundance and biomass, three main crustacean assemblages were detected: the first corresponding to shallow bottoms (10, 25 meters), the second corresponding to intermediate waters (75 meters) and the last to deeper waters (125, 200 meters). Depth was the main factor governing the distribution of megistobenthic crustacean in the area. Besides the depth, the structure of the sediment is the most important factor in determining the crustacean assemblage. Therefore, all factors governing the crustacean distribution were found to be related to the bottom depth. The population of Parapenaeus longirostris in the Antalya Gulf showed significant differences in depth. It was found that females dominated the population of the study area (65.11%), and were significantly larger than males for each cohort identified. The size-weight relationships revealed a slight negative allometry in growth, a bit more pronounced in females than in males.

Assessment of the resilience of coral reefs to natural and human disturbances by means of recruitment panels in Indo-Pacific

Indo-Pacific region encompasses about 75% of world's coral reefs, but hard coral cover in this region experienced a 32% region-wide decline since 1970s. This great change is primarily ascribable to natural and anthropogenic pressures, including climate change and human activities effects. Coral reef conservation requires management strategies oriented to maintain their diversity and the capacity to provide ecosystem goods and services. Coral reef resilience, i.e. the capacity to recover after disturbances, is critical to their long-term persistence. The aims of the present study were to design and to test field experiments intended to measure changes in recruitment processes, as a fundamental aspect of the coral reef resilience. Recruitment experiments, using artificial panels suspended in the water column, were carried out in two Indo-Pacific locations affected by different disturbances: a new mine in Bangka Island (Indonesia), and the increased sedimentation due to coastal dynamics in Vavvaru Island (Maldives). One (or more) putatively disturbed site(s) was selected to be tested against 3 randomly selected control sites. Panels’ arrangement simulates 2 proximities to living corals, i.e. the sources of propagules: few centimetres and 2 meters over. Panels were deployed simultaneously at each site and left submerged for about five months. Recruits were identified to the lowest possible taxonomic level and recruited assemblages were analysed in terms of percent cover. In general it was not possible to detect significant differences between the benthic assemblages recruited in disturbed and control sites. The high variability observed in recruits assemblages structure among control sites may be so large to mask the possible disturbance effects. Only few taxa showed possible effects of the disturb they undergo. The field tests have highlighted strengths and weaknesses of the proposed approach and, based on these results, some possible improvements were suggested.

Impact of irrigation in the Po valley (Italy) on meteorological parameters through numerical simulation

Recent studies found that soil-atmosphere coupling features, through soil moisture, have been crucial to simulate well heat waves amplitude, duration and intensity. Moreover, it was found that soil moisture depletion both in Winter and Spring anticipates strong heat waves during the Summer. Irrigation in geophysical studies can be intended as an anthropogenic forcing to the soil-moisture, besides changes in land proprieties. In this study, the irrigation was add to a LAM hydrostatic model (BOLAM) and coupled with the soil. The response of the model to irrigation perturbation is analyzed during a dry Summer season. To identify a dry Summer, with overall positive temperature anomalies, an extensive climatological characterization of 2015 was done. The method included a statistical validation on the reference period distribution used to calculate the anomalies. Drought conditions were observed during Summer 2015 and previous seasons, both on the analyzed region and the Alps. Moreover July was characterized as an extreme event for the referred distribution. The numerical simulation consisted on the summer season of 2015 and two run: a control run (CTR), with the soil coupling and a perturbed run (IPR). The perturbation consists on a mask of land use created from the Cropland FAO dataset, where an irrigation water flux of 3 mm/day was applied from 6 A.M. to 9 A.M. every day. The results show that differences between CTR and IPR has a strong daily cycle. The main modifications are on the air masses proprieties, not on to the dynamics. However, changes in the circulation at the boundaries of the Po Valley are observed, and a diagnostic spatial correlation of variable differences shows that soil moisture perturbation explains well the variation observed in the 2 meters height temperature and in the latent heat fluxes.On the other hand, does not explain the spatial shift up and downslope observed during different periods of the day. Given the results, irrigation process affects the atmospheric proprieties on a larger scale than the irrigation, therefore it is important in daily forecast, particularly during hot and dry periods.

Phenology for tropoclimatological surveys and large-scale mapping

Biotic and abiotic phenological observations can be collected from continental to local spatial scale. Plant phenological observations may only be recorded wherever there is vegetation. Fog, snow and ice are available as phenological para-meters wherever they appear. The singularity of phenological observations is the possibility of spatial intensification to a microclimatic scale where the equipment of meteorological measurements is too expensive for intensive campaigning. The omnipresence of region-specific phenological parameters allows monitoring for a spatially much more detailed assessment of climate change than with weather data. We demonstrate this concept with phenological observations with the use of a special network in the Canton of Berne, Switzerland, with up to 600 observations sites (more than 1 to 10 km² of the inhabited area). Classic cartography, gridding, the integration into a Geographic Information System GIS and large-scale analysis are the steps to a detailed knowledge of topoclimatic conditions of a mountainous area. Examples of urban phenology provide other types of spatially detailed applications. Large potential in phenological mapping in future analyses lies in combining traditionally observed species-specific phenology with remotely sensed and modelled phenology that provide strong spatial information. This is a long history from cartographic intuition to algorithm-based representations of phenology.

A remotely interrogated all-optical Rb-87 magnetometer

Atomic magnetometry was performed at Earth's magnetic field over a free-space distance of ten meters. Two laser beams aimed at a distant alkali-vapor cell excited and detected the Rb-87 magnetic resonance, allowing the magnetic field within the cell to be interrogated remotely. Operated as a driven oscillator, the magnetometer measured the geomagnetic field with less than or similar to 3.5 pT precision in a similar to 2 s data acquisition; this precision was likely limited by ambient field fluctuations. The sensor was also operated in self-oscillating mode with a 5.3 pT root Hz noise floor. Further optimization will yield a high-bandwidth, fully remote magnetometer with sub-pT sensitivity. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4747206]

Assessment of three frequently used blood glucose monitoring devices in clinical routine

Self-monitoring of blood glucose plays an important role in the management of diabetes and has been shown to improve metabolic control. The use of blood glucose meters in clinical practice requires sufficient reliability to allow adequate treatment. Direct comparison of different blood glucose meters in clinical practice, independent of the manufactures is scarce. We, therefore, aimed to evaluate three frequently used blood glucose meters in daily clinical practice.

Evaluating the Impact of the Number of Access Points in Mobile Robots Localization Using Artificial Neural Networks

Localization is information of fundamental importance to carry out various tasks in the mobile robotic area. The exact degree of precision required in the localization depends on the nature of the task. The GPS provides global position estimation but is restricted to outdoor environments and has an inherent imprecision of a few meters. In indoor spaces, other sensors like lasers and cameras are commonly used for position estimation, but these require landmarks (or maps) in the environment and a fair amount of computation to process complex algorithms. These sensors also have a limited field of vision. Currently, Wireless Networks (WN) are widely available in indoor environments and can allow efficient global localization that requires relatively low computing resources. However, the inherent instability in the wireless signal prevents it from being used for very accurate position estimation. The growth in the number of Access Points (AP) increases the overlap signals areas and this could be a useful means of improving the precision of the localization. In this paper we evaluate the impact of the number of Access Points in mobile nodes localization using Artificial Neural Networks (ANN). We use three to eight APs as a source signal and show how the ANNs learn and generalize the data. Added to this, we evaluate the robustness of the ANNs and evaluate a heuristic to try to decrease the error in the localization. In order to validate our approach several ANNs topologies have been evaluated in experimental tests that were conducted with a mobile node in an indoor space.

Assessing the Validity of Brain Glucose Concentration Measurement Using Microdialysis

Brain functions, such as learning, orchestrating locomotion, memory recall, and processing information, all require glucose as a source of energy. During these functions, the glucose concentration decreases as the glucose is being consumed by brain cells. By measuring this drop in concentration, it is possible to determine which parts of the brain are used during specific functions and consequently, how much energy the brain requires to complete the function. One way to measure in vivo brain glucose levels is with a microdialysis probe. The drawback of this analytical procedure, as with many steadystate fluid flow systems, is that the probe fluid will not reach equilibrium with the brain fluid. Therefore, brain concentration is inferred by taking samples at multiple inlet glucose concentrations and finding a point of convergence. The goal of this thesis is to create a three-dimensional, time-dependent, finite element representation of the brainprobe system in COMSOL 4.2 that describes the diffusion and convection of glucose. Once validated with experimental results, this model can then be used to test parameters that experiments cannot access. When simulations were run using published values for physical constants (i.e. diffusivities, density and viscosity), the resulting glucose model concentrations were within the error of the experimental data. This verifies that the model is an accurate representation of the physical system. In addition to accurately describing the experimental brain-probe system, the model I created is able to show the validity of zero-net-flux for a given experiment. A useful discovery is that the slope of the zero-net-flux line is dependent on perfusate flow rate and diffusion coefficients, but it is independent of brain glucose concentrations. The model was simplified with the realization that the perfusate is at thermal equilibrium with the brain throughout the active region of the probe. This allowed for the assumption that all model parameters are temperature independent. The time to steady-state for the probe is approximately one minute. However, the signal degrades in the exit tubing due to Taylor dispersion, on the order of two minutes for two meters of tubing. Given an analytical instrument requiring a five μL aliquot, the smallest brain process measurable for this system is 13 minutes.

Using Micro-Gravity Techniques to Map Alluvium Thickness and Pleistocene Location of the West Branch of the Susquehanna River Near Muncy, Pennsylvania

Laurentide glaciation during the early Pleistocene (~970 ka) dammed the southeast-flowing West Branch of the Susquehanna River (WBSR), scouring bedrock and creating 100-km-long glacial Lake Lesley near the Great Bend at Muncy, Pennsylvania (Ramage et al., 1998). Local drill logs and well data indicate that subsequent paleo-outwash floods and modern fluvial processes have deposited as much as 30 meters of alluvium in this area, but little is known about the valley fill architecture and the bedrock-alluvium interface. By gaining a greater understanding of the bedrock-alluvium interface the project will not only supplement existing depth to bedrock information, but also provide information pertinent to the evolution of the Muncy Valley landscape. This project determined if variations in the thickness of the valley fill were detectable using micro-gravity techniques to map the bedrock-alluvium interface. The gravity method was deemed appropriate due to scale of the study area (~30 km2), ease of operation by a single person, and the available geophysical equipment. A LaCoste and Romberg Gravitron unit was used to collect gravitational field readings at 49 locations over 5 transects across the Muncy Creek and Susquehanna River valleys (approximately 30 km2), with at least two gravity base stations per transect. Precise latitude, longitude and ground surface elevation at each location were measured using an OPUS corrected Trimble RTK-GPS unit. Base stations were chosen based on ease of access due to the necessity of repeat measurements. Gravity measurement locations were selected and marked to provide easy access and repeat measurements. The gravimeter was returned to a base station within every two hours and a looping procedure was used to determine drift and maximize confidence in the gravity measurements. A two-minute calibration reading at each station was used to minimize any tares in the data. The Gravitron digitally recorded finite impulse response filtered gravity measurements every 20 seconds at each station. A measurement period of 15 minutes was used for each base station occupation and a minimum of 5 minutes at all other locations. Longer or multiple measurements were utilized at some sites if drift or other externalities (i.e. train or truck traffic) were effecting readings. Average, median, standard deviation and 95% confidence interval were calculated for each station. Tidal, drift, latitude, free-air, Bouguer and terrain corrections were then applied. The results show that the gravitational field decreases as alluvium thickness increases across the axes of the Susquehanna River and Muncy Creek valleys. However, the location of the gravity low does not correspond with the present-day location of the West Branch of the Susquehanna River (WBSR), suggesting that the WBSR may have been constrained along Bald Eagle Mountain by a glacial lobe originating from the Muncy Creek Valley to the northeast. Using a 3-D inversion model, the topography of the bedrock-alluvium interface was determined over the extent of the study area using a density contrast of -0.8 g/cm3. Our results are consistent with the bedrock geometry of the area, and provide a low-cost, non-invasive and efficient method for exploring the subsurface and for supplementing existing well data.

Design Simulation for Solar Energy Research and Education

Solar research is primarily conducted in regions with consistent sunlight, severely limiting research opportunities in many areas. Unfortunately, the unreliable weather in Lewisburg, PA, can prove difficult for such testing to be conducted. As such, a solar simulator was developed for educational purposes for the Mechanical Engineering department at Bucknell University. The objective of this work was to first develop a geometric model to evaluate a one sun solar simulator. This was intended to provide a simplified model that could be used without the necessity of expensive software. This model was originally intended to be validated experimentally, but instead was done using a proven ray tracing program, TracePro. Analyses with the geometrical model and TracePro demonstrated the influence the geometrical properties had results, specifically the reflector (aperture) diameter and the rim angle. Subsequently, the two were approaches were consistent with one another for aperture diameters 0.5 m and larger, and for rim angles larger than 45°. The constructed prototype, that is currently untested, was designed from information provided by the geometric model, includes a metal halide lamp with a 9.5 mm arc diameter and parabolic reflector with an aperture diameter of 0.631 meters. The maximum angular divergence from the geometrical model was predicted to be 30 mRadians. The average angular divergence in TraceProof the system was 19.5 mRadians, compared to the sun’s divergence of 9.2 mRadians. Flux mapping in TracePro showed an intensity of 1000 W/m2 over the target plane located 40 meters from the lamp. The error between spectrum of the metal halide lamp and the solar spectrum was 10.9%, which was found by comparing their respective Plank radiation distributions. The project did not satisfy the original goal of matching the angular divergence of sunlight, although the system could still to be used for optical testing. The geometric model indicated performance in this area could be improved by increasing the diameter of the reflector, as well as decreasing the source diameter. Although ray tracing software provides more information to analyze the simulator system, the geometrical model is adequate to provide enough information to design a system.

Seasonal variation in orthopedic health services utilization in Switzerland: the impact of winter sport tourism

BACKGROUND: Climate- or holiday-related seasonality in hospital admission rates is well known for many diseases. However, little research has addressed the impact of tourism on seasonality in admission rates. We therefore investigated the influence of tourism on emergency admission rates in Switzerland, where winter and summer leisure sport activities in large mountain regions can generate orthopedic injuries. METHODS: Using small area analysis, orthopedic hospital service areas (HSAo) were evaluated for seasonality in emergency admission rates. Winter sport areas were defined using guest bed accommodation rate patterns of guest houses and hotels located above 1000 meters altitude that show clear winter and summer peak seasons. Emergency admissions (years 2000-2002, n = 135'460) of local and nonlocal HSAo residents were evaluated. HSAo were grouped according to their area type (regular or winter sport area) and monthly analyses of admission rates were performed. RESULTS: Of HSAo within the defined winter sport areas 70.8% show a seasonal, summer-winter peak hospital admission rate pattern and only 1 HSAo outside the defined winter sport areas shows such a pattern. Seasonal hospital admission rates in HSAo in winter sport areas can be up to 4 times higher in winter than the intermediate seasons, and they are almost entirely due to admissions of nonlocal residents. These nonlocal residents are in general -and especially in winter- younger than local residents, and nonlocal residents have a shorter length of stay in winter sport than in regular areas. The overall geographic distribution of nonlocal residents admitted for emergencies shows highest rates during the winter as well as the summer in the winter sport areas. CONCLUSION: Small area analysis using orthopedic hospital service areas is a reliable method for the evaluation of seasonality in hospital admission rates. In Switzerland, HSAo defined as winter sport areas show a clear seasonal fluctuation in admission rates of only nonlocal residents, whereas HSAo defined as regular, non-winter sport areas do not show such seasonality. We conclude that leisure sport, and especially ski/snowboard tourism demands great flexibility in hospital beds, staff and resource planning in these areas.

Effect of mesh distortion on the accuracy of high order vorticity-velocity CFD approaches

The numerical solution of the incompressible Navier-Stokes equations offers an alternative to experimental analysis of fluid-structure interaction (FSI). We would save a lot of time and effort and help cut back on costs, if we are able to accurately model systems by these numerical solutions. These advantages are even more obvious when considering huge structures like bridges, high rise buildings or even wind turbine blades with diameters as large as 200 meters. The modeling of such processes, however, involves complex multiphysics problems along with complex geometries. This thesis focuses on a novel vorticity-velocity formulation called the Kinematic Laplacian Equation (KLE) to solve the incompressible Navier-stokes equations for such FSI problems. This scheme allows for the implementation of robust adaptive ordinary differential equations (ODE) time integration schemes, allowing us to tackle each problem as a separate module. The current algortihm for the KLE uses an unstructured quadrilateral mesh, formed by dividing each triangle of an unstructured triangular mesh into three quadrilaterals for spatial discretization. This research deals with determining a suitable measure of mesh quality based on the physics of the problems being tackled. This is followed by exploring methods to improve the quality of quadrilateral elements obtained from the triangles and thereby improving the overall mesh quality. A series of numerical experiments were designed and conducted for this purpose and the results obtained were tested on different geometries with varying degrees of mesh density.

Paleoenvironmental reconstructions from cave sediments of the Moravian Karst, Czech Republic / by Pavel Sroubek.

This thesis presents a paleoclimatic/paleoenvironmental study conducted on clastic cave sediments of the Moravian Karst, Czech Republic. The study is based on environmental magnetic techniques, yet a wide range of other scientific methods was used to obtain a clearer picture of the Quaternary climate. My thesis also presents an overview of the significance of cave deposits for paleoclimatic reconstructions, explains basic environmental magnetic techniques and offers background information on the study area – a famous karst region in Central Europe with a rich history. In Kulna Cave magnetic susceptibility variations and in particular variations in pedogenic susceptibility yield a detailed record of the palaeoenvironmental conditions during the Last Glacial Stage. The Kulna long-term climatic trends agree with the deep-sea SPECMAP record, while the short-term oscillations correlate with rapid changes in the North Atlantic sea surface temperatures. Kulna Cave sediments reflect the intensity of pedogenesis controlled by short-term warmer events and precipitation over the mid-continent and provide a link between continental European climate and sea surface temperatures in the North Atlantic during the Last Glacial Stage. Given the number of independent climate proxies determined from the entrance facies of the cave and their high resolution, Kulna is an extremely important site for studying Late Pleistocene climate. In the interior of Spiralka Cave, a five meter high section of fine grained sediments deposited during floods yields information on the climatic and environmental conditions of the last millenium. In the upper 1.5 meters of this profile, mineral magnetic and other non-magnetic data indicate that susceptibility variations are controlled by the concentration of magnetite and its magnetic grain size. Comparison of our susceptibility record to the instrumental record of winter temperature anomalies shows a remarkable correlation. This correlation is explained by coupling of the flooding events, cultivation of land and pedogenetic processes in the cave catchment area. A combination of mineral magnetic and geochemical proxies yields a detail picture of the rapidly evolving climate of the near past and tracks both natural and human induced environmental changes taking place in the broader region.