On this page, the current GeoSim students present their projects. Click the names in the matrix to navigate directly to each student's project description!
RA I, Solid Earth
Tobias Willenbockel
Education
Diploma in Mathematics at Freie Universität Berlin. Diploma Thesis Title: 'Clusteranalyse in Netzwerken – Am Beispiel des Variational Bayesian Expectation Maximization Algorithmus'.
Current project with GeoSim (since 01.10.11)
‘Modelling of Seismic Activity via Network and Time Series Analysis’ at Freie Universität Berlin, Department of Mathematics and Computer Science, Institute of Mathematics.
Supervisors: Christof Schütte (FU Mathematics), Frank Scherbaum (UP Geophysics).
The aim of my work is to achieve a better understanding of the dynamics of seismic activity by using a scalefree modelation of earthquakes via a network approach. For capturing the dynamics of seismic activity in such a network, it has to be coarse grained and analyzed mathematically for introducing, in a second step, fitting methods from time series analysis.
With this approach, a mathematical model which will allow probability based short term predictions of seismic activity is to be developed.
Olga Zakharova

Education
Bsc (2005) and MSc (2007) in Plasma Physics at Kazakh National University named after al-Farabi.
Current project with GeoSim (since 01.04.2011)
'Correlations of aftershock productivity according to the tectonic regime'
Supervisors: Sebastian Hainzl (GFZ 2.1), Matthias Holschneider (UP Mathematics)
Earthquake clustering is a first-order feature of seismicity with important consequences for seismic hazard estimation. The dependence of the total aftershock numbers can be described by the function of the magnitude M of the mainshock, N = K exp(alpha*M). We analyze the total aftershock numbers globally and on regional scales.
Conference abstracts
O. Zakharova (2012), Correlations of aftershock productivity and tectonic regime, in IUGG (conference), edited, Edinburgh, UK.
O. Zakharova (2012), Dependence of aftershock parameters on tectonic regime, in AGU Fall Meeting, edited, San Francisco, USA.
Natalia Zamora
Education
BSc Geology, Master in Project Management and MSc in Marine and Coastal Geology.
Current project with GeoSim (since 01.10.11)
'Probabilistic tsunami hazard assessment for Central America' at department 2.5, GFZ. I use statistical techniques and tsunami modeling to determine main tsunamigenic sources, and to quantify tsunami hazard along Middle America Trench.
Supervisors: Andrey Babeyko (GFZ 2.5), Matthias Holschneider (UP Mathematics)
Conference abstracts
N. Zamora (2012), Validation of Hyflux2 numerical code: Nias earthquake 2005, in Tsunami Symposium, edited, Ispra, Italy.
N. Zamora (2012), 1. Seismotectonic and geodynamic data for tsunami modeling in Central America. 2. Tsunami hazard along Nicoya Peninsula., in INQUA Active Tectonics Workshop, edited, Morelia, Mexico.
N. Zamora (2012), Tsunami potential along Nicoya peninsula and first results of tsunami risk perception in Tamarindo community., in Field work in Nicoya Peninsula Costa Rica and work at University of Costa Rica, edited, Along Nicoya Península and San José, Costa Rica.
RA I, Atmosphere-Hydrosphere
Walter Acevedo

Education
MSc in Physics at national University of Colombia (2006), Thesis title: Quantum transport in a spatially periodic magnetic billiard.
Current project with GeoSim (since 01.07.2011)
'Data Assimilation of Paleo-Proxy Records' at the Institute of Meteorology, Freie Universität Berlin.
Supervisors: Ulrich Cubasch (FU Meteorology), Katja Matthes (Uni Kiel) and Sebastian Reich (UP Mathematics).
Objective: Implementation and testing of techniques for assimilating proxy records, particularly tree-ring width chronologies, into global climate models with the aim of generating improved climate datasets for the last centuries.
Conference abstracts
W. Acevedo, U. Cubasch, K. Matthes, and S. Reich (2012), Data assimilation approach to paleoclimate reconstruction using tree-ring growth forward models and ensemble Kalman filter techniques, in Dendrosymposium TRACE 2012, Tree Rings in Archaeology, Climatology and Ecology, edited, p. 52, Potsdam and Eberswalde, Germany.
W. Acevedo, I. Dorado, U. Cubasch, K. Matthes, and S. Reich (2012), Process-based forward modelling of European tree-ring width chronologies using models of the Vaganov-Shashkin family, in 3rd International Conference on Earth System Modelling, edited, pp. 3ICESM-396, Copernicus Publications, Hamburg, Germany.

Nico Becker

