Tsunami simulation and detection using global navigation satellite system reflectometry

In this thesis, research for tsunami remote sensing using the Global Navigation Satellite System-Reflectometry (GNSS-R) delay-Doppler maps (DDMs) is presented. Firstly, a process for simulating GNSS-R DDMs of a tsunami-dominated sea sur- face is described. In this method, the bistatic scattering Zavorotny-Voronovich (Z-V) model, the sea surface mean square slope model of Cox and Munk, and the tsunami- induced wind perturbation model are employed. The feasibility of the Cox and Munk model under a tsunami scenario is examined by comparing the Cox and Munk model- based scattering coefficient with the Jason-1 measurement. A good consistency be- tween these two results is obtained with a correlation coefficient of 0.93. After con- firming the applicability of the Cox and Munk model for a tsunami-dominated sea, this work provides the simulations of the scattering coefficient distribution and the corresponding DDMs of a fixed region of interest before and during the tsunami. Fur- thermore, by subtracting the simulation results that are free of tsunami from those with presence of tsunami, the tsunami-induced variations in scattering coefficients and DDMs can be clearly observed. Secondly, a scheme to detect tsunamis and estimate tsunami parameters from such tsunami-dominant sea surface DDMs is developed. As a first step, a procedure to de- termine tsunami-induced sea surface height anomalies (SSHAs) from DDMs is demon- strated and a tsunami detection precept is proposed. Subsequently, the tsunami parameters (wave amplitude, direction and speed of propagation, wavelength, and the tsunami source location) are estimated based upon the detected tsunami-induced SSHAs. In application, the sea surface scattering coefficients are unambiguously re- trieved by employing the spatial integration approach (SIA) and the dual-antenna technique. Next, the effective wind speed distribution can be restored from the scat- tering coefficients. Assuming all DDMs are of a tsunami-dominated sea surface, the tsunami-induced SSHAs can be derived with the knowledge of background wind speed distribution. In addition, the SSHA distribution resulting from the tsunami-free DDM (which is supposed to be zero) is considered as an error map introduced during the overall retrieving stage and is utilized to mitigate such errors from influencing sub- sequent SSHA results. In particular, a tsunami detection procedure is conducted to judge the SSHAs to be truly tsunami-induced or not through a fitting process, which makes it possible to decrease the false alarm. After this step, tsunami parameter estimation is proceeded based upon the fitted results in the former tsunami detec- tion procedure. Moreover, an additional method is proposed for estimating tsunami propagation velocity and is believed to be more desirable in real-world scenarios. The above-mentioned tsunami-dominated sea surface DDM simulation, tsunami detection precept and parameter estimation have been tested with simulated data based on the 2004 Sumatra-Andaman tsunami event.

Forecasting the term structure of volatility of crude oil price changes

Peer reviewed

Boron isotope and B/Ca ratios of planktonic foraminifera

Here a new analytical methodology is described for measuring the isotopic composition of boron in foraminifera using multicollector inductively coupled plasma mass spectrometry (MC-ICPMS). This new approach is fast (~10 samples analysed in duplicate per analytical session) and accurate (to better than 0.25 per mil at 95% confidence) with acceptable sample size requirements (1-3 mg of carbonate). A core top calibration of several common planktic and two benthic species from geographically widespread localities shows a very close agreement between the isotopic composition measured by MC-ICPMS and the isotopic composition of B(OH)-4 in seawater (as predicted using the recently measured isotopic equilibrium factor of 1.0272) at the depth of habitat. A down core and core top investigation of boron concentration (B/Ca ratio) shows that the partition coefficient is influenced by [CO2-3] complicating the application of this proxy. Nevertheless, it is demonstrated that these two proxies can be used to fully constrain the carbonate system of surface water in the Caribbean Sea (ODP Site 999A) over the last 130 kyr. This reconstruction shows that during much of the Holocene and the last interglacial period surface water at Site 999A was in equilibrium with the atmosphere with respect to CO2. During the intervening colder periods although the surface water pCO2 was lower than the Holocene, it was a minor to significant source of CO2 to the atmosphere possibly due to either an expansion of the eastern equatorial Atlantic upwelling zone, or a more local expansion of coastal upwelling in the southern Caribbean. Such reorganisation of the oceanic carbonate system in favour of a larger source of CO2 to the atmosphere from the equatorial ocean may require mechanisms responsible for lowering atmospheric CO2 during glacial periods to be more efficient than previously supposed.

Privileged Words and Sturmian Words

