Towards the Use of Radar Winds in Numerical Weather Prediction

Numerical weather prediction (NWP) models provide the basis for weather forecasting by simulating the evolution of the atmospheric state. A good forecast requires that the initial state of the atmosphere is known accurately, and that the NWP model is a realistic representation of the atmosphere. Data assimilation methods are used to produce initial conditions for NWP models. The NWP model background field, typically a short-range forecast, is updated with observations in a statistically optimal way. The objective in this thesis has been to develope methods in order to allow data assimilation of Doppler radar radial wind observations. The work has been carried out in the High Resolution Limited Area Model (HIRLAM) 3-dimensional variational data assimilation framework. Observation modelling is a key element in exploiting indirect observations of the model variables. In the radar radial wind observation modelling, the vertical model wind profile is interpolated to the observation location, and the projection of the model wind vector on the radar pulse path is calculated. The vertical broadening of the radar pulse volume, and the bending of the radar pulse path due to atmospheric conditions are taken into account. Radar radial wind observations are modelled within observation errors which consist of instrumental, modelling, and representativeness errors. Systematic and random modelling errors can be minimized by accurate observation modelling. The impact of the random part of the instrumental and representativeness errors can be decreased by calculating spatial averages from the raw observations. Model experiments indicate that the spatial averaging clearly improves the fit of the radial wind observations to the model in terms of observation minus model background (OmB) standard deviation. Monitoring the quality of the observations is an important aspect, especially when a new observation type is introduced into a data assimilation system. Calculating the bias for radial wind observations in a conventional way can result in zero even in case there are systematic differences in the wind speed and/or direction. A bias estimation method designed for this observation type is introduced in the thesis. Doppler radar radial wind observation modelling, together with the bias estimation method, enables the exploitation of the radial wind observations also for NWP model validation. The one-month model experiments performed with the HIRLAM model versions differing only in a surface stress parameterization detail indicate that the use of radar wind observations in NWP model validation is very beneficial.

Exploiting ground-based measurements of the Global Positioning System for numerical weather prediction

Data assimilation provides an initial atmospheric state, called the analysis, for Numerical Weather Prediction (NWP). This analysis consists of pressure, temperature, wind, and humidity on a three-dimensional NWP model grid. Data assimilation blends meteorological observations with the NWP model in a statistically optimal way. The objective of this thesis is to describe methodological development carried out in order to allow data assimilation of ground-based measurements of the Global Positioning System (GPS) into the High Resolution Limited Area Model (HIRLAM) NWP system. Geodetic processing produces observations of tropospheric delay. These observations can be processed either for vertical columns at each GPS receiver station, or for the individual propagation paths of the microwave signals. These alternative processing methods result in Zenith Total Delay (ZTD) and Slant Delay (SD) observations, respectively. ZTD and SD observations are of use in the analysis of atmospheric humidity. A method is introduced for estimation of the horizontal error covariance of ZTD observations. The method makes use of observation minus model background (OmB) sequences of ZTD and conventional observations. It is demonstrated that the ZTD observation error covariance is relatively large in station separations shorter than 200 km, but non-zero covariances also appear at considerably larger station separations. The relatively low density of radiosonde observing stations limits the ability of the proposed estimation method to resolve the shortest length-scales of error covariance. SD observations are shown to contain a statistically significant signal on the asymmetry of the atmospheric humidity field. However, the asymmetric component of SD is found to be nearly always smaller than the standard deviation of the SD observation error. SD observation modelling is described in detail, and other issues relating to SD data assimilation are also discussed. These include the determination of error statistics, the tuning of observation quality control and allowing the taking into account of local observation error correlation. The experiments made show that the data assimilation system is able to retrieve the asymmetric information content of hypothetical SD observations at a single receiver station. Moreover, the impact of real SD observations on humidity analysis is comparable to that of other observing systems.

Finnish Dosimetric Practice for Epithermal Neutron Beam Dosimetry in Boron Neutron Capture Therapy

