Investigations into the degassing and eruption mechanisms of Nyamuragira volcano, Democratic Republic of the Congo (Africa)

One of two active volcanoes in the western branch of the East African Rift, Nyamuragira (1.408ºS, 29.20ºE; 3058 m) is located in the D.R. Congo. Nyamuragira emits large amounts of SO2 (up to ~1 Mt/day) and erupts low-silica, alkalic lavas, which achieve flow rates of up to ~20 km/hr. The source of the large SO2 emissions and pre-eruptive magma conditions were unknown prior to this study, and 1994-2010 lava volumes were only recently mapped via satellite imagery, mainly due to the region’s political instability. In this study, new olivine-hosted melt inclusion volatile (H2O, CO2, S, Cl, F) and major element data from five historic Nyamuragira eruptions (1912, 1938, 1948, 1986, 2006) are presented. Melt compositions derived from the 1986 and 2006 tephra samples best represent pre-eruptive volatile compositions because these samples contain naturally glassy inclusions that underwent less post-entrapment modification than crystallized inclusions. The total amount of SO2 released from the 1986 (0.04 Mt) and 2006 (0.06 Mt) eruptions are derived using the petrologic method, whereby S contents in melt inclusions are scaled to erupted lava volumes. These amounts are significantly less than satellite-based SO2 emissions for the same eruptions (1986 = ~1 Mt; 2006 = ~2 Mt). Potential explanations for this observation are: 1) accumulation of a vapor phase within the magmatic system that is only released during eruptions, and/or 2) syn-eruptive gas release from unerupted magma. Post-1994 Nyamuragira lava volumes were not available at the beginning of this study. These flows (along with others since 1967) are mapped with Landsat MSS, TM, and ETM+, Hyperion, and ALI satellite data and combined with published flow thicknesses to derive volumes. Satellite remote sensing data was also used to evaluate Nyamuragira SO2 emissions. These results show that the most recent Nyamuragira eruptions injected SO2 into the atmosphere between 15 km (2006 eruption) and 5 km (2010 eruption). This suggests that past effusive basaltic eruptions (e.g., Laki 1783) are capable of similar plume heights that reached the upper troposphere or tropopause, allowing SO2 and resultant aerosols to remain longer in the atmosphere, travel farther around the globe, and affect global climates.

Estimation of vertical groundwater fluxes into a streambed through continuous temperature profile monitoring and the relationship of groundwater fluxes to coaster brook trout spawning habitat

We hypothesized that the spatial distribution of groundwater inflows through river bottom sediments is a critical factor associated with the selection of coaster brook trout (a life history variant of Salvelinus fontinalis,) spawning sites. An 80-m reach of the Salmon Trout River, in the Huron Mountains of the upper peninsula of Michigan, was selected to test the hypothesis based on long-term documentation of coaster brook trout spawning at this site. Throughout this site, the river is relatively similar along its length with regard to stream channel and substrate features. A monitoring well system consisting of an array of 27 wells was installed to measure subsurface temperatures underneath the riverbed over a 13-month period. The monitoring well locations were separated into areas where spawning has and has not been observed. Over 200,000 total temperature measurements were collected from 5 depths within each of the 27 monitoring wells. Temperatures within the substrate at the spawning area were generally cooler and less variable than river temperatures. Substrate temperatures in the non-spawning area were generally warmer, more variable, and closely tracked temporal variations in river temperatures. Temperature data were inverted to obtain subsurface groundwater velocities using a numerical approximation of the heat transfer equation. Approximately 45,000 estimates of groundwater velocities were obtained. Estimated velocities in the spawning and non-spawning areas confirmed that groundwater velocities in the spawning area were primarily in the upward direction, and were generally greater in magnitude than velocities in the non-spawning area. In the non-spawning area there was a greater occurrence of velocities in the downward direction, and velocity estimates were generally lesser in magnitude than in the spawning area. Both the temperature and velocity results confirm the hypothesis that spawning sites correspond to areas of significant groundwater influx to the river bed.

Effects of fixed night shift work on 24 hour blood pressure regulation, state anxiety levels and total sleep time

