Stationary Eddies and Zonal Variations of the Global Hydrological Cycle in a Changing Climate

This thesis advances our physical understanding of the sensitivity of the hydrological cycle to global warming. Specifically, it focuses on changes in the longitudinal (zonal) variation of precipitation minus evaporation (P - E), which is predominantly controlled by planetary-scale stationary eddies. By studying idealized general circulation model (GCM) experiments with zonally varying boundary conditions, this thesis examines the mechanisms controlling the strength of stationary-eddy circulations and their role in the hydrological cycle. The overarching goal of this research is to understand the cause of changes in regional P - E with global warming. An understanding of such changes can be useful for impact studies focusing on water availability, ecosystem management, and flood risk.

Based on a moisture-budget analysis of ERA-Interim data, we establish an approximation for zonally anomalous P - E in terms of surface moisture content and stationary-eddy vertical motion in the lower troposphere. Part of the success of this approximation comes from our finding that transient-eddy moisture fluxes partially cancel the effect of stationary-eddy moisture advection, allowing divergent circulations to dominate the moisture budget. The lower-tropospheric vertical motion is related to horizontal motion in stationary eddies by Sverdrup and Ekman balance. These moisture- and vorticity-budget balances also hold in idealized and comprehensive GCM simulations across a range of climates.

By examining climate changes in the idealized and comprehensive GCM simulations, we are able to show the utility of the vertical motion P - E approximation for splitting changes in zonally anomalous P - E into thermodynamic and dynamic components. Shifts in divergent stationary-eddy circulations dominate changes in zonally anomalous P - E. This limits the local utility of the "wet gets wetter, dry gets drier” idea, where existing P - E patterns are amplified with warming by the increase in atmospheric moisture content, with atmospheric circulations held fixed. The increase in atmospheric moisture content manifests instead in an increase in the amplitude of the zonally anomalous hydrological cycle as measured by the zonal variance of P - E. However, dynamic changes, particularly the slowdown of divergent stationary-eddy circulations, limit the strengthening of the zonally anomalous hydrological cycle. In certain idealized cases, dynamic changes are even strong enough to reverse the tendency towards "wet gets wetter, dry gets drier” with warming.

Motivated by the importance of stationary-eddy vertical velocities in the moisture budget analysis, we examine controls on the amplitude of stationary eddies across a wide range of climates in an idealized GCM with simple topographic and ocean-heating zonal asymmetries. An analysis of the thermodynamic equation in the vicinity of topographic forcing reveals the importance of on-slope surface winds, the midlatitude isentropic slope, and latent heating in setting the amplitude of stationary waves. The response of stationary eddies to climate change is determined primarily by the strength of zonal surface winds hitting the mountain. The sensitivity of stationary-eddies to this surface forcing increases with climate change as the slope of midlatitude isentropes decreases. However, latent heating also plays an important role in damping the stationary-eddy response, and this damping becomes stronger with warming as the atmospheric moisture content increases. We find that the response of tropical overturning circulations forced by ocean heat-flux convergence is described by changes in the vertical structure of moist static energy and deep convection. This is used to derive simple scalings for the Walker circulation strength that capture the monotonic decrease with warming found in our idealized simulations.

Through the work of this thesis, the advances made in understanding the amplitude of stationary-waves in a changing climate can be directly applied to better understand and predict changes in the zonally anomalous hydrological cycle.

I. Efforts toward the synthesis of aliphatic iodonium salts. II. Flourine-19 nuclear magnetic resonance spectroscopy of cyclic and bicyclic systems

The synthesis of iodonium salts of the general formula [C₆H₅IR]⁺X^-, where R is an alkyl group and x- is a stabilizing anion, was attempted. For the choice of R three groups were selected, whose derivatives are known to be sluggish in S_N1 and S_N2 substitutions: cyclopropyl, 7, 7 -dimethyl-1-norbornyl, and 9 -triptycyl. The synthetic routes followed along classical lines which have been exploited in recent years by Beringer and students. Ultimately, the object of the present study was to study the reactions of the above salts with nucleophiles. In none of the three cases, however, was it possible to isolate a stable salt. A thermodynamic argument suggests that this must be due to kinetic instability rather than thermodynamic instability. Only iodocyclopropane and 1-iodoapocamphane formed isolable iododichlorides.