Education
BSc in Meteorology at FU Berlin (2008), Thesis: Vergleich verschiedener Konvektionsparametrisierungen anhand von Niederschlägen in Ostafrika.
MSc in Meteorology at FU Berlin (2011), Thesis: Extreme Windereignisse im Mittelmeerraum in regionalen Simulationen mit verschiedener Auflösung.
Current project with GeoSim (since 01.04.2011)
'Error propagation in regional climate models' at the Institute of Meteorology, FU Berlin.
Supervisors: Uwe Ulbrich (FU Meteorology), Rupert Klein (FU Mathematics)
Regional climate modeling has evolved into an important part of climate research. The concept of adding regional details to large scale data - also called dynamical downscaling - implies certain sources of errors, for example the different resolutions at the model boundaries. These errors have to be revealed and quantified, to improve regional model output and its scientific interpretation.
Conference abstracts
N. Becker, (2012), Differences in wind speed and circulation between CCLM simulations and ECHAM5 driving data, in COSMO-CLM User Seminar, edited, Offenbach (Main), Germany.
N. Becker (2012), Wind speed differences between a regional and global climate model simulation and their relation to the circulation, in EGU General Assembly, edited, Vienna.
N. Becker (2012), Wind speed differences between a regional and global climate model simulation and their relation to the circulation, in E2SCMS, edited, Kos, Greece.
RA I, Hydrological Flow
No projects in this domain so far.
RA II, Solid Earth
Øystein Thorden Haug
Education
MSc in physics at the University Of Oslo, Norway (2012), Thesis title: Quantitative experimental modeling of fragmentation during phreatomagmatic and phreatic eruptions. Supervisors: Olivier Galland and Galen Gisler
Current project with GeoSim (since 01.04.12)
'Scale transitions in landslide and other crustal fragmentation processes' at Section 3.1, GFZ Potsdam
Supervisors: Matthias Rosenau, Onno Oncken (GFZ 3.1), Christof Schuette (FU Mathematics)
Landslides are one of the agents that put at risk the infrastructure and the people in mountainous regions. This thesis aims at understanding fundamental aspects of landslide dynamics.
Field observations have shown that the runout distance for small to normal sized landslides (Volume <10^6 m3) can be explained by a simple Coulomb friction law, i.e. their runout can be considered independent of gravity and volume. Large volume landslides (>10^6 m3), on the other hand, show a dependence between runout and volume, suggesting fundamentally different dynamics of small and large landslides. Several mechanisms have been suggested to explain this scale dependent behavior, as of yet, however, no general consensus has been reached.
We study landslides in a controlled laboratory environment, where landslides are modeled by a brittle material that is continuously fragmenting, while moving down a chute. An analogue material has been developed, where the material properties are defined from scaling arguments. The experiments are monitored with a digital video camera and particle image velocimetry (PIV) is used to study the surface dynamics. Additionally, a measure for the fragmentation, such as the fragment size distribution, will be measured.
From this information we aim to quantify the effect the fragmentation process has on the runout of the experimental slides. We hypothesis that the behavior of the avalanche is highly dependent on the amount of fragmentation, meaning that a non-fragmented material (a solid block) can be explained by the simple Coulomb friction law, while a highly fragmented material will have a behavior more similar to a liquid. This change in dynamic behavior occurs during the landslide event either continuously or episodic and is controlled by the fragmentation process. We finally aim to derive scaling laws that links the fragmentation and the flow behavior of the material. Such scaling laws, if validated against natural observations, might be an integral part of hazard analysis as it might be possible to better predict the runout behavior of large landslides.
Conference abstracts
O. Thorden, (2012), Quantitative experimental modeling of fragmentation during phreatic and phreatomagmatic eruptions, in Geomod, edited, Luasanne.
O. Thorden (2012), Quantitative experimental modeling of fragmentation during phreatic and phreatomagmatic eruptions, in EGU General Assembly, edited, Vienna.
Elias Pipping
Education
Diploma in Mathematics (2010) at Free University Berlin.
Current project with GeoSim (since 01.01.2012)
'Variational methods for rate- and state-dependent friction models' at Free University Berlin.
Supervisors: Ralf Kornhuber (FU Mathematics), Onno Oncken (GFZ 3.1).
Rate- and state-dependent friction laws can be used to model the tribological behaviour of a number of materials. They have, so far, been met with little interest from the mathematical community. The aim of this project is to embed problems of rate- and state-dependent friction in a variational framework, in order to open them to modern analytical tools as well as fast and robust numerical algorithms.
Conference abstract
Pipping, E. (2012), Variational methods for rate- and state-dependent friction models, in GeoMod, edited, Lausanne, Switzerland.
RA II, Atmosphere-Hydrosphere
since 01.04.2012: Nawinda Chutsagulprom at the Institute of Mathematics, Potsdam University. Supervisors: Sebastian Reich (UP Mathematics), Peter Nevir (FU Meteorology)
Christoph Glawe
Education
Diplom in Physics (2011) at RWTH Aachen University. Diploma thesis at Department Nonlinear Dynamics of Laser Processing, NLD, and Fraunhofer Institute for Laser Technology, ILT. Thesis title: 'Analysis of Gas Flows in Laser Cutting Processes with Discontinuous-Galerkin-Methods'.
Current project with GeoSim (since 01.10.2011)
'Numerical modeling and analysis of breaking gravity waves using the One-Dimensional Turbulence model' at the 'Juniorprofessur for flow modeling', BTU Cottbus and the Institute of Mathematics, Department of Mathematics and Computer Sciences, Freie Universtät Berlin.
Supervisors: Heiko Schmidt (BTU Cottbus), Rupert Klein (FU Mathematics)
The energy and momentum deposition through breaking atmospheric gravity waves need to be described in weather and climate models. A 3D turbulence model called ODTLES, based on the stochastic One-Dimensional Turbulence (ODT) model is able to describe the fully developed turbulent spectrum from gravity wave length down to Kolmogorov scale with moderate computing demand. The goal is to use these models stand-alone and coupled with other solvers.
Conference abstracts
C. Glawe, (2012), Stochastic Methods in Fluid Mechanics, in Workshop at "Centre International des Sciences Mecaniques", edited, Udine.
C. Glawe, (2012), Iterative Solvers and Parallelisation, in Workshop at "Hochleistungsrechenzentrum Stuttgart" (High Performance Computing Centre Stuttgart), edited, Stuttgart.
C. Glawe, R. Klein, K. Alan, and H. Schmidt (2012), Towards the simulation of gravity waves using the One-Dimensional Turbulence model, in European Geosciences Union General Assembly, edited.
Verena Molina