Boron neutron capture therapy (BNCT) is a form of chemically targeted radiotherapy that utilises the high neutron capture cross-section of boron-10 isotope to achieve a preferential dose increase in the tumour. The BNCT dosimetry poses a special challenge as the radiation dose absorbed by the irradiated tissues consists of several dose different components. Dosimetry is important as the effect of the radiation on the tissue is correlated with the radiation dose. Consistent and reliable radiation dose delivery and dosimetry are thus basic requirements for radiotherapy. The international recommendations for are not directly applicable to BNCT dosimetry. The existing dosimetry guidance for BNCT provides recommendations but also calls for investigating for complementary methods for comparison and improved accuracy. In this thesis the quality assurance and stability measurements of the neutron beam monitors used in dose delivery are presented. The beam monitors were found not to be affected by the presence of a phantom in the beam and that the effect of the reactor core power distribution was less than 1%. The weekly stability test with activation detectors has been generally reproducible within the recommended tolerance value of 2%. An established toolkit for epithermal neutron beams for determination of the dose components is presented and applied in an international dosimetric intercomparison. The measured quantities (neutron flux, fast neutron and photon dose) by the groups in the intercomparison were generally in agreement within the stated uncertainties. However, the uncertainties were large, ranging from 3-30% (1 standard deviation), emphasising the importance of dosimetric intercomparisons if clinical data is to be compared between different centers. Measurements with the Exradin type 2M ionisation chamber have been repeated in the epithermal neutron beam in the same measurement configuration over the course of 10 years. The presented results exclude severe sensitivity changes to thermal neutrons that have been reported for this type of chamber. Microdosimetry and polymer gel dosimetry as complementary methods for epithermal neutron beam dosimetry are studied. For microdosimetry the comparison of results with ionisation chambers and computer simulation showed that the photon dose measured with microdosimetry was lower than with the two other methods. The disagreement was within the uncertainties. For neutron dose the simulation and microdosimetry results agreed within 10% while the ionisation chamber technique gave 10-30% lower neutron dose rates than the two other methods. The response of the BANG-3 gel was found to be linear for both photon and epithermal neutron beam irradiation. The dose distribution normalised to dose maximum measured by MAGIC polymer gel was found to agree well with the simulated result near the dose maximum while the spatial difference between measured and simulated 30% isodose line was more than 1 cm. In both the BANG-3 and MAGIC gel studies, the interpretation of the results was complicated by the presence of high-LET radiation.

The effect of ageing on microstructure and nanoindentation behaviour of dc magnetron sputter deposited nickel rich NiTi films

Nickel rich NiTi films were sputter deposited on p-doped Si left angle bracket1 0 0right-pointing angle bracket substrates maintained at 300 °C. The films were subsequently solution treated at 700 °C for 30 min followed by ageing at 400 and 500 °C for 5 h. The microstructure of the films was examined by TEM and these studies revealed that the NiTi films were mostly amorphous in the as-deposited condition. The subsequent solution treatment and ageing resulted in crystallization of the films with the film aged at 400 °C exhibiting nanocrystalline grains and three phases viz. B2 (austenite), R and Ni3Ti2 whereas the film aged at 500 °C shows micron sized grains and two phases viz. R and Ni3Ti2. Nanoindentation studies revealed that the nature of the load versus indentation depth response for the films aged at 400 and 500 °C was different. For the same load, the indenter penetrated to a much greater depth for the film aged at 400 °C as compared to the film aged at 500 °C. Also the ratio of the residual indentation depth (hf) to maximum indentation depth (hmax) is lower for the film aged at 400 °C as compared to the film aged at 500 °C. This was attributed to the occurrence of stress induced martensitic transformation of the B2 phase present in the film aged at 400 °C during indentation loading which results in a transformation strain in addition to the normal elastic and plastic strains and its subsequent recovery on unloading. The hardness and elastic modulus measured using the Oliver and Pharr analysis was also found to be lower for the film aged at 400 °C as compared to the film aged at 500 °C which was also primarily attributed to the same effect.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Thermodynamic Properties and Phase Diagram for the System MoO2–TiO2

The activity of molybdenum dioxide (MoO2) in the MoO2–TiO2 solid solutions was measured at 1600 K using a solid-state cell incorporating yttria-doped thoria as the electrolyte. For two compositions, the emf was also measured as a function of temperature. The cell was designed such that the emf is directly related to the activity of MoO2 in the solid solution. The results show monotonic variation of activity with composition, suggesting a complete range of solid solutions between the end members and the occurrence of MoO2 with a tetragonal structure at 1600 K. A large positive deviation from Raoult's law was found. Excess Gibbs energy of mixing is an asymmetric function of composition and can be represented by the subregular solution model of Hardy as follows.The temperature dependence of the emf for two compositions is reasonably consistent with ideal entropy of mixing. A miscibility gap is indicated at a lower temperature with the critical point characterized by Tc (K)=1560 and . Recent studies indicate that MoO2 undergoes a transition from a monoclinic to tetragonal structure at 1533 K with a transition entropy of 9.91 J·(mol·K)−1. The solid solubility of TiO2 with rutile structure in MoO2 with a monoclinic structure is negligible. These features give rise to a eutectoid reaction at 1412 K. The topology of the computed phase diagram differs significantly from that suggested by Pejryd.

Thermodynamic activities in the Pb(Zr1-XTiX)O-3 solid solution at 1373 K