Nearly 22 million Americans operate as shift workers, and shift work has been linked to the development of cardiovascular disease (CVD). This study is aimed at identifying pivotal risk factors of CVD by assessing 24 hour ambulatory blood pressure, state anxiety levels and sleep patterns in 12 hour fixed shift workers. We hypothesized that night shift work would negatively affect blood pressure regulation, anxiety levels and sleep patterns. A total of 28 subjects (ages 22-60) were divided into two groups: 12 hour fixed night shift workers (n=15) and 12 hour fixed day shift workers (n=13). 24 hour ambulatory blood pressure measurements (Space Labs 90207) were taken twice: once during a regular work day and once on a non-work day. State anxiety levels were assessed on both test days using the Speilberger’s State Trait Anxiety Inventory. Total sleep time (TST) was determined using self recorded sleep diary. Night shift workers demonstrated increases in 24 hour systolic (122 ± 2 to 126 ± 2 mmHg, P=0.012); diastolic (75 ± 1 to 79 ± 2 mmHg, P=0.001); and mean arterial pressures (90 ± 2 to 94 ± 2mmHg, P<0.001) during work days compared to off days. In contrast, 24 hour blood pressures were similar during work and off days in day shift workers. Night shift workers reported less TST on work days versus off days (345 ± 16 vs. 552 ± 30 min; P<0.001), whereas day shift workers reported similar TST during work and off days (475 ± 16 minutes to 437 ± 20 minutes; P=0.231). State anxiety scores did not differ between the groups or testing days (time*group interaction P=0.248), suggesting increased 24 hour blood pressure during night shift work is related to decreased TST, not short term anxiety. Our findings suggest that fixed night shift work causes disruption of the normal sleep-wake cycle negatively affecting acute blood pressure regulation, which may increase the long-term risk for CVD.

Evolution of the Southern Kenya Rift from Miocene to present with a focus on the Magadi area

The Kenya (a.k.a., Gregory) Rift is a geologically active area located within the eastern branch of the larger East African Rift System (EARS). The study area is located in the southern Kenya Rift between 1° South and the Kenya-Tanzania border (covering approximately 1.5 square degrees, semi-centered on Lake Magadi) and is predominantly filled with extrusive igneous rocks (mostly basalts, phonolites and trachytes) of Miocene age or younger. Sediments are thin, less than 1.5Ma, and are confined to small grabens. The EARS can serve both as an analogue for ancient continental rifting and as a modern laboratory to observe the geologic processes responsible for rifting. This study demonstrates that vintage (as in older, quality maps published by the Kenya Geological Survey, that may be outdated based on newer findings) quarter-degree maps can be successfully combined with recently published data, and used to interpret satellite (mainly Landsat 7) images to produce versatile, updated digital maps. The study area has been remapped using this procedure and although it covers a large area, the mapping retains a quadrangle level of detail. Additionally, all geologic mapping elements (formations, faults, etc.) have been correlated across older map boundaries so that geologic units don't end artificially at degree boundaries within the study area. These elements have also been saved as individual digital files to facilitate future analysis. A series of maps showing the evolution of the southern Kenya rift from the Miocene to the present was created by combining the updated geologic map with age dates for geologic formations and fault displacements. Over 200 age dates covering the entire length of the Kenya Rift have been compiled for this study, and 6 paleo-maps were constructed to demonstrate the evolution of the area, starting with the eruption of the Kishalduga and Lisudwa melanephelinites onto the metamorphic basement around 15Ma. These eruptions occurred before the initial rift faulting and were followed by a massive eruption of phonolites between 13-10 Ma that covered most of the Kenya dome. This was followed by a period of relative quiescence, until the initial faulting defined the western boundary of the rift around 7Ma. The resulting graben was asymmetrical until corresponding faults to the east developed around 3Ma. The rift valley was flooded by basalts and trachytes between 3Ma and 700ka, after which the volcanic activity slowed to a near halt. Since 700ka most of the deposition has been comprised of sediments, mainly from lakes occupying the various basins in the area. The main results of this study are, in addition to a detailed interpretation of the rift development, a new geologic map that correlates dozens of formations across old map boundaries and a compilation of over 300 age dates. Specific products include paleomaps, tables of fault timing and displacement, and volume estimates of volcanic formations. The study concludes with a generalization of the present environment at Magadi including discussions of lagoon chemistry, mantle gases in relation to the trona deposit, and biology of the hot springs. Several biologic samples were collected during the 2006 field season in an attempt to characterize the organisms that are commonly seen in the present Lake Magadi environment. Samples were selected to represent the different, distinctive forms that are found in the hotsprings. Each sample had it own distinctive growth habit, and analysis showed that each was formed by a different cyanobacterial. Actual algae was rare in the collected samples, and represented by a few scattered diatoms.