Education
Diplom in Mathematics (2010), Institute of Mathematics, Department of Mathematics and Computer Sciences, Freie Universität Berlin. Diploma thesis title: "Multiple-Scales Analysis of Deep Convection in a Squall Line"
Current project with GeoSim (since 01.04.2011)
'Asymptotic and Numerical Modeling of Cloud Processes in the Atmosphere' at the Institute of Mathematics, Department of Mathematics and Computer Sciences, Freie Universtät Berlin.
Supervisors: Rupert Klein (FU Mathematics), Uwe Ulbrich (FU Meteorology)
Derivation and validation of consistent mathematical models of dynamical and microphysical processes in earth's atmosphere.
Annette Müller

Education
Diplom in Mathematics (2011) at TU Berlin
Current project with GeoSim (since 01.04.2012)
'A transition scenario from three to two-dimensional regimes of turbulence based on the energy-vorticity theory of fluid mechanics using an infinite cascade of helicity related constraints' at the Institute of Meteorology, FU Berlin.
Supervisors: Peter Nevir (FU Meteorology), Sebastian Reich (UP Mathematics)
The aim of this project is to combine different physical theories of fluid mechanics focusing on the energy-vorticity theory (EVT) to model a transition scenario between two-dimensional and three-dimensional turbulence theories.
Especially from the meteorological point of view, the challenge is to combine the different 2D-and 3D-turbulence models. A tornado is an example of a meteorological event appearing in the 3D convective scale. Those events are characteristic for three-dimensional turbulence and can be modelled as vortex filaments. Two-dimensional turbulence theory contains larger-scale processes. For example low-pressure and high-pressure areas. These synoptic processes can be described by point vortices.
The dynamics of point vortices and vortex filaments can be described by the EVT. The EVT focuses on two conserved quantities, equally important for the nonlinear evolution of fluids. The first conserved quantity in two and three-dimensional turbulence theory is the kinetic energy. The second one, e.g. enstrophy (2D) and helicity (3D), determines the different topological behaviour of vortex motion.
The main part of the project is to understand the topological constraints of vorticity conserved quantities and to apply the results to model atmospheric motion.
From the meteorological point of view, modelling the transition scenario will connect the 3D convective scale with the 2D synoptic scale.
Lisa Schielicke

Education
Diplom in Meteorology at Freie Universität Berlin. Thesis title: On the theory of frequency distributions of tornadoes, hurricanes and low pressure systems of the mid-latitudes.
Current project with GeoSim (since 01.10.2011)
'Scale-depending analysis of vortex structures in numerical simulations and observations under different climate conditions' at the Institute of Meteorology, FU Berlin.
Supervisors: Uwe Ulbrich (FU Meteorology), Peter Nevir (FU Meteorology), Rupert Klein (FU Mathematics)
Atmospheric vortices exist on different scales in time and space and can vary considerably in their intensity. This diversity in the properties of vortices mirror the complexity of the flow field and the underlying processes. Modified climate conditions change the conditions for the genesis of vortices. This should be reflected in a change of the properties of the whole vortex ensemble. In my GeoSim-Project I will study the influence of modified forcing-dissipation settings on the vortex properties and on their frequency as well as their effect on the interactions between vortex structures on different scales.
Publications
L. Schielicke and P. Névir (2011), Introduction of an atmospheric moment combining Eulerian and Lagrangian aspects of vortices: Application to tornadoes, Atm. Res, 100, 357-365.
L. Schielicke and P. Névir (2013), Comprehensive analysis of tornado statistics in comparison to earthquakes: intensity and temporal behaviour, Nonlin. Processes Geophys., 20, 47-57., doi:10.5194/npg-20-47-2013
Conference abstracts
L. Schielicke and P. Névir (2011), Power laws in the atmosphere, in 6th European Conference on Severe Storms, edited, Palma de Mallorca, Spain.
L. Schielicke (2012), Testing a different method for identifying and tracking atmospheric vortices, in DAMES 2012: Data analysis and modeling in Earth sciences, edited, Potsdam, Germany.
L. Schielicke, P. Névir, U. Ulbrich, and R. Klein (2012), Testing a different method for identifying and tracking atmospheric vortices, in 3rd International Conference on Data Analysis and Modeling in Earth Sciences (DAMES 2012), edited, Potsdam, Germany.
RA II, Hydrological Flow
Matthias Munz