Although Pb(Zr1-XTiX)O-3 solid solution is the cornerstone of the piezoelectric ceramics, there is no information in the literature on thermodynamic activities of the component phases in the solid solution. Using inter-crystalline ion exchange equilibria between Pb(Zr1-XTiX)O-3 solid solution with cubic perovskite structure and (Zr1-YTiY)O-2 solid solutions with monoclinic and tetragonal structures, activities of PbTiO3 and PbZrO3 in the perovskite solid solution have been derived at 1373 K using the modified Gibbs-Duhem integration technique of Jacob and Jeffes. Tie-lines from the cubic solid solution are skewed towards the ZrO2 corner. Activities in the zirconia-rich (Zr1-YTiY)02 solid solutions are taken from a recent emf study. The results for the perovskite solid solution at 1373 K can be represented by a sub-regular solution model: Delta G(E.M) (J mol(-1)) = X-PbTiO3 X-PbZrO3(5280X(PbTiO3) - 1980X(PbZrO3)) where Delta G(E.M) is the excess Gibbs energy of mixing of the cubic solid solution and Xi represents the mole fraction of component i. There is a significant positive deviation from ideality for PbTiO3-rich compositions and mild negative deviation near the PbZrO3 corner. The cubic solid solution is intrinsically stable against composition fluctuations at temperatures down to 840 K. The results contrast sharply with the recent calorimetric data on enthalpy of mixing which signal instability of the cubic perovskite solid solution. (C) 2007 Elsevier B.V. All rights reserved.

A study of the nanostructure of the cell wall of the tracheids of conifer xylem by x-ray scattering

The purpose of this study was to develop practical and reliable x-ray scattering methods to study the nanostructure of the wood cell wall and to use these methods to systematically study the nanostructure of Norway spruce and Scots pine grown in Finland and Sweden. Methods to determine the microfibril angle (MFA) distribution, the crystallinity of wood, and the average size of cellulose crystallites using wide-angle x-ray scattering were developed and these parameters were determined as a function of the number of the year ring. The mean MFA in Norway spruce decreases rapidly as a function of the number of the year ring and after the 7th year ring it varies between 6° and 10°. The mean MFA of Scots pine behaves the same way as the mean MFA of Norway spruce. The thickness of cellulose crystallites for Norway spruce and Scots pine appears to be constant as a function of the number of the year ring. The obtained mean values are 32 Å for Norway spruce and 31 Å for Scots pine. The length of the cellulose crystallites was also quite constant as a function of the year ring. The mean length of the crystallites for Norway spruce was 364 Å, while the standard deviation was 27 Å. The mass fraction of crystalline cellulose in wood is the crystallinity of wood and the intrinsic crystallinity of cellulose is the crystallinity of cellulose. The crystallinity of wood increases from the 2nd year ring to the 10th year ring from the pith and is constant after the 10th year ring. The crystallinity of cellulose obtained using nuclear magnetic resonance spectroscopy was 52% for both species. The crystallinity of wood and the crystallinity of cellulose behave the same way in Norway spruce and Scots pine. The methods were also applied to studies on thermally modified Scots pine wood grown in Finland. Wood is modified thermally by heating and steaming in order to improve its properties such as biological resistance and dimensional stability. Modification temperatures varied from 100 °C to 240 °C. The thermal modification increases the crystallinity of wood and the thickness of cellulose crystallites but does not influence the MFA distribution. When the modification temperature was 230 °C and time 4 h, the thickness of the cellulose crystallites increased from 31 Å to 34 Å.

A Fair Contract for Managing Water Scarcity

In public utilities, under supply constraints, fairness considerations lead to a market failure. This paper characterizes a two-period principal-agent contract for demand management, that mitigates this market failure in urban water systems. The contract is designed as an extensive form mechanism using subgame perfect Nash equilibrium (SPNE) as the solution concept. The contract is fair; and is shown to be economically efficient if, in case of deviation by the agent, the gain to the agent and the loss to the principal are small. It is shown that the assumption can be avoided in an infinite horizon contract.

Metal-spinel-corundum three phase equilibria in the system Ni-Cr-Al-O at 1373 K

The tie-lines representing the inter-crystalline ion exchange equilibria between the NiCr2O4-NiAl2O4 spinet solid solution and Cr2O3-Al2O3 corundum solid solution are determined by electron microprobe andEDAX pointcountanalysis of the oxide phases equilibrated with metallic Ni at 1373 K. The component activities in the spinet solid solution are derived from the tie-lines and thermodynamic data for Cr2O3-Al2O3 solid solution available in the literature. The Gibbs energy of mixing of the spinet solid solution calculated from the experimental data is discussed in relation to the values derived from the cation distribution models which assume random mixing of cations on both tetrahedral and octahedral sites. Positive deviation from the models is observed indicating significant positive enthalpy contribution arising form the size mismatch between Al+3 and Ni+2 ions on the tetrahedral site and Al+3, Ni+2 and Cr+3 on the octahedral site. Variation of the oxygen potential for threephase equilibrium involving metallic nickel, spinet solid solution and corundum solid solution is computed as a function of composition of the solid solutions at 1373 K. The oxygen potential exhibits a minimum at aluminum cationic fraction eta(Al)/(eta(Al) + eta(Cr)) = 0.524 in the oxide solid solutions.