This dissertation has two almost unrelated themes: privileged words and Sturmian words. Privileged words are a new class of words introduced recently. A word is privileged if it is a complete first return to a shorter privileged word, the shortest privileged words being letters and the empty word. Here we give and prove almost all results on privileged words known to date. On the other hand, the study of Sturmian words is a well-established topic in combinatorics on words. In this dissertation, we focus on questions concerning repetitions in Sturmian words, reproving old results and giving new ones, and on establishing completely new research directions. The study of privileged words presented in this dissertation aims to derive their basic properties and to answer basic questions regarding them. We explore a connection between privileged words and palindromes and seek out answers to questions on context-freeness, computability, and enumeration. It turns out that the language of privileged words is not context-free, but privileged words are recognizable by a linear-time algorithm. A lower bound on the number of binary privileged words of given length is proven. The main interest, however, lies in the privileged complexity functions of the Thue-Morse word and Sturmian words. We derive recurrences for computing the privileged complexity function of the Thue-Morse word, and we prove that Sturmian words are characterized by their privileged complexity function. As a slightly separate topic, we give an overview of a certain method of automated theorem-proving and show how it can be applied to study privileged factors of automatic words. The second part of this dissertation is devoted to Sturmian words. We extensively exploit the interpretation of Sturmian words as irrational rotation words. The essential tools are continued fractions and elementary, but powerful, results of Diophantine approximation theory. With these tools at our disposal, we reprove old results on powers occurring in Sturmian words with emphasis on the fractional index of a Sturmian word. Further, we consider abelian powers and abelian repetitions and characterize the maximum exponents of abelian powers with given period occurring in a Sturmian word in terms of the continued fraction expansion of its slope. We define the notion of abelian critical exponent for Sturmian words and explore its connection to the Lagrange spectrum of irrational numbers. The results obtained are often specialized for the Fibonacci word; for instance, we show that the minimum abelian period of a factor of the Fibonacci word is a Fibonacci number. In addition, we propose a completely new research topic: the square root map. We prove that the square root map preserves the language of any Sturmian word. Moreover, we construct a family of non-Sturmian optimal squareful words whose language the square root map also preserves.This construction yields examples of aperiodic infinite words whose square roots are periodic.

Release currents of IP₃ receptor channel clusters and concentration profiles

We simulate currents and concentration profiles generated by Ca2+ release from the endoplasmic reticulum (ER) to the cytosol through IP3 receptor channel clusters. Clusters are described as conducting pores in the lumenal membrane with a diameter from 6 nm to 36 nm. The endoplasmic reticulum is modeled as a disc with a radius of 1–12 mm and an inner height of 28 nm. We adapt the dependence of the currents on the trans Ca2+ concentration (intralumenal) measured in lipid bilayer experiments to the cellular geometry. Simulated currents are compared with signal mass measurements in Xenopus oocytes. We find that release currents depend linearly on the concentration of free Ca2+ in the lumen. The release current is approximately proportional to the square root of the number of open channels in a cluster. Cytosolic concentrations at the location of the cluster range from 25 μM to 170 μM. Concentration increase due to puffs in a distance of a few micrometers from the puff site is found to be in the nanomolar range. Release currents decay biexponentially with timescales of < 1 s and a few seconds. Concentration profiles decay with timescales of 0.125–0.250 s upon termination of release.

Periodic emission of droplets from an oscillating electrified meniscus of a low-viscosity, highly conductive liquid