Several methylated 2, 2-difluoronorbornanes were prepared with the intent of correlating fluorine -19 chemical shifts with geometric features in a rigid system. The effect of a methyl group on the shielding of a β -fluorine is dependent upon the dihedral angle; the maximum effect (an upfield shift of the resonance) occurs at 0° and 180°, whereas almost no effect is felt at a dihedral angle of 120°. The effect of a methyl group on a γ -fluorine is to strongly shift the resonance downfield when fluorine and methyl group are in a 1, 3 - diaxial-like relationship. Molecular orbital calculations of fluorine shielding in a variety of molecules were carried out using the formalism developed by Pople; the results are, at best, in modest agreement with experiment.

Dynamical models for the Earth's geoid

The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.

Correlations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.

Temperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.

Thermodynamic and kinetic behavior of an argon plasma jet. Decomposition of nitric oxide between 1300 - 1750 degrees K in an argon plasma

In the first part of the study, an RF coupled, atmospheric pressure, laminar plasma jet of argon was investigated for thermodynamic equilibrium and some rate processes.

Improved values of transition probabilities for 17 lines of argon I were developed from known values for 7 lines. The effect of inhomogeneity of the source was pointed out.

The temperatures, T, and the electron densities, n_e , were determined spectroscopically from the population densities of the higher excited states assuming the Saha-Boltzmann relationship to be valid for these states. The axial velocities, v_z, were measured by tracing the paths of particles of boron nitride using a three-dimentional mapping technique. The above quantities varied in the following ranges: 10¹² ˂ n_e ˂ 10¹⁵ particles/cm³, 3500 ˂ T ˂ 11000 °K, and 200 ˂ v_z ˂ 1200 cm/sec.

The absence of excitation equilibrium for the lower excitation population including the ground state under certain conditions of T and n_e was established and the departure from equilibrium was examined quantitatively. The ground state was shown to be highly underpopulated for the decaying plasma.

Rates of recombination between electrons and ions were obtained by solving the steady-state equation of continuity for electrons. The observed rates were consistent with a dissociative-molecular ion mechanism with a steady-state assumption for the molecular ions.

In the second part of the study, decomposition of NO was studied in the plasma at lower temperatures. The mole fractions of NO denoted by x_NO were determined gas-chromatographically and varied between 0.0012 ˂ x_NO ˂ 0.0055. The temperatures were measured pyrometrically and varied between 1300 ˂ T ˂ 1750°K. The observed rates of decomposition were orders of magnitude greater than those obtained by the previous workers under purely thermal reaction conditions. The overall activation energy was about 9 kcal/g mol which was considerably lower than the value under thermal conditions. The effect of excess nitrogen was to reduce the rate of decomposition of NO and to increase the order of the reaction with respect to NO from 1.33 to 1.85. The observed rates were consistent with a chain mechanism in which atomic nitrogen and oxygen act as chain carriers. The increased rates of decomposition and the reduced activation energy in the presence of the plasma could be explained on the basis of the observed large amount of atomic nitrogen which was probably formed as the result of reactions between excited atoms and ions of argon and the molecular nitrogen.

Photorefractive performance of a novel multifunctional inorganic-organic hybridized nanocomposite sensitized by CdS nanoparticles