Education
Geoecology at Potsdam University. Majors: Modelling in Geoecology und Hydrology. Diploma thesis (2011): "Sand box experiments to evaluate the influence of subsurface temperature probe installation on temperature based water flux calculation"
Current project with GeoSim (since 01.04.2011)
'Coupled water and heat transport modelling at the surface water-groundwater interface - from a highly instrumented stream stretch to the entire stream' at the Research Group “Water and Matter Transport in Complex Landscapes”, Institute of Earth and Environmental Science, Potsdam University.
Supervisors: Sascha Oswald (UP Hydrology), Ralf Kornhuber (FU Mathematics)
Surface waters and groundwater are the interconnected parts of stream catchments. They are coupled by variable water, solute and heat exchange processes through streambed sediments. In consequence there is an increasing multi-disciplinary interest in investigating hydrological, biogeochemical and ecological processes at the direct interface between river and aquifer (hyporheic zone), because they affect ecosystem functioning in natural and human impacted water resources systems and water quality as well as water discharges.
The goal of this research project will be to simulate hyporheic zone temperature patterns to advance the understanding of spatially and temporally highly variable thermal exchange between rivers and groundwater, using existing numerical models -HydroGeoSphere. Received subsurface temperature patterns will be used to delineate river aquifer interaction zones and could be linked to microbial activity attached to sediments and to nutrient cycles in the subsurface as these ecological and biogeochemical processes are temperature sensitive.
Publications
M. Munz, S. Krause, C. Tecklenburg, A. Binley, 2011. Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns. Hydrological Processes. Published online: 12 May 2011. DOI: 10.1002/hyp.8080.
Conference abstracts
M. Munz, C. Schmidt, J. Fleckenstein, and S. Oswald (2012), Coupled water and heat transport modelling to evaluate hyporheic temperatures and exchange fluxes., in European Geosciences Union - General Assembly, edited, Wien, Austria.
M. Munz, C. Schmidt, J. Fleckenstein, and S. Oswald (2012), Modelling spatial and temporal variability of surface water - groundwater heat exchange along a lowland river reach., in AGU Fall Meeting, edited, San Francisco, USA.

Ulrich Wilbrandt
Education
MSc in mathemathics / scientific computing at FU Berlin.
Current project with GeoSim (since 01.10.2011)
We study coupled Stokes and Darcy equations which for example model surface water / ground water interactions. We will particularly focus on implementation, analysis, inf-sup stable versus stabilized discretizations. Furthermore we will apply these methods to real world problems from geophysics. This includes simulations in three space dimensions and time dependencies. We hope to compare our results to existing commercial codes available.
Supervisors: Volker John (WIAS), Sascha Oswald (UP Meteorology)
RA III, Solid Earth
Camilla Cattania
Education
BA. and MSc. (2011) in Experimental and Theoretical Physics at the University of Cambridge.
Current project with GeoSim (since 01.11.2011)
'Improvement of seismicity models based on Coulomb stress interactions and rate-state dependent friction'
Supervisors: Sebastian Hainzl (GFZ 2.1), Frank Roth (GFZ 2.1)
Physics based aftershock models consider changes in the stress field caused by a mainshock (Coulomb stress), combined with the frictional response of the surrounding faults (Rate-and-State friction). Several sources of uncertainties are present in Coulomb stress calculations, such as errors in the slip models and the choice of receiver faults: one goal of my project is to improve forecast models by correctly accounting for these uncertainties. While Coulomb stress is widely recognized as playing a role in aftershock triggering, the relative effect of co-seismic and post-seismic slip is subject of debate: the model is therefore extended to include an arbitrary (time dependent) stressing history, and the impact of afterslip on the forecasted seismicity is analyzed.
Conference abstracts
C. Cattania and S. Hainzl (2012), Uncertainties in Coulomb Rate-and-State seismicity models: effects of receiver fault orientation, in IUGG conference, edited, Edimburgh.
C. Cattania, S. Hainzl, and F. Roth (2012), Towards more realistic Coulomb-Rate-and-State seismicity models: effect of receiver fault orientation and afterslip in Fall Meeting, AGU, edited, pp. S54B-03, San Francisco, California.
Juliane Dannberg