The generation of identical droplets of controllable size in the micrometer range is a problem of much interest owing to the numerous technological applications of such droplets. This work reports an investigation of the regime of periodic emission of droplets from an electrified oscillating meniscus of a liquid of low viscosity and high electrical conductivity attached to the end of a capillary tube, which may be used to produce droplets more than ten times smaller than the diameter of the tube. To attain this periodic microdripping regime, termed axial spray mode II by Juraschek and Röllgen [R. Juraschek and F. W. Röllgen, Int. J. Mass Spectrom. 177, 1 (1998)], liquid is continuously supplied through the tube at a given constant flow rate, while a dc voltage is applied between the tube and a nearby counter electrode. The resulting electric field induces a stress at the surface of the liquid that stretches the meniscus until, in certain ranges of voltage and flow rate, it develops a ligament that eventually detaches, forming a single droplet, in a process that repeats itself periodically. While it is being stretched, the ligament develops a conical tip that emits ultrafine droplets, but the total mass emitted is practically contained in the main droplet. In the parametrical domain studied, we find that the process depends on two main dimensionless parameters, the flow rate nondimensionalized with the diameter of the tube and the capillary time, q, and the electric Bond number BE, which is a nondimensional measure of the square of the applied voltage. The meniscus oscillation frequency made nondimensional with the capillary time, f, is of order unity for very small flow rates and tends to decrease as the inverse of the square root of q for larger values of this parameter. The product of the meniscus mean volume times the oscillation frequency is nearly constant. The characteristic length and width of the liquid ligament immediately before its detachment approximately scale as powers of the flow rate and depend only weakly on the applied voltage. The diameter of the main droplets nondimensionalized with the diameter of the tube satisfies dd≈(6/π)1/3(q/f)1/3, from mass conservation, while the electric charge of these droplets is about 1/4 of the Rayleigh charge. At the minimum flow rate compatible with the periodic regimen, the dimensionless diameter of the droplets is smaller than one-tenth, which presents a way to use electrohydrodynamic atomization to generate droplets of highly conducting liquids in the micron-size range, in marked contrast with the cone-jet electrospray whose typical droplet size is in the nanometric regime for these liquids. In contrast with other microdripping regimes where the mass is emitted upon the periodic formation of a narrow capillary jet, the present regime gives one single droplet per oscillation, except for the almost massless fine aerosol emitted in the form of an electrospray.

Photocurrents in P3MeT Schottky barrier diodes

The photocurrent action spectrum of a Schottky diode formed from electrodeposited poly(3-methylthiophene) is shown to follow closely the UV-visible absorption spectrum. At low forward bias, the peak photocurrent asymptotes to the expected square-root dependence on total potential, V-t, across the depletion region. At high reverse bias the superlinear dependence of the photocurrent on V-t suggests that internal photoemission from the rectifying aluminium electrode may be the dominant process.

Photocurrents in P3MeT Schottky barrier diodes

The photocurrent action spectrum of a Schottky diode formed from electrodeposited poly(3-methylthiophene) is shown to follow closely the UV-visible absorption spectrum. At low forward bias, the peak photocurrent asymptotes to the expected square-root dependence on total potential, V-t, across the depletion region. At high reverse bias the superlinear dependence of the photocurrent on V-t suggests that internal photoemission from the rectifying aluminium electrode may be the dominant process.

Mathematical pedagogical content knowledge in early childhood education: A study in initial teacher education in Portugal

The study aims to explore the specificity of mathematics Pedagogical Content Knowledge in Early Childhood Education Pedagogy. The pedagogy of ECE (Siraj-Blatchford, 2010) and the didactics of ECE (Pramling & Pramling-Samuelsson, 2011) suggest dimensions of knowledge that require strong content and PC knowledge of teachers. Recent studies about PCK of ECE teachers highlight similar specific dimensions: organization of educational environment and interactions with children (Lee, 2010, McCray, 2008, Rojas, 2008). The current framework for ECE Teacher Education in Portugal (since 2007) focuses both content knowledge and subject didactics. PCK has been labelled the 'great unknown' in ECE (Rojas, 2008) in traditions where the child's development is considered as the main knowledge base for ECE (Chen & McNamee, 2006, Cullen, 2005, Hedges & Cullen, 2005). We studied the perspectives of 27 initial teacher education students about knowledge for teaching and about ECE Pedagogy. We used one open-ended questionnaire and students' analysis of episodes focusing children's answers or discourse relevant for mathematics (about high numbers and square root). The questionnaire was anonymous and students’ permission to use the answers was obtained. In the questionnaire, interactions with children (62%) and organization of the educational environment (38%) are highlighted as the most important focus for the teacher. Students suggested tasks that were adult planned and oriented to further the situations presented in the episodes. Very few references to children's exploratory actions (Bonawitz et al., 2011) were made. The specificity of ECE (child initiated activities, e.g.) needs to be further developed in initial teacher education.

Manipulating the structural and electronic properties of epitaxial NaNbO3 films via strain and stoichiometry

