Extracting S-matrix poles for resonances from numerical scattering data: Type-II Pade reconstruction

We present a FORTRAN 77 code for evaluation of resonance pole positions and residues of a numerical scattering matrix element in the complex energy (CE) as well as in the complex angular momentum (CAM) planes. Analytical continuation of the S-matrix element is performed by constructing a type-II Pade approximant from given physical values (Bessis et al. (1994) [421: Vrinceanu et al. (2000) [24]; Sokolovski and Msezane (2004) [23]). The algorithm involves iterative 'preconditioning' of the numerical data by extracting its rapidly oscillating potential phase component. The code has the capability of adding non-analytical noise to the numerical data in order to select 'true' physical poles, investigate their stability and evaluate the accuracy of the reconstruction. It has an option of employing multiple-precision (MPFUN) package (Bailey (1993) [451) developed by D.H. Bailey wherever double precision calculations fail due to a large number of input partial waves (energies) involved. The code has been successfully tested on several models, as well as the F + H-2 -> HE + H, F + HD : HE + D, Cl + HCI CIH + Cl and H + D-2 -> HD + D reactions. Some detailed examples are given in the text.

Dynamics of the Indonesian seas circulation. Part II – The role of pressure head

The aim of this paper is to analyze the role of the pressure head, i.e., the difference of total pressure forces acting on the Indonesian seas waters from the western Pacific and the eastern Indian Ocean, in driving the Indonesian Throughflow (ITF) and in determining the total transport of the ITF. These questions have been discussed in the literature but no consensus has been reached. A regional model of the Indonesian seas circulation has been developed that properly resolves all major topographic features in the region. The results of model runs have been used to calculate all components of the overall momentum balance. The estimates disclose that the dynamical balance is primarily between the volume integrated Coriolis acceleration, pressure gradient and the area integral of local wind stress. It is shown that consideration of components of momentum balance in the direction of the outflow through the Indian Ocean port leads to the formulation of a diagnostic relation between total inflow transports due to the Mindanao and New Guinea Coastal Currents and the external pressure head, internal pressure head, bottom form stress, and area integrated wind stress. Based on this relation, it is concluded that the external pressure head is not the major driving force of the ITF, which is why there is no unique relation between the total transport of the ITF and the external pressure head. However, Wyrtki's suggestion to monitor the variability of the total transport of the ITF by measurement of the sea-surface-height difference between the western Pacific and the eastern Indian Ocean is validated.

Signal and variability within a Holocene peat bog - Chronological uncertainties of pollen, macrofossil and fungal proxies

A single raised bog from the eastern Netherlands has been repeatedly analysed and ¹⁴C dated over the past few decades. Here we assess the within-site variability of fossil proxy data through comparing the regional
pollen, macrofossils and non-pollen palynomorphs of four of these profiles. High-resolution chronologies were obtained using ¹⁴C dating and Bayesian age-depth modelling. Where chronologies of profiles overlap, proxy curves are compared between the profiles using greyscale graphs that visualise chronological uncertainties. Even at this small spatial scale, there is considerable variability of the fossil proxy curves. Implications regarding signal (climate) and noise (internal dynamics) of the different types of fossil proxies are discussed. Single cores are of limited value for reconstructing centennial-scale climate change, and only by combining multiple cores and proxies can we obtain a reliable understanding of past environmental change and possible forcing factors (e.g., solar variability).

Ionically-Mediated Electromechanical Hysteresis in Transition Metal Oxides

Nanoscale electromechanical activity, remanent polarization states, and hysteresis loops in paraelectric TiO2 and SrTiO3 thin films are observed using scanning probe microscopy. The coupling between the ionic dynamics and incipient ferroelectricity in these materials is analyzed using extended Landau-Ginzburg-Devonshire (LGD) theory. The possible origins of electromechanical coupling including ionic dynamics, surface-charge induced electrostriction, and ionically induced ferroelectricity are identified. For the latter, the ionic contribution can change the sign of first order LGD expansion coefficient, rendering material effectively ferroelectric. The lifetime of these ionically induced ferroelectric states is then controlled by the transport time of the mobile ionic species and well above that of polarization switching. These studies provide possible explanation for ferroelectric-like behavior in centrosymmetric transition metal oxides.

Assessment of a zoomed global model for the North Sea by comparison with a conventional nested regional model

Assessment of marine downscaling of global model simulations to the regional scale is a prerequisite for understanding ocean feedback to the atmosphere in regional climate downscaling. Major difficulties arise from the coarse grid resolution of global models, which cannot provide sufficiently accurate boundary values for the regional model. In this study, we first setup a stretched global model (MPIOM) to focus on the North Sea by shifting poles. Second, a regional model (HAMSOM) was performed with higher resolution, while the open boundary values were provided by the stretched global model. In general, the sea surface temperatures (SSTs) in the two experiments are similar. Major SST differences are found in coastal regions (root mean square difference of SST is reaching up to 2°C). The higher sea surface salinity in coastal regions in the global model indicates the general limitation of this global model and its configuration (surface layer thickness is 16 m). By comparison, the advantage of the absence of open lateral boundaries in the global model can be demonstrated, in particular for the transition region between the North Sea and Baltic Sea. On long timescales, the North Atlantic Current (NAC) inflow through the northern boundary correlates well between both model simulations (R~0.9). After downscaling with HAMSOM, the NAC inflow through the northern boundary decreases by ~10%, but the circulation in the Skagerrak is stronger in HAMSOM. The circulation patterns of both models are similar in the northern North Sea. The comparison suggests that the stretched global model system is a suitable tool for long-term free climate model simulations, and the only limitations occur in coastal regions. Regarding the regional studies focusing on the coastal zone, nested regional model can be a helpful alternative.