Using water balance models to approximate the effects of climate change on spring catchment discharge : Mt. Hanang, Tanzania

This project addresses the potential impacts of changing climate on dry-season water storage and discharge from a small, mountain catchment in Tanzania. Villagers and water managers around the catchment have experienced worsening water scarcity and attribute it to increasing population and demand, but very little has been done to understand the physical characteristics and hydrological behavior of the spring catchment. The physical nature of the aquifer was characterized and water balance models were calibrated to discharge observations so as to be able to explore relative changes in aquifer storage resulting from climate changes. To characterize the shallow aquifer supplying water to the Jandu spring, water quality and geochemistry data were analyzed, discharge recession analysis was performed, and two water balance models were developed and tested. Jandu geochemistry suggests a shallow, meteorically-recharged aquifer system with short circulation times. Baseflow recession analysis showed that the catchment behavior could be represented by a linear storage model with an average recession constant of 0.151/month from 2004-2010. Two modified Thornthwaite-Mather Water Balance (TMWB) models were calibrated using historic rainfall and discharge data and shown to reproduce dry-season flows with Nash-Sutcliffe efficiencies between 0.86 and 0.91. The modified TMWB models were then used to examine the impacts of nineteen, perturbed climate scenarios to test the potential impacts of regional climate change on catchment storage during the dry season. Forcing the models with realistic scenarios for average monthly temperature, annual precipitation, and seasonal rainfall distribution demonstrated that even small climate changes might adversely impact aquifer storage conditions at the onset of the dry season. The scale of the change was dependent on the direction (increasing vs. decreasing) and magnitude of climate change (temperature and precipitation). This study demonstrates that small, mountain aquifer characterization is possible using simple water quality parameters, recession analysis can be integrated into modeling aquifer storage parameters, and water balance models can accurately reproduce dry-season discharges and might be useful tools to assess climate change impacts. However, uncertainty in current climate projections and lack of data for testing the predictive capabilities of the model beyond the present data set, make the forecasts of changes in discharge also uncertain. The hydrologic tools used herein offer promise for future research in understanding small, shallow, mountainous aquifers and could potentially be developed and used by water resource professionals to assess climatic influences on local hydrologic systems.

MIMO waveform design using continuous phase modulation

Hybrid MIMO Phased-Array Radar (HMPAR) is an emerging technology that combines MIMO (multiple-in, multiple-out) radar technology with phased-array radar technology. The new technology is in its infancy, but much of the theoretical work for this specific project has already been completed and is explored in great depth in [1]. A brief overview of phased-array radar systems, MIMO radar systems, and the HMPAR paradigm are explored in this paper. This report is the culmination of an effort to support research in MIMO and HMPAR utilizing a concept called intrapulse beamscan. Using intrapulse beamscan, arbitrary spatial coverage can be achieved within one MIMO beam pulse. Therefore, this report focuses on designing waveforms for MIMO radar systems with arbitrary spatial coverage using that phenomenon. With intrapulse beamscan, scanning is done through phase-modulated signal design within one pulse rather than phase-shifters in the phased array over multiple pulses. In addition to using this idea, continuous phase modulation (CPM) signals are considered for their desirable peak-to-average ratio property as well as their low spectral leakage. These MIMO waveforms are designed with three goals in mind. The first goal is to achieve flexible spatial coverage while utilizing intrapulse beamscan. As with almost any radar system, we wish to have flexibility in where we send our signal energy. The second goal is to maintain a peak-to-average ratio close to 1 on the envelope of these waveforms, ensuring a signal that is close to constant modulus. It is desired to have a radar system transmit at the highest available power; not doing so would further diminish the already very small return signals. The third goal is to ensure low spectral leakage using various techniques to limit the bandwidth of the designed signals. Spectral containment is important to avoid interference with systems that utilize nearby frequencies in the electromagnetic spectrum. These three goals are realized allowing for limitations of real radar systems. In addition to flexible spatial coverage, the report examines the spectral properties of utilizing various space-filling techniques for desired spatial areas. The space-filling techniques examined include Hilbert/Peano curves and standard raster scans.

Modeling particulate organic matter diagenesis with SED2K