Education
Diploma in Geophysics (2012) at Friedrich-Schiller-University of Jena, Germany. Thesis Title: “Modelling of thermochemical mantle plumes and their influence on surface topography”
Current Project with GeoSim (since 01.04.2012)
“Dynamics of thermo-chemical mantle plumes and plume-lithosphere interaction” at Section 2.5 Geodynamic Modelling, GFZ German Research Centre for Geosciences
Supervisors: Stephan V. Sobolev (GFZ) and Volker John (WIAS)
Geodynamic processes occur on various scales both temporally and spatially. Mantle plumes are a link between these scales: They form a part of large-scale mantle convection as well as effecting surface magmatism. According to widely recognised models, large igneous provinces (LIPs) develop as a result of plumes ascending from the core-mantle boundary and the associated massive melting when the plume head reaches the base of the lithosphere. This triggers the eruption of an immense amount of magma in less than a million years. While classical models of purely thermal plumes can explain most of the observations in a first-order approximation, recent studies have shown that a thermo-chemical plume, which contains a fraction of recycled oceanic crust, fits the observations much better. The objective of this project is to investigate the dynamics of ascending thermo-chemical mantle plumes and their interaction with the lithosphere, including their impact on the Earth’s surface such as melting events and their influence on surface topography. Both the development of the plume from a thermo-chemical boundary layer and its movement in the whole mantle employing a temperature- and depth-dependent rheology as well as a realistic chemical density contrast will be considered. For this purpose a three-dimensional numerical model will be developed. Challenges include the implementation of adaptive mesh refinement, melt generation and transport as well as a precise method to calculate the surface topography.
Conference abstracts
J. Dannberg and S. V. Sobolev (2012), Dynamics of thermo-chemical mantle plumes and their influence on surface topography, in EGU General Assembly, edited, Vienna.
J. Dannberg and S. V. Sobolev (2012), Dynamics of thermo-chemical mantle plumes and their impact on surface topography in Geomod, edited.
J. Dannberg and S. V. Sobolev (2012), Dynamics of low-buoyancy mantle plumes, in Deutscher Geodynamik-Workshop, edited.
Annabel Händel

Education
Diploma in Geophysics at the Institute of Earth- and Environmental Sciences of the University of Potsdam (2012), Thesis title: “Data driven selection and ranking of ground-motion prediction equations for subduction zones”.
Current project with GeoSim (since 01.04.2012)
“Determination of attenuation properties from ambient vibrations” at the Institute of Earth- and Environmental Sciences, University of Potsdam.
Supervisor: Matthias Ohrnberger (UP Geophysics & Seismology), Frank Scherbaum (UP Geophysics & Seismology), Volker John (WIAS)
Understanding how earthquake waves are attenuated in the ground is crucial for any seismic hazard analysis. Especially the high-frequency filtering effect of the uppermost layers (described by the parameter kappa) can have a considerably influence on the proper prediction of future ground motions and needs to be assessed.
The main goal of this reserch project is the evolution of a non-destructive method to determine attenuation from noise data. We will invesitgate how seismic interferometry can be used to extract amplitude and attenuation information from ambient vibrations. We will further simulate a homogeneous scattering medium using finite differences to validate the obtained results and to better understand the mechanisms of seismic wave attenuation in the ground.
Conference abstracts
C. Beauval, F. Cotton, N. Abrahamson, N. Theodulidis, E. Delavaud, L. Rodriguez, F. Scherbaum, and A. Haendel (2012), Regional differences in subduction ground motions, in World Conference on Earthquake Engineering (WCEE), edited, Lisboa, Portugal.
A. Haendel, N. Kuehn, and F. Scherbaum (2012), Generation of a mixture model ground-motion prediction equation for Northern Chile, in AGU Fall Meeting, edited, San Francisco, USA.
Shaoyang Li