A novel multifunctional inorganic-organic photorefractive (PR) poly(N-vinyl)-3-[p-nitrophenylazolcarbazolyl-CdS nanocomposites with different molar ratios of US to poly(N-vinyl)-3-[p-nitrophenylazo]carbazolyl (PVNPAK) were synthesized via a postazo-coupling reaction and chemically hybridized approach, respectively. The nanocomposites are highly soluble and could be obtained as film-forming materials with appreciably high molecular weights and low glass transition temperature (T,) due to the flexible spacers. The PVNPAK matrix possesses a highest-occupied molecular orbital value of about -5.36 eV determined from cyclic voltammetry. Second harmonic generation (SHG) could be observed in PVNPAK film without any poling procedure and 4.7 pm/V of effective second-order nonlinear optical susceptibility is obtained. The US particles as photosensitizers had a nanoscale size in PVNPAK adopting transmission electron microscopy. The improvement of interface quality between US and polymer matrix is responsible for efficient photoinduced charge generation efficiency in the nanocomposites. An asymmetric optical energy exchange between two beams on the polymer composites PVNPAK-CdS/ECZ has been found even without an external field in two-beam coupling (TBC) experiment, and the TBC gain and diffraction efficiency of 14.26 cm(-1) and 3.4% for PVNPAK-5-CdS/ECZ, 16.43 cm(-1) and 4.4% for PVNPAK-15-CdS/ECZ were measured at a 647.1 nm wavelength, respectively.

Miniature pulse compressor of deep-etched gratings

We propose a miniature pulse compressor that can be used to compensate the group velocity dispersion that is produced by a commercial femtosecond laser cavity. The compressor is composed of two identical highly efficient deep-etched transmissive gratings. Compared with prism pairs, highly efficient deep-etched transmissive grating pairs are lightweight and small. With an optimized groove depth and a duty cycle, 98% diffraction efficiency of the -1 transmissive order can be achieved at a wavelength of 800 nm under Littrow conditions. The deep-etched gratings are fabricated in fused silica by inductively coupled plasma etching. With a pair of the fabricated gratings, the input positively chirped 73.9 fs pulses are neatly compressed into the nearly Fourier transform-limited 43.2 fs pulses. The miniature deep-etched grating-based pulse compressor should be of interest for practical applications. (c) 2008 Optical Society of America

I. The structure and properties of closed circular duplex DNA. II. Chemically interacting systems at equilibrium in a buoyant density gradient

I. The binding of the intercalating dye ethidium bromide to closed circular SV 40 DNA causes an unwinding of the duplex structure and a simultaneous and quantitatively equivalent unwinding of the superhelices. The buoyant densities and sedimentation velocities of both intact (I) and singly nicked (II) SV 40 DNAs were measured as a function of free dye concentration. The buoyant density data were used to determine the binding isotherms over a dye concentration range extending from 0 to 600 µg/m1 in 5.8 M CsCl. At high dye concentrations all of the binding sites in II, but not in I, are saturated. At free dye concentrations less than 5.4 µg/ml, I has a greater affinity for dye than II. At a critical amount of dye bound I and II have equal affinities, and at higher dye concentration I has a lower affinity than II. The number of superhelical turns, τ, present in I is calculated at each dye concentration using Fuller and Waring's (1964) estimate of the angle of duplex unwinding per intercalation. The results reveal that SV 40 DNA I contains about -13 superhelical turns in concentrated salt solutions.

The free energy of superhelix formation is calculated as a function of τ from a consideration of the effect of the superhelical turns upon the binding isotherm of ethidium bromide to SV 40 DNA I. The value of the free energy is about 100 kcal/mole DNA in the native molecule. The free energy estimates are used to calculate the pitch and radius of the superhelix as a function of the number of superhelical turns. The pitch and radius of the native I superhelix are 430 Å and 135 Å, respectively.

A buoyant density method for the isolation and detection of closed circular DNA is described. The method is based upon the reduced binding of the intercalating dye, ethidium bromide, by closed circular DNA. In an application of this method it is found that HeLa cells contain in addition to closed circular mitochondrial DNA of mean length 4.81 microns, a heterogeneous group of smaller DNA molecules which vary in size from 0.2 to 3.5 microns and a paucidisperse group of multiples of the mitochondrial length.

II. The general theory is presented for the sedimentation equilibrium of a macromolecule in a concentrated binary solvent in the presence of an additional reacting small molecule. Equations are derived for the calculation of the buoyant density of the complex and for the determination of the binding isotherm of the reagent to the macrospecies. The standard buoyant density, a thermodynamic function, is defined and the density gradients which characterize the four component system are derived. The theory is applied to the specific cases of the binding of ethidium bromide to SV 40 DNA and of the binding of mercury and silver to DNA.

Rydberg states of diatomic and polyatomic molecules using model potentials

A simple model potential is used to calculate Rydberg series for the molecules: nitrogen, oxygen, nitric oxide, carbon monoxide, carbon dioxide, nitrogen dioxide, nitrous oxide, acetylene, formaldehyde, formic acid, diazomethane, ketene, ethylene, allene, acetaldehyde, propyne, acrolein, dimethyl ether, 1, 3-butadiene, 2-butene, and benzene. The model potential for a molecule is taken as the sum of atomic potentials, which are calibrated to atomic data and contain no further parameters. Our results agree with experimentally measured values to within 5-10% in all cases. The results of these calculations are applied to many unresolved problems connected with the above molecules. Some of the more notable of these problems are the reassignment of states in carbon monoxide, the first ionization potential of nitrogen dioxide, the interpretation of the V state in ethylene, and the mystery bands in substituted ethylenes, the identification of the R and R’ series in benzene and the determination of the orbital scheme in benzene from electron impact data.

Bicyclo[3.2.0]hepta-1,4,6-triene: synthesis, thermal rearrangement, and anion formation

In order to determine the properties of the bicycloheptatrienyl anion (Ia) (predicted to be conjugatively stabilized by Hückel Molecular Orbital Theory) the neutral precursor, bicyclo[3. 2. 0] hepta-1, 4, 6-triene (I) was prepared by the following route.

Reaction of I with potassium-t-butoxide, potassium, or lithium dicyclohexylamide gave anion Ia in very low yield. Reprotonation of I was found to occur solely at the 1 or 5 position to give triene II, isolated as to its dimers.

A study of the acidity of I and of other conjugated hydrocarbons by means of ion cyclotron resonance spectroscopy resulted in determination of the following order of relative acidities:

H₂S ˃ C₅H₆ ˃ CH₃NO₂ ˃ 1, 4- C₅H₈ ˃ I ˃ C₂H₅OH ˃ H₂O; cyclo-C₇H₈ ˃ C₂ H₅OH; фCH₃ ˃ CH₃OH

In addition, limits for the proton affinities of the conjugate bases were determined:

350 kcal/mole ˂ PA(C₅ H₅^-) ˂ 360 kcal/mole

362 kcal/mole ˂ PA(C₅H₇^-, Ia, cyclo-C₇H₇^-) ˂ 377 kcal/mole PA(фCH₂^-) ˂ 385 kcal/mole

Gas phase kinetics of the trans-XVIII to I transformation gave the following activation parameters: E_a = 43.0 kcal/mole, log A = 15.53 and ∆S^ǂ (220°) = 9.6 cu. The results were interpreted as indicating initial 1,2 bond cleavage to give the 1,3-diradical which closed to I. Similar studies on cis-XVIII gave results consistent with a surface component to the reaction (E_a = 22.7 kcal/mole; log A = 9.23, ∆S^ǂ (119°) = -18.9 eu).

The low pressure (0.01 to 1 torr) pyrolysis of trans-XVIII gave in addition to I, fulvenallene (LV), ethynylcyclopentadiene (LVI) and heptafulvalene (LVII). The relative ratios of the C₇H₆ isomers were found to be dependent upon temperature and pressure, higher relative pressure and lower temperatures favoring formation of I. The results were found to be consistent with the intermediacy of vibrationally excited I and subsequent reaction to give LV and LVI.

Finite-Element Simulations of Earthquakes

This thesis discusses simulations of earthquake ground motions using prescribed ruptures and dynamic failure. Introducing sliding degrees of freedom led to an innovative technique for numerical modeling of earthquake sources. This technique allows efficient implementation of both prescribed ruptures and dynamic failure on an arbitrarily oriented fault surface. Off the fault surface the solution of the three-dimensional, dynamic elasticity equation uses well known finite-element techniques. We employ parallel processing to efficiently compute the ground motions in domains containing millions of degrees of freedom.

Using prescribed ruptures we study the sensitivity of long-period near-source ground motions to five earthquake source parameters for hypothetical events on a strike-slip fault (M_w 7.0 to 7.1) and a thrust fault (M_w 6.6 to 7.0). The directivity of the ruptures creates large displacement and velocity pulses in the ground motions in the forward direction. We found a good match between the severity of the shaking and the shape of the near-source factor from the 1997 Uniform Building Code for strike-slip faults and thrust faults with surface rupture. However, for blind thrust faults the peak displacement and velocities occur up-dip from the region with the peak near-source factor. We assert that a simple modification to the formulation of the near-source factor improves the match between the severity of the ground motion and the shape of the near-source factor.

For simulations with dynamic failure on a strike-slip fault or a thrust fault, we examine what constraints must be imposed on the coefficient of friction to produce realistic ruptures under the application of reasonable shear and normal stress distributions with depth. We found that variation of the coefficient of friction with the shear modulus and the depth produces realistic rupture behavior in both homogeneous and layered half-spaces. Furthermore, we observed a dependence of the rupture speed on the direction of propagation and fluctuations in the rupture speed and slip rate as the rupture encountered changes in the stress field. Including such behavior in prescribed ruptures would yield more realistic ground motions.

Microscopic Origin of Macroscopic Strength in Granular Media: A Numerical and Analytical Approach

Constitutive modeling in granular materials has historically been based on macroscopic experimental observations that, while being usually effective at predicting the bulk behavior of these type of materials, suffer important limitations when it comes to understanding the physics behind grain-to-grain interactions that induce the material to macroscopically behave in a given way when subjected to certain boundary conditions.

The advent of the discrete element method (DEM) in the late 1970s helped scientists and engineers to gain a deeper insight into some of the most fundamental mechanisms furnishing the grain scale. However, one of the most critical limitations of classical DEM schemes has been their inability to account for complex grain morphologies. Instead, simplified geometries such as discs, spheres, and polyhedra have typically been used. Fortunately, in the last fifteen years, there has been an increasing development of new computational as well as experimental techniques, such as non-uniform rational basis splines (NURBS) and 3D X-ray Computed Tomography (3DXRCT), which are contributing to create new tools that enable the inclusion of complex grain morphologies into DEM schemes.

Yet, as the scientific community is still developing these new tools, there is still a gap in thoroughly understanding the physical relations connecting grain and continuum scales as well as in the development of discrete techniques that can predict the emergent behavior of granular materials without resorting to phenomenology, but rather can directly unravel the micro-mechanical origin of macroscopic behavior.

In order to contribute towards closing the aforementioned gap, we have developed a micro-mechanical analysis of macroscopic peak strength, critical state, and residual strength in two-dimensional non-cohesive granular media, where typical continuum constitutive quantities such as frictional strength and dilation angle are explicitly related to their corresponding grain-scale counterparts (e.g., inter-particle contact forces, fabric, particle displacements, and velocities), providing an across-the-scale basis for better understanding and modeling granular media.

In the same way, we utilize a new DEM scheme (LS-DEM) that takes advantage of a mathematical technique called level set (LS) to enable the inclusion of real grain shapes into a classical discrete element method. After calibrating LS-DEM with respect to real experimental results, we exploit part of its potential to study the dependency of critical state (CS) parameters such as the critical state line (CSL) slope, CSL intercept, and CS friction angle on the grain's morphology, i.e., sphericity, roundness, and regularity.

Finally, we introduce a first computational algorithm to ``clone'' the grain morphologies of a sample of real digital grains. This cloning algorithm allows us to generate an arbitrary number of cloned grains that satisfy the same morphological features (e.g., roundness and aspect ratio) displayed by their real parents and can be included into a DEM simulation of a given mechanical phenomenon. In turn, this will help with the development of discrete techniques that can directly predict the engineering scale behavior of granular media without resorting to phenomenology.

Part I: Total yield measurements for the ^(21)Ne(α,n)^)(24)Mg reaction. Part II: A study of the ^(12)C(He,p)^(14)N reaction