Diagenesis of particulate organic matter in lake sediments consumes and produces chemical species that have significant effects on water quality, e.g. oxygen and nitrate depletion and attendant mediation of nutrient and metal recycling. A mechanistic, mass balance model (SED2K) is applied here in quantifying the time course and magnitude of sediment response to reductions in depositional fluxes of organic matter. In applying the model, direct, site-specific measurements of the sedimentation and POM rates in Onondaga Lake are used, leaving only the diagenesis coefficient (solubilization) for estimation by fit to downcore POM profiles. Model calibration is constrained by the dual requirement that both POM profiles and the time series of efflux of the products of diagenesis must be matched. Simulations point to the existence of POM preservation processes at depth, a phenomenon that may enhance the timing and magnitude of lake recovery.

Optical hysteresis of ridge waveguide magnetic garnet films and cavities in two-dimensional magneto-photonic crystal slabs

Magnetic iron garnets as well as magnetic photonic crystals are of great interests in magneto-optic applications such as isolators, current captors, circulators, TE-TM mode conversion, wavelength accordable filters, optical sensors and switches, all of which provide a promising platform for future integrated optical circuits. In the present work, two topics are studied based on magnetic iron garnet films. In the first part, the characteristics of the magnetization are investigated for ridge waveguides fabricated on (100) oriented iron garnet thin films. The magnetic response in magneto-optic waveguides patterned on epitaxial magnetic garnet films depends on the crystallographic orientation of the waveguides and the magnetic anisotropy of the material. These can be studied by polarization rotation hysteresis loops, which are related to the component of magnetization parallel to the light propagation direction and the linear birefringence. Polarization rotation hysteresis loops for low birefringence waveguides with different orientations are experimentally investigated. Asymmetric stepped curves are obtained from waveguides along, due to the large magnetocrystalline anisotropy in the plane. A model based on the free energy density is developed to demonstrate the motion of the magnetization and can be used in the design of magneto-optic devices. The second part of this thesis focuses on the design and fabrication of high-Q cavities in two-dimensional magneto-photonic crystal slabs. The device consists of a layer of silicon and a layer of iron garnet thin film. Triangular lattice elliptical air holes are patterned in the slab. The fundamental TM band gap overlaps with the first-order TE band gap from 0374~0.431(a/λ) showing that both TE and TM polarization light can be confined in the photonic crystals. A nanocavity is designed to obtain both TE and TM defect modes in the band gaps. Additional work is needed to overlap the TE and TM defect modes and obtain a high-Q cavity so as to develop miniaturized Faraday rotators.

Synthesis of novel solid materials from carbon dioxide and carbon monoxide

The increase of atmospheric CO2 has been identified as the primary cause for the observed global warming over the past century. The geological and oceanic sequestration of CO2 has issues, such as cost and leakage as well as effects on sea biota. The ideal solution should be the conversion of CO2 into useful materials. However, most processes require high energy input. Therefore, it is necessary to explore novel processes with low energy demands to convert CO2 to useful solid materials. Amorphous carbon nitride and graphone received much attention due to their unusual structures and properties as well as their potential applications. However, to date there has been no attempt to synthesize those solid materials from CO2. Lithium nitride (Li3N) and lithium imide (Li2NH) are important hydrogen storage materials. However, their optical properties and reactivity has not yet studied. This dissertation research is aimed at the synthesis of carbon nitrides and graphone from CO2 and CO via their reaction with Li3N and Li2NH. The research was focused on (1) the evaluation of Li3N and Li2NH properties, (2) thermodynamic analysis of conversion of carbon dioxide and carbon monoxide into carbon nitride and other solid materials, (3) synthesis of carbon nitride from carbon dioxide, and (4) synthesis of graphone from carbon monoxide. First, the properties of Li3N, Li2NH, and LiNH2 were investigated. The X-ray diffraction measurements revealed that heat-treatment at 500°C introduce a phase transformation of β-Li3N to α-Li3N. Furthermore, the UV-visible absorption evaluation showed that the energy gaps of α-Li3N and β-Li3N are 1.81 and 2.14 eV, respectively. The UV-visible absorption measurements also revealed that energy gaps are 3.92 eV for Li2NH and 3.93 eV for LiNH2. This thermodynamic analysis was performed to predict the reactions. It was demonstrated that the reaction between carbon dioxide and lithium nitride is thermodynamically favorable and exothermic, which can generate carbon nitride and lithium cyanamide. Furthermore, the thermodynamic calculation indicated that the reaction between carbon monoxide and lithium imide can produce graphone and lithium cyanamide along with releasing heat. Based on the above thermodynamic analysis, the experiment of CO2 and Li3N reaction and CO and Li2NH were carried out. It was found that the reaction between CO2 and Li3N is very fast and exothermic. The XRD and element analysis revealed that the products are crystal lithium cyanamide and amorphous carbon nitrides with Li2O and Li2CO3. Furthermore, TEM images showed that carbon nitrides possess layer-structure, namely, it is graphene-structured carbon nitride. It was found that the reaction between Li2NH and CO was also exothermic, which produced graphone instead of carbon nitride. The composition and structures of graphone were evaluated by XRD, element analysis, TEM observation, and Raman spectra.