Education:
MSc in Cartography and GIS (2012) at Beijing Normal University (Beijing, China); MSc in Applied Geophysics (2011) at Université Pierre et Marie Curie (Paris VI, France)
Current project with GeoSim (since 01.09.2012):
“Geomechanical modeling of earthquake cycles in the north Chilean seismic gap”
Supervisors: Marcos Moreno (GFZ 3.1) Onno Oncken (GFZ 3.1)
Focusing on North Chilean subduction seismic gap, my proposed project is composed of two main parts, one that establishes the kinematics at the North Chilean plate margin using inversion of GPS data and other information (seismology, tectonics, InSAR, etc.) with Finite Element Models and the second part is dynamic modeling of the system behavior in the seismic cycles. The latter involves analogue experiments as well as numerical simulation, both constraint by the then established kinematics of the North Chilean plate boundary with realistic geometries and laboratory-derived constitutive fault friction rheologies.
Christian Molkenthin
Education
Geoecology and Geophysics at University of Potsdam. Diploma of Geoecology in 2011.
Current project with GeoSim (since 1.10.2011)
'Sensitivity studies in probabilistic seismic hazard analysis (PSHA).' Institute of Earth and Environmental Sciences, University of Potsdam.
Supervisors: Frank Scherbaum (UP Geophysics), Prof. Dr. Sebastian Reich (UP Mathematics)
Seismic hazard analysis is an important task for reducing the consequences due to devastating earthquakes in vulnerable regions worldwide. The appropriate method for the assessment is probabilistic seismic hazard analysis (PSHA).
PSHA treats ground motion as a random variable and takes into account all possible earthquake occurrences and ground motion scenarios. Therefore, the probabilistic approach enables a complete quantitative consideration of uncertainty. The outcome is a hazard curve quantifying the annual rate of exceedance (or its reciprocal: return period) for different levels of a selected ground motion parameter, such as peak ground acceleration (PGA) or spectral acceleration (SA).
PSHA is based on input assumptions about several parameters, for example characterizing the relative frequencies of occurrence of different earthquakes (e.g. a-value, b-value), the distribution of epicenters or hypocenters (e.g. zonation) and the attenuation and scaling of the amplitude of ground motions (e.g. fc, kappa). It is not always obvious how and to what extent these parameters influence the shape of the hazard curve. In order to determine the impact of individual parameters and also to identify those uncertainties that can be reduced, each PSHA should be accompanied by a sensitivity study. Due to the lack of appropriate tools, this is currently not the case.
The aim of this dissertation is to approach sensitivity analyses in PSHA in a theoretically sound and systematic way, providing quantitative estimates of sensitivity with respect to every relevant parameter.
The approach in favor is known as 'adjoint model' which has been widely used within complex meteorological and oceanographic simulation models to asses the parameter driven uncertainty in model predictions. The goal of the thesis is to implement the 'adjoint model' approach for PSHA and explore its potential for large scale seismic hazard problems.
Elvira Mulyukova

Education
MSc in Physics (2011), University of Oslo, Norway. Thesis title: Thermomechanical Modeling of stability of Large Low Shear Velocity Provinces. Supervisors: Marcin Dabrowski, Trond H. Torsvik, Daniel W. Schmid.
Current project with GeoSim (since 01.10.11)
'Thermochemical Modeling: Stability of the Large Low Shear Velocity Provinces' at Section 2.5, Geodynamic Modelling, GFZ.
Supervisors: Bernhard Steinberger (GFZ 2.5), Volker John (WIAS)
Tectonic plate motions result in many important geodynamic processes, including formation of mountain belts, earthquakes and volcanic eruptions. Mechanisms that drive plate tectonics remain uncertain, but the convective dynamics of the mantle have been suggested to play an important role. The focus of my PhD project is to study processes that take place in the deeper mantle.
One of the most robust results from tomographic studies is the existence of two antipodally located Large Low Shear Velocity Provinces (LLSVPs) at the base of the mantle. Reconstruction of the eruption sites of Large Igneous Provinces and hotspot volcanoes of the last 200 Ma has shown that these project radially downward to the margins of the LLSVPs. This has led to inferences that the LLSVPs are stable, long-lived and impose the planform of flow in the mantle and of plate tectonics at the surface.
The objective of my M.Sc. project was to study the inferred gravitational and lateral stability of the LLSVPs and their relation to the generation of plumes. I developed a thermomechanical FEM code for simulating convection of a viscous Newtonian fluid in a two-dimensional rectangular box. The model included a heterogeneous basal layer, which density and viscosity were varied in order to study their effects on the gravitational stability and topography of the layer.
One of the objectives of my PhD project is to further develop the thermomechanical code for simulating vigorous convection of the mantle with introduced heterogeneities. The depth-dependence of properties such as viscosity, thermal expansivity, and conductivity will be included, as well as the discontinuities associated with chemical or phase variations.
The geodynamic model will be constrained by results from several disciplines: seismology, experimental mineral physics and field observations. Numerical challenges include implementation of realistic boundary conditions, developing a robust non-diffusive advection solver, and devising an automatic grid refinement strategy.
Conference abstracts
E. Mulyukova, M. Dabrowski, and B. Steinberger (2012), Numerical Modeling of the Deep Mantle Convection: Accuracy of the Tracer Methods. Abstract, EGU general assembly.
E. Mulyukova (2012), Links between Plumes and Large Igneous, in Course on Mineralogy, Chemistry and Dynamics of Earth and Other Planets Provinces., edited, Oslo, Norway.
E. Mulyukova (2012), Numerical Modelling of Deep Mantle Flow: Accuracy of the Marker-Methods, in EGU, edited, Vienna, Austria.
E. Mulyukova (2012), Numerical Modelling of Deep Mantle Flow:Thermochemical Evolution., in Geomod, edited, Lausanne, Switzerland.
E. Mulyukova (2012), Numerical modelling of Deep Mantle Flow, in Geodynamic workshop, edited, Wandlitz, Germany.
Antonia Oelke