Due to their intriguing dielectric, pyroelectric, elasto-electric, or opto-electric properties, oxide ferroelectrics are vital candidates for the fabrication of most electronics. However, these extraordinary properties exist mainly in the temperature regime around the ferroelectric phase transition, which is usually several hundreds of K away from room temperature. Therefore, the manipulation of oxide ferroelectrics, especially moving the ferroelectric transition towards room temperature, is of great interest for application and also basic research. In this thesis, we demonstrate this using examples of NaNbO3 films. We show that the transition temperature of these films can be modified via plastic strain caused by epitaxial film growth on a structurally mismatched substrate, and this strain can be fixed by controlling the stoichiometry. The structural and electronic properties of Na1+xNbO3+δ thin films are carefully examined by among others XRD (e.g. RSM) and TEM and cryoelectronic measurements. Especially the electronic features are carefully analyzed via specially developed interdigitated electrodes in combination with integrated temperature sensor and heater. The electronic data are interpreted using existing as well as novel theories and models, they are proved to be closely correlated to the structural characteristics. The major results are: -Na1+xNbO3+δ thin films can be grown epitaxially on (110)NdGaO3 with a thickness up to 140 nm (thicker films have not been studied). Plastic relaxation of the compressive strain sets in when the thickness of the film exceeds approximately 10 – 15 nm. Films with excess Na are mainly composed of NaNbO3 with minor contribution of Na3NbO4. The latter phase seems to form nanoprecipitates that are homogeneously distributed in the NaNbO3 film which helps to stabilize the film and reduce the relaxation of the strain. -For the nominally stoichiometric films, the compressive strain leads to a broad and frequency-dispersive phase transition at lower temperature (125 – 147 K). This could be either a new transition or a shift in temperature of a known transition. Considering the broadness and frequency dispersion of the transition, this is actually a transition from the dielectric state at high temperature to a relaxor-type ferroelectric state at low temperature. The latter is based on the formation of polar nano-regions (PNRs). Using the electric field dependence of the freezing temperature, allows a direct estimation of the volume (70 to 270 nm3) and diameter (5.2 to 8 nm, spherical approximation) of the PNRs. The values confirm with literature values which were measured by other technologies. -In case of the off-stoichiometric samples, we observe again the classical ferroelectric behavior. However, the thermally hysteretic phase transition which is observed around 620 – 660 K for unstrained material is shifted to room temperature due to the compressive strain. Beside to the temperature shift, the temperature dependence of the permittivity is nearly identical for strained and unstrained materials. -The last but not least, in all cases, a significant anisotropy in the electronic and structural properties is observed which arises automatically from the anisotropic strain caused by the orthorhombic structure of the substrate. However, this anisotropy cannot be explained by the classical model which tries to fit an orthorhombic film onto an orthorhombic substrate. A novel “square lattice” model in which the films adapt a “square” shaped lattice in the plane of the film during the epitaxial growth at elevated temperature (~1000 K) nicely explains the experimental results. In this thesis we sketch a way to manipulate the ferroelectricity of NaNbO3 films via strain and stoichiometry. The results indicate that compressive strain which is generated by the epitaxial growth of the film on mismatched substrate is able to reduce the ferroelectric transition temperature or induce a phase transition at low temperature. Moreover, by adding Na in the NaNbO3 film a secondary phase Na3NbO4 is formed which seems to stabilize the main phase NaNbO3 and the strain and, thus, is able to engineer the ferroelectric behavior from the expected classical ferroelectric for perfect stoichiometry to relaxor-type ferroelectric for slightly off-stoichiometry, back to classical ferroelectric for larger off-stoichiometry. Both strain and stoichiometry are proven as perfect methods to optimize the ferroelectric properties of oxide films.

Atomic Charge Transfer-counter Polarization Effects Determine Infrared Ch Intensities Of Hydrocarbons: A Quantum Theory Of Atoms In Molecules Model.

Atomic charge transfer-counter polarization effects determine most of the infrared fundamental CH intensities of simple hydrocarbons, methane, ethylene, ethane, propyne, cyclopropane and allene. The quantum theory of atoms in molecules/charge-charge flux-dipole flux model predicted the values of 30 CH intensities ranging from 0 to 123 km mol(-1) with a root mean square (rms) error of only 4.2 km mol(-1) without including a specific equilibrium atomic charge term. Sums of the contributions from terms involving charge flux and/or dipole flux averaged 20.3 km mol(-1), about ten times larger than the average charge contribution of 2.0 km mol(-1). The only notable exceptions are the CH stretching and bending intensities of acetylene and two of the propyne vibrations for hydrogens bound to sp hybridized carbon atoms. Calculations were carried out at four quantum levels, MP2/6-311++G(3d,3p), MP2/cc-pVTZ, QCISD/6-311++G(3d,3p) and QCISD/cc-pVTZ. The results calculated at the QCISD level are the most accurate among the four with root mean square errors of 4.7 and 5.0 km mol(-1) for the 6-311++G(3d,3p) and cc-pVTZ basis sets. These values are close to the estimated aggregate experimental error of the hydrocarbon intensities, 4.0 km mol(-1). The atomic charge transfer-counter polarization effect is much larger than the charge effect for the results of all four quantum levels. Charge transfer-counter polarization effects are expected to also be important in vibrations of more polar molecules for which equilibrium charge contributions can be large.