A REGIONAL ASSESSMENT OF GROUNDWATER AVAILABILITY AS CONSTRAINED BY LOCAL HYDROGEOLOGY AND ENVIRONMENTAL LIMITS TO STREAMFLOW DEPLETION

Groundwater pumping from aquifers in hydraulic connection with nearby streams is known to cause adverse impacts by decreasing flows to levels below those necessary to maintain aquatic ecosystems. The recent passage of the Great Lakes--St. Lawrence River Basin Water Resources Compact has brought attention to this issue in the Great Lakes region. In particular, the legislation requires the Great Lakes states to enact measures for limiting water withdrawals that can cause adverse ecosystem impacts. This study explores how both hydrogeologic and environmental flow limitations constrain groundwater availability in the Great Lakes Basin. A methodology for calculating maximum allowable pumping rates is presented. Groundwater availability across the basin is shown to be constrained by a combination of hydrogeologic yield and environmental flow limitations varying over both local and regional scales. The results are sensitive to factors such as pumping time and streamflow depletion limits as well as streambed conductance. Understanding how these restrictions constrain groundwater usage and which hydrogeologic characteristics and spatial variables have the most influence on potential streamflow depletions has important water resources policy and management implications.

Paleomagnetic Study of the Portage Lake Volcanics exposed in the Quincy Mine

A detailed paleomagnetic and rock-magnetic investigation was conducted on thirty six basaltic flows of the ~1095 Ma Portage Lake Volcanics. The flows were sampled along the East Adit of the Quincy Mine (Hancock, MI). Thirty two flows yielded well-defined primary magnetization directions carried by magnetite. A secondary magnetization component carried by hematite was also found in twenty nine flows. After correction for serial correlation between the flows, nineteen independent mean directions were calculated. The corresponding paleomagnetic pole is located at 25.5 °N, 182.1 °W (A95 = 3.5°). The new pole overlaps with the pole from the ~1087 Ma Lake Shore Traps suggesting a standstill of the North American plate during that time period. The low angular dispersion of virtual geomagnetic poles (S = 7.9°) suggests that the flows were erupted within a short time period, or that the strength of geomagnetic secular variation was lower than that of the recent field.

Multimetallic complexes based on phosphine- and phosphine oxide- appended p-hydroquinones

Presented here, is the work done with a series of binucleating ligands based on phosphine and phosphine oxide appended p-hydroquinones and their reactions towards various metals sources. The long term goal of the project was to produce coordination polymers that would have novel electronic, magnetic, and optical properties which would be of use in the field of molecular electronics. Binucleating ligands contained a p-hydroquinone motif in which various phosphine- and phosphine oxide substituents have been placed in the ortho position relative to each of the hydroxy position were synthesized. A previously published synthetic method for such lugands utilized n-BuLi to form a phenyl lithium intermediate which was quenched with chlorodiphenylphosphine. This technique was also used to produce a ligand with diisopropylphosphine groups. Phosphine ligands, containing the same structural motif, were also generated using LDA as the lithiating agent. This technique was found to be higher yielding. Phosphine chalcogenide ligands were accessed by further oxidizing the low valent phosphorous centers with either hydrogen peroxide or with elemental sulfur. These ligands were characterized using multinuclear NMR, low and high resolution mass spectroscopy, FTIR, and single crystal X-ray diffraction. Their electrochemical properties were explored with cyclic voltammetry. The phosphine appended ligands were used in the synthesis of a several bimetallic complexes. It was found that the ligands readily reacted with NiCp2 and NiCp*2, displacing one of the cyclopentadiene (Cp) or pentamethylcyclopentadiene (Cp*) rings. A cyclopentadiene complexes, containing diisopropylphine, was readily oxidized by[FeCp2]PF6 to give a NMR silent mixed valence complex. Cyclic voltammetry of these complexes showed a number of reversible waves with a large potential separation. The mixed valence compounds also showed a large absorbance band in the NIR region which was assigned to be an intervalence charge transfer. The cyclic voltammetry and NIR spectroscopy suggest that these systems are very capable of efficient metal-to-metal charge transfer. These complexes were characterized by multinuclear NMR, single crystal X-ray diffraction, UV/VIS-NIR spectroscopy and elemental analysis. The phosphine oxide ligands were reacted with a variety of different metal sources but limited success was gained in obtaining single crystals, allowing structural characterization of these compounds. Single crystals were obtained from products generated by reacting the diphenylphosphine oxide ligand with (Bipy)Cu(NO3)2 and Cu(NO3)2. In all cases the ligand had been further oxidized to a 2,5-dihydroxy-1,4-benzoquinone motif. In the reaction between the diphenylphosphine oxide ligand and (Bipy)Cu(NO3)2 it was found that the phosphine oxide moiety was involved with intermolecular coordination leading to the formation of a one-dimensional polymer composed of a series of bimetallic complexes tethered together. When NaSbF6 was present in the reaction with (Bipy)Cu(NO3)2 a unique tetrametallic complex was formed. Here the phospine oxide moiety was oriented so that two bimetallic complexes were bound together. If only Cu(NO3)2 was present, a two-dimensional polymeric sheet was formed where the ligand was present in two different coordination modes. The electronic properties of these complexes remained to be assessed.