Education
Diplom in Physics (2010) at Freie Universität Berlin. MSc thesis at Department of Geophysics under supervision of Dr. Stefan Buske (now at Freiberg University). Thesis title: 'Seismic Imaging of the Whataroa'98 Reflection Profile across the Alpine Fault, South Island, New Zealand '.
Current project with GeoSim (since 01.04.2011)
'Reflection coefficients at hydraulic fractures' at the Department of Geophysics, Freie Universität Berlin.
Supervisors: Serge Shapiro (FU Geophysics) and Ralf Kornhuber (FU Mathematics)
Conference abstracts
A. Oelke, S. Buske, and S. Bannister (2010), Seismic Imaging of the Alpine Fault, New Zealand, at Whataroa River., in the 14th International Symposium on "Deep Seismic Profiling of the Continents and their Margins", edited, Cairns, Australia.
A. Oelke., D. Alexandrov, I. Abakumov, V. Troyan, B. Kashtan, and S. Shapiro (2011), Reflection Coefficients at a Thin Fluid Layer as a Model of a Hydraulic Fracture. , in 73rd EAGE Conference & Exhibition incorporating SPE EUROPEC 2011 edited, p. 353, Extended abstract, Vienna.
A. Oelke, S. Bannister, and B. S. (2011), Seismic Imaging of the Alpine Fault at Whataroa (NZ), in Fall Meeting, AGU, edited, pp. T13G-06, San Francisco, California.
A. Oelke O. Krüger, S. Shapiro, and I. Abakumov (2012), Modeling Reflection Coefficients at a Thin Fluid Layer Representing a Hydraulic Fracture, in SEG Technical Program Expanded, edited, pp. 1-5, Abstracts 2012. link to abstract
A. Oelke (2012), Reflection coefficients at a thin fluid layer as a model of a hydraulic fracture (& Approximate reflection coefficients at a fracture for microseismic needs), in 5th Saint-Petersburg International Conference & Exhibition 2012, edited, Saint Petersburg, Russia.
A. Oelke (2012), in SEG Annual Meeting, edited, Las Vegas, Nevada, USA.
A. Oelke, D. Alexandrov, I. Abakumov, R.Shigapov, N. Troyan, B. Kashtan and S.Shapiro (2012), Reflection coefficients at a thin fluid layer as a model of a hydraulic fracture. 5th Saint-Petersburg International Conference & Exhibition 2012 Saint-Petersburg, Extended Abstract, P069.
Shigapov, R., A. Oelke, N. Troyan, B. Kashtan, and S. Shapiro (2012), Approximate reflection coefficients at a fracture for microseismic needs, in 5th Saint-Petersburg International Conference & Exhibition edited, p. P068, Saint-Petersburg.
Sanjay Singh Bora
Education
M.Tech (Computational Seismology) Indian Institute of Technology Kharagpur, India. Title of the master thesis: “Deterministic Seismic Hazard Assessment in India & Adjoining Regions and A Probabilistic Pilot Study in Gujarat, Western India.”
Current project with GeoSim (since 1.10.2011)
‘Regionally adaptable Ground Motion Prediction Equations’ at the Institute of Earth and Environmental Sciences, University of Potsdam.
Supervisors: Frank Scherbaum (UP Geophysics), Sebastian Reich (UP Mathematics)
The precise prediction of ground motion for future earthquake scenarios is a key issue in seismic hazard analysis. In this context empirical Ground Motion Prediction Equations (GMPEs) are used. The GMPEs are derived by regression analysis from observed ground motion data such as peak ground acceleration (PGA) or response spectral values against a chosen set of predictor values. There are many regions in the world where the number of the earthquake recordings are not sufficient in quantity and quality to generate endemic GMPEs particularly in intra-plate regions e.g. Gujarat western India. As a result GMPEs have to be imported from other regions. That requires modification to accommodate changes in source and path characteristics such as stress drop, quality factor, and kappa values between the host region and the target region. The idea for my PhD thesis is to generate GMPEs with a different approach using the linear system theory. I will use the Random Vibration Theory (RVT) to calculate the response spectrum from Fourier spectrum and duration model of ground motion. The changes in stress drop, kappa, and quality factor values will be accommodated by adding an appropriate frequency response function in Fourier spectrum model. Consequently the Fourier spectrum model will be combined with the duration model to give the modified response spectrum model. Since the Fourier spectrum model is considered as representing, the actual physics of the process, pertaining to some assumptions. Therefore the proposed research work is believed to add new insights in seismic hazard analysis by developing a response spectrum model through Fourier spectrum model.
Conference abstracts
S. Singh Bora, (2012), OralGround Motion Prediction: A Fourier Spectrum Approach, GEOSIM Fall Workshop, Golm
S. Singh Bora, F. Scherbaum, N. Kuehn and P. Stafford (2013), Fourier Spectral- and Duration Models for the Generation of Response Spectra Adjustable to Different Source-, Path- and Site conditions, SSA Annual Meeting, Salt Lake City, USA. Winner of best student poster award (top 10%)
RA III, Atmosphere-Hydrosphere
Linda Michalk
Education
Diplom in mathematics at Freie Universität Berlin. Thesis title: "Numerical methods for the linear advection equation: Plateaus vs. Extrema"
Current project with GeoSim: (since 01.10.2011)
"Advanced finite volume schemes numerics for advection dominated problems" at the Institute of Mathematics, Department of Mathematics and Computer Sciences, Freie Universtät Berlin
Supervisors: Rupert Klein (FU Mathematics), Uwe Ulbrich (FU Meteorology)
The goal of my work is to compare existing numerical methods, find improvements or develope new ones, that yield satisfying results for solving geophysical problems - in particular advection dominated problems.
Jana de Wiljes