A reverse Stroop task with mouse tracking

In a reverse Stroop task, observers respond to the meaning of a color word irrespective of the color in which the word is printed—for example, the word red may be printed in the congruent color (red), an incongruent color (e.g., blue), or a neutral color (e.g., white). Although reading of color words in this task is often thought to be neither facilitated by congruent print colors nor interfered with incongruent print colors, this interference has been detected by using a response method that does not give any bias in favor of processing of word meanings or processing of print colors. On the other hand, evidence for the presence of facilitation in this task has been scarce, even though this facilitation is theoretically possible. By modifying the task such that participants respond to a stimulus color word by pointing to a corresponding response word on a computer screen with a mouse, the present study investigated the possibility that not only interference but also facilitation would take place in a reverse Stroop task. Importantly, in this study, participants’ responses were dynamically tracked by recording the entire trajectories of the mouse. Arguably, this method provided richer information about participants’ performance than traditional measures such as reaction time and accuracy, allowing for more detailed (and thus potentially more sensitive) investigation of facilitation and interference in the reverse Stroop task. These trajectories showed that the mouse’s approach toward correct response words was significantly delayed by incongruent print colors but not affected by congruent print colors, demonstrating that only interference, not facilitation, was present in the current task. Implications of these findings are discussed within a theoretical framework in which the strength of association between a task and its response method plays a critical role in determining how word meanings and print colors interact in reverse Stroop tasks.

(2-Chlorobenzo[h]quinolin-3-yl)methanol

In the title molecule, C14H10ClNO, all non-H atoms are coplanar (r.m.s deviation = 0.0266 angstrom). In the crystal, symmetry-related molecules are hydrogen bonded via intermolecular O-H center dot center dot center dot O interactions, forming chains along the b axis.

1-[(2-Chloro-7-methyl-3-quinolyl)methyl]pyridin-2(1H)-one

In the title compound, C16H13ClN2O, the quinoline ring system is essentially planar, with a maximum deviation of 0.021 (2) angstrom. The pyridone ring is oriented at a dihedral angle of 85.93 (6)degrees with respect to the quinoline ring system. In the crystal structure, intermolecular C-H center dot center dot center dot O hydrogen bonds link the molecules along the b axis. Weak pi-pi stacking interactions [centroid-centroid distances = 3.7218 (9) and 3.6083 (9) angstrom] are also observed.

1-[(2-Chloro-7,8-dimethylquinolin-3-yl)methyl]pyridin-2(1H)-one

In the title compound, C17H15ClN2O, the quinoline ring system is nearly planar, with a maximum deviation from the mean plane of 0.074 (2) angstrom, and makes a dihedral angle of 81.03 (7)degrees with the pyridone ring. The crystal packing is stabilized by pi-pi stacking interactions between the pyridone and benzene rings of the quinoline ring system [centroid-centroid distance = 3.6754 (10) angstrom]. Furthermore, weak intermolecular C-H center dot center dot center dot O hydrogen bonding links molecules into supramolecular chains along [001].

Monoclinic polymorph of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe(anhydrous). Parallel packing of 3(10)-/alpha-helices and a transition of helix type

The structures of two crystal forms of Boc-Trp-Ile-Ala-Aib-Ile-Val-Aib-Leu-Aib-Pro-OMe have been determined. The triclinic form (P1, Z = 1) from DMSO/H2O crystallizes as a dihydrate (Karle, Sukumar & Balaram (1986) Proc, Natl, Acad. Sci. USA 83, 9284-9288). The monoclinic form (P2(1), Z = 2) crystallized from dioxane is anhydrous. The conformation of the peptide is essentially the same in both crystal system, but small changes in conformational angles are associated with a shift of the helix from a predominantly alpha-type to a predominantly 3(10)-type. The r.m.s. deviation of 33 atoms in the backbone and C beta positions of residues 2-8 is only 0.29 A between molecules in the two polymorphs. In both space groups, the helical molecules pack in a parallel fashion, rather than antiparallel. The only intermolecular hydrogen bonding is head-to-tail between helices. There are no lateral hydrogen bonds. In the P2(1) cell, a = 9.422(2) A, b = 36.392(11) A, c = 10.548(2) A, beta = 111.31(2) degrees and V = 3369.3 A for 2 molecules of C60H97N11O13 per cell.