Regulation of woody plant development by gibberellin catabolic and signaling genes

Gibberellin (GA) is a growth promoting hormone implicated in regulating a diversity of plant processes. This dissertation examines the role of GA metabolic and signaling genes in woody plant growth and development. Transgenic modifications, expression analysis, physiological/biochemical assays, biometric measurements and histological analysis were used to understand the regulatory roles these genes play in the model woody plant, Populus. Our results highlight the importance of GA regulatory genes in woody perennial growth, including: phenology, wood formation, phenotypic plasticity, and growth/survival under field conditions. We characterize two putative Populus orthologs of the SHORT INTERNODES (SHI) gene from Arabidopsis, a negative regulator of GA signaling. RNAi-mediated suppression of Populus SHI-like genes increased several growth-related traits, including extent of xylem proliferation, in a dose-dependent manner. Three Populus genes, sharing sequence homology to the positive regulator of GA signaling gene PHOTOPERIOD-RESPONSIVE 1 (PHOR1) from Solanum, are up-regulated in GA-deficient and insensitive plants suggesting a conserved role in GA signaling. We demonstrate that Populus PHOR1-like genes have overlapping and divergent function(s). Two PHOR1-like genes are highly expressed in roots, predominantly affect root growth (e.g., morphology, starch quantity and gravitropism), and induced by short-days (SD). The other PHOR1-like gene is ubiquitously expressed with a generalized function in root and shoot development. The effects of GA catabolic and signaling genes on important traits (e.g., adaptive and productivity traits) were studied in a multi-year field trial. Transgenics overexpressing GA 2-oxidase (GA2ox) and DELLA genes showed tremendous variation in growth, form, foliage, and phenology (i.e., vegetative and reproductive). Observed gradients in trait modifications were correlated to transgene expression levels, in a manner suggesting a dose-dependent relationship. We explore GA2ox and DELLA genes involvement in mediating growth responses to immediate short-term drought stress, and SD photoperiods, signaling prolonged periods of stress (e.g., winter bud dormancy). GA2ox and DELLA genes show substantial up-regulation in response to drought and SDs. Transgenics overexpressing homologs of these genes subjected to drought and SD photoperiods show hypersensitive growth restraint and increased stress resistances. These results suggest growth cessation (i.e., dormancy) in response to adverse conditions is mediated by GA regulatory genes.

Problems in the classical calculus of variations

This dissertation concerns convergence analysis for nonparametric problems in the calculus of variations and sufficient conditions for weak local minimizer of a functional for both nonparametric and parametric problems. Newton's method in infinite-dimensional space is proved to be well-defined and converges quadratically to a weak local minimizer of a functional subject to certain boundary conditions. Sufficient conditions for global converges are proposed and a well-defined algorithm based on those conditions is presented and proved to converge. Finite element discretization is employed to achieve an implementable line-search-based quasi-Newton algorithm and a proof of convergence of the discretization of the algorithm is included. This work also proposes sufficient conditions for weak local minimizer without using the language of conjugate points. The form of new conditions is consistent with the ones in finite-dimensional case. It is believed that the new form of sufficient conditions will lead to simpler approaches to verify an extremal as local minimizer for well-known problems in calculus of variations.