Education
Diplom in Mathematics at Freie Universität Berlin. Thesis title: Adopting a Bayesian framework to multidimensional cluster modeling
Current project with GeoSim (since 01.04.2011)
‘Discrete-continuous non-stationary methods of time series analysis of coupled economics and climate processes’ at the Institute of Mathematics, Department of Mathematics and Computer Sciences, Freie Universtät Berlin
Supervisors: Rupert Klein (FU Mathematics), Illia Horenko (Uni Lugano, CH), Maik Thomas (GFZ 1.3)
My current work aims at extending the existent methods of time series analysis towards the hierarchical non-stationary and non-homogenous space-time clustering of geophysical data, with a special emphasis on the sea ice coverage data. The main mathematical target is to develop a data-based percolation model for the arctic ice dynamics using the information gained from the analysis of spatiotemporal arctic satellite data. The main geophysical target of my research will be to use this stochastic model to statistically investigate the influence of various global and local climate factors on the dynamics of the ice cover.
RA III, Hydrological Flow
Catherine Abon
Education
MSc in Geology at the University of the Philippines Diliman, Quezon City
Current project with GeoSim (since 01.10.11)
‘Hydrologic simulations using rainfall radar data for flood forecasting in Marikina River Basin, Philippines’ at the Hydrology and Climatology group, Institute of Earth and Environmental Science, Potsdam University.
Supervisors: Axel Bronstert (UP Hydrology), Ralf Kornhuber (FU Berlin)
Improving flood preparedness through early warning systems is a necessity for countries like the Philippines where flooding is frequent. Particularly for river basins that host huge settlements such as the Marikina River Basin (MRB), an early flood forecast can potentially reduce the hazards associated with floods. The MRB has an area of 535 km2, and for huge typhoons such as Ketsana in 2009 the reaction time is four to five hours. For an early warning system, a hydrological model will be implemented which is capable of producing flood forecasts that can give a lead within the reaction time of the basin. The model uses a flexible modeling framework, the framework Eco-Hydrological Simulation Environment (ECHSE), in order to consider specific hydrological processes in the catchment. Beyond, I will evaluate the potential of rainfall radar data as precipitation forcing—a novel application in the Philippines as a new network of C- and S-band radars is just being implemented. For this purpose, the new Open-Source radar processing library wradlib will be employed.
Conference abstracts
C. Abon, (2012), Hydrologic simulations using rainfall radar data for flood forecasting in Marikina River Basin, Philippines, in EGU 2012, edited, Vienna, Austria.
Klaus Vormoor
Education
Diploma in Geography at the University of Bonn in 2009. Thesis: "Regionalization of runoff data in Norway."
Current project with GeoSim (since 01.04.2012)
“Uncertainties in climate change impacts on hydrological extremes”, Hydrology and Climatology Group, Institute for Earth and Environmental Science, University of Potsdam
Supervisors: Axel Bronstert (UP Hydrology), Sebastian Reich (UP Mathematics), Deborah Lawrence (external, Norwegian Water Resources and Energy Directorate)
Projections on the hydrologic impacts of climate change at the catchment scale are highly uncertain – especially for extremes. Uncertainties emerge from various sources and accumulate throughout the model chain, beginning with large scale Global Circulation Models estimating changes in climate variables based on emission scenarios, over dynamical or statistical downscaling techniques intending to adjust this information to smaller scales, and finally to the hydrological models which introduce uncertainties due to model parameterization and model structure. This study will simulate changes in flood frequency in selected Norwegian catchments due to climate change and focus on the uncertainties that are affiliated with these projections. A special emphasis will be on the uncertainties that arise from statistical downscaling techniques and hydrologic model parameterizations.
Conference abstract
K. Vormoor and T. Skaugen (2012), Temporal disaggregation of daily meteorological grid data, in EGU General Assembly, edited, Vienna.

News
5th GeoSim application round for 10 PhD positions has now been closed. We will open one more application round this fall.
Upcoming seminar on Friday, 7 June, by dr. S Hergarten (U. Freiburg): "Complexity in geosciences". Agenda.
Sanjay Singh Bora received an award for his poster presentation (best 10%) during the SSA annual meeting in Salt Lake City, 17-19 April 2013