Critical behavior of the susceptible-infected-recovered model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c(0)=0.176 500 5(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda(c)=(1-c(0))/c(0)=4.665 71(3) and a net transmissibility of (1-c(0))/(1+3c(0))=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

Statistics of opinion domains of the majority-vote model on a square lattice

The existence of juxtaposed regions of distinct cultures in spite of the fact that people's beliefs have a tendency to become more similar to each other's as the individuals interact repeatedly is a puzzling phenomenon in the social sciences. Here we study an extreme version of the frequency-dependent bias model of social influence in which an individual adopts the opinion shared by the majority of the members of its extended neighborhood, which includes the individual itself. This is a variant of the majority-vote model in which the individual retains its opinion in case there is a tie among the neighbors' opinions. We assume that the individuals are fixed in the sites of a square lattice of linear size L and that they interact with their nearest neighbors only. Within a mean-field framework, we derive the equations of motion for the density of individuals adopting a particular opinion in the single-site and pair approximations. Although the single-site approximation predicts a single opinion domain that takes over the entire lattice, the pair approximation yields a qualitatively correct picture with the coexistence of different opinion domains and a strong dependence on the initial conditions. Extensive Monte Carlo simulations indicate the existence of a rich distribution of opinion domains or clusters, the number of which grows with L(2) whereas the size of the largest cluster grows with ln L(2). The analysis of the sizes of the opinion domains shows that they obey a power-law distribution for not too large sizes but that they are exponentially distributed in the limit of very large clusters. In addition, similarly to other well-known social influence model-Axelrod's model-we found that these opinion domains are unstable to the effect of a thermal-like noise.

Macroscopic root water uptake distribution using a matric flux potential approach

Hydrological models featuring root water uptake usually do not include compensation mechanisms such that reductions in uptake from dry layers are compensated by an increase in uptake from wetter layers. We developed a physically based root water uptake model with an implicit compensation mechanism. Based on an expression for the matric flux potential (M) as a function of the distance to the root, and assuming a depth-independent value of M at the root surface, uptake per layer is shown to be a function of layer bulk M, root surface M, and a weighting factor that depends on root length density and root radius. Actual transpiration can be calculated from the sum of layer uptake rates. The proposed reduction function (PRF) was built into the SWAP model, and predictions were compared to those made with the Feddes reduction function (FRF). Simulation results were tested against data from Canada (continuous spring wheat [(Triticum aestivum L.]) and Germany (spring wheat, winter barley [Hordeum vulgare L.], sugarbeet [Beta vulgaris L.], winter wheat rotation). For the Canadian data, the root mean square error of prediction (RMSEP) for water content in the upper soil layers was very similar for FRF and PRF; for the deeper layers, RMSEP was smaller for PRF. For the German data, RMSEP was lower for PRF in the upper layers and was similar for both models in the deeper layers. In conclusion, but dependent on the properties of the data sets available for testing,the incorporation of the new reduction function into SWAP was successful, providing new capabilities for simulating compensated root water uptake without increasing the number of input parameters or degrading model performance.

A split-pot experiment with sorghum to test a root water uptake partitioning model

Correct modeling of root water uptake partitioning over depth is an important issue in hydrological and crop growth models. Recently a physically based model to describe root water uptake was developed at single root scale and upscaled to the root system scale considering a homogeneous distribution of roots per soil layer. Root water uptake partitioning is calculated over soil layers or compartments as a function of respective soil hydraulic conditions, specifically the soil matric flux potential, root characteristics and a root system efficiency factor to compensate for within-layer root system heterogeneities. The performance of this model was tested in an experiment performed in two-compartment split-pot lysimeters with sorghum plants. The compartments were submitted to different irrigation cycles resulting in contrasting water contents over time. The root system efficiency factor was determined to be about 0.05. Release of water from roots to soil was predicted and observed on several occasions during the experiment; however, model predictions suggested root water release to occur more often and at a higher rate than observed. This may be due to not considering internal root system resistances, thus overestimating the ease with which roots can act as conductors of water. Excluding these erroneous predictions from the dataset, statistical indices show model performance to be of good quality.