PART I

The total cross-section for the reaction ²¹Ne(α, n)²⁴Mg has been measured in the energy range 1.49 Mev ≤ E_cm ≤ 2.6 Mev. The cross-section factor, S(O), for this reaction has been determined, by means of an optical model calculation, to be in the range 1.52 x 10¹² mb-Mev to 2.67 x 10¹² mb-Mev, for interaction radii in the range 5.0 fm to 6.6 fm. With S(O) ≈ 2 x 10¹² mb-Mev, the reaction ²¹Ne(α, n)²⁴Mg can produce a large enough neutron flux to be a significant astrophysical source of neutrons.

PART II

The reaction¹²C(³He, p)¹⁴N has been studied over the energy range 12 Mev ≤ E_lab ≤ 18 Mev. Angular distributions of the proton groups leading to the lowest seven levels in ¹⁴N were obtained.

Distorted wave calculations, based on two-nucleon transfer theory, were performed, and were found to be reliable for obtaining the value of the orbital angular momentum transferred. The present work shows that such calculations do not yield unambiguous values for the spectroscopic factors.

Motion simulation of annular polishing

Based on the Coulomb friction model, the frictional motion model of workpiece relating to the polishing pad was presented in annular polishing. By the dynamic analysis software, the model was simulated and analysed. The conclusions from the results were that the workpiece did not rotate steadily. When the angular velocity of ring and the direction were the same as that of the polishing pad, the angular velocity of workpiece hoicked at the beginning and at the later stage were the same as that of the polishing pad before contacting with the ring. The angular velocity of workpiece vibrated at the moment of contacting with the ring. After that the angular velocity of workpiece increased gradually and fluctuated at a given value, while the angular velocity of ring decreased gradually and also fluctuated at a given value. Since the contact between the workpiece and the ring was linear, their linear velocities and directions should be the same. But the angular velocity of workpiece was larger than that of the polishing pad on the condition that the radius of the workpiece was less than that of the ring. This did not agree with the pure translation principle and the workpiece surface could not be flat, either. Consequently, it needed to be controlled with the angular velocity of ring and the radii of the ring and the workpiece, besides friction to make the angular velocity of workpiece equal to that of the polishing pad for obtaining fine surface flatness of the workpiece. Copyright © 2007 Inderscience Enterprises Ltd.}

Optical properties of a square-lattice photonic crystal within the partial bandgap

Optical properties of a two-dimensional square-lattice photonic crystal are systematically investigated within the partial bandgap through anisotropic characteristics analysis and numerical simulation of field pattern. Using the plane-wave expansion method and Hellmann-Feynman theorem, the relationships between the incident and refracted angles for both phase and group velocities are calculated to analyze light propagation from air to photonic crystals. Three kinds of flat slab focusing are summarized and demonstrated by numerical simulations using the multiple scattering method. (c) 2007 Optical Society of America

Theoretical research on noncollinear match conditions of the type I optical parametric process

Broad bandwidth group match conditions are reported for a noncollinear type I optical parametric process. The theoretical calculations corresponding to two special situations in practice were made, respectively, which are in accordance with the published experimental results. Furthermore, we provide a method to not only achieve maximal parametric bandwidth output but also match the group velocities between three waves. (c) 2006 Optical Society of America.